%利用蒙特卡罗算法计算pi P=rand(10000,2); x=P(:,1); y=P(:,2); s=find(x.^2+y.^2<1); length(s)/10000*4 plot(x(s),y(s),'g.')</p>
阅读全文>>
%利用蒙特卡罗算法计算pi %思路:在0<x<1和0<y<1中随机出很多的点,然后计算出现在1/4圆内的圆的个数,得到出现在1/4圆的概率即1/4圆与1*1正方形的面积比，得到pi n=input('a big num:'); count=0; for i=1:n x(i)=rand(1); y(i)=rand(1); if x(i)^2+y(i)^2<1 count=count+1; end end 4*count/n