viete.m

algotithm to calculate pi

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)