近似求PI

#include<stdio.h>
double fact(int n)
{
    int i;
    double num=1;
    for(i=1;i<=n;i++)
    {
        num*=i;
    }
    return num;
}
double fact1(int n)
{
    int i;
    double num=1;
    for(i=1;i<=n;i++)
    {
        num*=2*i+1;
    }
    return num;
}
int main()
{
    double pi=2,eps,a=1;
    int i=0;
    scanf("%le",&eps);
    while(fact(i)/fact1(i)>=eps)
    {
        i++;
        pi+=(fact(i)/fact1(i))*2;
 
    }
	printf("PI = %.5lf",pi);
}