Plot a histogram and fit a normal pdf in python

python_normal_pdf_fit.py

plt.hist(data),normed=True,stacked=True,facecolor='0.5');
from scipy.stats import norm
mu, std = norm.fit(data)
std = float(std)
xmin, xmax = plt.xlim()
x = np.linspace(xmin, xmax, 100)
p = norm.pdf(x, mu, std)
plt.plot(x, p, 'k', linewidth=2);