I'm finding that hist has problems computing on 2d arrays.
mu, sigma = 2, 0.5
v = numpy.random.normal(mu,sigma,160000)
pylab.hist(v, bins=1000, normed=1)
This works without any problems. But if you try this:
pylab.hist(w, bins=1000, normed=1)
it doesn't come back on my machine until all of memory is used up. However:
n,bins = numpy.histogram(w,bins=1000,normed=1)
works just fine.