实验6-1 近似求PI (15分)

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

​2
​
​π
​​ =1+
​3
​
​1!
​​ +
​3×5
​
​2!
​​ +
​3×5×7
​
​3!
​​ +⋯+
​3×5×⋯×(2×i+1)
​
​i!
​​ +⋯

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

scanf("%le", &eps);
输出格式：
在一行内，按照以下格式输出π的近似值（保留小数点后5位）：

PI = 近似值
输入样例：
1E-5
输出样例：
PI = 3.14158

#include<stdio.h>
#include<ctype.h>
#include<string.h>
#include<math.h>
#include<stdlib.h>
double fact(double n){
    double result = 1, i;
    if ( n == 0 || n == 1)
        return 1;
    for ( i = 2; i <= n; i++)
        result *= i;
    return result;
}
double ffact(double n){
    double result = 1, i;
    for ( i = 1; i <= n; i++)
        result *= (2*i + 1);
    return result;
}
int main()
{	
	double sum = 1, eps, n = 1, item = 1;		//要设置初始时的最后一项为1，不然后面过不去第三个测试点。
	scanf("%le", &eps);
    
    while(item > eps){
        item = fact(n)/ffact(n);
        sum += item;
        n++;
    }
    printf("PI = %.5f", 2*sum);
	return 0;
}