simplefixedpoint

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 import numpy as np def fx(x): return np.exp(-1*x) - x def gx(x): return np.exp(-1*x) def simple(x0, imax, es): iter = 1 ea = 0 xi = x0 print(" i xi ea") while True: xold = xi xi = gx(xold) iter += 1 ea = abs((xi - xold)/xi)*100 print("%4d %.6f %.2f"%(iter, xi, ea)) if iter >= imax or ea