Метод Симпсона и трапеций

import sympy
 
x = sympy.Symbol('x')
 
f_x = x / (x**4 + 4)
 
 
def sympson(left, right, n, function):
    h = (right - left) / (2 * n)
 
    tmp_sum = float(function.subs({x: left})) +\
        float(function.subs({x: right}))
 
    for step in range(1, 2 * n):
        if step % 2 != 0:
            tmp_sum += 4 * float(function.subs({x: left + step * h}))
        else:
            tmp_sum += 2 * float(function.subs({x: left + step * h}))
 
    return tmp_sum * h / 3
 
 
print(sympson(0, 5, 5, f_x))

def f(x):
    return x / (x**4 + 4)
 
 
def sympson(left, right, n, function):
    h = (right - left) / (2 * n)
 
    tmp_sum = function(left) + function(right)
 
    for step in range(1, 2 * n):
        if step % 2 != 0:
            tmp_sum += 4 * function(left + step * h)
        else:
            tmp_sum += 2 * function(left + step * h)
 
    return tmp_sum * h / 3
 
 
def trapezian(left, right, n, function):
    h = (right - left) / (n)
 
    return (function(left) + function(right) +
            sum(function(left + step * h) for step in range(1, n))) * h
 
 
print(sympson(0, 5, 15, f))
print(trapezian(0, 5, 30, f))