x, y = 3/2, 5/3 for i in range(20): z = y y = 2003 - 6002/y + 4000/(x*y) x = z print(y)