#include int main() { const int max_it = 13; const int max_it_j = 10; double a1 = 1.0, a2 = 0.0, d1 = 3.2; std::cout << " i d\n"; for (int i = 2; i <= max_it; ++i) { double a = a1 + (a1 - a2) / d1; for (int j = 1; j <= max_it_j; ++j) { double x = 0.0; double y = 0.0; for (int k = 1; k <= 1 << i; ++k) { y = 1.0 - 2.0*y*x; x = a - x * x; } a -= x / y; } double d = (a1 - a2) / (a - a1); printf("%2d %.8f\n", i, d); d1 = d; a2 = a1; a1 = a; } return 0; }