StoneHaos

math9

Dec 25th, 2020
878
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
  1. import math
  2.  
  3. f = lambda x: 1/x
  4.  
  5. a = -2
  6. b = -1
  7. n = 10
  8.  
  9. s = 0
  10. h = (b - a) / n
  11. x = a
  12. for i in range(n + 1):
  13.     print("x = {0}, y = {1}".format(round(x, 4), round(f(x), 4)))
  14.     s += f(x)
  15.     x += h
  16.  
  17. Sn = h * (s - f(b))
  18. Si = h * (s - f(a))
  19. print("Sнед = {0}*({1}) = {2}, Sизб = {0}*({3}) = {4}".format(round(h, 4), round(s - f(b), 4), round(Sn, 4), round(s - f(a), 4), round(Si, 4)))
  20. St = h * (s - f(a) - f(b) + (f(a) + f(b)) / 2)
  21. print("Sтрап = {0}".format(round(St, 4)))
  22.  
  23. s1 = 0
  24. s2 = 0
  25. x = a + h
  26. for i in range(1, n):
  27.     if (i % 2 == 0):
  28.         s2 += f(x)
  29.     else:
  30.         s1 += f(x)
  31.     x += h
  32.  
  33. Ss = (b - a) / (3 * n) * (f(a) + f(b) + 4 * s1 + 2 * s2)
  34. print("{0}({1}+{2}+4({3})+2({4}) = {0}({5}+{6}+{7}) = {0}*{8} = {9}".format(round((b - a) / (3 * n), 4), round(f(a), 4), round(f(b), 4), round(s1, 4), round(s2, 4), round(f(a) + f(b), 4), round(4 * s1, 4), round(2 * s2, 4), round(f(a) + f(b) + 4 * s1 + 2 * s2, 4), round(Ss, 4)))
RAW Paste Data