Hi David,

How about:

import pylab as P

dist = P.rand(100)
x=P.arange(0.,1.,0.01)
fx=10.*x**4
   
P.hist(dist)
P.plot(x,fx,'r-')
P.axis([0.,1.,0.,20.])
P.show()