odeint для матриц

Question

Я хотел бы закодировать в python связанную систему дифференциальных уравнений, такую ​​как:

$dF_{ij}/dt=A_{ij}$ where $A_{ij}$ является функцией матрицы $F$.

Значение системы Дифференциальные уравнения 1-го порядка.

Для той же системы, но с векторами: $dF_{i}/dt=A_{i}$ Я использую scipy.integrate.odeint.

Но я пытался, и это не работает для матриц.

tmin, tmax, tstep = (0., 200., 1)
t_test=np.arange(tmin, tmax, tstep) #time vector

dydt_testm=np.array([[0.,1.],[2.,3.]])
Y0_test=np.array([[0,1],[0,1]])

def dydt_test(y,t):
    return dydt_testm

result = si.odeint(dydt_test, #Computes the derivative of y at t (voir ligne 82)
                                     Y0_test,               #Initial condition on y (can be a vector).
                                     t_test)

ValueError: Начальное условие y0 должно быть одномерным.

Как мне это исправить?

Answer 1

Как прокомментировал Уоррен Векессер в комментариях, odeintw делает свою работу.

from odeintw import odeintw
import numpy as np

Y0_test=np.array([[0,1],[0,1]])
tmin, tmax, tstep = (0., 200., 1)
t_test=np.arange(tmin, tmax, tstep) #time vector

dydt_testm=np.array([[0.,1.],[2.,3.]])

def dydt_test(y,t):
    return dydt_testm

result = odeintw(dydt_test, #Computes the derivative of y at t 
                                     Y0_test,               #Initial condition on y (can be a vector).
                                     t_test)
plt.plot(t_test,result[:,0,1])
plt.show()