在不使用符号计算软件的情况下,如何通过C++或其他类似的编程语言计算一般非线性多项式系统的Jacobian?as symbolic Jacobian- only dependent on C++ or similar programming languages;
如果我不得不用有限差分法来获得近似雅可比,所选择的步长将如何影响最终解的精度?如何确定步长,以便在相同的计算精度水平下
为此,我需要能够访问ODE求解器中的步长,以计算此时间步长。我使用的是scipy.integrate.odeint --有没有办法在定义odes的函数中访问时间步长?问题是,在第一个时间步长上,如果时间步长是1,那么总体将变得太低,但由于时间步长将很小,最初您可以从x获取所有内容。matplotlib.pyplot as plt
tspan = np.linspace(0.0, tend, 1000)
B