import math


def f(x, n):
    return x * x - n


def fp(x):
        return 2 * x


def pierwiastek(n):
    epsilon = 0.001
    xk = n - f(n, n)/fp(n)
    xk1 = xk - f(xk, n) / fp(xk)

    while math.fabs(xk1 - xk) > epsilon:
        xk = xk1
        xk1 = xk - f(xk, n) / fp(xk)

    return xk1


def odleglosc(x, y, z):
    return pierwiastek(x*x + y*y + z*z)


n = 4
punkty = [
    [ 5.1,  4.8,    2.1 ],
    [ -2.1, 2.5,    42.31 ],
    [ -4,   2.1,    -3.33 ],
    [ -4.9, -2.6,   -5.4 ]
]


for punkt in punkty:
    wynik = int(odleglosc(punkt[0], punkt[1], punkt[2]) * 100)
    przyblizenie = wynik/100
    print(przyblizenie)