Scipy.integrate.solve_ivp Unclear How To Solve Implicit ODE Of Form 0=F(t, Y(t), Y'(t))

Currently, I do use assimulos solver suite to solve an implicit differential equation of the form 0 = F(t, y(t), y'(t)). I would like to use a solver which comes with a native scip

Solution 1:

Such an equation is called a differential algebraic equation (DAE). As far as I know, SciPy doesn't have a built-in DAE solver.

This answer to this similar question provides an example of using fsolve to solve the algebraic part of the problem in combination with odeint. This might be applicable to your problem.

You might also be interested in the ODES package, which provides additional solvers as an extension for SciPy, including a DAE solver.