R中没有循环的蒙特卡罗模拟概率练习

Question

我需要使用Monte Carlo Simulation在R中进行此练习，而无需使用循环，这是我苦苦挣扎的原因。这是练习：让Z（n）为标准正态分布中n个观测值的最大值。估计n必须取P（Z（n）> 4）= 0.25的值。

这是我得到的所有代码

n = 10000

x = replicate(n, max(rnorm(?)))


prob = (sum(x>=4))/n


prob must be = 0.25