我想得到一个非线性方程组的数值jacobian矩阵。我在python 3.4上使用模块numdifftools 0.9.16。
我建立了一个有趣的功能＆＃39;采用数组X并返回数组Y.
我使用命令jac=nd.Jacobian(fun)(X)获得jacobian并且工作正常。
我还使用命令jac=nd.Jacobian(fun,step=1e-1)(X)控制步长，也没关系。
现在，我希望根据等式有两个不同的步长：jac=nd.Jacobian(fun,step=[1,1,1e-1,1e-1])(X)，假设有4个方程式。我想在前两个方程中使用1，而在其他方程中使用1e-1。
不幸的是，这个命令不起作用。我无法找到任何强制不同步长的例子，但阅读手册应该是可行的。没人有经验吗？感谢。