# Método de Monte Carlo para estimar Pi
monte_carlo_pi <- function(iterations) {
  inside_circle <- 0
  for (i in 1:iterations) {
    x <- runif(1)
    y <- runif(1)
    if (x^2 + y^2 <= 1) {
      inside_circle <- inside_circle + 1
    }
  }
  return(4 * inside_circle / iterations)
}
 
# Número de iteraciones (puedes modificar este valor para más o menos datos)
iterations <- 10^6
 
# Temporizador de ejecución
start_time <- Sys.time()
pi_estimate <- monte_carlo_pi(iterations)
end_time <- Sys.time()
 
print(paste("Estimación de Pi:", pi_estimate))
print(paste("Tiempo de ejecución:", end_time - start_time))
 

IyBNw6l0b2RvIGRlIE1vbnRlIENhcmxvIHBhcmEgZXN0aW1hciBQaQptb250ZV9jYXJsb19waSA8LSBmdW5jdGlvbihpdGVyYXRpb25zKSB7CiAgaW5zaWRlX2NpcmNsZSA8LSAwCiAgZm9yIChpIGluIDE6aXRlcmF0aW9ucykgewogICAgeCA8LSBydW5pZigxKQogICAgeSA8LSBydW5pZigxKQogICAgaWYgKHheMiArIHleMiA8PSAxKSB7CiAgICAgIGluc2lkZV9jaXJjbGUgPC0gaW5zaWRlX2NpcmNsZSArIDEKICAgIH0KICB9CiAgcmV0dXJuKDQgKiBpbnNpZGVfY2lyY2xlIC8gaXRlcmF0aW9ucykKfQoKIyBOw7ptZXJvIGRlIGl0ZXJhY2lvbmVzIChwdWVkZXMgbW9kaWZpY2FyIGVzdGUgdmFsb3IgcGFyYSBtw6FzIG8gbWVub3MgZGF0b3MpCml0ZXJhdGlvbnMgPC0gMTBeNgoKIyBUZW1wb3JpemFkb3IgZGUgZWplY3VjacOzbgpzdGFydF90aW1lIDwtIFN5cy50aW1lKCkKcGlfZXN0aW1hdGUgPC0gbW9udGVfY2FybG9fcGkoaXRlcmF0aW9ucykKZW5kX3RpbWUgPC0gU3lzLnRpbWUoKQoKcHJpbnQocGFzdGUoIkVzdGltYWNpw7NuIGRlIFBpOiIsIHBpX2VzdGltYXRlKSkKcHJpbnQocGFzdGUoIlRpZW1wbyBkZSBlamVjdWNpw7NuOiIsIGVuZF90aW1lIC0gc3RhcnRfdGltZSkpCg==