simulation

最近、アジカンのワールド ワールド ワールドをエンドレスで聞いてます。
特に「惑星」が良いです。

で、モデルですが。
いつまでもMathematicaのNDSolveで解いているわけにはいかないので、C++で書いてみました。
オイラー法って、テイラー展開の第一項と第二項で近似しているみたいなもんでしょうか？
より精度を求めるなら、ルンゲクッタ法らしい。
けど、そんなに重くないので、時間の刻み幅を小さくしてなんとかなりそうです。

今日の一枚は、我が家の前に舞い降りたコクワガタのメス？