蒙特·卡罗方法计算圆周率近似值

wrj001

假设有一块边长为 2 的正方形木板，上面画一个单位圆，然后随意往木板上扔飞镖，落点坐标(x, y)必然在木板上（更多的时候是落在单位圆内），如果扔的次数足够多，那么落在单位圆内的次数除以总次数再乘以 4，这个数字会无限逼近圆周率的值。这就是蒙特·卡罗发明的用于计算圆周率近似值的方法

times := 100000; 设定投掷次数
hits := 0; 落在圆内次数
for i from 1 to 100000 do
x := rand(0 .. 100);
y := rand(0 .. 100);
a := x()/100;
b := y()/100;
if a*a+b*b <= 1 then
hits := hits+1
end if
end do;
4.0*hits/times