浙大版《C语言程序设计实验与习题指导（第4版）》题目集 实验6-1 近似求PI

题目：

本题要求编写程序，根据下式求π的近似值，直到最后一项小于给定精度eps。

2π​=1+31!​+3×52!​+3×5×73!​+⋯+3×5×⋯×(2×i+1)i!​+⋯

---

问题描述

输入格式：

输入在一行中给出精度eps，可以使用以下语句来读输入：

scanf("%le", &eps);

输出格式：

在一行内，按照以下格式输出π的近似值（保留小数点后5位）：

PI = 近似值

---

样例分析：

输入样例：

1E-5

输出样例：

PI = 3.14158

---

代码展示：

#include<stdio.h>
#include<math.h>
double fact(int n){
    double f=1*1.0;//f可能很大很大，所以不能用int；
    for(int i=1;i<=n;i++){
        f=f*i;
    }
    return f;
}
int main(){
    double eps;
    int  i=1;    
    scanf("%le", &eps);
    double t=1.0;
    double sum=0.0;
    double sum1=1.0;//sum1可能很大很大，所以不能用int；
    while(t>=eps){
        sum=sum+t;
        sum1=sum1*(2*i+1);
        t=fact(i)*1.0/sum1;
        i++;
    }
    sum=sum+t;
    double g=sum*2;
    printf("PI = %.5f\n", g);
    return 0;
}