Babylonian Method of Calculating Square Roots Using Python

babylonian.py

number = abs(float(raw_input("Calculate square root of? ")))
guess = abs(float(raw_input("Initial guess? ")))

epsilon = 0.001

while True:
    difference = guess**2 - number
    print('{} : {}'.format(round(guess, 4), round(difference, 4)))
    if abs(difference) <= epsilon:
        break
    guess = (guess + number / guess) / 2.0