import numpy as np
import matplotlib.pyplot as plt
mu, sigma = 100, 15
np_list = np.random.normal(mu, sigma, 1000)
plt.title('$\mu=100$,$\sigma=15$',color='r')
count, bins, _ = plt.hist(np_list, 50, normed=True, color='g')
plt.plot(bins, 1./(np.sqrt(2*np.pi)*sigma)*np.exp(-(bins-mu)**2/(2*sigma**2)), lw=1, c='r')
plt.show()