 Solution

Write $k=GM$ and $r=|\mathbf{r}|$. The gravitational central force has zero torque, so the specific angular momentumis constant. For a nonradial orbit $h=|\mathbf{h}|>0$, motion lies in the plane normal to $\mathbf{h}$, and plane polar coordinates give $r^2\dot{\theta }=h$. With $u=1/r$ and primes denoting $\theta$ derivatives,The radial equation $\ddot{r}-r{\dot{\theta }}^2=-k/r^2$ is therefore the Binet equationIts general solution, after choosing the polar origin of angle, isThe origin is a focus of this conic section. To identify which conic section occurs, use the conserved specific energyA bound, nonradial orbit has $\epsilon <0$, hence $0\le e<1$ and is an ellipse. For example, putting $x=r\cos \theta$, $y=r\sin \theta$ in $r+ex=p$ givesan explicit ellipse equation. This is the intended low, bound Kepler orbit. The inverse-square equation alone does not imply an ellipse: an initially tangential speed larger than $\sqrt{2GM/r}$ gives $\epsilon >0$, $e>1$ and a hyperbolic Kepler orbit. Zero angular momentum instead gives radial motion. Those qualifications are necessary for the elliptic conclusion.

In a circular orbit, centripetal acceleration balances the gravitational acceleration:Here $\Omega$ is the positive angular speed; either sense of circulation has this period.

With the stated drag force, use the full angular momentum $\mathbf{L}=m\mathbf{r}\times \dot{\mathbf{r}}$. Gravitational torque is still zero, while the drag torque isTherefore $\boxed{\mathbf{L}(t)=\mathbf{L}(0)e^{-At/m}}$ for constant $A>0$. Its direction remains fixed, and its magnitude has exponential decay with time constant $m/A$; the orbit remains planar although its shape is no longer a fixed ellipse.