Calculation of the definition of fundamental integrals

numeric_lab4/rectangles.py

@@ -0,0 +1,32 @@
+import math
+
+
+def f(x):
+    return x ** 2 / math.sqrt(x ** 2 + 3)
+
+
+def F(x):
+    return x / 2 * math.sqrt(x ** 2 + 3) + 3 / 2 * math.log(x + math.sqrt(x ** 2 + 3))
+
+
+def check(a, b):
+    print(abs(F(b) - F(a)))
+
+
+def compute_right_rectangles(interactive=False):
+    integral = 0
+    a = int(input("Enter lower bound: ")) if interactive else 1
+    b = int(input("Enter upper bound: ")) if interactive else 2
+    n = int(input("Enter number of divisions: ")) if interactive else 100
+    h = (b - a) / n
+    x = a
+    for i in range(1, n):
+        integral += f(x)
+        x += h
+    integral *= h
+    print(integral)
+    check(a, b)
+
+
+if __name__ == '__main__':
+    compute_right_rectangles()