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()