Advertisement
Guest User

Untitled

a guest
Jun 28th, 2017
51
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.00 KB | None | 0 0
  1. def poly(func,x):
  2. answer=0
  3. for coex in func:
  4. polyterm=coex[0]*x**coex[1]
  5. answer= answer+polyterm
  6. return answer
  7.  
  8. def f_of(x):
  9. return x**3+5*x**2+x-6
  10.  
  11.  
  12.  
  13.  
  14. def root_a(f, left, right):
  15. left = float(left)
  16. right = float(right)
  17.  
  18. if left > right :
  19. tmp = a
  20. left = right
  21. right = tmp
  22.  
  23. mid = (right+left)/2
  24.  
  25. while right-left > 0.000001 :
  26. if poly.poly(f, mid)*poly.poly(f, left) < 0 :
  27. right = mid
  28. else:
  29. left = mid
  30. mid = (right+left)/2
  31. print mid
  32. return mid
  33.  
  34.  
  35. def root_b(left, right):
  36. left = float(left)
  37. right = float(right)
  38.  
  39. if left > right :
  40. tmp = left
  41. left = right
  42. right = tmp
  43.  
  44. mid = (right+left)/2
  45.  
  46. while right-left > 0.000001 :
  47. if f_of(mid)*f_of(left) < 0 :
  48. right = mid
  49. else:
  50. left = mid
  51. mid = (right+left)/2
  52. print mid
  53. return mid
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement