Untitled

 avatar
unknown
plain_text
5 days ago
301 B
6
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)
Leave a Comment