changed time from 5000 to 50000 to more adequately cover dynamics
removed some now-obsolete scaffolding
changed from random package to numpy.random, tweaked some plots, and added kolmogorov-smirnov test with plot; takes 20 min to run