Advertisement
Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- def poly(func,x):
- answer=0
- for coex in func:
- polyterm=coex[0]*x**coex[1]
- answer= answer+polyterm
- return answer
- def f_of(x):
- return x**3+5*x**2+x-6
- def root_a(f, left, right):
- left = float(left)
- right = float(right)
- if left > right :
- tmp = a
- left = right
- right = tmp
- mid = (right+left)/2
- while right-left > 0.000001 :
- if poly.poly(f, mid)*poly.poly(f, left) < 0 :
- right = mid
- else:
- left = mid
- mid = (right+left)/2
- print mid
- return mid
- def root_b(left, right):
- left = float(left)
- right = float(right)
- if left > right :
- tmp = left
- left = right
- right = tmp
- mid = (right+left)/2
- while right-left > 0.000001 :
- if f_of(mid)*f_of(left) < 0 :
- right = mid
- else:
- left = mid
- mid = (right+left)/2
- print mid
- return mid
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement