标签： c++ numeric solver odeint
我想用一个隐式求解器解决一个（可能是刚性的）数值问题。我实现了ODE的右侧，并且能够使用odeint的显式求解器。 但是我对如何向隐式求解器提供雅可比矩阵感到困惑。由于ODE太复杂，我也无法计算。 不过：有没有一种方法可以根据给定的ODE计算Jacobian？还是有任何隐式求解器库在没有提供Jacobian的情况下工作？
谢谢！