1

我想使用 numjac 计算我的系统的数值雅可比行列式。

它适用于常规语法:

 [Jac, FAC] = numjac(@(t,y)model(t,y,X,VAR),ttot,ytot',dydt,jac_tol,FAC);

但是我不需要完整的雅可比,而只需要关于某些 y 的导数。我可以计算完整的然后选择所需的行,但这是浪费计算时间。如何仅在需要的地方计算雅可比?

4

0 回答 0