% function M=MC_PI(N)
%
% Calcola il valore di pi-greco utilizzando il metodo montecarlo su N
% campioni.

 function m = mc_pi(n) 

     % estrae le coordinate
     X = uniform_rnd(-1,1,1,n);
     Y = uniform_rnd(-1,1,1,n);

% se il punto (X,Y) dista da (0,0) piu' di uno, allora cade fuori dal cerchio di raggio 1

D = sqrt(X.^2+Y.^2); % distanza da (0,0)
p = length(find(D<=1)); % per quanti campioni D e' minore di 1 
media_camp = p/n;

% il valore di pi-greco e' un quattro volte il valore atteso della media campionaria
m = media_camp*4;

endfunction
