创建标准正态分布随机变量的函数（Box-Muller）

Question

我想编写一个模拟偶数个标准正态分布随机变量的函数（使用Box-Muller变换）。 我是一个初学者，并试图编写下面的功能。目前我不知道如何继续得到我想要的方式。

我目前的代码：

snorm <- function(){
  u1 <- runif(1)
  u2 <- runif(1)
  z1 <- ((((-2)*log(u1))^(1/2))*cos(2*u2*pi))
  z2 <- ((((-2)*log(u1))^(1/2))*sin(2*u2*pi))
  z <- list(z1,z2)
  return(z)
}
snorm()