用c语言实现割圆法计算圆周率
“割圆法”,是用圆内接正多边形的面积去无限逼近圆面积并以此求取圆周率的方法。计算的关键是,从当前的多边形如何计算加倍后的多边形的边长。
1.变量说明
变量i:第几次切割
变量s:多边形的边数
变量k:圆周率近似值
变量h:临时变量
2.程序代码
#include <stdio.h>
#include <math.h>
void pai(){
int i=0,s=6; //i:第几次切割 s:多边形的边数
double k=3.0,len=1.0,h=0;
while(i<=50){
printf("第 %2d次切割,为正 %10d边形,PI=%.15f\n",i,s,k); //k:圆周率祭祀值
h=1-sqrt(1-pow(len/2,2)); //临时变量
len=sqrt(pow(len/2,2)+pow(h,2));
i++;
s*=2;
k=s*len/2;
if((k-3.0)<1e-16){
break;
}
}
}
int main(){
pai(); //调用pai()函数
return 0;
}
3.运行结果
4.分析结论
大概循环到24次,精确度就可以准确到小数点后15位,效率较高。“割圆法”,是用圆内接正多边形的面积去无限逼近圆面积并以此求取圆周率。根据圆的面积公式,当圆的半径为1时,圆的面积等于π,算半径为1 的正多边形的面积,即可得到近似的π。通过观察图像,S圆>S正24边形>S正12边形>S正6边形,计算结果的精确度随多边形的边数的增加而提高,但始终小于S圆,即始终小于π。