（C语言）圆周率 PI算法

#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