Вариант 17 / 4 / 4t

import math

def f(x):
    return 5 * (math.sin(x)) ** 3 + (math.cos(x)) ** 3

def trap(a, b):
    e = 1e-4
    n = 1
    h = (b - a)
    s = (f(a) + f(b)) * h / 2
    s1 = 1e-15
    while abs(s - s1) / 3 >= e:
        n = 2 * n
        h = (b - a)/ n
        s1 = s
        s = f(a) + f(b)
        for i in range(1, n - 1):
            s = s + 2*f(a + i*h)
        s = s*h/2
    print(f'Приближённое значение интеграла: {s}')

trap(1, 2)