Question

使用程序Matlab，割圆术求圆周率π

function calpi(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
有没有人知道这个m文件是什么意思？ 解释一下也好

Answer 1

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 %显示最终结果

追问

= = 我想知道的是为什么这么计算，是周长还是面积计算 也就是计算的方法是什么？