【matlab 圆周率计算】matlab 求圆周率的两种算法实现比较

%author:laidefa
%data:2014-09-19
%丘德诺夫斯基公式求圆周率
function mpi=qdnfsj(m)
i=m;
s=13591409;
for n=1:i
A=(factorial(6*n)(13591409+54514013*n))/(factorial(3*n)*factorial(n)^3(-640320)^(3*n));
s=s+A;
end
mpi=(426880*sqrt(10005))/s;
end

% % %调用例子
% clc,clear
% mpi=qdnfsj(20);
% vpa(mpi,50)

ans =

3.1415926535897411575604110112180933356285095214844

%author:laidefa
%data:2014-09-19
%function:蒙特卡罗方法求圆周率
function [api]=montcpi(n)
format long;
m=0;
x=2*rand(n,2)-1;
for i=1:n
if x(i,1)^2+x(i,2)^2<=1
m=m+1;
end
end
api=4*m/n;

% %调用例子
% clc,clear
% [api]=montcpi(19999999);
%vpa(api,50)

ans =

3.1417683570884178401172448502620682120323181152344