Achievements Unlocked

Today is a big day for akickinthediscovery.com! Not only is today π\pi day, basically the holiday that celebrates math(s), but this π\pi day article will be the 50th50^{th} article written on akickinthediscovery.com. Also, just as significant, it’s Einstein‘s birthday! Last π\pi day feels like yesterday, when we went through one of Euler’s solutions to the Basel problem. Today, we will keep building on our knowledge of π\pi and prove that π\pi is an irrational number. We will follow a brilliant proof devised by Ivan Niven (also, there is a typo in his paper, so we fixed it!), but before we show that, we will show 2\sqrt{2} and ee (Euler’s number) are irrational as well. This is so that there is one short resource that has proof of all three!

Setting the Stage

Before we move on, just to make sure that we are all on the same page, let’s recall what a rational number is.

(Incomplete) Definition (Rational Numbers \mathbb{Q}): We call the set of integer fractions, denoted ,\mathbb{Q},the set of rational numbers. That is,

={ab:a,b&b0}.\mathbb{Q} = \left\{\frac{a}{b}\;:\;a,b\in \Z \;\;\And \;\;b\neq 0\right\}.

Remark: The reason why the above definition is incomplete has to do with the fact that we are not giving structure to \mathbb{Q} that it ought to have. For example, we know that 12=36=816=4386=,\frac{1}{2} = \frac{3}{6} = \frac{-8\;\;}{16} = \frac{43}{86} = \cdots, however, this is not deducible from our “definition” of .\mathbb{Q}. But, for our purposes, we will think of \mathbb{Q} as the set of integer fractions and assume that all the properties and rules that we know about manipulating fractions hold without issue.

Now that we know what rational numbers are, we can define irrational numbers in a similarly incomplete way.

(Incomplete) Definition (Irrational Numbers /\R/\mathbb{Q}): We call real numbers that are not rational numbers irrational numbers. That is,

/={r:r}.\R/\mathbb{Q} = \left\{r\in \R \;:\; r\notin \mathbb{Q}\right\}.

Remark: This definition is similarly incomplete. For instance, we are assuming that we know that the real numbers are. But this discussion will take us too much into the weeds. So, for those interested, you can search for the construction of the real numbers on Google. The two classic ways are to use Dedekind cuts or to use equivalence classes of Cauchy sequences of rational numbers.

We are now set up to show three things: (i) 2\sqrt{2} is irrational; (ii) ee is irrational; and the coup de grâce (iii) π\pi is irrational.

2\sqrt{2} Using Fermat’s Infinite Descent

Theorem: The square root of 2 is not a rational number,

2.\sqrt{2}\notin \mathbb{Q}.
Proof: (Click in the Discovery)

Key Idea: We will see a trend in all the proofs presented here; in particular, every proof is by contradiction. It’s usually very difficult to prove that something is impossible, and in our case, it’s very difficult to show it’s impossible to find integers whose ratio equals 2\sqrt{2}, i.e., that it’s impossible to find a,b,a,b\in \Z, such that ab=2.\frac{a}{b} =\sqrt{2}. It’s far easier to assume there are two integers a,ba,b\in \Z such that ab=2,\frac{a}{b} =\sqrt{2}, and then deduce a contradiction. This will be our strategy.

Beginning of proof: Assume for the hope of a contradiction that 2\sqrt{2} is a rational number. It follows that we can find two integers a,ba,b\in \Z such that ab=2.\frac{a}{b} =\sqrt{2}. Note that we can safely assume that aa and bb are both positive. Since working with square roots and fractions can be annoying, let’s square and rearrange everything to deduce,

a2=2b2.a^2 = 2b^2.

Intuitively, we know that if a2a^2 is even, then aa must also be even. Rigorously, this is a consequence of Euclid’s lemma (or Euler’s Lemma). (Or, you can do a mini proof by contrapositive: if aa is not even (odd), then a2a^2 is not even (odd).) In any case, we conclude that aa is even, and therefore a=2a0,a = 2a_0, for some a0.a_0\in \Z. Let’s plug this into our equation relating aa and b,b,

(2a0)2=4a02=2b2(2a_0)^2 = 4a_0^2 = 2b^2

Canceling a factor of 2 gives,

2a02=b2.2a_0^2 = b^2.

Similarly, we see that b2b^2 is even, and hence bb is even. Therefore, b=2b0b = 2b_0 for some b0.b_0\in \Z. It follows,

2a02=(2b0)2=4b02.2a_0^2 = (2b_0)^2 = 4b_0^2.

Again, canceling a factor of 2 gives,

a02=2b02.a_0^2 = 2b_0^2.

Observe that the above equation implies that 2=a0b0.\sqrt{2}=\frac{a_0}{b_0} . This means that we have found a new way to express the square root of 2 as a fraction of positive integers:

2=ab=a0b0,where0<a0<aand0<b0<b.\sqrt{2}=\frac{a}{b} =\frac{a_0}{b_0}, \qquad\qquad \mathrm{where} \;\; 0<a_0<a\;\; \mathrm{and} \;\; 0<b_0<b.

Note that 0<a0<a0<a_0<a and 0<b0<b0<b_0<b follow from a=2a0a = 2a_0 and b=2b0.b = 2b_0. Interesting…

Taking a look back at a02=2b02,a_0^2 = 2b_0^2, you might notice that we are in the exact same situation that we started, except we are using the letters a0a_0 and b0b_0 instead of aa and b.b. Thus, by going through the same argument, we can find two more positive integers a1a_1 and b1b_1 such that: (i) a12=2b12;a_1^2 = 2b_1^2;(ii) 2=a1b1;\sqrt{2}=\frac{a_1}{b_1}; and (iii) a0=2a1a_0 = 2a_1 and b0=2b1.b_0 = 2b_1. In particular, (ii) and (iii) imply that we have found yet another way to express 2\sqrt{2} as a fraction of smaller positive integers: 0<a1<a0<a0<a_1<a_0<a and 0<b1<b0<b.0<b_1<b_0<b. And by starting with (i), we can go through the same reasoning yet again… and then again, and again indefinitely, always finding smaller and smaller numbers an,bna_n,b_n\in \N to express 2\sqrt{2} as a fraction of integers.

2=ab=a0b0=a1b1=anbn,where0<<an<<a1<a0<aand0<<bn<<b1<b0<b.\sqrt{2}=\frac{a}{b} =\frac{a_0}{b_0} =\frac{a_1}{b_1} \cdots =\frac{a_n}{b_n}\cdots, \qquad\qquad \begin{matrix}\mathrm{where} \;\; 0<\cdots< a_n<\cdots<a_1<a_0<a\;\;\\ \;\,\mathrm{and} \;\;\; 0<\cdots< b_n<\cdots<b_1<b_0<b. \end{matrix}

But, this process must stop at some point since there are only a finite number of positive integers that ana_n and bnb_n could be. In particular, the numerator ana_n can only be one of the numbers 1,2˙,3˙,,a11,\.2,\.3,\cdots, \,a-1 and the denominator bnb_n can only be one of the numbers 1,2,3,,b˙1.1,\,2,\,3,\,\cdots, \.b-1. This is our contradiction and hence 2.\sqrt{2}\notin \mathbb{Q}.

Remark:  Fermat’s infinite descent is the particular type of contradiction that we reached. That is, assuming 2\sqrt{2} is rational, let us find an infinite sequence of fractions to represent 2,\sqrt{2}, where each new representation is constructed from the previous one and uses smaller positive integers. However, there is not an infinite number of positive integers less than aa or b.b. Thus, a contradiction.

\square

Euler’s Number ee is not rational

Theorem: Define Euler’s number e,e, to be equal to,

e=n=01n!.e = \sum_{n=0}^{\infty} \frac{1}{n!}.

Then, Euler’s number is an irrational number.

This proof is a classic!

One quick fact that we will use is that 2<e<3.2<e< 3. This follows from the following two observations: e>10!+11!=2e>\frac{1}{0!}+ \frac{1}{1!} =2 and,

e=n=01n!=10!+n=11n!\;\;\;\;\;\;\;\;\;e = \sum_{n=0}^{\infty} \frac{1}{n!} = \frac{1}{0!}+\sum_{n=1}^{\infty} \frac{1}{n!}
<1+n=112n1\lt 1 +\sum_{n=1}^{\infty} \frac{1}{2^{n-1}}
=1+n=012n= 1 + \sum_{n=0}^{\infty} \frac{1}{2^n} \;\;\;\;
=1+111/2= 1+ \frac{1}{1-1/2} \;\;\;
=3.= 3 .\qquad\qquad\;\;\;\;

Where we used that 211=1!,2^{1-1} =1!, and 221=2!,2^{2-1} =2!, and 2n1<n!2^{n-1}\lt n! for all n3n\geq 3 in going from the first line to the second. (I challenge you to prove this inequality, have fun!)

Proof: (Click in the Discovery)

Key Idea: Factorials grow super-duper fast. So fast that the series that defines the number ee converges quickly. So quick that, in some sense, ee cannot be a rational number. Let’s make this idea more precise.

