Untitled
N <- 1e6 triangle_exists <- logical(N) #pre-allocation for(i in 1:N) { x <- runif(1, 0,1) y <- runif(1, 0,1) a <- min(x,y) b <- max(x,y) cond1 <- a<0.5 cond2 <- b>0.5 cond3 <- (b-a)<0.5 triangle_exists[i] <- (cond1) & (cond2) & (cond3) } appr_p <- mean(triangle_exists) print(appr_p)
Leave a Comment