pi_monte_carlo.r

estimatePi <- function(numSims){
    x <- runif(numSims)
    y <- runif(numSims)
    inUnitCircle <- as.integer(x^2 + y^2 <= 1)
    return(4 * sum(inUnitCircle) / numSims)
}

set.seed(1)
estimatePi(3000)
estimatePi(30000)