symulacje

##########################################################################################
###                                                                                Realizacja algorytmu Boxa-Muellera w R                                                           ###
##########################################################################################

s=2
while (s>=1)
{

u1=runif(1)
u2=runif(1)
v1=2*u1-1
v2=2*u2-1
s=v1^2+v2^2

}
x1=v1*sqrt(-2*log(s)/s)
x2=v2*sqrt(-2*log(s)/s)
print(x1)
print(x2)