コンピューター物理実験は誤差と時間との闘い・・・・

ルンゲ・クッタ法で３体問題をある初期条件で解くと、きれいな三ツ輪もようになります。
プロットの都合で偏重してるけど、同じ質量の質点を質量中心周りの角運動量が等しくなる接線速度で、互いの相互作用ハミルトニアンのみが働く系で計算させてみた。（はじめの図）