mygauss.m

来自「algotithm to calculate pi」· M 代码 · 共 15 行

function y=mygauss(n)
a=1;
b=1/sqrt(2);
t=1/4;
x=1;
for k=1:n
    l=a;
    a=(a+b)/2;
    b=sqrt(b.*l);
    t=t-x*(l-a).^2;
    x=2*x;
end
temp=(a+b).^2/(4*t);
vpa(temp)