如何将Jacobian馈送到scikits.odes.odeint

Question

我正在使用scikits.odes.odeint来解决相对较大的ODE系统（40个变量，100个动力学参数）。我想通过向求解器提供一个jacobian来加快此过程。但是，我不知道如何为此目的最好地定义jacobian并将其提供给求解器。

Answer 1

这是一个非常基本的符号Jacobion，但适用于大多数ODE求解器。

:goose

还有更深入的指南here。