编写程序使用 Box-Muller 方法生成 1000 个𝑁(5,4)正态分布随机数，并绘制 正态 QQ 图检验随机数的正态性（可使用 R 函数：qqnorm）。

# 生成1000个N(5,4)正态分布随机数
set.seed(123)
u1 <- runif(1000)
u2 <- runif(1000)
z1 <- sqrt(-2*log(u1))*cos(2*pi*u2)
z2 <- sqrt(-2*log(u1))*sin(2*pi*u2)
x <- 5 + sqrt(4)*z1

# 绘制正态 QQ 图
qqnorm(x)
qqline(x)

结果图如下：


可以看出，随机数的分布与正态分布非常接近，符合正态性。