#include <iostream>
#include <cmath>
 
int main() {
    double c = 15.0;      // liczba, z której liczymy pierwiastek
    double eps = 0.001;   // dokładność
    double x = c;         // początkowe przybliżenie
    double x_next;
 
    do {
        x_next = 0.5 * (x + c / x);
        x = x_next;
    } while (fabs(x * x - c) > eps);
 
    std::cout << "Przyblizona wartosc pierwiastka kwadratowego z "
              << c << " wynosi: " << x << std::endl;
 
    return 0;
}
 

I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+CgppbnQgbWFpbigpIHsKICAgIGRvdWJsZSBjID0gMTUuMDsgICAgICAvLyBsaWN6YmEsIHoga3TDs3JlaiBsaWN6eW15IHBpZXJ3aWFzdGVrCiAgICBkb3VibGUgZXBzID0gMC4wMDE7ICAgLy8gZG9rxYJhZG5vxZvEhwogICAgZG91YmxlIHggPSBjOyAgICAgICAgIC8vIHBvY3rEhXRrb3dlIHByenlibGnFvGVuaWUKICAgIGRvdWJsZSB4X25leHQ7CgogICAgZG8gewogICAgICAgIHhfbmV4dCA9IDAuNSAqICh4ICsgYyAvIHgpOwogICAgICAgIHggPSB4X25leHQ7CiAgICB9IHdoaWxlIChmYWJzKHggKiB4IC0gYykgPiBlcHMpOwoKICAgIHN0ZDo6Y291dCA8PCAiUHJ6eWJsaXpvbmEgd2FydG9zYyBwaWVyd2lhc3RrYSBrd2FkcmF0b3dlZ28geiAiCiAgICAgICAgICAgICAgPDwgYyA8PCAiIHd5bm9zaTogIiA8PCB4IDw8IHN0ZDo6ZW5kbDsKCiAgICByZXR1cm4gMDsKfQo=