x = np.linspace(xmin, xmax, 100)
p = norm.pdf(x, mu, std)
plt.plot(x, p, 'k', linewidth=2)
title = 'Fit results: mu = %.2f,  std = %.2f' % (mu, std)
plt.title(title)
plt.show()
Fitting a Normal distribution to 1D data