In this series, I’m discussing how ideas from calculus and precalculus (with a touch of differential equations) can predict the precession in Mercury’s orbit and thus confirm Einstein’s theory of general relativity. The origins of this series came from a class project that I assigned to my Differential Equations students maybe 20 years ago.
In this part of the series, we will show that if the motion of a planet around the Sun is expressed in polar coordinates , with the Sun at the origin, then under Newtonian mechanics (i.e., without general relativity) the motion of the planet follows the differential equation
,
where and
is a certain constant. Deriving this governing differential equation will require some principles from physics. If you’d rather skip the physics and get to the mathematics, we’ll get to solving this differential equations in the next post.
Part of the derivation of this governing differential equation will involve Newton’s Second Law
,
where is the mass of the planet and the force
and the acceleration
are vectors. In usual rectangular coordinates, the acceleration vector would be expressed as
,
where the components of the acceleration in the and
directors are
and
, and the unit vectors
and
are perpendicular, pointing in the positive
and positive
directions.
Unfortunately, our problem involves polar coordinates, and rewriting the acceleration vector in polar coordinates, instead of rectangular coordinates, is going to take some work.
Suppose that the position of the planet is in polar coordinates, so that the position in rectangular coordinates is
. This may be rewritten as
,
where
is a unit vector that points away from the origin. We see that this is a unit vector since
.
We also define
to be a unit vector that is perpendicular to ; it turns out that
points in the direction of increasing
. To see that
and
are perpendicular, we observe
.
Computing the velocity and acceleration vectors in polar coordinates will have a twist that’s not experienced with rectangular coordinates since both and
are functions of
. Indeed, we have
.
Furthermore,
.
These two equations will be needed in the derivation below.
We are now in position to express the velocity and acceleration of the orbiting planet in polar coordinates. Clearly, the position of the planet is , or a distance
from the origin in the direction of
. Therefore, by the Product Rule, the velocity of the planet is
We now apply the Chain Rule to the second term:
.
Differentiating a second time with respect to time, and again using the Chain Rule, we find
.
This will be needed in the next post, when we use both Newton’s Second Law and Newton’s Law of Gravitation, expressed in polar coordinates.


