Bonjour,

J'ai besoin de résoudre un couple d'équation différentielles:

dxdt = a - b*x*z
dydt = c*y - y*exp(-E)
dzdt = dxdt - sum(dydt*dE)

avec a, b et c des constantes et E, un tableau.

Comme vous le voyez, je dois procéder à une intégration sur les E pour déterminer dzdt.
Cela implique que dydt soit un tableau également.

J'ai d'abord essayé d'utiliser odeint mais on ne peut pas initialiser le système en mélangeant scalaires et tableau:
x=odeint(fun, initialValues, t) ---> initialValues doit être du type [0.0 0.0 0.0] alors qu'il me faudrait quelque chose comme [0.0 zeros(len(E)) 0.0].

Pourriez-vous m'indiquer la manière de résoudre mon système d'équations?

Merci beaucoup par avance