计算圆周率

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

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

输入格式：

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

输出格式：

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

输入样例：

0.01

输出样例：

3.132157

#include<stdio.h>
double my(int x)
{
    if (x == 0)
        return 1;
    return x * my(x - 1);
}
double My(int x)
{
    if (x == 0)
        return 1;
    return (2 * x + 1) * My(x - 1);
}
int main()
{
    double k = 1;
    double n = 0;
    scanf("%lf", &n);
    double pi = 0;
    int i = 0;
    double sum = 1.0;
    while (k>n)
    {
        i++;
        k = my(i) / My(i);
        sum += k;
    }
    printf("%lf", 2*sum);
    return 0;
}