# -*- coding: utf-8 -*- def f(x): return x**2-x-1 def g(x): return 1+1/x epsilon = 0.001 x = 2 x_before = 1 i = 2 while (abs(x - x_before) > epsilon): x_before = x x = g(x) print(i, "{0:.6f}".format(x)) i=i+1