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.
Under general relativity, the motion of a planet around the Sun —in polar coordinates , with the Sun at the origin — satisfies the initial-value problem
,
,
,
where , , , is the gravitational constant of the universe, is the mass of the planet, is the mass of the Sun, is the constant angular momentum of the planet, is the speed of light, and is the smallest distance of the planet from the Sun during its orbit (i.e., at perihelion).
We now take the perspective of a student who is taking a first-semester course in differential equations. There are two standard techniques for solving a second-order non-homogeneous differential equations with constant coefficients. One of these is the method of variation of parameters. First, we solve the associated homogeneous differential equation
.
The characteristic equation of this differential equation is , which clearly has the two imaginary roots . Therefore, two linearly independent solutions of the associated homogeneous equation are and .
(As an aside, this is one answer to the common question, “What are complex numbers good for?” The answer is naturally above the heads of Algebra II students when they first encounter the mysterious number , but complex numbers provide a way of solving the differential equations that model multiple problems in statics and dynamics.)
According to the method of variation of parameters, the general solution of the original nonhomogeneous differential equation
is
,
where
,
,
and is the Wronskian of and , defined by the determinant
.
Well, that’s a mouthful.
The only good news is that is easy to compute. Since and , we have
from the usual Pythagorean trigonometric identity. Therefore, the denominators in the integrals for and essentially disappear.
Unfortunately, computing and , using
,
is a beast, requiring the creative use of multiple trigonometric identities. We begin with , using the substitution :
,
where we use for the constant of integration instead of the usual . Second,
.
Unfortunately, this is not easily simplified with a substitution, so we have to expand the integrand:
,
using for the constant of integration. Therefore, by variation of parameters, the general solution of the nonhomogeneous differential equation is
,
where is another arbitrary constant.
Next, we use the initial conditions to find the constants and . From the initial condition , we obtain
,
so that
.
Next, we compute and use the initial condition :
.
Substituting these values for and , we finally arrive at the solution
.