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)
# wyświetlenie dwóch wartości losowych
# z rozkładu normalnego standardowego
print(x1)
print(x2)