import numpy as np
import scipy.integrate as si

def f(x):
    return x ** 2
x_vals = np.arange(-20, 21, 1)
y_vals = si.odeint(lambda y, x: f(x), 0, x_vals)
Scipy: эффективно генерировать серию интеграции (интегральная функция)