数值分析--matlab刘徽割圆术

作业要求：

利⽤刘徽割圆术思想，编写计算圆周率的程序，并与pi⽐较精度。

---

代码：

n = 6;
sin = sqrt(3) / 2;
i = 1;
while (i <= 10)
	cos = sqrt(1 - sin * sin);
	sin = sqrt(1 / 2 * (1 - cos));
	n = n * 2;
	i = i + 1;
end
vpa(pi, 10)
vpa(sin * n / 2, 10)

运行结果