Confirming Einstein’s Theory of General Relativity With Calculus, Part 6c: Solving New Differential Equation with Variation of Parameters

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 (r,\theta), with the Sun at the origin — satisfies the initial-value problem

u''(\theta) + u(\theta) = \displaystyle \frac{1}{\alpha} + \delta \left( \frac{1 + \epsilon \cos \theta}{\alpha} \right)^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).

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

u''(\theta) + u(\theta) = 0.

The characteristic equation of this differential equation is r^2 + 1 = 0, which clearly has the two imaginary roots r = \pm i. Therefore, two linearly independent solutions of the associated homogeneous equation are u_1(\theta) = \cos \theta and u_2(\theta) = \sin \theta.

(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 i, 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

u''(\theta) + u(\theta) = g(\theta)

is

u(\theta) = f_1(\theta) u_1(\theta) + f_2(\theta) u_2(\theta),

where

f_1(\theta) = -\displaystyle \int \frac{u_2(\theta) g(\theta)}{W(\theta)} d\theta ,

f_2(\theta) = \displaystyle \int \frac{u_1(\theta) g(\theta)}{W(\theta)} d\theta ,

and W(\theta) is the Wronskian of u_1(\theta) and u_2(\theta), defined by the determinant

W(\theta) = \displaystyle \begin{vmatrix} u_1(\theta) & u_2(\theta) \\ u_1'(\theta) & u_2'(\theta) \end{vmatrix}  = u_1(\theta) u_2'(\theta) - u_1'(\theta) u_2(\theta).

Well, that’s a mouthful.

The only good news is that W is easy to compute. Since u_1(\theta) = \cos \theta and u_2(\theta) = \sin \theta, we have

W(\theta) = (\cos \theta)(\cos \theta) - (\sin \theta)(-\sin \theta) = \cos^2 \theta + \sin^2 \theta = 1

from the usual Pythagorean trigonometric identity. Therefore, the denominators in the integrals for f_1(\theta) and f_2(\theta) essentially disappear.

Unfortunately, computing f_1(\theta) and f_2(\theta), using

g(\theta) = \displaystyle \frac{1}{\alpha} + \delta \left( \frac{1 + \epsilon \cos \theta}{\alpha} \right)^2,

is a beast, requiring the creative use of multiple trigonometric identities. We begin with f_1(\theta), using the substitution t = \cos \theta:

f_1(\theta) = -\displaystyle \int \left[ \frac{1}{\alpha} + \delta \left( \frac{1 + \epsilon \cos \theta}{\alpha} \right)^2 \right] \sin \theta \, d\theta

= \displaystyle \int \left[ \frac{1}{\alpha} + \delta \left( \frac{1 + \epsilon t}{\alpha} \right)^2 \right] \, dt

= \displaystyle \int \left[ \frac{1}{\alpha} + \delta \left( \frac{1 + \epsilon t}{\alpha} \right)^2 \right] \, dt

= \displaystyle \frac{t}{\alpha} +  \frac{\alpha\delta}{3\epsilon} \left( \frac{1 + \epsilon t}{\alpha} \right)^3 + a

= \displaystyle \frac{\cos \theta}{\alpha} +  \frac{\alpha\delta}{3\epsilon} \left( \frac{1 + \epsilon \cos \theta}{\alpha} \right)^3 + a,

where we use +a for the constant of integration instead of the usual +C. Second,

f_2(\theta) = \displaystyle \int \left[ \frac{1}{\alpha} + \delta \left( \frac{1 + \epsilon \cos \theta}{\alpha} \right)^2 \right] \cos\theta \, d\theta.

Unfortunately, this is not easily simplified with a substitution, so we have to expand the integrand:

f_2(\theta) = \displaystyle \int \left[ \frac{\cos \theta}{\alpha} + \frac{\delta \cos \theta}{\alpha^2} + \frac{2 \delta \epsilon \cos^2 \theta}{\alpha^2} + \frac{\delta \epsilon^2 \cos^3 \theta}{\alpha^2} \right] \, d\theta

= \displaystyle \int \left[ \frac{\cos \theta}{\alpha} + \frac{\delta \cos \theta}{\alpha^2} + \frac{\delta \epsilon (1 + \cos 2 \theta)}{\alpha^2} + \frac{\delta \epsilon^2 \cos \theta \cos^2 \theta}{\alpha^2} \right] \, d\theta

= \displaystyle \int \left[ \frac{\cos \theta}{\alpha} + \frac{\delta \cos \theta}{\alpha^2} + \frac{\delta \epsilon}{\alpha^2} + \frac{\delta \epsilon \cos 2 \theta}{\alpha^2}+ \frac{\delta \epsilon^2 \cos \theta (1- \sin^2 \theta)}{\alpha^2} \right] \, d\theta

= \displaystyle \int \left[ \frac{\cos \theta}{\alpha} + \frac{\delta (1 + \epsilon^2) \cos \theta}{\alpha^2}  + \frac{\delta \epsilon}{\alpha^2} + \frac{\delta \epsilon \cos 2 \theta}{\alpha^2} - \frac{\delta \epsilon^2 \cos \theta \sin^2 \theta}{\alpha^2} \right] \, d\theta

= \displaystyle \frac{\sin\theta}{\alpha} + \frac{\delta (1 + \epsilon^2) \sin\theta}{\alpha^2}  + \frac{\delta \epsilon \theta}{\alpha^2} + \frac{\delta \epsilon \sin 2 \theta}{2\alpha^2} - \frac{\delta \epsilon^2 \sin^3 \theta}{3\alpha^2} + b,

using +b for the constant of integration. Therefore, by variation of parameters, the general solution of the nonhomogeneous differential equation is

u(\theta) = f_1(\theta) u_1(\theta) + f_2(\theta) u_2(\theta)

= \displaystyle \left( \frac{\cos \theta}{\alpha} +  \frac{\delta}{3\alpha^2\epsilon} \left( 1 + \epsilon \cos \theta \right)^3 + a \right) \cos \theta

+ \displaystyle \left(\frac{\sin\theta}{\alpha} + \frac{\delta (1 + \epsilon^2) \sin\theta}{\alpha^2}  + \frac{\delta \epsilon \theta}{\alpha^2} + \frac{\delta \epsilon \sin 2 \theta}{2\alpha^2} - \frac{\delta \epsilon^2 \sin^3 \theta}{3\alpha^2} + b  \right) \sin \theta

= a\cos \theta + b \sin \theta + \displaystyle \frac{\cos^2 \theta + \sin^2 \theta}{\alpha} +  \frac{\delta \cos \theta}{3\alpha^2\epsilon} +  \frac{\delta \cos^2 \theta}{\alpha^2} +  \frac{\delta \epsilon \cos^3 \theta}{\alpha^2}+  \frac{\delta\epsilon^2 \cos^4 \theta}{3\alpha^2}

+ \displaystyle \frac{\delta \sin^2 \theta}{\alpha^2} + \frac{\delta \epsilon^2 \sin^2\theta}{\alpha^2}  + \frac{\delta \epsilon \theta \sin \theta}{\alpha^2} + \frac{\delta \epsilon \sin \theta\sin 2 \theta}{2\alpha^2} - \frac{\delta \epsilon^2 \sin^4 \theta}{3\alpha^2}

= a\cos \theta + \displaystyle \frac{\delta \cos \theta}{3\alpha^2\epsilon} + b \sin \theta + \displaystyle \frac{1}{\alpha} +  \frac{\delta (\cos^2 \theta+\sin^2 \theta)}{\alpha^2} +  \frac{\delta \epsilon \cos^3 \theta}{\alpha^2}+  \frac{\delta\epsilon^2 \cos^4 \theta}{3\alpha^2}

\displaystyle + \frac{\delta \epsilon^2 \sin^2\theta}{\alpha^2}  + \frac{\delta \epsilon \theta \sin \theta}{\alpha^2} + \frac{\delta \epsilon \sin \theta\sin 2 \theta}{2\alpha^2} - \frac{\delta \epsilon^2 \sin^4 \theta}{3\alpha^2}

= a\cos \theta+ \displaystyle \frac{\delta \cos \theta}{3\alpha^2\epsilon}  + b \sin \theta + \displaystyle \frac{1}{\alpha} +  \frac{\delta}{\alpha^2} + \frac{\delta \epsilon \theta \sin \theta}{\alpha^2} +  \frac{\delta \epsilon \cos^3 \theta}{\alpha^2} + \frac{\delta \epsilon \sin \theta\sin 2 \theta}{2\alpha^2}

\displaystyle + \frac{\delta \epsilon^2 \sin^2\theta}{\alpha^2}  + \frac{\delta \epsilon^2 (\cos^4 \theta -\sin^4 \theta)}{3\alpha^2}

= a\cos \theta + \displaystyle \frac{\delta \cos \theta}{3\alpha^2\epsilon} + b \sin \theta + \displaystyle \frac{1}{\alpha} +  \frac{\delta}{\alpha^2} + \frac{\delta \epsilon \theta \sin \theta}{\alpha^2}  +  \frac{\delta \epsilon \cos^3 \theta}{\alpha^2} + \frac{\delta \epsilon \sin \theta \cdot 2 \sin \theta \cos \theta}{2\alpha^2}

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

= a\cos \theta+ \displaystyle \frac{\delta \cos \theta}{3\alpha^2\epsilon} + b \sin \theta + \displaystyle \frac{1}{\alpha} +  \frac{\delta}{\alpha^2} + \frac{\delta \epsilon \theta \sin \theta}{\alpha^2}  +  \frac{\delta \epsilon \cos^3 \theta}{\alpha^2} + \frac{\delta \epsilon \sin^2 \theta\cos 2 \theta}{\alpha^2}

\displaystyle + \frac{\delta \epsilon^2}{2\alpha^2} -\frac{\delta \epsilon^2 \cos 2\theta}{2\alpha^2}  + \frac{\delta \epsilon^2 \cos 2 \theta}{3\alpha^2}

= a\cos \theta + \displaystyle \frac{\delta \cos \theta}{3\alpha^2\epsilon}+ b \sin \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{\delta \epsilon \cos^3 \theta}{\alpha^2} + \frac{\delta \epsilon \sin^2 \theta\cos \theta}{\alpha^2}

= a\cos \theta + \displaystyle \frac{\delta \cos \theta}{3\alpha^2\epsilon}+ b \sin \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{\delta \epsilon \cos \theta (\cos^2 \theta + \sin^2 \theta)}{\alpha^2}

= a\cos \theta + \displaystyle \frac{\delta \cos \theta}{3\alpha^2\epsilon}+ b \sin \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} +  \frac{\delta \epsilon \cos \theta}{\alpha^2}

