I read your proof, and I don’t quite get it. Or rather, I get most of the steps in a hand-waving manner, but some are fairly obscure and the others I can’t quite fully justify to myself. This is just a reflection of my own limits and not your description, of course.
That said, Niven’s proof on the Wikipedia page was understandable to the point that I can describe it from memory. I’ll repeat it in my own words since I did have to think through a few things that confused me a bit initially.
This is a proof by contradiction, where we start by assuming pi=a/b (a and b being positive integers). We examine this function:
f(x) = x[sup]n[/sup](a - bx)[sup]n[/sup] / n!
Technically this should be f[sub]n/sub, since n is a parameter, but I’ll leave it off and we’ll just keep in mind that we can change n to whatever we want.
First, note that f(0) = f(pi). If you substitute x => (a/b - x) and do a little algebra, you get the same equation (i.e., f(x) = f(pi - x)). So whatever we conclude about f(0) will apply to f(pi) later.
Let’s examine the derivatives: we can see that the lowest order term (of x) is x[sup]n[/sup], which means that f(0) is equal to 0 for the first n-1 derivatives. Note that as we’re doing this, we’re progressively canceling out the n! term at the bottom. At the nth derivative, there’s finally a constant term (not multiplied by a power of x), but nothing left in the denominator. After that, no fractions are introduced in further derivatives. At the (2n)th derivative, there’s only a constant term, and after that nothing.
The big implication here is that f(0), f’(0), etc. are all integers. We don’t know (or care) what their values are, just that there are no fractions in there. Furthermore, f(pi) (and its derivatives) is an integer since it’s equal to f(0).
The next step involves this formula:
F(x) = f(x) - f’'(x) + f[sup]4/sup - … -1[sup]n[/sup]f[sup]2n/sup
It’s pretty straightforward to see that this holds:
F’’ + F = f
Since the terms alternate in sign, F’’ mostly cancels the terms in F, with the exception of f(x) and the final term. But that final term is the 2n+2 derivative, which we know to be 0. So it’s really just f.
The next step is to show this is true:
(F’sin - Fcos)’ = f*sin
A little application of the product rule (and linearity of derivatives) get you this:
F’'sin + Fsin = f*sin
And then:
(F’’ + F)sin = f*sin
So we see that the formula works. Going back, we can also turn this (va fundamental theorem of calculus):
(F’sin - Fcos)’ = f*sin
Into this:
Int[0…pi] f*sin = [0…pi] (F’sin - Fcos)
Looking at the right hand side, we see that we can ignore the F’sin term because sin is 0 at 0 and pi. From before, F(0) and F(pi) are integers, and cos(0)=1 and cos(pi)=-1. From all of this, we can conclude that:
Int[0…pi] f*sin = [0…pi] (F’sin - Fcos) = an integer
But let’s try to compute an upper bound for:
Int[0…pi] f*sin
We’ll so this by looking at the upper bounds of the different components. sin(x) of course has a maximum of 1.0 in that range. f is a little trickier, but let’s do a little transformation:
x[sup]n[/sup](a - bx)[sup]n[/sup] / n! = (ax - bx[sup]2[/sup])[sup]n[/sup] / n!
That bx[sup]2[/sup] term is negative… so it can only make the value smaller. We’ll ignore it. Therefore, one upper bound is:
(ax)[sup]n[/sup] / n!
Put together, this is an upper bound for the integral:
pi * (ax)[sup]n[/sup] / n!
But look: that n! grows much faster than (ax)[sup]n[/sup]. No matter what a is, and even with the maximum value of pi for x, n! can be bigger yet. So we can always choose an n where our integral is less than 1.
However, we already decided that the integral must evaluate to an integer. If it can’t even be as large as 1, it must be zero.
Well, it can’t be that either: f(x) and sin(x) are all positive for 0<x<pi. Therefore the integral is positive. Small, maybe, but not zero.
There’s no integer between 0 and 1, so we have our contradiction. Pi can’t be expressed as a/b and is irrational.