def arrel_quadrada(n, eps):
    x1 = n
    x2 = 0.5 * (x1+ n/x1)
    while abs(x2-x1) >= eps:
        x1 = x2
        x2 = 0.5 * (x1+ n/x1)
    return x2