目前，我确实使用 assimulos 求解器套件来求解形式为 0 = F(t， y(t(， y'(t(( 的隐式微分方程。我想使用本机 scipy 安装附带的求解器，并查看了scipy.integrate.solve_ivp的文档。文档说求解器求解了一个形式的系统：
求解常微分方程系统的初值问题。 该函数在数值上积分一个常微分方程组，给定一个初始值
dy / dt = f(t, y)y(t0) = y0
我确实怀念解决形式dy / dt = f(t, y, y')系统的能力。这可以用scipy.integrate.solve_ivp吗？
这样的方程称为微分代数方程(DAE(。据我所知，SciPy 没有内置的 DAE 求解器。
这个类似问题的答案提供了一个使用fsolve结合odeint来解决问题的代数部分的示例。这可能适用于您的问题。
您可能还对 ODES 包感兴趣，该包提供了其他求解器作为 SciPy 的扩展，包括 DAE 求解器。