function calpi(n) %定义函数，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);%求最终结果Pi=y %显示最终结果
