Deriving π Using Area of a Circle © Laura Shears, 2010
As a student in the 1980's I wanted to prove that the area of a circle is indeed a constant times the radius squared, but what I ended up doing was deriving that constant. I started off with the assumption that the circumference of a circle is proportional to its diameter. This seems reasonable since the perimeter of an n-sided regular polygon is proportional to its diameter and a circle is the limit as n goes towards infinity of an n-sided regular polygon. Let's go ahead and call that constant π.
Break up the circle into n cocentric circles with adjacent circles having a difference of radius of dx. Since we have n of these circles, we know dx=r/n. Now we calculate the area of circles by cutting along the circular lines and unrolling each piece. As n goes to infinity the top and bottom of each unrolled piece will get closer in size to each other thus making the unrolled pieces approximate a rectangle. The radius of the ith circle is
Next, take the limit as n goes to infinity of the above expression to get A = πr2. So the original objective of showing that the area of a circle is proportional to the square of its radius has been achieved. That means that the area of a circle is a constant times its radius squared and we can call that constant π. The next object is to derive π. We can accomplish this by using another method for calculating the area of the circle and then setting the two methods equal to each other.
We could also calculate the area of a circle by taking the limit of the Reimann Sum of the rectangles in the first quadrant and multiplying that area by 4. If we center the circle at the origin, the equation of the circle is:
Setting this equals to πr2 gives us: . When I first came up with this formula, I owned a TI-58C calculator. I program it with the formula and left it plugged in for several days after putting a value in for n large enough to get just 5 decimal place accuracy. Later, I learned about the Trapezoid Rule and Simpson's Approximation and I rederived the formula for π using those techniques. With the Trapezoid Rule it took a smaller value of n to get 5 decimal place accuracy and my calculator was able to do it in a few hours instead of days. Simpon's Rule took only about 10 minutes to achieve the same level of accuracy.
Class Discussion Problems:
Try calculating pi with various values of n using the formula above and then try it again with the formulas that come from the Trapezoid Rule and Simpon's rule.
For some reason there is a limit to the number n that can make this program work of only 46340. However when I hard code a value for n into the program, n can be as large approximately 37,000,000, thus creating a bit more accuracy.