这是一个简单的例子函数计算π的值通过生成均匀分布的点在一个正方形的边长1和确定这些点的分数落在圆。四次,分数将π的近似值。
函数mypi = approxpi (n)
输入:% n =数量的点来生成
%默认是n = 1 e6
% n值的增大应该表现的更好
如果输入参数个数< 1
n = 1 e6;
结束
%产生均匀分布的点
% x [0, 1] [0, 1]
xy =兰德(n, 2);
%计算距离(0.5,0.5)
r =√(xy (: 1) -0.5)。^ 2 + (xy (:, 2) -0.5)。^ 2);
%计算分数1/2单位内的点(0.5,0.5)
压裂=总和(r < = 0.5) / n;
%因为广场边1,圆半径(1/2)
%和应该的π* (1/2)^ 2
%压裂大约是π/ 4π大约是4 *压裂
mypi = 4 *压裂;
另一个例子使用MATLAB蒙特卡罗模拟可以发现:
