Вариант 4 / 4 / simpson

import math

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

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

simp(1, 2)