normal = function(mu,sigma){
 z = (-4000:4000)/1000
 x = mu + sigma*z
 Density = 1/(sigma*sqrt(2*pi)) * exp(-z^2/2)
 plot(x,Density,type="l")
}