我一直在寻找一种方法来找到微分方程改变其符号的点,它可以是零交叉或求根。
有人告诉我使用Brent's Method
,但这并不准确。例如,图形在0.12
处切零,但是该方法获得了0.19
之类的值或更多,这是不可接受的。
现在,我发现了该算法scipy.integrate.solve_ivp,在这种方法中,它表示为The solver looks for a sign change over each step, so if multiple zero crossings occur within one step, events may be missed.
我不知道该怎么做。
考虑我们有一个微分方程
dE/dz = f(E,z)
我们知道dE/dz
的答案,我们知道起始点E0
,但是我们想知道在哪个z
中答案为dE/dz
。因此我们应该获取z
的值,从而得出dE/dz
请帮忙吗?
谢谢