from math import * def f(x): return x - cos(x) def fp(x): return 1.0 + sin(x) x = 1.0 n = 5 x = float(x) print "x0", " = ", x for i in range( n ): x = x - f(x) / fp(x) print "x" + str(i+1), " = ", x