（C语言）圆周率 PI算法

最新推荐文章于 2024-03-12 20:52:11 发布

原创最新推荐文章于 2024-03-12 20:52:11 发布 · 1.3w 阅读

20 ·

CC 4.0 BY-SA版权

文章标签：

#PI算法

PI算法研究专栏收录该内容

1 篇文章

订阅专栏

本文提供了一种使用特定级数来计算圆周率π的方法，通过迭代更新算法实现了π值的高精度计算，并展示了部分计算结果。

#include<stdio.h>

#include<stdlib.h>

int main()

{
double x=2, z=2;
int a=1, b=3;
while(z>1e-15)
   {
   z = z*a/b;
   x += z;
   a++;
   b+=2;
   }

printf("Pi=%.13f", x);

return 0;

}

//应用第三个级数进行计算，级数具体内容请参考高等数学，E-mail:chisir2000@163.com

//输出结果：Pi=3.1415926535898

#include <stdio.h>

int main(){
const int ARRSIZE=1010, DISPCNT=1000; //定义数组大小，显示位数
//const ARRSIZE=10100, DISPCNT=10000;
char x[ARRSIZE], z[ARRSIZE]; //x[0] x[1] . x[2] x[3] x[4] .... x[ARRSIZE-1]
int a=1, b=3, c, d, Run=1, Cnt=0,i;
for (i=0;i<ARRSIZE;i++){
x[i]=0;z[i]=0;
}
//memset(x,0,ARRSIZE);
//memset(z,0,ARRSIZE);
x[1] = 2;
z[1] = 2;
while(Run && (++Cnt<200000000)){
//z*=a;
d = 0;
for(i=ARRSIZE-1; i>0; i--)
{
c = z[i]*a + d;
z[i] = c % 10;
d = c / 10;
}
//z/=b;
d = 0;
for(i=0; i<ARRSIZE; i++){
c = z[i]+d*10;
z[i] = c / b;
d = c % b;
}
//x+=z;
Run = 0;
for(i=ARRSIZE-1; i>0; i--){
c = x[i] + z[i];
x[i] = c%10;
x[i-1] += c/10;
Run |= z[i];
}
a++;
b+=2;
//printf("计算了 %d 次\n",Cnt);
}
printf("计算了 %d 次\r\n",Cnt);
printf("Pi=%d%d.\r\n", x[0],x[1]);
for(i=0; i<DISPCNT; i++){
if(i && ((i%100)==0))
printf("\r\n");
printf("%d",(int)x[i+2]);
}
return 0;

}

//应用第三个级数进行计算，级数具体内容请参考高等数学，E-mail:chisir2000@163.com