# Box Muller commands
u1=runif(1000)
u2=runif(1000)
x=sqrt(-2*log(u1))*cos(2*pi*u2)
y=sqrt(-2*log(u1))*sin(2*pi*u2)
par(mfrow=c(2,2))
hist(x)
hist(y)
qqnorm(x)
qqnorm(y)
cor(cbind(x,y))