Untitled
unknown
plain_text
9 months ago
301 B
9
Indexable
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)Editor is loading...
Leave a Comment