#
#  Quarto di cerchio:
#   valutazioni area e pi greco mediante campionamento
                                        #

N = 10000

x <- runif(N)
y <- runif(N)
plot(x, y, pch=19, cex=0.2, col='blue', asp=1)

IN <- (x^2 + y^2) <= 1
points(x[IN], y[IN], pch=19, cex=0.25, col='red')

n.IN <- sum(IN)

area.stimata =  n.IN/N
cat( sprintf("Area stimata:  %f\n", area.stimata) )

pi.stimato = n.IN/N * 4
cat( sprintf("Pi greco stimato:  %f\n", pi.stimato) )