Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- def f(x):
- return x**5 - x**3 + 3*x - 5
- LIMIT = 50
- EPSILON = 10**-5
- a = 1
- b = 2
- c = (a+b)/2
- qt = 1
- while (qt < LIMIT and f(a)*f(b) < 0 and not(-EPSILON <= f(c) <= EPSILON)):
- c = (a+b)/2
- qt += 1
- if (f(a) * f(c) > 0):
- a = c
- else:
- b = c
- print(c)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement