cutcircle.m

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

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);