#include <iostream>
#include <cmath>
using namespace std;
 
int main() {
    double c = 15.0;
    double eps = 0.001;
 
    double x = c;   // przybliżenie początkowe
    double x_next;
    int iter = 0;
 
    do {
        x_next = 0.5 * (x + c / x);
        iter++;
        x = x_next;
    } while (fabs(x * x - c) > eps);
 
    cout << "Pierwiastek (Newton-Raphson): " << x << endl;
    cout << "Liczba iteracji: " << iter << endl;
 
    return 0;
}
 

I2luY2x1ZGUgPGlvc3RyZWFtPgojaW5jbHVkZSA8Y21hdGg+CnVzaW5nIG5hbWVzcGFjZSBzdGQ7CgppbnQgbWFpbigpIHsKICAgIGRvdWJsZSBjID0gMTUuMDsKICAgIGRvdWJsZSBlcHMgPSAwLjAwMTsKCiAgICBkb3VibGUgeCA9IGM7ICAgLy8gcHJ6eWJsacW8ZW5pZSBwb2N6xIV0a293ZQogICAgZG91YmxlIHhfbmV4dDsKICAgIGludCBpdGVyID0gMDsKCiAgICBkbyB7CiAgICAgICAgeF9uZXh0ID0gMC41ICogKHggKyBjIC8geCk7CiAgICAgICAgaXRlcisrOwogICAgICAgIHggPSB4X25leHQ7CiAgICB9IHdoaWxlIChmYWJzKHggKiB4IC0gYykgPiBlcHMpOwoKICAgIGNvdXQgPDwgIlBpZXJ3aWFzdGVrIChOZXd0b24tUmFwaHNvbik6ICIgPDwgeCA8PCBlbmRsOwogICAgY291dCA8PCAiTGljemJhIGl0ZXJhY2ppOiAiIDw8IGl0ZXIgPDwgZW5kbDsKCiAgICByZXR1cm4gMDsKfQo=