%《MATLAB及其在理工課程中的應用指南》(十一五規劃版)第五章例5-6-5程序
% 蒙特卡羅方法求圓周率
% 西安電子科技大學出版社出版  陳懷琛編著 2007年7月
%
N=input('取的總點數= ');
x=2*(rand(1,N)-0.5);	% 生成均值為零，數值在-1~1均勻分布的隨機數
y=2*(rand(1,N)-0.5);   % 同上
k=0; 						%  k為落入圓內部的點數，初始值為零
for i=1:N
    if x(i)*x(i)+y(i)*y(i)<=1 k=k+1; end  % 計算落入圓中點的個數
end
pihat=4*k/N			% 圓面積＝正方形面積×k/N