mon_car.m

来自「generate a normal sequence from rayleigh」· M 代码 · 共 13 行

%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;