#-------------------------------------------------
# Script per valutare pi greco tramite Monte Carlo
#
# GdA 29/2/2024
#-------------------------------------------------

N=10000; x<-runif(N); y<-runif(N)
plot(x,y,col='cyan',asp=1)
sel <- x^2+y^2 <= 1
points(x[sel],y[sel],col='red')
pi.MC <- 4*sum(sel)/N
cat(sprintf("pi greco valutato mediante MC usando %d punti: %.3f\n",
            N, pi.MC))

# eseguire lo script più volte, anche cambiando N
#
# Nota: per help sulle funzioni usate, ad es. runif(), 
# usare il comando "?funzione", ad es ?runif