求PI

esybq

#include <stdio.h>
#include <conio.h>
#include <math.h>
#define JD 1e-6
/*

pi = 4/1 -4/3 + 4/5 - 4/7.......


pi/2 = (2*2)/(1*3) + (4*4)/(3*5) + (6*6)/(5*7) +.........

精度 10^-6  : (1e-6) :
通项的 绝对值 大于 JD

*/

int main()
{
    float pi = 0;
    float i = 1;// 通项序号
    float sign = 1; //符号位
    float term = sign*4/(2*i -1); //通项

    while( fabs(term) > JD )
    {
        pi += term;
        //为下一次循环做准备 ： 求下个通项
        i++;
        sign = -sign;
        term = sign*4/(2*i -1);
    }

    printf("pi = %.6f",pi);
    return 0;
}