#numpy, matplotlib import numpy      from matplotlib import mlab from matplotlib import pyplot mu = 100      sigma = 15   x = mu + sigma * numpy.random.randn(1000)     num_bins = 50     n, bins, patches = pyplot.hist(x, num_bins, normed=1, facecolor='green', alpha=0.5)     y = mlab.normpdf(bins, mu, sigma)     pyplot.plot(bins, y, 'r--')     pyplot.xlabel('x-axis')     pyplot.ylabel('y-axis')     pyplot.title('$\mu=100$, $\sigma=15$')     pyplot.subplots_adjust(left=0.15)     pyplot.show()