# Número de puntos a simular
n <- 10000
# Generar n puntos para x y y
x <- runif(n, 0, 1)
y <- runif(n, 0, 1)
# Calcular la distancia al cuadrado desde el centro (0.5,0.5)
distancias <- (x - 0.5)^2 + (y - 0.5)^2
# Determinar si cada punto está dentro del círculo (distancia al cuadrado < 0.25)
dentro_circulo <- distancias < 0.25
# Contar cuántos puntos están dentro del círculo
puntos_dentro <- sum(dentro_circulo)
# Estimar π (4 veces la fracción de puntos dentro del círculo)
pi_estimado <- 4 * (puntos_dentro / n)
# Mostrar la estimación de π
print(pi_estimado)## [1] 3.1268