"Пример scipy.integrate.romberg.\nСреда: встроенный MicroPython FEMM 5.\nУточнение трапеций экстраполяцией Ромберга.\nРезультат: float — оценка интеграла."

import scipy as sp
def f(i):
    # Учебное потокосцепление psi=2I; интеграл по току — коэнергия, Дж.
    return 2*i
r=sp.integrate.romberg(f,0.0,2.0,eps=1e-8)
assert abs(r-4)<1e-6
print(r)