Advertisement
Guest User

Untitled

a guest
Sep 25th, 2017
69
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.28 KB | None | 0 0
  1. import math
  2. from pathlib import _Accessor
  3.  
  4. min_precision = 0.000001
  5.  
  6. def my_function(x):
  7. return x*x*x*x+(7/2)*(x*x*x)-(29/2)*(x*x)-(55/2)*x+(75/2)
  8.  
  9. def go_deeper(x, intervalas, depth = 0):
  10. answer = my_function(x)
  11. if abs(answer) < min_precision:
  12. print("Gylis: ", depth)
  13. return x
  14. positive = is_positive(answer)
  15. if (bool(is_positive(my_function(x+(intervalas/2)))) ^ bool(positive)):
  16. return go_deeper(x,intervalas/2, depth+1)
  17. else:
  18. return go_deeper(x+intervalas/2, intervalas/2, depth+1)
  19.  
  20. def is_positive(n):
  21. if n > 0:
  22. return True
  23. else:
  24. return False
  25.  
  26.  
  27. B = (55)/2
  28. k = 4 - 2
  29. R = 1 + pow(B/1,1/k)
  30. print("B =", B, "k =", k, "Rteig =", R)
  31.  
  32. B=(29)/2
  33. k = 4 - 3
  34. Rneig = 1 + pow(B/1, 1/k)
  35. print("B =", B, "k =", k, "Rneig =", Rneig)
  36. print("Imame nuo",-R,",iki",R,". Nes R teig yra mažesnis")
  37. intervalas =0.001
  38. print("intervalas yra:", intervalas)
  39. i = -R
  40. positive = False
  41. while i < R:
  42. answer = my_function(i)
  43. positive = is_positive(answer)
  44. i += intervalas
  45.  
  46.  
  47. # if abs(answer) < 0.1:
  48. # print(answer, i)
  49. if (is_positive(my_function(i)) is not positive):
  50. answer = go_deeper(i-intervalas, intervalas)
  51. print("Ats: ", answer, " Value:", my_function(answer))
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement