DOP853(ルンゲクッタ8次)のExcel VBA移植版でピタゴラスの三体問題を計算

久々にDOP853(Dormand-Prince法)で計算シリーズ。

ピタゴラスの三体問題は、質量mが3,4,5の三体が辺の長さが3,4,5の直角三角形の頂点におかれている場合にどのような動きをするか、というもの。(Burrauさんが問題提起して、Szebehelyさんが数値的に解決したということです）。

詳細はこちら。

http://www.ucolick.org/~laugh/oxide/projects/burrau.html

で計算結果はこちら。

リンク先の図を再現してる。が、これを描くにはDOP853の誤差指定を相当小さくしないとだめだった。そりゃこんな動きをするならそうだろうが、最初気付かず全然再現しなかった。。。だめだめ。