Integrating DAEs/ implicit ODE models discovered by SINDy

[lucafusarbassini]

Hi! I'm running SINDy on a 2D system (x0, x1 variables) with a SINDy-PI library because from biological theory we expect functional nonlinearities. anyway, starting from the Michaelis-Menten tutorial I got to a number of models discovered by SINDy, for example:

x0x0_t = 0.195 1 + 0.247 x0 + -0.045 x0x1 + -0.001 x1x1 + 0.035 x1x0_t + 0.013 x0x1x0_t + -0.136 x0x1x1_t + -0.218 x1x1x0_t

now, however, the tutorial uses odeint to integrate the discovered system because it's a trivial 1D system, which is not my case here. I've tried integrating it by diffeqpy and fipy with very poor results. something coarse with time discretization would be more than enough for my purposes.

thank you very much,
Luca