---

def my_poly_integral(f, x_var):

    if f.diff(x_var) == 0:
        return f * x_var
    
    coeffs = f.coefficients(sparse=False)
    
    integral_result = 0
    for power, coeff in enumerate(coeffs):

        new_power = power + 1
        integral_result += (coeff / new_power) * (x_var^new_power)
        
    return integral_result


var('x')
p = 3*x^2 + 2*x + 5
print(f"Интеграл от {p} равен: {my_poly_integral(p, x)}")