numerical experiment.
report, orbital evolution of solar systems.
Eccentricity = 0.0 (t=0).
Time evolution of the magnitude of the eccentricity vector.
N=32.
https://youtu.be/Zd7TeECvfuI?si=wWU6xSPNXU_w81rl
https://youtu.be/eDsN9gEhS4w?si=ldrdr02TGHg0cQhq
- Makino J. (2001), 重力多体系の数値計算
- Kepler問題におけるKustaanheimo-Stiefel変換
- 谷川清隆 訳、Prasenjit Saha and Scott Tremaine (1992), Symplectic integrators for solar system dynamics
- 正則化法と全保存型差分法を用いた重力N体問題のシミュレーション
- Makino J. , 6 計算法 --- 時間領域
- 重力多体系の天体物理
- E.ハイラー, S.P.ネルセット, G.ヴァンナー著「常微分方程式の数値解法 I」
- 天文学特別講義IV
- Sverre J. Aarseth, "Gravitational N-Body Simulations: Tools and Algorithms (Cambridge Monographs on Mathematical Physics)"