Kepler's First Law of Planetary Motion

[Brunetto]

I'm reviewing my old calculus textbook and I stumbled upon a proof of Kepler's First Law of Planetary Motion which uses vector valued functions along with all of the operations to demonstrate the material. I understand the math and how to to DO it but what I am struggling with is why.

It goes from:

v X h = GM u + c

and then does the dot product of both sides with r.


Why is the above equation "dotted" with r.

This is apparently found in most calculus textbooks but if it is not:
v is the velocity
u is the unit vector of direction
h is a constant vector of rXv
c is a constant vector of integrating u'

 

[D H]

The typical approach is to start by showing that the specific angular momentum \vec h = \vec r \times \vec v and the eccentricity vector \vec e = \vec v \times \vec h / \mu - \hat r are constants of motion, where \mu=GM.

That \vec h is constant is a direct consequence of central motion, any kind of central motion. That the eccentricity vector is constant is only true for inverse square law central motion, \vec a = -\mu/r^2 \hat r.

To show that the eccentricity vector is constant, look at the time derivative of \vec v \times \vec h:
\begin{aligned}<br /> \frac{d}{dt}(\vec v \times \vec h)<br /> &amp;= \frac{d\vec v}{dt}\times h \\<br /> &amp;= -\frac{\mu}{r^2}\hat r \times \vec h \\<br /> &amp;= -\frac{\mu h}{r^2} \hat r \times \hat h \\<br /> &amp;= \mu\dot{\theta}\hat{\theta} \\<br /> &amp;= \mu\frac{d\hat r}{dt}<br /> \end{aligned}
The penultimate step relies upon h=r^2\dot{\theta} and \hat h \times \hat r = \hat{\theta}. The final step relies upon the definition of the circular unit vectors \hat r and \hat{\theta}.

From the above, d/dt(\vec v \times \vec h) = \mu d/dt(\hat r), or d/dt(\vec v \times \vec h - \mu\hat r) = 0. Thus \mu\vec e = \vec v \times \vec h - \mu\hat r is a constant of motion.

Now we're at the point where we can take the inner product with the position vector.
\begin{aligned}<br /> \mu\vec e \cdot \vec r<br /> &amp;= (\vec v \times \vec h)\cdot \vec r - \mu \hat r \cdot \vec r \\<br /> &amp;= (\vec r \times \vec v)\cdot \vec h - \mu r \\<br /> &amp;= h^2 - \mu r<br /> \end{aligned}

Another way to express the above inner product is in terms of the magnitudes of the two vectors and the angle f between them: \mu\vec e \cdot \vec r = \mu e r \cos f. Equating these approaches yields
\mu e r \cos f = h^2 - \mu r
from which
r = \frac{h^2/\mu}{1+e\cos f}
Defining p=h^2/\mu yields r=p/(1+e\cos f), which is the equation of a conic section with semi latus rectum p, eccentricity e, and the origin at one of the foci.