关于余弦函数的约等计算

Question

下面是一个余弦函数的大致计算程序：但是仍然有两个测试点无法通过，不知道是为什么啊？

double fact(int n)
{
    int i;
    double a=1;
    for (i=1;i<=n;i++)
    {
        a*=i;
    }
    return a;
}

double MyCos(double x,double epsilon)
{
    int i;
    int flag=-1;
    double sum=0;
    for (i=2;;i+=2)
    {
        if (pow(x,i)/fact(i)<epsilon)
            return sum+1.0;
        else
            sum+=flag*pow(x,i)/fact(i);
            flag=-flag;
    }
}

Answer 1

这个是算法问题（我以前遇过），当-40>=x>=40时,算法的结果就会不对了
试试在程序中加

    double sum=0; //原程序
    while(x>2*M_PI)  //增加
      x-=2*M_PI;//增加
    for (i=2;;i+=2) //原程序