viete.m

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

function y=viete(n)
pi=1;a(1)=sqrt(2)/2;
for i=1:n-1
a(i+1)=sqrt(0.5+(1/2)*a(i));
pi=pi*a(i);
end;
vpa(2/pi)