回複：小推廣的Runge-Kutta方法正確嗎？

我以前用 Runge-Kutta 方法算過 chaos 的問題，記得不太清楚了，隻能說個大概的印象。

你的天體軌道應該是二維平麵問題， 有兩個 coupled differential equations， 時間是變量參數，the increment （step size h） 應該總是 > 0。

或者，如果你能用 polar coordinate 表達軌道方程， 以角度為參數，那麽角度的 increment 也是大於零的。

沒有時間細看你的程序，不一定說得對，僅供參考。