from numpy import *
from matplotlib.mlab import *
from matplotlib.pyplot import *

mu, sigma = 100, 15
x = mu + sigma*random.randn(10000)
n, bins, patches = hist(x, 50, normed=1, \
   facecolor='green', alpha=0.75)

# Linie mit 'best fit' einzeichnen
y = normpdf( bins, mu, sigma)
plot(bins, y, 'r--', linewidth=1)

axis([40, 160, 0, 0.03])
show()