半径为1的圆（假定圆心为坐标原点）面积为（以下用pi表示），作圆外接正方形（各边与坐标轴平行），则正方形面积为4。设(X,Y)服从该正方形中的二维均匀分布，则随机分布点位于圆内的概率P=pi/4，于是pi=4*P，对于概率P可用蒙特卡洛方法模拟计算，下面给出一个简例。

n=input('请输入随机数个数');
x=2*(rand(1,n)-0.5);
y=2*(rand(1,n)-0.5);
r=sqrt(x.^2+y.^2);
r(r<=1)=1;
r(r>1)=0;
pai=4*sum(r)/n;