Fig. 1: Estimação Monte Carlo para π

x1 <- runif(1,min=0,max=1)
x2 <- runif(1,min=0,max=1)

e verificando-se a proporção dos que se encontram no interior do círculo de diâmetro unitário (raio 0,5 e área π/4),

if ( sqrt( (x1-0.5)^2 + (x2-0.5)^2 ) <= 0.5 ) {
    pontos.no.interior <- pontos.no.interior + 1
}
else {
    pontos.no.exterior <- pontos.no.exterior + 1
}

pode-se estimar o valor de π pela expressão:

pi.est <- 4 * (pontos.no.interior)/(pontos.no.interior + pontos.no.exterior)

• uma função em R/S-PLUS para estimar o valor de π, sendo o número de pontos gerados n um parâmetro de entrada;

• uma verifição do comportamento do estimador para π em função de n.