ルンゲ＝クッタ法
The two ordinary differential equations are solved using Runge–Kutta methods of orders 1, 3, and 5 with the same time step, to show the effects of computational accuracy in the results.