MATLAB 蒙特卡罗法求π

n = [1:10000];
PI=[];
for i = n
    m=0;
    for j = 1:i
	x = rand;
	y = rand;
    
    if (x^2 + y^2) < 1   
		m = m + 1;
    end
    
    end
    p=4 * m/i;
    PI=[PI,p];
end
plot(PI)
fprintf('圆周率值是：%d',4 * m/n);  

思路：

1 蒙特卡罗法：使用类似抛针法，统计投入圆内的数量

2 建立一个抛掷次数的list，储存次数(1~10000)，用于用图像表明随着实验次数的增加，π值趋近于真实值

3 嵌套循环：对于实验次数n，我们建立一个子循环，通过if语句判断每次产生的实验结果，最后在母循环里进行统计