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

N = 10000      # nr di punti

x <- runif(N)
y <- runif(N)

IN <- (x^2 + y^2) <= 1

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) )