a = rnormal(1000,1,1,0.5)