用matlab怎么编程求pi啊 ？

求用MATLAB编写求pi的m文件 急需~

一、复化梯形法
clear;
f=inline('4./(1+x.*x)');
a=0;b=1;n=1;
h=(b-a)/n;
t1=h/2*(f(a)+f(b));
er=1;k=1;
while er>1.0e-5;
s=0;
for i=1:n
s=s+f(a+(i-1/2)*h);
end
t2=(t1+h*s)/2;
er=abs(t2-t1);
fprintf('n=%.0f,p=%.6f,r=%.6f
',k,t2,er);
n=2*n;h=h/2;t1=t2;
k=k+1;
end
二、Ramanujan公式法
clear;
n=0;
r=0;
while n<=5;
r=r+(factorial(4*n)/(factorial(n))^4)*(1103+26390*n)/396^(4*n);
pi_zhi=9801/(2*sqrt(2)*r);
vpa(pi_zhi,60)
n=n+1
end
三、泰勒级数法
clear;
n=0;
r=1;
p=0;
k=-1;
a=1;
b=1;
while r>=1.0e-9
n=n+1;
k=k*(-1);
a=4*a;b=9*b;
pl=p+k/(2*n-1)*(2/a+3/b);
r=abs(4*(pl-p));
fprintf('n=%.0f,p=%.10f
',n,4*pl);
p=pl;
end
。。。

fun=@(x,y)y.*log10(y)-x创建关于x,y的匿名函数
log10(x)由于Matlab中自然对数的符号是log(与我们数学上的不同)，以10为底的对数符号是log10，同理有log2，但是没有log3等哦
vpa将得到的符号解转化为指定精度的小数，比如vpa('sqrt(2)'，10)将sqrt(2)表示10位精度的小数
fzero是Matlab求解一元函数零点的内置函数
fsolve是Matlab中求解方程组解(可以求解复数)的内置函数，功能强于fzero
fval是返回fun(x)的值，也就是说fzero和fsolve返回方程的解x，那么此时的fun(x)就是fval

这个方法比较多 可以数值也可以解析求解

1.解析求解

>> solve('y*log10(y)=3','y')

ans =

3*log(10)/lambertw(3*log(10))


>> vpa(ans)

ans =

4.5555357051951280236911974082484

2.数值解法这个可以使用fsolve和fzero，前者的功能更加强大，但是对于这个问题fzero就可以了

>>x=3;
>>fun=@(y)y*log10(y)-x;
%使用fsolve
>> [x,fval]=fsolve(fun,rand)
Optimization terminated: first-order optimality is less than options.TolFun.

x =

4.5555


fval =

3.9620e-008

%使用fzero
>> [x,fval]=fzero(fun,5)

x =

4.5555


fval =

-8.8818e-016

你啥教育水平？如果学过高等数学的话，有若干种方法求Pi的。还记得那个欧拉证明的pi的公式嘛？就是pi^2/6=1+1/2^2+1/3^2+......
就用一个for循环就遍出来了，你循环里面的n越大，pi的值就越精确。不过不要把n改太大了，免得死机。

程序是：
n=1000;
sigma=zeros(1,n);
for i=1:n
sigma(i+1)=sigma(i)+(1/(i^2));
end
pi=max(sqrt(6*sigma));%最后pi的值就是π啦

你可以看到随着n增加，pi的值增加的过程。
不过貌似要用到vpa命令修改小数点的位数，你自己改改吧

不用编，内置的就有这个常量；名称就是

pi

下面是源代码：

function [varargout] = pi(varargin)
%PI 3.1415926535897....
% PI = 4*atan(1) = imag(log(-1)) = 3.1415926535897....

% Copyright 1984-2003 The MathWorks, Inc.
% $Revision: 5.7.4.2 $ $Date: 2004/04/16 22:05:16 $
% Built-in function.

if nargout == 0
builtin('pi', varargin{:});
else
[varargout{1:nargout}] = builtin('pi', varargin{:});
end

---

用matlab怎么编程求pi啊 ？视频