Processing math: 100%

Monday, November 25, 2013

Trigonometry and Nested Radicals

Early last month, I was chatting with one of my officemates about a curious problem I had studied in high school. I hadn't written any of the results down, so much of the discussion involved me rediscovering the results and proving them with much more powerful tools than I once possessed.

Before writing about the problem I had played around with, I want to give a brief motivation. For as long as humans have been doing mathematics, finding values of π has been deemed worthwhile (or every generation has just found it worthwhile to waste time computing digits).

One such way the Greeks (particularly Archmides) computed π was by approximating a circle by a regular polygon and letting the number of sides grow large enough so that the error between the area of the unit circle (π12) and the area of the polygon would be smaller than some fixed threshold. Usually these thresholds were picked to ensure that the first k digits were fully accurate (for some appropriate value of k).

In many introductory Calculus courses, this problem is introduced exactly when the limit is introduced and students are forced to think about the area problem in the regular polygon:
Given N sides, the area is NTN where TN is the area of each individual triangle given by one side of the polygon and the circumcenter.

Call one such triangle ΔABC and let BC be the side that is also a side of the polygon while the other sides have |AB|=|AC|=1 since the polygon is inscribed in a unit circle. The angle BAC=2πN since each of the triangles has the same internal angle and there are N of them. If we can find the perpendicular height h from AB to C, the area will be 12h|AB|=h2. But we also know that
sin(BAC)=h|AC|h=sin(2πN). Combining all of these, we can approximate π with the area:
πN2sin(2πN)=πsin(2πN)2πN. As I've shown my Math 1A students, we see that
limNπsin(2πN)2πN=πlimx0sin(x)x=π so these are indeed good approximations.

Theory is Nice, But I Thought We Were Computing Something

Unfortunately for us (and Archimedes), computing sin(2πN) is not quite as simple as dividing by N, so often special values of N were chosen. In fact, starting from N and then using 2N, the areas could be computed via a special way of averaging the previous areas. Lucky for us, such a method is equivalent to the trusty half angle identities (courtesy of Abraham De Moivre). To keep track of these polygons with a power of two as the number of sides, we call An=2n2sin(2π2n).

Starting out with the simplest polygon, the square with N=22 sides, we have
A2=2sin(π2)=2. Jumping to the octagon (no not that "The Octagon"), we have
A3=4sin(π4)=422=22. So far, the toughest thing we've had to deal with is a 45 angle and haven't yet had to lean on Abraham (himnot him) for help. The hexadecagon wants to change that:
A4=8sin(π8)=81cos(π4)2=8224=422.
To really drill home the point (and motivate my next post) we'll compute this for the 32-gon (past the point where polygons have worthwhile names):
A5=16sin(π16)=161cos(π8)2. Before, we could rely on the fact that we know that a 454590 triangle looked like, but now, we come across cos(π8), a value which we haven't seen before. Luckily, Abraham has help here as well:
cos(π8)=1+cos(π4)2=2+24=122+2 which lets us compute
A5=161122+22=822+2.

So why have I put you through all this? If we wave our hands like a magician, we can see this pattern continues and for the general n
An=2n222+2++2
where there are n3 nested radicals with the sign and only one minus sign at the beginning.

This motivates us to study two questions, what is the limiting behavior of such a nested radical:
2+s12+s2 as the signs s1,s2, takes values in {1,1}. Recasting in terms of the discussion above, we want to know how close we are to π as we increase the number of sides.

When I was in high school, I just loved to nerd out on any and all math problems, so I studied this just for fun. Having heard about the unfathomable brain of Ramanujan and the fun work he had done with infinitely nested radicals, I wanted to examine which sequences of signs (s1,s2,) produced an infinite radical that converged and what the convergence behavior was.

I'm fairly certain my original questions came from an Illinois Council of Teachers of Mathematics (ICTM) contest problem along the lines of
Find the value of the infinite nested radical 2+2+ or maybe the slightly more difficult Find the value of the infinite nested radical 22+22+. Armed with my TI-83, I set out to do some hardcore programming and figure it out. It took me around a month of off-and-on tinkering. This second time around as a mathematical grown-up, it took me the first half of a plane ride from SFO to Dallas.

In the next few weeks/months, I hope to write a few blog posts, including math, proofs and some real code on what answers I came up with and what other questions I have.