solving 7 coupled odes with different delays

[atefehneu]

Hi Zulko,

firstly, thank you for writing ddeint!. I want to solve 7 coupled differential equations with different delays (e.g.,
z1 = (Jd1d1 * Y1(t-d1) + Jd1d2 * Y2(t-d2) + Jd1fsi * Y3(t-d3) + Jd1TA * Y4(t-d4) + Jd1TI * Y5(t-d5) + Jd1ctx * (2. * np.sin(q(t-d6)) + 2.5))
z2 = (Jd2d1 * Y1(t-d21) + Jd2d2 * Y2(t-d22) + Jd2fsi * Y3 + Jd2TA * Y4 + Jd2TI * Y5(t-d25) + Jd2ctx * (2. * np.sin(q(t-d26)) + 2.5))
z3 = (Jfsifsi * Y3(t-d33) + JfsiTA * Y4(t-d34) + JfsiTI * Y5 + Jfsictx * (2. * np.sin(q(t-d36)) + 2.5))
, ...
)
and i'm afraid if ddeint can solve equations like what I want.

Thank you so much,
Atefeh

[Menginventor]

I implemented dde solver base on ddeint with scipy.integrate.solve_ivp which is newer version.
https://github.com/Menginventor/dde_ivp
It support rk23, rk45 adaptive step size. vectorized variable, dense output, and event detection
I would be more efficient and accurate the existing ddeint. Feel free to try it out.