用Scipy求解微分代数方程

Question

我希望（在Python中）为定义状态变量x约束的函数g:R^n->R^m求解形式为x'(t) = f(x(t),y(t))的微分代数方程。我认为不需要使用fsolve从而手动建立一个数值方法就可以做到这一点

例如，使用sci.integrate.solve_ivp或类似的方法似乎是合理的。你能建议一些具体有效的选择吗