回复：小推广的Runge-Kutta方法正确吗？

我以前用 Runge-Kutta 方法算过 chaos 的问题，记得不太清楚了，只能说个大概的印象。

你的天体轨道应该是二维平面问题， 有两个 coupled differential equations， 时间是变量参数，the increment （step size h） 应该总是 > 0。

或者，如果你能用 polar coordinate 表达轨道方程， 以角度为参数，那么角度的 increment 也是大于零的。

没有时间细看你的程序，不一定说得对，仅供参考。