7-15 计算圆周率 (15 分)

根据下面关系式，求圆周率的值，直到最后一项的值小于给定阈值。

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

输入格式：

输入在一行中给出小于1的阈值。

输出格式：

在一行中输出满足阈值条件的近似圆周率，输出到小数点后6位。

输入样例：

0.01

输出样例：

3.132157

代码如下：

#include <stdio.h>
double fac(int n);
double fun(int n);
int main()
{
    double min, pi = 0;
    scanf("%lf", &min);
    int i = 0;
    for (i = 0;; i++)
    {
        pi += fac(i) / fun(i);
        if (fac(i) / fun(i) < min)
            break;
    }
    printf("%.6lf", 2.0 * pi);

    return 0;
}
double fac(int n)
{
    double s = 1;
    int i;
    if (n == 0)
        return 1;
    for (i = 1; i <= n; i++)
        s *= i;
    return s;
}
double fun(int n)
{
    double s = 1;
    int i;
    for (i = 3; i <= 2 * n + 1; i += 2)
        s *= i;
    return s;
}