l = 1;
n = 50000;
m = 0;
for k = 1: n
x = unifrnd(0,1);
y = unifrnd(0,1);
if y < sqrt(1-x^2)
m = m+1;
end
end
m/n
matlab 蒙特卡洛
最新推荐文章于 2024-04-18 20:52:58 发布
l = 1;
n = 50000;
m = 0;
for k = 1: n
x = unifrnd(0,1);
y = unifrnd(0,1);
if y < sqrt(1-x^2)
m = m+1;
end
end
m/n