changed error of epsilon from uniform to normal

in simulation function I had np.rand(). Should have been np.randn().
The plot doesn’t look different, but the convergence look weird. In
some sense it doesn’t really matter, because one thing is the true
process and other thing is what the firm thinks is the true process
Added verbose option for simulate