Beginning of proof: Assume for the hope of a contradiction that ee were rational. It follows that we can find two integers a,ba,b\in \Z such that ab=e.\frac{a}{b} =e. First, note that we can assume that aa and bb are both positive. Second, note that b2b\geq2 by observing that 2<e<32<e< 3 and there is no way to express a number greater than 2 and less than 3 as a fraction of positive integers with the denominator equal to 1. So, we have ab=e\frac{a}{b} =e for a,ba,b\in \N and b2.b\geq2 .

Now, consider e(b!).e \cdot (b!).

eb!=n=0b!n!=b!0!+b!1!+b!2!++b!b!+b!(b+1)!+.e\cdot b! = \sum_{n=0}^{\infty} \frac{b!}{n!} = \frac{b!}{0!}+\frac{b!}{1!}+ \frac{b!}{2!} + \cdots + \frac{b!}{b!} + \frac{b!}{(b+1)!} + \cdots.

Observe that the denominator is less than the numerator for the first b+1.b+1. Consequently, the first b+1b+1 terms are integers, and hence the sum of the first b+1b+1 terms is an integer, which we denote by N.N. To highlight this observation, let’s group the integer terms,

eb!=(b!0!+b!1!+b!2!++b!b!)=N+b!(b+1)!+b!(b+2)!+.e\cdot b! = \underbrace{\Bigg( \frac{b!}{0!}+\frac{b!}{1!}+ \frac{b!}{2!} + \cdots + \frac{b!}{b!}\Bigg)}_{=\;N \in \Z} + \frac{b!}{(b+1)!} + \frac{b!}{(b+2)!} + \cdots.

After some rearrangement, we have

eb!N=1b+1+1(b+1)(b+2)+.e\cdot b! – N= \frac{1}{b+1} + \frac{1}{(b+1)(b+2)} + \cdots.

Moreover, eb!e\cdot b! is an integer as well. Indeed, since e=abe = \frac{a}{b} it follows eb!=a(b1)!.e\cdot b! = a\cdot(b-1)!\in \Z. Thus, the left-hand side is an integer, and we conclude

Integer=1b+1+1(b+1)(b+2)+.\mathrm{Integer}\;=\; \frac{1}{b+1} + \frac{1}{(b+1)(b+2)} + \cdots.

However, we claim that the right-hand side is not an integer, and once we show this, we will have found our contradiction. Observe that,

0<1b+1+1(b+1)(b+2)+<1b+1b2+=n=11bn=1/b11/b=1b11forb2.0<\frac{1}{b+1} + \frac{1}{(b+1)(b+2)} + \cdots \lt \frac{1}{b} + \frac{1}{b^2} + \cdots = \sum_{n=1}^{\infty} \frac{1}{b^n} = \frac{1/b}{1-1/b} = \frac{1}{b-1}\leq 1 \qquad \mathrm{for} \; b\geq 2.

So that the tail of the sum is greater than 0 but less than 1. We have found our contradiction! We conclude that ee is irrational.

\square

The Coup de Grâce π\pi is not rational.

Theorem: π\pi is irrational.

Again, we follow a brilliant proof by Ivan Niven.

Proof: (Click in the Discovery)

Key Idea: The idea behind this proof is to find a quantity that we can show to be an integer under the assumption that π is rational, which must be greater than 0 but less than 1. The difficulty and brilliance come from actually finding that quantity.

Beginning of proof: Assume for the hope of a contradiction that π\pi is a rational number. It follows that we can find two integers a,ba,b\in \Z such that ab=π.\frac{a}{b} =\pi. We can safely assume that aa and bb are both positive. Now consider the polynomials,

f(x)=xn(abx)nn!,f(x) = \frac{x^n(a-bx)^n}{n!} ,

and

F(x)=k=0n(1)kf(2k)(x).F(x) = \sum_{k=0}^{n}(-1)^{k}f^{(2k)}(x) .

Our goal is to find some nn\in \N that leads to a contradiction. With this goal in mind, we make some observations regarding f(x)f(x) and F(x).F(x).


8 Key Observations regarding f(x)f(x) and F(x)F(x):

\bigstarObservation (1): f(x)n!f(x) \cdot n! is a polynomial with integer coefficients, i.e., f(x)n![x].f(x) \cdot n! \in \Z[x].

\bigstarObservation (2): f(x)>0f(x) \gt 0 for all x(0,π).x\in (0,\pi). Indeed, since x<π=abx<\pi=\frac{a}{b} implies a>bx.a>bx.

