c – 用Chudnovsky算法计算Pi数

#include <iostream>
#include <cmath>
#include <iomanip>

int fact(int digit)
{
    int result = 1;
    for (int i=digit; i>1; i--)
    {
        result *= i;
    }
    return result;
}

int main()
{
    long double pi = 0.0;
    int precision = 10;
    for ( int k=0; k<=precision; k++)
    {
        pi += ( pow(-1,k)*fact((6*k))*(13591409+545140134*k))/(fact(3*k)*pow(fact(k),3)*pow(640320,(3*k+1.5)));
    }
    pi = 1.0/(12*pi);
    std::cout<<std::setprecision(80)<<pi<<'\n'<<M_PI<<'\n';
}

问题是这段代码只返回了15位正确的PI编号,另一部分是垃圾…当精度变量设置为12或更大时,它返回-nan.任何人都可以帮助我吗？ TIA
编辑：Ubuntu 12.10 x64,gcc 4.7.2