# 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))