---

def f(x):
    return 1/(x**3+1)**0.5
a = float(input("a: "))
b = float(input("b: "))
n = int(input("n: "))
h = (b-a)/n
s = 0

x=a+h
while x<=b-h:
    s+=f(x)
    x+=h
    print("x: ",x)

sl = h*(f(a)+s)
sp = h(s+f(b))
print(f"Левый прямоугольник: {sl}")
print(f"Правый прямоугольник: {sp}")