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