\bigstarObservation (3): f(x)<(πa)nn!f(x) \lt \frac{(\pi a)^n}{n!} for all x(0,π).x\in (0,\pi). This is seen by noting that x<πx \lt \pi for x(0,π)x\in (0,\pi) and abx<aa-bx \lt a for x(0,π).x\in (0,\pi).

\bigstarObservation (4): the lowest power of xx that shows up in f(x)f(x) is xn.x^n.

\bigstarObservation (5): f(x)=bnf(πx).f(x) =b^n f(\pi -x). Indeed, since ab=π\frac{a}{b} =\pi we have,

bnf(πx)=bnf(abx)=bn(abx)n(ab(abx))nn!b^nf(\pi -x) =b^nf\left(\frac{a}{b} -x\right) = b^n\frac{\left(\frac{a}{b}-x \right)^n\left(a-b\left(\frac{a}{b}-x\right)\right)^n}{n!}
=(abx)n(aa+x)nn!=xn(abx)nn!\qquad= \frac{ \left(a-bx \right)^n \left(a-a+x\right)^n}{n!} = \frac{x^n(a-bx)^n}{n!}
=f(x).= f(x) .\qquad\qquad\qquad\qquad\qquad\qquad\;\;\;\;

A corollary of this result is the following: f(k)(x)=bnf(k)(πx).f^{(k)}(x) =b^n f^{(k)}(\pi -x). Which we will make use of.

\bigstarObservation (6): f(k)(0)=f(k)(π)=0f^{(k)}(0) =f^{(k)}(\pi) =0 for all 0k<n.0\leq k\lt n. This follows from observations (4) and (5). Since the lowest power of xx in f(x)f(x) is xn,x^n, the kthk^{th} derivative for 0k<n0\leq k\lt n will not yield a constant term, that is, it leaves an xx in every term in f(k)(x).f^{(k)}(x) . Hence, f(k)(0)=0.f^{(k)}(0) =0. Now, using (5) we conclude that f(k)(0)=bnf(k)(π)=0.f^{(k)}(0) =b^nf^{(k)}(\pi) =0.

\bigstarObservation (7): f(k)(0)=bnf(k)(π)f^{(k)}(0) =b^nf^{(k)}(\pi) \in \Z for all nk2n.n\leq k\leq 2n. This observation requires more justification.

Note that when we expand the product in f(x),f(x) , we get

f(x)==0ncn!xn+,f(x) = \sum_{\ell=0}^n \frac{c_{\ell}}{n!} \cdot x^{n+\ell},

where c.c_{\ell}\in \Z. We are interested in determining f(k)(0).f^{(k)}(0) . This means we can focus on the constant term in f(k)(0)f^{(k)}(0) because all other terms will be 0. Observe that the constant term in f(k)(0)f^{(k)}(0) will be the term with xkx^{k} in f(x).f(x) . Then, by taking the kthk^{th} derivative, we bring down a factor of k!k! and we deduce that the constant term in f(k)(0)f^{(k)}(0) is equal to k!cknn!.\frac{k!\cdot c_{k-n}}{n!}. And, since nk2nn\leq k\leq 2n we see that k!cknn!.\frac{k!\cdot c_{k-n}}{n!}\in \Z.

Now, using observation (5) we deduce f(k)(0)=bnf(k)(π).f^{(k)}(0) =b^nf^{(k)}(\pi) \in \Z.

\bigstarObservation (8): F(0)F(0) and F(π)F(\pi) are integers. This follows from the previous items (6) and (7). In particular, the sum F(0)+F(π)F(0) + F(\pi) is an integer.


Now, let’s set up a differential equation in F(x),F(x) ,

