Confirming Einstein’s Theory of General Relativity With Calculus, Part 6f: Rationale for Method of Undetermined Coefficients III

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.

We have shown that the motion of a planet around the Sun, expressed in polar coordinates $(r,\theta)$ with the Sun at the origin, under general relativity follows the initial-value problem

$u''(\theta) + u(\theta) = \displaystyle \frac{1}{\alpha} + \frac{\delta}{\alpha^2} + \frac{\delta \epsilon^2}{2\alpha^2} + \frac{2\delta \epsilon \cos \theta}{\alpha^2} + \frac{\delta \epsilon^2 \cos 2\theta}{2\alpha^2}$ ,

$u(0) = \displaystyle \frac{1}{P}$ ,

$u'(0) = 0$ ,

where $u = \displaystyle \frac{1}{r}$ , $\displaystyle \frac{1}{\alpha} = \frac{GMm^2}{\ell^2}$ , $\delta = \displaystyle \frac{3GM}{c^2}$ , $G$ is the gravitational constant of the universe, $m$ is the mass of the planet, $M$ is the mass of the Sun, $\ell$ is the constant angular momentum of the planet, $c$ is the speed of light, and $P$ is the smallest distance of the planet from the Sun during its orbit (i.e., at perihelion).

In the previous post, I used a standard technique from differential equations to find the general solution of

$u^{(4)} + 2u'' + u = 0$ .

to be

$u(theta) = c_1 \cos \theta + c_2 \sin \theta + c_3 \theta \cos \theta + c_4 \theta \sin \theta$ .

However, as much as possible in this series, I want to take the perspective of a talented calculus student who has not yet taken differential equations — so that the conclusion above is far from obvious. How could this be reasonable coaxed out of such a student?

To begin, we observe that the characteristic equation is

$r^4 + 2r^2 + 1 = 0$ ,

$(r^2 + 1)^2 = 0$ .

Clearly this has the same roots as the simpler equation $r^2 + 1 = 0$ , which corresponds to the second-order differential equation $u'' + u = 0$ . We’ve already seen that $u_1(\theta) = \cos \theta$ and $u_2(\theta) = \sin \theta$ are solutions of this differential equation; perhaps they might also be solutions of the more complicated differential equation also? The answer, of course, is yes:

$u_1^{(4)} + 2 u_1'' + u_1 = \cos \theta - 2 \cos \theta + \cos \theta = 0$

and

$u_2^{(4)} + 2u_2'' + u_2 = \sin \theta - 2 \sin \theta + \sin \theta = 0$ .

The far trickier part is finding the two additional solutions. To find these, we use a standard trick/technique called reduction of order. In this technique, we guess that any additional solutions much have the form of either

$u(\theta) = f(\theta) \cos \theta \qquad \hbox{or} \qquad u(\theta) = f(\theta) \sin \theta$ ,

where $f(\theta)$ is some unknown function that we’re multiplying by the solutions we already have. We then substitute this into the differential equation $u^{(4)} + 2u'' + u = 0$ to form a new differential equation for the unknown $f$ , which we can (hopefully) solve.