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