なおこのシミュレーションは数値積分を計算してるわけなのですが ここで示したのはオイラー法と言って一番単純な方法であり、理解しやすくて計算も速いというメリットがあります。 反面、精度はよくないです。 この精度を上げるためにいろいろ工夫されており、ニュートン法、ルンゲ=クッタ法、シンプレクティク法、ニューマークβ法等があります。 興味があれば調べてみてください。