数值分析--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)

运行结果

在这里插入图片描述


版权声明:本文为weixin_45850271原创文章,遵循CC 4.0 BY-SA版权协议,转载请附上原文出处链接和本声明。