mon_car.m

generate a normal sequence from rayleigh distribution

%monte carlo simulation
 
function p = mon_car(N)
N = 10000;
hit = 0;
for i=1:N,
x = 2*rand-1;
y = 2*rand-1;
if x^2+y^2<1,
    hit = hit+1;
end
end
p = hit/N;