#include <cstdio>

int main()
{
    float x = 3;
    float y = 0.25;
    std::printf("x = %g, y0 = %g ~ 1/%g (0 < y0 < 2/%g)\n", x, y, x, x);
    for (int i = 0; 5 != i; ++i) {
        std::printf("y%d = %g dy = %g (dy = y%d - x * (y%d)^2)\n", i, y, y - x * y * y, i, i);
        y = y * (2.f - x * y);
    }
    return 0;
}