你自己再考虑一下 Runge-cotta 方法的应用范围

它是为了解决已有解析式微分方程的非解析解问题，(iteration method), 所以 step size 取负值是没有意义的。

当然，如果你原方程中的函数在变量到一定值时在所涵盖的范围内发生了变化就是另一回事。

所有跟帖： 

• 假如ODE的解是椭圆，解上半圆h大于0，解下半圆h小于0， -9$- ♂ (116 bytes) () 12/13/2012 postreply 09:31:01

• 这个用 polar coordinate （r，theta）的牛顿方程 -孩子长了翅膀- ♀ (29 bytes) () 12/13/2012 postreply 10:10:24

• 用polar coordinate解也是上半圆h大于0，解下半圆h小于0，一样滴，嘿嘿 -参谋总长- ♂ (0 bytes) () 12/13/2012 postreply 15:03:30

• 打转转不会想不通的，就是x有点想不通，走到头了回不来 ：） -孩子长了翅膀- ♀ (0 bytes) () 12/13/2012 postreply 15:09:22

• 我支持9美刀，坚决不打转转，就是政府好玩，就是要玩政府 ：） -参谋总长- ♂ (0 bytes) () 12/13/2012 postreply 15:14:31

• 不是玩政府吧，是玩新 Kepler？ -孩子长了翅膀- ♀ (0 bytes) () 12/13/2012 postreply 15:21:39

• 是玩正负吧，跟玩政府一样，我看没什么不同，嘿嘿 -参谋总长- ♂ (0 bytes) () 12/13/2012 postreply 15:58:55