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

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 this post, we will use the guesses

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

that arose from the technique/trick of reduction of order, where $f(\theta)$ is some unknown function, to find the general solution of the differential equation

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

To do this, we will need to use the Product Rule for higher-order derivatives that was derived in the previous post:

$(fg)'' = f'' g + 2 f' g' + f g''$

and

$(fg)^{(4)} = f^{(4)} g + 4 f''' g' + 6 f'' g'' + 4f' g''' + f g^{(4)}$ .

In these formulas, Pascal’s triangle makes a somewhat surprising appearance; indeed, this pattern can be proven with mathematical induction.

We begin with $u(\theta) = f(\theta) \cos \theta$ . If $g(\theta) = \cos \theta$ , then

$g'(\theta) = - \sin \theta$ ,

$g''(\theta) = -\cos \theta$ ,

$g'''(\theta) = \sin \theta$ ,

$g^{(4)}(\theta) = \cos \theta$ .

Substituting into the fourth-order differential equation, we find the differential equation becomes

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

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

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

$f^{(4)} \cos \theta - 4 f''' \sin \theta - 4 f'' \cos \theta = 0$

The important observation is that the terms containing $f$ and $f'$ cancelled each other. This new differential equation doesn’t look like much of an improvement over the original fourth-order differential equation, but we can make a key observation: if $f'' = 0$ , then differentiating twice more trivially yields $f''' = 0$ and $f^{(4)} = 0$ . Said another way: if $f'' = 0$ , then $u(\theta) = f(\theta) \cos \theta$ will be a solution of the original differential equation.

Integrating twice, we can find $f$ :

$f''(\theta) = 0$

$f'(\theta) = c_1$

$f(\theta) = c_1 \theta + c_2$ .

Therefore, a solution of the original differential equation will be

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

We now repeat the logic for $u(\theta) = f(\theta) \sin \theta$ :

$(f \sin \theta)^{(4)} + 2 (f \sin \theta)'' + f \sin \theta = 0$

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

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

$f^{(4)} \sin\theta - 4 f''' \cos\theta - 4 f'' \sin\theta = 0$ .

Once again, a solution of this new differential equation will be $f(\theta) = c_3 \theta + c_4$ , so that $f'' = f''' = f^{(4)} = 0$ . Therefore, another solution of the original differential equation will be

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

Adding these provides the general solution of the differential equation:

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

Except for the order of the constants, this matches the solution that was presented earlier by using techniques taught in a proper course in differential equations.

One thought on “Confirming Einstein’s Theory of General Relativity With Calculus, Part 6g: Rationale for Method of Undetermined Coefficients IV”

Pingback: Confirming Einstein’s General Theory of Relativity with Calculus: Index – Mean Green Math

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

Published by John Quintanilla

One thought on “Confirming Einstein’s Theory of General Relativity With Calculus, Part 6g: Rationale for Method of Undetermined Coefficients IV”

Leave a comment Cancel reply

Share this:

Related

Published by John Quintanilla

One thought on “Confirming Einstein’s Theory of General Relativity With Calculus, Part 6g: Rationale for Method of Undetermined Coefficients IV”

Leave a comment Cancel reply