import sys

n = float(sys.argv[1])
EPSILON = sys.float_info.epsilon

t = n
while abs(t - n/t) > (EPSILON*t):
    t = (n / t + t) / 2

print(f"sqrt({n}) ~= {t}")