#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;
}