generatePoint <- function () {
  pct <- c(runif(1), runif(1))
  return (pct)
}

getDistance <- function(p1, p2)
{
  ans <- sqrt(sum((p1 - p2) ^ 2))

  return (ans)
}


getPPoints <- function (n=1000) {
  center = c(1/2,1/2)
  p <- 0
  for (i in 1:n){
    point <- generatePoint()
    distance <- getDistance(center, point)
    if (distance <= 0.5) {
      p <- p + 1
    }
  }
  return (p)
}

n <- 1000000
p <- getPPoints(n)

my_pi <- (p*4) / n
my_pi

pi # real pi
#lucky https://i.imgur.com/2zCZrUp.png