虚数が実数から現れる時（惑星の軌道編） - Qiita

前回に引き続き、Pythonを用いて、実数のみを用いた簡単な関係式をシミュレートすると不思議な結果が出てくること及び、虚数を使うとそれがうまく説明できるようなケースを紹介することで、より皆様に虚数に親しみを持ってもらおうと思います。 惑星の運動方程式 突然ですが、万有引力による惑星の運動方程式は以下のようなものになります。 \displaylines{ \vec{F} = m\vec{a} = -\frac{GMm}{r^3}\vec{r} \\ \vec{F}:力 \\ m:惑星の質量 \\ \vec{a}:惑星の加速度 \\ G:万有引力定数 \\ M:恒星の質量 \\ m:惑星の質量 \\ \vec{r}:惑星の位置 \\ r:\vec{r}の長さ} 解析的に解くには2回積分をする必要があり、かなりめんどくさいのですが、そこはPythonに任せて力任せにシミュレートすればなんともない