c – 用Chudnovsky算法计算Pi数

概述所以这里有一些根据Chudnovsky算法计算Pi的C代码： #include <iostream>#include <cmath>#include <iomanip>int fact(int digit){ int result = 1; for (int i=digit; i>1; i--) { result *= i; } r 所以这里有一些根据Chudnovsky算法计算Pi的C代码：

#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

解决方法 在fact()中使用double作为fact()返回类型以及“result”变量. 总结

以上是内存溢出为你收集整理的c – 用Chudnovsky算法计算Pi数全部内容，希望文章能够帮你解决c – 用Chudnovsky算法计算Pi数所遇到的程序开发问题。

如果觉得内存溢出网站内容还不错，欢迎将内存溢出网站推荐给程序员好友。