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