数値積分にはルンゲ・クッタ法を使用しています。 「おかし！力学的エネルギーが保存されねえじゃねえか！」という人は変数Hの値を小さくすると幸せになれるかも。