The circular motion

The two graph above I made is to show the unstability of EulerForward method during long period of computing.
上面是一个模拟万有引力控制的有心运动，用向前Euler格式造成的误差使能量不守恒，变成了向心运动。
图2用了Runge-Kutta方法，其精度有所提高，但仍不能从根本上避免能量减少。
其实最好的办法是采用辛几何，保证Hamilton量守恒，不过暂时不会，呵呵。