r - 创建标准正态分布随机变量的函数 (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()