Advertisement
Guest User

Untitled

a guest
Apr 24th, 2017
43
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 0.34 KB | None | 0 0
  1. def f1(x):
  2. return x**3-(2.*x)-5.
  3.  
  4. def derivative_f1(x):
  5. return (3*x**2)-2.
  6.  
  7. def Newton(f, df, x, tol):
  8.  
  9. while True:
  10. x1 = x - (f(x)/df(x)) # it was f(x) - (f(x)/df(x))
  11.  
  12. t = abs(x1-x)
  13.  
  14. if t < tol:
  15. break
  16. x = x1
  17.  
  18. return x
  19.  
  20. initial = 2
  21.  
  22. print(Newton(f1,derivative_f1,initial,0.000001))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement