Almkvist-Giullera formula for pi: Difference between revisions
m
Minor edit to C++ code
(Added C++ solution) |
m (Minor edit to C++ code) |
||
Line 198:
<< '\n';
rational sum = 0, prev = 0;▼
▲ rational prec(1, pow(big_int(10), 70));
for (int n = 0;; ++n) {
rational term(almkvist_giullera(n), pow(big_int(10), 6 * n + 3));
sum += term;
if (abs(pi - prev) < epsilon)
break;
prev =
}
std::cout << "\nPi to 70 decimal places is:\n"
▲ big_float pi2(1 / sum);
<< std::
}</lang>
|