cutcircle.m

algotithm to calculate pi

function y=cutcircle(n)
syms a;
for i=1:n
    a=sqrt(2-sqrt(4-a^2));
end
a=subs(a,'a','1');
y=3*2^n*vpa(a,n+5);