#include <iostream>
#include <cmath>
using namespace std;
int main() {
double c = 15.0;
double eps = 0.001;
double x = c; // początkowe przybliżenie
double x_prev;
do {
x_prev = x;
x = 0.5 * (x_prev + c / x_prev);
} while (fabs(x - x_prev) > eps);
cout << "Pierwiastek kwadratowy z " << c << " wynosi: " << x << endl;
return 0;
}
I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+Cgp1c2luZyBuYW1lc3BhY2Ugc3RkOwoKaW50IG1haW4oKSB7CiAgICBkb3VibGUgYyA9IDE1LjA7CiAgICBkb3VibGUgZXBzID0gMC4wMDE7CiAgICBkb3VibGUgeCA9IGM7ICAgICAgICAvLyBwb2N6xIV0a293ZSBwcnp5YmxpxbxlbmllCiAgICBkb3VibGUgeF9wcmV2OwoKICAgIGRvIHsKICAgICAgICB4X3ByZXYgPSB4OwogICAgICAgIHggPSAwLjUgKiAoeF9wcmV2ICsgYyAvIHhfcHJldik7CiAgICB9IHdoaWxlIChmYWJzKHggLSB4X3ByZXYpID4gZXBzKTsKCiAgICBjb3V0IDw8ICJQaWVyd2lhc3RlayBrd2FkcmF0b3d5IHogIiA8PCBjIDw8ICIgd3lub3NpOiAiIDw8IHggPDwgZW5kbDsKCiAgICByZXR1cm4gMDsKfQo=