为什么这个操作不能给我更高的精度?

时间:2013-10-11 23:48:26

标签: c++ floating-point

我正在使用以下公式估算Pi的值:

enter image description here

使用以下C ++代码:

double sub = 0;
int prec = 1000; //How many iterations to use in estimate.

for(int i = 1; i <= prec; i++){
    double frac = 1/((3+(2*(i-1))) * pow(3, i));
    sub += (i == 1) ? 1-frac : (i%2) ? -frac : frac;
}

double pi = sqrt(12)*sub;

cout << "Pi estimated with precision of " << prec << " iterations is " << pi << ".\n";

我的问题是,即使在1000(或100000这个问题)迭代,我得到的最高精度是3.14159。我尝试在计算中的每个数字上使用static_cast<double>(),但仍然得到相同的结果。我在这里做错了什么,或者这是这种方法产生的最大精度?我是C ++的新手,但不是编程。

1 个答案:

答案 0 :(得分:5)

问题是你不打印所有精度。你需要打电话,

std::cout << std::setprecision(10) << ...