title: norm fit producing incorrect fit
body: <p>Why does the following produce an incorrect output for muf and stdf?</p>
  
  <pre><code>import numpy as np
  from scipy.stats import norm
  x=np.linspace(-50,50,100)
  sig=10
  mu=0
  y=1/np.sqrt(2*sig*sig*np.pi)*np.exp(-(x-mu)*(x-mu)/(2*sig*sig))
  muf, stdf = norm.fit(y)
  print muf,stdf
  </code></pre>
  
  <p>This prints
  <code>0.00989999568634 0.0134634293279</code> </p>
  
  <p>Thanks.</p>