#-------------------------------------------------
# Script per eseguire l'integrale tramite Monte Carlo
# della funzione sin^2(x)

# GdA 29/2/2024
#-------------------------------------------------

N=10000; x<-runif(N,0,pi); y<-runif(N)
plot(x,y,col='cyan')
sel <- y <= sin(x)^2
points(x[sel],y[sel],col='red')
int.MC <- sum(sel)/N * (1*pi)
cat(sprintf("Integrale valutato mediante MC usando %d punti: %.3f\n",
            N, int.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