= \displaystyle \left(a +  \displaystyle \frac{\delta}{3\alpha^2\epsilon} +  \frac{\delta \epsilon}{\alpha^2} \right) \cos \theta + b \sin \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 A \cos \theta + b \sin \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} ,

where A is another arbitrary constant.

Next, we use the initial conditions to find the constants A and b. From the initial condition u(0) = \displaystyle \frac{1}{P} = \frac{1+\epsilon}{\alpha}, we obtain

u(0) = \displaystyle A \cos 0 + b \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} = A + \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} = A + \displaystyle \frac{3\delta +\delta \epsilon^2}{3\alpha^2},

so that

A = \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 -A \sin \theta + b \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 -A \sin 0 + b \cos 0 + \frac{\delta \epsilon}{\alpha^2} (\sin 0 + 0  \cos 0) + \frac{\delta \epsilon^2 \sin 0}{3\alpha^2}

0 = b.

Substituting these values for A and b, 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.

Confirming Einstein’s Theory of General Relativity With Calculus, Part 5d: Deriving Orbits under Newtonian Mechanics Using Variation of Parameters

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 previously showed that if the motion of a planet around the Sun is expressed in polar coordinates (r,theta), with the Sun at the origin, then under Newtonian mechanics (i.e., without general relativity) the motion of the planet follows the differential equation