ddx(F(x)sin(x)F(x)cos(x))=F(x)sin(x)+F(x)cos(x)F(x)cos(x)+F(x)sin(x)\frac{d}{dx} \Bigg(F'(x) \;\sin(x) – F(x) \;\cos(x) \Bigg) = F”(x) \;\sin(x) + F'(x) \;\cos(x)- F'(x) \;\cos(x) +F(x) \;\sin(x)
=F(x)sin(x)+F(x)sin(x). = F”(x) \;\sin(x) + F(x) \;\sin(x) .
=f(x)sin(x).=f(x) \;\sin(x).\qquad\qquad\qquad\;\;\,

Therefore,

0πf(x)sin(x)dx=(F(x)sin(x)F(x)cos(x))|0π=F(π)+F(0).\int_0^{\pi}f(x) \,\sin(x)\;dx =\Big(F'(x) \;\sin(x) – F(x) \;\cos(x)\Big){\Huge |}_0^{\pi} = F(\pi) + F(0)\in \Z.

We will use the fact that the integral above is an integer to find a contradiction. Our strategy is to use our observations about f(x)f(x) to bound the integrand, then use this bound to find an nn\in \N such that the integral cannot be an integer.

Observe that the integrand, f(x)sin(x),f(x) \sin(x), is strictly bounded by,

0<f(x)sin(x)<(πa)nn!0\lt\;f(x) \sin(x)\lt\; \frac{(\pi a)^n}{n!}

for 0<x<π.0\lt x\lt \pi. These follow from Observation (2) and Observation (3). It follows that the integral, 0πf(x)sin(x)dx,\int_0^{\pi}f(x) \,\sin(x)\;dx , is positive and bounded by

0<0πf(x)sin(x)dxsupx[0,π](f(x)sin(x))0πdxπ(πa)nn!.0<\int_0^{\pi}f(x) \,\sin(x)\;dx \leq \sup_{x\in [0,\pi]} {(f(x) \,\sin(x))} \cdot \int_0^{\pi} \;dx \leq\pi\frac{(\pi a)^n}{n!} .

for any choice of n.n. However, we can see from the upper bound above that the integral can be made to be strictly less than by choosing a sufficiently large n.n.

0<0πf(x)sin(x)dx<1forsufficientlylargen.0<\int_0^{\pi}f(x) \,\sin(x)\;dx <1 \qquad \mathrm{for}\;\mathrm{sufficiently }\;\mathrm{large} \;n.

This is our contradiction! The integral cannot both be an integer and be in the interval (0,1).(0,1).

Thus, π\pi is irrational!

\square

Wow…

The number π is incredibly fascinating to people who study mathematics. There are an infinity of reasons why. It might be because π is defined to be equal to the ratio of the circumference of a circle to its diameter, and yet somehow π shows up mysteriously all over the place! For instance, π shows up in the sums,

n=11n2=π26,\sum_{n=1}^{\infty} \frac{1}{n^2} = \frac{\pi^2}{6},
n=1(1)α(n)n=π&n=1(1)β(n)n=π2,\sum_{n=1}^{\infty} \frac{(-1)^{\alpha(n)}}{n} = \pi\qquad\And \qquad \sum_{n=1}^{\infty} \frac{(-1)^{\beta(n)}}{n} = \frac{\pi}{2},

where α(n)=#{primefactorsofncongruentto1(mod4)}\alpha(n) = \#\{\mathrm{prime}\;\mathrm{factors}\;\mathrm{of}\;n\;\mathrm{congruent}\;\mathrm{to}\;1 \;\;(\mathrm{mod}\;{4})\}, and β(n)=#{primefactorsofncongruentto3(mod4)}.\beta(n) = \#\{\mathrm{prime}\;\mathrm{factors}\;\mathrm{of}\;n\;\mathrm{congruent}\;\mathrm{to}\; 3\;\;(\mathrm{mod}\;{4})\}. Or, π shows up in,

(p1(mod4)pp1)(p3(mod4)pp+1)=π4,\Bigg(\prod_{p\equiv 1 \;(\mathrm{mod}\;4)}\frac{p}{p-1} \Bigg)\cdot \Bigg(\prod_{p\equiv 3 \;(\mathrm{mod}\;4)}\frac{p}{p+1} \Bigg) = \frac{\pi}{4},

And those are formulas only discovered by Euler! We also have the Wallis product

n=1(2n)2(2n1)(2n+1)=(2123)(4345)(6567)=π2\prod_{n=1}^{\infty} \frac{(2n)^2}{(2n-1)(2n+1)} = \Big(\frac{2}{1}\cdot\frac{2}{3}\Big) \Big(\frac{4}{3}\cdot\frac{4}{5}\Big) \Big(\frac{6}{5}\cdot\frac{6}{7}\Big)\cdots = \frac{\pi}{2}

If sums and products aren’t enough,

ex2dx=π.\int_{-\infty}^{\infty}e^{-x^2}\;dx = \sqrt{\pi}.

And last but not least,

eiπ+1=0.e^{i\pi}+1 =0.

And, other than that last identity (Euler’s Identity), it’s not at all obvious why π, a number defined in terms of circles, would show up! Furthermore, π also has some mysteries left for us to solve. For instance, we don’t know if π is a normal number. Or, whether ππππ\pi^{{\pi}^{\pi^{\pi}}} is an integer! How amazing! Maybe you will solve these problems!

Be Kind. Be Curious. Be Compassionate. Be Creative.

And Have Fun!

Leave a comment