Confirming Einstein’s Theory of General Relativity With Calculus, Part 6k: Solving New Differential Equation with Method of Undetermined Coefficients

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 recent posts, we used the method of undetermined coefficients to show that the general solution of the differential equation is

$u(\theta) = c_1 \cos \theta + c_2 \sin \theta + \displaystyle \frac{1}{\alpha} + \frac{\delta}{\alpha^2} + \frac{\delta \epsilon^2}{2\alpha^2} +\frac{\delta \epsilon}{\alpha^2} \theta \sin \theta- \frac{\delta \epsilon^2}{6\alpha^2} \cos 2\theta$ .

We now use the initial conditions to find the constants $c_1$ and $c_2$ . (We did this earlier when we solved the differential equation via variation of parameters, but we repeat the argument here for completeness.) From the initial condition $u(0) = \displaystyle \frac{1}{P} = \frac{1+\epsilon}{\alpha}$ , we obtain

$u(0) = \displaystyle c_1 \cos 0 + c_2 \sin 0 + \displaystyle \frac{1}{\alpha} + \frac{\delta}{\alpha^2} + \frac{\delta \epsilon^2}{2\alpha^2} + \frac{\delta \epsilon \cdot 0 \cdot \sin 0}{\alpha^2} -\frac{\delta \epsilon^2 \cos 0}{6\alpha^2}$

$\displaystyle \frac{1+\epsilon}{\alpha} = c_1 + \displaystyle \frac{1}{\alpha} + \frac{\delta}{\alpha^2} + \frac{\delta \epsilon^2}{2\alpha^2} -\frac{\delta \epsilon^2}{6\alpha^2}$

$\displaystyle \frac{\epsilon}{\alpha} = c_1 + \displaystyle \frac{3\delta +\delta \epsilon^2}{3\alpha^2}$ ,

so that

$c_1 = \displaystyle \frac{\epsilon}{\alpha} - \frac{\delta(3 + \epsilon^2)}{3\alpha^2}$ .

Next, we compute $u'(\theta)$ and use the initial condition $u'(0) = 0$ :

$u'(\theta) = \displaystyle -c_1 \sin \theta + c_2 \cos \theta + \frac{\delta \epsilon}{\alpha^2} (\sin \theta + \theta \cos \theta) + \frac{\delta \epsilon^2 \sin 2\theta}{3\alpha^2}$

$u'(0) = \displaystyle -c_1 \sin 0 + c_2 \cos 0 + \frac{\delta \epsilon}{\alpha^2} (\sin 0 + 0 \cos 0) + \frac{\delta \epsilon^2 \sin 0}{3\alpha^2}$

$0 = c_2$ .

Substituting these values for $c_1$ and $c_2$ , we finally arrive at the solution

$u(\theta) = \displaystyle \left(\frac{\epsilon}{\alpha} - \frac{\delta(3 + \epsilon^2)}{3\alpha^2} \right) \cos \theta + \displaystyle \frac{1}{\alpha} + \frac{\delta}{\alpha^2} + \frac{\delta \epsilon^2}{2\alpha^2} + \frac{\delta \epsilon \theta \sin \theta}{\alpha^2} -\frac{\delta \epsilon^2 \cos 2\theta}{6\alpha^2}$

$= \displaystyle \frac{1 + \epsilon \cos \theta}{\alpha} + \frac{\delta}{\alpha^2} + \frac{\delta \epsilon^2}{2\alpha^2} + \frac{ \delta\epsilon}{\alpha^2} \theta \sin \theta - \frac{ \delta \epsilon^2}{6\alpha^2} \cos 2\theta - \frac{\delta(3+\epsilon^2)}{3\alpha^2} \cos \theta$ .

One thought on “Confirming Einstein’s Theory of General Relativity With Calculus, Part 6k: Solving New Differential Equation with Method of Undetermined Coefficients”

Confirming Einstein’s Theory of General Relativity With Calculus, Part 6k: Solving New Differential Equation with Method of Undetermined Coefficients

Published by John Quintanilla

One thought on “Confirming Einstein’s Theory of General Relativity With Calculus, Part 6k: Solving New Differential Equation with Method of Undetermined Coefficients”

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 6k: Solving New Differential Equation with Method of Undetermined Coefficients”

Leave a comment Cancel reply