El objetivo de este proceso es estimar el valor de π utilizando la probabilidad de que un punto esté dentro del círculo inscrito en un cuadrado. Esto se logra contando cuántos puntos caen dentro del círculo y calculando la proporción de puntos dentro del círculo en relación con el total de puntos generados. Esta proporción es igual a π/4 debido a la relación de áreas entre el círculo y el cuadrado y que es la probabilidad de que un punto esté dentro del circulo. Posteriormente, esta estimación se iguala a π/4 y se despeja π, lo que proporciona una aproximación del valor de π. A continuación, se presenta una descripción detallada de los pasos del código.

###Para facilitar la generación de estimaciones de π con diferentes valores de n, se ha encapsulado el proceso en una función que se puede invocar con diferentes valores de n y ademas se inicializa una variable puntos en cero para contar los puntos que caerán dentro del circulo. 

estimacion <- function(n) {
  puntos <- 0

### Lo primero que se hara es generar n puntos aleatorios con coordenadas aleatorias X, Y para cada punto en el cuadrado [0,1] por medio de una distribución uniforme

  for (i in 1:n) {
    x <- runif(1, 0, 1)
    y <- runif(1, 0, 1)

### Despues se calcula la distancia que tiene cada punto desde el centro del circulo (0.5,0.5) hasta el punto con el objetivo de comprobar si el punto está dentro del circulo (distancia <= 0.5, teniendo que en cuenta que si esta en el borde tambien se considerara dentro del circulo

    distancia <- sqrt((x-0.5)^2 + (y-0.5)^2)
    
### Finalmente, al contar el número de puntos dentro del círculo y dividirlo por el número total de puntos generados (n), se está estimando la probabilidad de que un punto aleatorio caiga dentro del círculo. Esta estimación es π/4 y a partir de esta igualdad se despeja la estimación del valor de π.

    if (distancia <= 0.5) {
      puntos <- puntos + 1
    }
  }
  
  proporcion <- puntos / n
  pi_estimacion <- proporcion * 4
  
  return(pi_estimacion)
}

Ahora probamos la funcion de estimacion para los siguientes valores de n: 1000, 10000

# Número de puntos a generar 
n <- 10000

estimacion <- estimacion(n)

cat("Estimación de π:", estimacion, "\n")
## Estimación de π: 3.1508