というわけで

惑星の運動の微分方程式を Mathematica で解いて、 軌道が円錐曲線 (典型的な場合として楕円) になることを試しにやってみませんか？