u''(\theta) + u(\theta) = \displaystyle \frac{1}{\alpha},

where u = 1/r and \alpha is a certain constant. We will also impose the initial condition that the planet is at perihelion (i.e., is closest to the sun), at a distance of P, when \theta = 0. This means that u obtains its maximum value of 1/P when \theta = 0. This leads to the two initial conditions

u(0) = \displaystyle \frac{1}{P} \qquad \hbox{and} \qquad u'(0) = 0;

the second equation arises since u has a local extremum at \theta = 0.

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

u''(\theta) + u(\theta) = 0.

The characteristic equation of this differential equation is r^2 + 1 = 0, which clearly has the two imaginary roots r = \pm i. Therefore, two linearly independent solutions of the associated homogeneous equation are u_1(\theta) = \cos \theta and u_2(\theta) = \sin \theta.

(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 i, 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

u''(\theta) + u(\theta) = g(\theta)

is

u(\theta) = f_1(\theta) u_1(\theta) + f_2(\theta) u_2(\theta),

where

f_1(\theta) = -\displaystyle \int \frac{u_2(\theta) g(\theta)}{W(\theta)} d\theta ,

f_2(\theta) = \displaystyle \int \frac{u_1(\theta) g(\theta)}{W(\theta)} d\theta ,

and W(\theta) is the Wronskian of u_1(\theta) and u_2(\theta), defined by the determinant

W(\theta) = \displaystyle \begin{vmatrix} u_1(\theta) & u_2(\theta) \\ u_1'(\theta) & u_2'(\theta) \end{vmatrix}  = u_1(\theta) u_2'(\theta) - u_1'(\theta) u_2(\theta).

Well, that’s a mouthful.

Fortunately, for the example at hand, these computations are pretty easy. First, since u_1(\theta) = \cos \theta and u_2(\theta) = \sin \theta, we have

W(\theta) = (\cos \theta)(\cos \theta) - (\sin \theta)(-\sin \theta) = \cos^2 \theta + \sin^2 \theta = 1

from the usual Pythagorean trigonometric identity. Therefore, the denominators in the integrals for f_1(\theta) and f_2(\theta) essentially disappear.

Since g(\theta) = \displaystyle \frac{1}{\alpha}, the integrals for f_1(\theta) and f_2(\theta) are straightforward to compute:

f_1(\theta) = -\displaystyle \int u_2(\theta) \frac{1}{\alpha} d\theta = -\displaystyle \frac{1}{\alpha} \int \sin \theta \, d\theta = \displaystyle \frac{1}{\alpha}\cos \theta + a,

where we use +a for the constant of integration instead of the usual +C. Second,

f_2(\theta) = \displaystyle \int u_1(\theta)  \frac{1}{\alpha} d\theta = \displaystyle \frac{1}{\alpha} \int \cos \theta \, d\theta = \displaystyle \frac{1}{\alpha}\sin \theta + b,

using +b for the constant of integration. Therefore, by variation of parameters, the general solution of the nonhomogeneous differential equation is

u(\theta) = f_1(\theta) u_1(\theta) + f_2(\theta) u_2(\theta)

= \left( \displaystyle \frac{1}{\alpha}\cos \theta + a \right) \cos \theta + \left( \displaystyle \frac{1}{\alpha}\sin\theta + b \right) \sin \theta

= a \cos \theta + b\sin \theta + \displaystyle \frac{\cos^2 \theta + \sin^2 \theta}{\alpha}

= a \cos \theta + b \sin \theta + \displaystyle \frac{1}{\alpha}.

Unsurprisingly, this matches the answer in the previous post that was found by the method of undetermined coefficients.

For the sake of completeness, I repeat the argument used in the previous two posts to determine a and b. This is require using the initial conditions u(0) = \displaystyle \frac{1}{P} and u'(0) = 0. From the first initial condition,

u(0) = a \cos 0 + b \sin 0 + \displaystyle \frac{1}{\alpha}

\displaystyle \frac{1}{P} = a + \frac{1}{\alpha}

\displaystyle \frac{1}{P} - \frac{1}{\alpha} = a

\displaystyle \frac{\alpha - P}{\alpha P} = a

From the second initial condition,

u'(\theta) = -a \sin \theta + b \cos \theta

u'(0) = -a \sin 0 + b \cos 0

0 = b.

From these two constants, we obtain

u(\theta) = \displaystyle \frac{\alpha - P}{\alpha P}  \cos \theta + 0 \sin \theta + \displaystyle \frac{1}{\alpha}

= \displaystyle \frac{1}{\alpha} \left(  1 + \frac{\alpha-P}{P} \cos \theta \right)

= \displaystyle \frac{1 + \epsilon \cos \theta}{\alpha},

where \epsilon = \displaystyle \frac{\alpha - P}{P}.

Finally, since r = 1/u, we see that the planet’s orbit satisfies

r = \displaystyle \frac{\alpha}{1 + \epsilon \cos \theta},

so that, as shown earlier in this series, the orbit is an ellipse with eccentricity \epsilon.

Engaging students: Deriving the double angle formulas for sine, cosine, and tangent

In my capstone class for future secondary math teachers, I ask my students to come up with ideas for engaging their students with different topics in the secondary mathematics curriculum. In other words, the point of the assignment was not to devise a full-blown lesson plan on this topic. Instead, I asked my students to think about three different ways of getting their students interested in the topic in the first place. I plan to share some of the best of these ideas on this blog (after asking my students’ permission, of course). This student submission comes from my former student Morgan Mayfield. His topic, from Precalculus: deriving the double angle formulas for sine, cosine, and tangent. green line How could you as a teacher create an activity or project that involves your topic? I want to provide some variety for opportunities to make this an engaging opportunity for Precalculus students and some Calculus students. Here are my three thoughts: IDEA 1: For precalculus students in a regular or advanced class, have them derive this formula in groups. After students are familiar with the Pythagorean identities and with angle sum identities, group students and ask them to derive a formula for double angles Sin(2θ), Cos(2θ), Tan(2θ). Let them struggle a bit, and if needed give them some hints such as useful formulas and ways to represent multiplication so that it looks like other operations. From here, encourage students to simplify when they can and challenge students to find the other formulas of Cos(2θ). Ask students to speculate instances when each formula for Cos(2θ) would be advantageous. This gives students confidence in their own abilities and show how math is interconnected and not just a bunch of trivial formulas. Lastly, to challenge students, have them come up with an alternative way to prove Tan(2θ), notably Sin(2θ)/Cos(2θ). This would make an appropriate activity for students while having them continue practicing proving trigonometric identities. IDEA 2: This next idea should be implemented for an advanced Precal class, and only when there is some time to spare. Euler was an intelligent man and left us with the Euler’s Formula: e^{ix}=\cos x + i \sin x. Have Precalculus students suspend their questions about where it comes from and what it is used for. This is not something they would use in their class. Reassure them that for what they will do, all they need to understand is imaginary numbers, multiplying imaginary numbers, and laws of exponents. Have them plug in x = A + B and simplify the right-hand side of the equation so that we get: \cos(A+B)+i\sin(A+B)= a + bi where a and b are two real numbers. The goal here is to get \cos(A+B)+i\sin(A+B)= \cos \theta \cos \theta - \sin \theta \sin \theta + (\sin \theta \cos \theta + \cos \theta \sin \theta)i. All the steps to get to this point is Algebra, nothing out of their grasp. Now, the next part is to really get their brains going about what meaning we can make of this. If they are struggling, have them think about the implications of two imaginary numbers being equal; the coefficient of the real parts and imaginary parts must be equal to each other. Lastly, ask them if these equations seem familiar, where are they from, and what are they called…the angle sum formulas. From here, this can lead into what if x=2A? Students will either brute force the formula again, and others will realize x = A + A and plug it in to the equation they just derived and simplify. This idea is a 2-in-1 steal for the angle sum formulas and double angle formulas. It’s biggest downside is this is for Sin(2θ) and Cos(2θ).   IDEA 3: Take IDEA 2, and put it in a Calculus 2 class. Everything that the precalculus class remains, but now have the paired students prove the Euler’s Formula using Taylor Series. Guide them through using the Taylor Series to figure out a Taylor Series representation of e^x, sin x,  and cos x. Then ask students to find an expanded Taylor Series of to 12 terms with ellipses, no need to evaluate each term, just the precise term. Give hints such as i^2= -1 and to consider i^3=i^2 \cdot i = -i and other similar cases. Lastly, ask students to separate the extended series in a way that mimics a + bi using ellipses to shows the series goes to infinity. What they should find is something like this:
Look familiar? Well it is the addition of two Taylor Series that represent Sin(x) and Cos(x). This is the last connection students need to make. Give hints to look through their notes to see why the “a” and “b” in the imaginary number look so familiar. This, is just one way to prove Euler’s Formula, then you can continue with IDEA 2 until your students prove the angle sum formulas and double angle formulas.green line How does this topic extend what your students should have learned in previous courses? Students in Texas will typically be exposed to the Pythagorean Theorem in 8th grade. At this stage, students use a^2+b^2=c^2 to find a missing side length. Students may also be exposed to Pythagorean triples at this stage. Then at the Geometry level or in a Trigonometry section, students will be exposed to the Pythagorean Identity. The Identity is \sin^2 \theta + \cos^2 \theta = 1.  I think that this is not fair for students to just learn this identity without connecting it to the Pythagorean Theorem. I think it would be a nice challenge student to solve for this identity by using a right triangle with hypotenuse c so that Sin (θ) = b/c and cos (θ) = a/c, one could then show either c^2 \sin^2 \theta + c^2 \cos^2 \theta = c^2 and thus c^2(\sin^2 \theta + \cos^2 \theta) = c^2 or one could show (a/c)^2 + (b/c)^2 = (c/c)^2 = 1 (using the Pythagorean theorem). From here, students learn about the angle addition and subtraction formulas in Precalculus. This is all that they need to derive the double angle formulas.

\sin(\alpha + \beta) = \sin \alpha \cos \beta + \cos \alpha \sin \beta

\sin(\alpha - \beta) = \sin \alpha \cos \beta - \cos \alpha \sin \beta

\cos(\alpha + \beta) = \cos \alpha \cos \beta - \sin \alpha \sin \beta

\cos(\alpha - \beta) = \cos \alpha \cos \beta + \sin \alpha \sin \beta

\tan(\alpha + \beta) = \displaystyle \frac{\tan \alpha + \tan \beta}{1 - \tan \alpha \tan \beta}

\tan(\alpha - \beta) = \displaystyle \frac{\tan \alpha - \tan \beta}{1 + \tan \alpha \tan \beta}

This would be a good challenge exercise for students to do in pairs. Sin(2θ) = Sin(θ + θ), Cos(2 θ) = Cos(θ + θ), Tan(2θ) = Tan(θ + θ). Now we can apply the angle sum formula where both angles are equal: Sin(2θ) = sin(θ)cos(θ) + cos(θ)sin(θ) = 2sin(θ)cos(θ) Cos(2θ) = cos(θ)cos(θ) – sin(θ)sin(θ) =  (We use a Pythagorean Identity here) Tan(2θ) = \displaystyle \frac{\tan \theta + \tan \theta}{1 - \tan^2 \theta} = \frac{2 \tan \theta}{1-\tan^2 \theta} Bonus challenge, use Sin(2θ) and Cos(2θ) to get Tan(2θ). Well, if \tan \theta = \displaystyle \frac{\sin \theta}{\cos \theta}, then

\tan 2\theta = \displaystyle \frac{\sin 2\theta}{\cos 2\theta}

= \displaystyle \frac{2 \sin \theta \cos \theta}{\cos^2 \theta - \sin^2 \theta}

= \displaystyle \frac{ \frac{2 \sin \theta \cos \theta}{\cos^2 \theta} }{ \frac{\cos^2 \theta - \sin^2 \theta}{\cos^2 \theta} }

= \displaystyle \frac{2 \tan \theta}{1 - \tan^2 \theta}

The derivations are straight forward, and I believe that many students get off the hook by not being exposed to deriving many trigonometric identities and taking them as facts. This is in the grasp of an average 10th to 12th grader. green line What are the contributions of various cultures to this topic? I have included four links that talk about the history of Trigonometry. It seemed that ancient societies would need to know about the Pythagorean Identities and the angles sum formulas to know the double angle formulas. Here is our problem, it’s hard to know who “did it first?” and when “did they know it?”. Mathematical proofs and history were not kept as neatly written record but as oral traditions, entertainment, hobbies, and professions. The truth is that from my reading, many cultures understood the double angle formula to some extent independently of each other, even if there was no formal proof or record of it. Looking back at my answer to B2, it seems that the double angle formula is almost like a corollary to knowing the angle sum formulas, and thus to understand one could imply knowledge of the other. Perhaps, it was just not deemed important to put the double angle formula into a category of its own. Many of the people who figured out these identities were doing it because they were astronomers, navigators, or carpenters (construction). Triangles and circles are very important to these professions. Knowledge of the angle sum formula was known in Ancient China, Ancient India, Egypt, Greece (originally in the form of broken chords theorem by Archimedes), and the wider “Medieval Islamic World”. Do note that that Egypt, Greece, and the Medieval Islamic World were heavily intertwined as being on the east side of the Mediterranean and being important centers of knowledge (i.e. Library of Alexandria.) Here is the thing, their knowledge was not always demonstrated in the same way as we know it today. Some cultures did have functions similar to the modern trigonometric functions today, and an Indian mathematician, Mādhava of Sangamagrāma, figured out the Taylor Series approximations of those functions in the 1400’s. Greece and China for example relayed heavily on displaying knowledge of trigonometry in ideas of the length of lines (rods) as manifestations of variables and numbers. Ancient peoples didn’t have calculators, and they may have defined trigonometric functions in a way that would be correct such as the “law of sines” or a “Taylor series”, but still relied on physical “sine tables” to find a numerical representation of sine to n numbers after the decimal point. How we think of Geometry and Trigonometry today may have come from Descartes’ invention of the Cartesian plane as a convenient way to bridge Algebra and Geometry. References: https://www.mathpages.com/home/kmath205/kmath205.htm https://en.wikipedia.org/wiki/History_of_trigonometry https://www.ima.umn.edu/press-room/mumford-and-pythagoras-theorem

My Favorite One-Liners: Part 120

I used these shirts as props when teaching Precalculus this week, and they worked like a charm.

After deriving the three Pythagorean identities from trigonometry, I told my class that I got these hand-made his-and-hers T-shirts for my wife’s birthday a couple of years ago. If you can’t see from the picture, one says \sin^2 \theta and the other \cos^2 \theta.

After holding up the shirts, I then asked the class what mathematical message was being communicated.

After a few seconds, someone ventured a guess: “We add up to 1?”

I answered, “That’s right. Together, we’re one.”

Whereupon the class spontaneously reacted with a loud “Awwwwwwwwww.”

I was exceedingly happy.

Engaging students: Deriving the double angle formulas for sine, cosine, and tangent

In my capstone class for future secondary math teachers, I ask my students to come up with ideas for engaging their students with different topics in the secondary mathematics curriculum. In other words, the point of the assignment was not to devise a full-blown lesson plan on this topic. Instead, I asked my students to think about three different ways of getting their students interested in the topic in the first place.

I plan to share some of the best of these ideas on this blog (after asking my students’ permission, of course).

This student submission comes from my former student Daniel Adkins. His topic, from Precalculus: deriving the double angle formulas for sine, cosine, and tangent.

green line

How does this topic extend what your students should have already learned?

A major factor that simplifies deriving the double angle formulas is recalling the trigonometric identities that help students “skip steps.” This is true especially for the Sum formulas, so a brief review of these formulas in any fashion would help students possibly derive the equations on their own in some cases. Listed below are the formulas that can lead directly to the double angle formulas.

A list of the formulas that students can benefit from recalling:

  • Sum Formulas:
    • sin(a+b) = sin(a)cos(b) + cos(a)sin(b)
    • cos(a+b) = cos(a)cos(b) – sin(a)sin(b)
    • tan(a+b) = [tan(a) +tan(b)] / [1-tan(a)tan(b)]

 

  • Pythagorean Identity:
    • Sin2 (a) + Cos2(a) = 1

 

This leads to the next topic, an activity for students to attempt the equation on their own.

 

 

green line

How could you as a teacher create an activity or project that involves your topic?

I’m a firm believer that the more often a student can learn something of their own accord, the better off they are. Providing the skeletal structure of the proofs for the double angle formulas of sine, cosine, and tangent might be enough to help students reach the formulas themselves. The major benefit of this is that, even though these are simple proofs, they have a lot of variance on how they may be presented to students and how “hands on” the activity can be.

I have an example worksheet demonstrating this with the first two double angle formulas attached below. This is in extremely hands on format that can be given to students with the formulas needed in the top right corner and the general position where these should be inserted. If needed the instructor could take this a step further and have the different Pythagorean Identities already listed out (I.e. Cos2(a) = 1 – Sin2(a), Sin2(a) = 1 – Cos2(a)) to emphasize that different formats could be needed. This is an extreme that wouldn’t take students any time to reach the conclusions desired. Of course a lot of this information could be dropped to increase the effort needed to reach the conclusion.

A major benefit with this also is that even though they’re simple, students will still feel extremely rewarded from succeeding on this paper on their own, and thus would be more intrinsically motivated towards learning trig identities.

 

 

green line

How can Technology be used to effectively engage students with this topic?

When it comes to technology in the classroom, I tend to lean more on the careful side. I know me as a person/instructor, and I know I can get carried away and make a mess of things because there was so much excitement over a new toy to play with. I also know that the technology can often detract from the actual math itself, but when it comes to trigonometry, and basically any form of geometric mathematics, it’s absolutely necessary to have a visual aid, and this is where technology excels.

The Wolfram Company has provided hundreds of widgets for this exact purpose, and below, you’ll find one attached that demonstrates that sin(2a) appears to be equal to its identity 2cos(a)sin(a). This is clearly not a rigorous proof, but it will help students visualize how these formulas interact with each other and how they may be similar. The fact that it isn’t rigorous may even convince students to try to debunk it. If you can make a student just irritated enough that they spend a few minutes trying to find a way to show you that you’re wrong, then you’ve done your job in that you’ve convinced them to try mathematics for a purpose.

After all, at the end of the day, it doesn’t matter how you begin your classroom, or how you engage your students, what matters is that they are engaged, and are willing to learn.

Wolfram does have a free cdf reader for its demonstrations on this website: http://demonstrations.wolfram.com/AVisualProofOfTheDoubleAngleFormulaForSine/

 

References

Different Ways of Solving a Contest Problem: Index

I’m doing something that I should have done a long time ago: collecting a series of posts into one single post. The following links comprised my series on different ways of solving the contest problem “If 3 \sin \theta = \cos \theta, what is \sin \theta \cos \theta?”

Part 1: Drawing the angle \theta

Part 2: A first attempt using a Pythagorean identity.

Part 3: A second attempt using a Pythagorean identity and the original hypothesis for \theta.

 

 

Different ways of computing a limit (Part 3)

One of my colleagues placed the following problem on an exam for his Calculus II course…

\displaystyle \lim_{x \to \infty} \frac{\sqrt{x^2+1}}{x}

and was impressed by the variety of correct responses that he received. I thought it would be fun to discuss some of the different ways that this limit can be computed.

Method #3. A trigonometric identity. When we see \sqrt{x^2+1} inside of an integral, one kneejerk reaction is to try the trigonometric substitution x = \tan \theta. So let’s use this here. Also, since x \to \infty, we can change the limit to be \theta \to \pi/2:

\displaystyle \lim_{x \to \infty} \frac{\sqrt{x^2+1}}{x} = \displaystyle \lim_{\theta \to \pi/2} \frac{\sqrt{\tan^2 \theta+1}}{\tan \theta}

= \displaystyle \lim_{\theta \to \pi/2} \frac{\sqrt{\sec^2 \theta}}{\tan \theta}

= \displaystyle \lim_{\theta \to \pi/2} \frac{ \sec \theta}{\tan \theta}

= \displaystyle \lim_{\theta \to \pi/2} \frac{ ~~\displaystyle \frac{1}{\cos \theta} ~~}{ ~~ \displaystyle \frac{\sin \theta}{\cos \theta} ~~ }

= \displaystyle \lim_{\theta \to \pi/2} \frac{ 1}{\sin \theta}

= 1.

Different ways of solving a contest problem (Part 3)

The following problem appeared on the American High School Mathematics Examination (now called the AMC 12) in 1988:

If 3 \sin \theta = \cos \theta, what is \sin \theta \cos \theta?

When I presented this problem to a group of students, I was pleasantly surprised by the amount of creativity shown when solving this problem.

Yesterday, I presented a solution using a Pythagorean identity, but I was unable to be certain if the final answer was a positive or negative without drawing a picture. Here’s a third solution that also use a Pythagorean trig identity but avoids this difficulty. Again, I begin by squaring both sides.

9 \sin^2 \theta = \cos^2 \theta

9 (1 - \cos^2 \theta) = \cos^2 \theta

9 - 9 \cos^2 \theta = \cos^2 \theta

9 = 10 \cos^2 \theta

\displaystyle \frac{9}{10} = \cos^2 \theta

\displaystyle \pm \frac{3}{\sqrt{10}} = \cos \theta

Yesterday, I used the Pythagorean identity again to find \sin \theta. Today, I’ll instead plug back into the original equation 3 \sin \theta = \cos \theta:

3 \sin \theta = \cos \theta

3 \sin \theta = \displaystyle \frac{3}{\sqrt{10}}

\sin \theta = \displaystyle \pm \frac{1}{\sqrt{10}}

Unlike the example yesterday, the signs of \sin \theta and \cos \theta must agree. That is, if \cos \theta = \displaystyle \frac{3}{\sqrt{10}}, then \sin \theta = \displaystyle \frac{1}{\sqrt{10}} must also be positive. On the other hand, if \cos \theta = \displaystyle -\frac{3}{\sqrt{10}}, then \sin \theta = \displaystyle -\frac{1}{\sqrt{10}} must also be negative.

If they’re both positive, then

\sin \theta \cos \theta = \displaystyle \left( \frac{1}{\sqrt{10}} \right) \left( \frac{3}{\sqrt{10}} \right) =\displaystyle \frac{3}{10},

and if they’re both negative, then

\sin \theta \cos \theta = \displaystyle \left( -\frac{1}{\sqrt{10}} \right) \left( -\frac{3}{\sqrt{10}} \right) = \displaystyle \frac{3}{10}.

Either way, the answer must be \displaystyle \frac{3}{10}.

This is definitely superior to the solution provided in yesterday’s post, as there’s absolutely no doubt that the product \sin \theta \cos \theta must be positive.

Different ways of solving a contest problem (Part 2)

The following problem appeared on the American High School Mathematics Examination (now called the AMC 12) in 1988:

If 3 \sin \theta = \cos \theta, what is \sin \theta \cos \theta?

When I presented this problem to a group of students, I was pleasantly surprised by the amount of creativity shown when solving this problem.

Yesterday, I presented a solution using triangles. Here’s a second solution that I received: begin by squaring both sides and using a Pythagorean trig identity.

9 \sin^2 \theta = \cos^2 \theta

9 (1 - \cos^2 \theta) = \cos^2 \theta

9 - 9 \cos^2 \theta = \cos^2 \theta

9 = 10 \cos^2 \theta

\displaystyle \frac{9}{10} = \cos^2 \theta

\displaystyle \pm \frac{3}{\sqrt{10}} = \cos \theta

We use the Pythagorean identity again to find \sin \theta:

\displaystyle \frac{9}{10} = \cos^2 \theta

\displaystyle \frac{9}{10} = 1 - \sin^2 \theta

\sin^2 \theta = \displaystyle \frac{1}{10}

\sin \theta = \displaystyle \pm \frac{1}{\sqrt{10}}

Therefore, we know that

\sin \theta \cos \theta = \displaystyle \left( \pm \frac{1}{\sqrt{10}} \right) \left( \pm \frac{3}{\sqrt{10}} \right) = \displaystyle \pm \displaystyle \frac{3}{10},

so the answer is either \displaystyle \frac{3}{10} or \displaystyle -\frac{3}{10}. However, this was a multiple-choice contest problem and \displaystyle -\frac{3}{10} was not listed as a possible answer, and so the answer must be \displaystyle \frac{3}{10}.

green lineFor a contest problem, the above logic makes perfect sense. However, the last step definitely plays to the fact that this was a multiple-choice problem, and the concluding step would not have been possible had \displaystyle -\frac{3}{10} been given as an option.