Python bindings to SUNDIALS differential algebraic equation solvers
ODE system solver using dG(q), time-discontinuous Galerkin with Lobatto basis
Extend scipy.integrate with various methods for solve_ivp
