exn565.m

来自「MATLAB及在电子信息课程中的应用.感觉很有用,本来也不花钱,有想用的来拿吧」· M 代码 · 共 13 行

%《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