Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- def NRroot(x,f,dfdx):
- x_val = [0]
- x_diff = 1
- counter = 0
- while True:
- x = x-f(x)/dfdx(x)
- x_val.append(x)
- x_diff = x_val[len(x_val)-1] - x_val[len(x_val)-2]
- if x_diff > 10**(-6):
- return x_val[:-1]
- print([(counter,x,x_diff)])
- counter += 1
- return x_val
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement