import math def babylon(root,guess): """The Babylonian algorithm""" oldGuess = guess*1.0 while True: newGuess = (oldGuess+root/oldGuess)/2 print oldGuess, newGuess if abs(oldGuess - newGuess) < 10 ** -10: break else: oldGuess = newGuess babylon(2,1)