Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- set.seed(1447)
- stunif_d <- function(N=10000000){
- U1 <- runif(N)
- U2 <- runif(N)
- U3 <- runif(N)
- return(sqrt(U1) + sqrt(U2) +sqrt(U3))
- }
- U <- stunif_d(100000000)
- mean(U[U>=0.8])
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement