import numpy as np import matplotlib.pyplot as plt import scipy.stats as st mu, sigma = 100, 15 s = st.norm(mu, sigma).rvs(1000) s_fit = np.linspace(s.min(), s.max()) plt.title('$\mu=100$,$\sigma=15$',color='r') count, bins, _ = plt.hist(s, 30, normed=True, color='g') plt.plot(s_fit, st.norm(mu, sigma).pdf(s_fit), lw=1, c='r') plt.show()