Proof that π is irrational

The number π (pi) has been studied since ancient times, and so has the concept of irrational numbers. An irrational number is any real number that cannot be expressed as a fraction a/b, where a is an integer and b is a non-zero integer.

It was not until the 18th century that Johann Heinrich Lambert proved that π is irrational. In the 19th century, Charles Hermite found a proof that requires no prerequisite knowledge beyond basic calculus. A simplification of Hermite's proof is due to Mary Cartwright. Two other such proofs are due to Ivan Niven and to Miklós Laczkovich.

In 1882, Ferdinand von Lindemann proved that π is not just irrational, but transcendental as well.[1]

Lambert's proof

Scan of formula on page 288 of Lambert's "Mémoires sur quelques propriétés remarquables des quantités transcendantes, circulaires et logarithmiques", Mémoires de l'Académie royale des sciences de Berlin (1768), 265322.

In 1761, Lambert proved that π is irrational by first showing that this continued fraction expansion holds:

\tan(x) = \cfrac{x}{1 - \cfrac{x^2}{3 - \cfrac{x^2}{5 - \cfrac{x^2}{7 - {}\ddots}}}}.

Then Lambert proved that if x is non-zero and rational then this expression must be irrational. Since tan(π/4) = 1, it follows that π/4 is irrational and therefore that π is irrational.[2] A simplification of Lambert's proof is given below. This result can also be proved using even more basic tools of calculus (integrals instead of series).[3][4]

Hermite's proof

This proof uses the characterization of π as the smallest positive number whose half is a zero of the cosine function and it actually proves that π2 is irrational.[5][6] As in many proofs of irrationality, the argument proceeds by reductio ad absurdum.

Consider the sequences (An)n  0 and (Un)n  0 of functions from R into R thus defined:

  1. A_0(x)=\sin(x);\,
  2. (\forall n\in\mathbb{Z}_+):A_{n+1}(x)=\int_0^xyA_n(y)\,dy;
  3. U_0(x)=\frac{\sin(x)}x;
  4. (\forall n\in\mathbb{Z}_+):U_{n+1}(x)=-\frac{U_n'(x)}x.

It can be proven by induction that

(\forall n\in\mathbb{Z}_+):A_n(x)=\frac{x^{2n+1}}{(2n+1)!!}-\frac{x^{2n+3}}{2\times(2n+3)!!}+\frac{x^{2n+5}}{2\times4\times(2n+5)!!}\mp\cdots

and that

(\forall n\in\mathbb{Z}_+):U_n(x)=\frac1{(2n+1)!!}-\frac{x^2}{2\times(2n+3)!!}+\frac{x^4}{2\times4\times(2n+5)!!}\mp\cdots

and therefore that

U_n(x)=\frac{A_n(x)}{x^{2n+1}}.

So

\frac{A_{n+1}(x)}{x^{2n+3}}=U_{n+1}(x)=-\frac{U_n'(x)}x=-\frac1x\frac d{dx}\left(\frac{A_n(x)}{x^{2n+1}}\right),

which is equivalent to

A_{n+1}(x)=(2n+1)A_n(x)-xA_n'(x)=(2n+1)A_n(x)-x^2A_{n-1}(x).\,

It follows by induction from this, together with the fact that A0(x) = sin(x) and that A1(x) = x cos(x) + sin(x), that An(x) can be written as P_n(x^2)\sin(x)+xQ_n(x^2)\cos(x), where Pn and Qn are polynomial functions with integer coefficients and where the degree of Pn is smaller than or equal to ⌊n2⌋. In particular,

A_n\left(\frac\pi2\right)=P_n\left(\frac{\pi^2}4\right).

Hermite also gave a closed expression for the function A_n, namely

A_n(x)=\frac{x^{2n+1}}{2^nn!}\int_0^1(1-z^2)^n\cos(xz)\,dz.

He did not justify this assertion, but it can be proved easily. First of all, this assertion is equivalent to

\frac{1}{2^nn!}\int_0^1(1-z^2)^n\cos(xz)\,dz=\frac{A_n(x)}{x^{2n+1}}=U_n(x).

Proceeding by induction, take n = 0.

\int_0^1\cos(xz)\,dz=\frac{\sin(x)}x=U_0(x)

and, for the inductive step, consider any n  Z+. If

\frac{1}{2^nn!}\int_0^1(1-z^2)^n\cos(xz)\,dz=U_n(x),

then, using integration by parts and Leibniz's rule, one gets

\begin{align}
& {}\quad \frac{1}{2^{n+1}(n+1)!}\int_0^1(1-z^2)^{n+1}\cos(xz)\,dz \\
& =\frac{1}{2^{n+1}(n+1)!}\Biggl(\overbrace{\left.(1-z^2)^{n+1}\frac{\sin(xz)}x\right|_{z=0}^{z=1}}^{=\,0} + \int_0^12(n+1)(1-z^2)^nz\frac{\sin(xz)}x\,dz\Biggr)\\[8pt]
&=\frac1x\cdot\frac1{2^nn!}\int_0^1(1-z^2)^nz\sin(xz)\,dz\\[8pt]
&=-\frac1x\cdot\frac d{dx}\left(\frac1{2^nn!}\int_0^1(1-z^2)^n\cos(xz)\,dz\right) \\[8pt]
& =-\frac{U_n'(x)}x = U_{n+1}(x).
\end{align}

If π2/4 = p/q, with p and q in N, then, since the coefficients of Pn are integers and its degree is smaller than or equal to ⌊n2⌋, qn/2⌋Pn(π2/4) is some integer N. In other words,

\begin{align}N&=q^{\left\lfloor\frac n2\right\rfloor}A_n\left(\frac\pi2\right)\\&=q^{\left\lfloor\frac n2\right\rfloor}\frac{\left(\frac pq\right)^{n+\frac 12}}{2^nn!}\int_0^1(1-z^2)\cos\left(\frac\pi2z\right)\,dz.\end{align}

But this number is clearly greater than 0; therefore, N  N. On the other hand, the integral that appears here is not greater than 1 and

\lim_{n\in\mathbb{N}}q^{\left\lfloor\frac n2\right\rfloor}\frac{\left(\frac pq\right)^{n+\frac 12}}{2^nn!}=0.

So, if n is large enough, N < 1. Thereby, a contradiction is reached.

Hermite did not present his proof as an end in itself but as an afterthought within his search for a proof of the transcendence of π. He discussed the differential-recurrent relations to motivate and to obtain the convenient integral representation. Once the integral is obtained, there are various ways to present a succinct and self-contained proof starting from the integral (as in Cartwright's or Niven's presentations), which Hermite could easily see (as he did in his proof of the transcendence of e[7]).

Moreover, Hermite's proof is closer to Lambert's proof than it seems. In fact, An(x) is the "residue" (or "remainder") of Lambert's continued fraction for tan(x).[4]

Cartwright's proof

Harold Jeffreys wrote that this proof was set as an example in an exam at Cambridge University in 1945 by Mary Cartwright, but that she had not traced its origin.[8]

Consider the integrals

I_n(x)=\int_{-1}^1(1 - z^2)^n\cos(xz)\,dz.

Two integrations by parts give the recurrence relation

\text{for }n>1\text{, } x^2I_n(x)=2n(2n-1)I_{n-1}(x)-4n(n-1)I_{n-2}(x).

If

J_n(x)=x^{2n+1}I_n(x),\,

then this becomes

J_n(x)=2n(2n-1)J_{n-1}(x)-4n(n-1)x^2J_{n-2}(x).\,

Also

J_0(x)=2\sin(x)\, and J_1(x)=-4x\cos(x)+4\sin(x).\,

Hence for all n  Z+,

J_n(x)=x^{2n+1}I_n(x)=n!\bigl(P_n(x)\sin(x)+Q_n(x)\cos(x)\bigr),\,

where Pn(x) and Qn(x) are polynomials of degree ≤ 2n, and with integer coefficients (depending on n).

Take x = π2, and suppose if possible that π2 = ba, where a and b are natural numbers (i.e., assume that π is rational). Then

 \frac{b^{2n+1}}{n!} I_n\left(\frac\pi2\right) = P_n\left(\frac\pi2\right)a^{2n+1}.

The right side is an integer. But 0 < In(π2) < 2 since the interval [1, 1] has length 2 and the function which is being integrated takes only values between 0 and 1. On the other hand,

 \frac{b^{2n+1}}{n!} \to 0\text{ as }n \to \infty.

Hence for sufficiently large n

 0 < \frac{b^{2n+1}I_n\left(\frac\pi2\right)}{n!} < 1,

that is, we could find an integer between 0 and 1. That is the contradiction that follows from the assumption that π is rational.

This proof is similar to Hermite's proof. Indeed,

\begin{align}J_n(x)&=x^{2n+1}\int_{-1}^1 (1 - z^2)^n \cos(xz)\,dz\\&=2x^{2n+1}\int_0^1 (1 - z^2)^n \cos(xz)\,dz\\&=2^{n+1}n!A_n(x).\end{align}

However, it is clearly simpler. This is achieved bypassing the inductive definition of the functions An and taking as a starting point their expression as an integral.

Niven's proof

This proof uses the characterization of π as the smallest positive zero of the sine function.[9]

Preparation: Suppose that π is rational, i.e. π = a/b for some integers a and b ≠ 0, which may be taken without loss of generality to be positive. Given any positive integer n, we define the polynomial function

 f(x) = \frac{x^n(a - bx)^n}{n!},\quad x\in\mathbb{R},\!

and denote by

 F(x) = f(x) + \cdots + (-1)^j f^{(2j)}(x) + \cdots + (-1)^n f^{(2n)}(x),\quad x\in\mathbb{R},\!

the alternating sum of f and its first n even derivatives.

Claim 1: F(0) + F(π) is an integer.

Proof: Expanding f as a sum of monomials, the coefficient of xk is a number of the form ck/n! where ck is an integer, which is 0 if k < n. Therefore, f(k)(0) is 0 when k < n and it is equal to (k!/n!)ck if nk ≤ 2n; in each case, f(k)(0) is an integer and therefore F(0) is an integer.

On the other hand, f(πx) = f(x) and so (–1)kf(k)(πx) = f(k)(x) for each non-negative integer k. In particular, (–1)kf(k)(π) = f(k)(0). Therefore, f(k)(π) is also an integer and so F(π) is an integer (in fact, it is easy to see that F(π) = F(0), but that is not relevant to the proof). Since F(0) and F(π) are integers, so is their sum.

Claim 2:

 \int_0^\pi f(x)\sin(x)\,dx=F(0)+F(\pi)

Proof: Since f(2n+2) is the zero polynomial, we have

 F'' + F = f.\,

The derivatives of the sine and cosine function are given by sin' = cos and cos' = sin. Hence the product rule implies

 (F'\cdot\sin - F\cdot\cos)' = f\cdot\sin\!

By the fundamental theorem of calculus

\int_0^\pi f(x)\sin(x)\,dx= \bigl(F'(x)\sin x - F(x)\cos x\bigr)\Big|_{0}^{\pi}.\!

Since sin 0 = sin π = 0 and cos 0 = –cos π = 1 (here we use the above mentioned characterization of π as a zero of the sine function), Claim 2 follows.

Conclusion: Since f(x) > 0 and sin x > 0 for 0 < x < π (because π is the smallest positive zero of the sine function), Claims 1 and 2 show that F(0) + F(π) is a positive integer. Since 0 ≤ x(abx) ≤ πa and 0 ≤ sin x ≤ 1 for 0 ≤ xπ, we have, by the original definition of f,

\int_0^\pi f(x)\sin(x)\,dx\le\pi\frac{(\pi a)^n}{n!}

which is smaller than 1 for large n, hence F(0) + F(π) < 1 for these n, by Claim 2. This is impossible for the positive integer F(0) + F(π).

The above proof is a polished version, which is kept as simple as possible concerning the prerequisites, of an analysis of the formula

\begin{align}
\int_0^\pi f(x)\sin(x)\,dx
&=\sum_{j=0}^n (-1)^j \bigl(f^{(2j)}(\pi)+f^{(2j)}(0)\bigr)\\
&\qquad+(-1)^{n+1}\int_0^\pi f^{(2n+2)}(x)\sin(x)\,dx,
\end{align}

which is obtained by 2n + 2 integrations by parts. Claim 2 essentially establishes this formula, where the use of F hides the iterated integration by parts. The last integral vanishes because f(2n+2) is the zero polynomial. Claim 1 shows that the remaining sum is an integer.

Niven's proof is closer to Cartwright's (and therefore Hermite's) proof than it appears at first sight.[4] In fact,

\begin{align}J_n(x)&=x^{2n+1}\int_{-1}^1(1-z^2)^n\cos(xz)\,dz\\&=\int_{-1}^1\bigl(x^2-(xz)^2\bigr)^nx\cos(xz)\,dz.\end{align}

Therefore, the substitution xz = y turns this integral into

\int_{-x}^x(x^2-y^2)^n\cos(y)\,dy.

In particular,

\begin{align}J_n\left(\frac\pi2\right)&=\int_{-\pi/2}^{\pi/2}\left(\frac{\pi^2}4-y^2\right)^n\cos(y)\,dy\\
&=\int_0^\pi\left(\frac{\pi^2}4-\left(y-\frac\pi2\right)^2\right)^n\cos\left(y-\frac\pi2\right)\,dy\\
&=\int_0^\pi y^n(\pi-y)^n\sin(y)\,dy\\&=\frac{n!}{b^n}\int_0^\pi f(x)\sin(x)\,dx.\end{align}

Another connection between the proofs lies in the fact that Hermite already mentions[5] that if f is a polynomial function and

F=f-f^{(2)}+f^{(4)}\mp\cdots,

then

\int f(x)\sin(x)\,dx=F'(x)\sin(x)-F(x)\cos(x),

from which it follows that

\int_0^\pi f(x)\sin(x)\,dx=F(\pi)+F(0).

Laczkovich's proof

Miklós Laczkovich's proof is a simplification of Lambert's original proof.[10] He considers the functions


\begin{align}
f_k(x) & = 1 - \frac{x^2}k+\frac{x^4}{2! k(k+1)}-\frac{x^6}{3! k(k+1)(k+2)} + \cdots \\
& {} \quad (k\notin\{0,-1,-2,\ldots\}).
\end{align}

These functions are clearly defined for all x  R. Besides

f_{1/2}(x)=\cos(2x)\text{ and }f_{3/2}(x)=\frac{\sin(2x)}{2x}.

Claim 1: The following recurrence relation holds:

 (\forall x\in\mathbb{R}):\frac{x^2}{k(k+1)}f_{k+2}(x)=f_{k+1}(x)-f_k(x).

Proof: This can be proved by comparing the coefficients of the powers of x.

Claim 2: For each x  R, \lim_{k\to+\infty}f_k(x)=1.

Proof: In fact, the sequence x2n/n! is bounded (since it converges to 0) and if C is an upper bound and if k > 1, then

\bigl|f_k(x)-1\bigr|\leqslant\sum_{n=1}^\infty\frac C{k^n}=C\frac{1/k}{1-1/k}=\frac C{k-1}.

Claim 3: If x  0 and if x2 is rational, then

(\forall k\in\mathbb{Q}\setminus\{0,-1,-2,\ldots\}):f_k(x)\neq0\text{ and }\frac{f_{k+1}(x)}{f_k(x)}\notin\mathbb{Q}.

Proof: Otherwise, there would be a number y  0 and integers a and b such that f_k(x)=ay\text{ and }f_{k+1}(x)=by. In order to see why, take y = fk + 1(x), a = 0 and b = 1 if fk(x) = 0; otherwise, choose integers a and b such that fk + 1(x)/fk(x) = b/a and define y = fk(x)/a = fk + 1(x)/b. In each case, y cannot be 0, because otherwise it would follow from claim 1 that each fk + n(x) (n  N) would be 0, which would contradict claim 2. Now, take a natural number c such that all three numbers bc/k, ck/x2 and c/x2 are integers and consider the sequence

g_n=\begin{cases}f_k(x)&\text{ if }n=0\\ \frac{c^n}{k(k+1)\cdots(k+n-1)}f_{k+n}(x)&\text{ otherwise.}\end{cases}

Then

g_0=f_k(x)=ay\in\mathbb{Z}y\text{ and }g_1=\frac ckf_{k+1}(x)=\frac{bc}ky\in\mathbb{Z}y.

On the other hand, it follows from claim 1 that

\begin{align}g_{n+2}&=\frac{c^{n+2}}{x^2k(k+1)\cdots(k+n-1)}\cdot\frac{x^2}{(k+n)(k+n+1)}f_{k+n+2}(x)\\
&=\frac{c^{n+2}}{x^2k(k+1)\cdots(k+n-1)}f_{k+n+1}(x)-\frac{c^{n+2}}{x^2k(k+1)\cdots(k+n-1)}f_{k+n}(x)\\
&=\frac{c(k+n)}{x^2}g_{n+1}-\frac{c^2}{x^2}g_n\\
&=\left(\frac{ck}{x^2}+\frac c{x^2}n\right)g_{n+1}-\frac{c^2}{x^2}g_n,\end{align}

which is a linear combination of g_{n+1} and g_n with integer coefficients. Therefore, each g_n is an integer multiple of y. Besides, it follows from claim 2 that each g_n>0 (and therefore that gn  |y|) if n is large enough and that the sequence of all g_n's converges to 0. But a sequence of numbers greater than or equal to |y| cannot converge to 0.

Since f1/2(π4) = cos(π2) = 0, it follows from claim 3 that π2/16 is irrational and therefore π is irrational.

On the other hand, since

\tan x=\frac{\sin x}{\cos x}=x\frac{f_{3/2}(x/2)}{f_{1/2}(x/2)},

another consequence of claim 3 is that, if x  Q\{0}, then tan x is irrational.

Laczkovich's proof is really about the hypergeometric function. In fact, fk(x) = 0F1(k; x2) and Gauss found a continued fraction expansion of the hypergeometric function using its functional equation.[11] This allowed Laczkovich to find a new and simpler proof of the fact that the tangent function has the continued fraction expansion that Lambert had discovered.

Laczkovich's result can also be expressed in Bessel functions of the first kind J_{\nu}(x). In fact, Γ(k)Jk  1(2x) = xk  1fk(x). So Laczkovich's result is equivalent to: If x  0 and if x2 is rational, then

(\forall k\in\mathbb{Q}\setminus\{0,-1,-2,\ldots\}):\frac{x J_k(x)}{J_{k-1}(x)}\notin\mathbb{Q}.

See also

References

  1. Lindemann, Ferdinand von (2004) [1882], "Ueber die Zahl π", in Berggren, Lennart; Borwein, Jonathan M.; Borwein, Peter B., Pi, a source book (3rd ed.), New York: Springer-Verlag, pp. 194–225, ISBN 0-387-20571-3
  2. Lambert, Johann Heinrich (2004) [1768], "Mémoire sur quelques propriétés remarquables des quantités transcendantes circulaires et logarithmiques", in Berggren, Lennart; Borwein, Jonathan M.; Borwein, Peter B., Pi, a source book (3rd ed.), New York: Springer-Verlag, pp. 129140, ISBN 0-387-20571-3
  3. Zhou, Li; Markov, Lubomir (2010), "Recurrent Proofs of the Irrationality of Certain Trigonometric Values", American Mathematical Monthly 117 (4): 360–362, arXiv:0911.1933, doi:10.4169/000298910x480838
  4. 4.0 4.1 4.2 Zhou, Li (2011), "Irrationality proofs a la Hermite", Math. Gazette (November), arXiv:0911.1929
  5. 5.0 5.1 Hermite, Charles (1873), "Extrait d'une lettre de Monsieur Ch. Hermite à Monsieur Paul Gordan", Journal für die reine und angewandte Mathematik (in French) 76: 303–311
  6. Hermite, Charles (1873), "Extrait d'une lettre de Mr. Ch. Hermite à Mr. Carl Borchardt", Journal für die reine und angewandte Mathematik (in French) 76: 342–344
  7. Hermite, Charles (1912) [1873], "Sur la fonction exponentielle", in Picard, Émile, Œuvres de Charles Hermite (in French) III, Gauthier-Villars, pp. 150–181
  8. Jeffreys, Harold (1973), Scientific Inference (3rd ed.), Cambridge University Press, p. 268, ISBN 0-521-08446-6
  9. Niven, Ivan (1947), "A simple proof that π is irrational" (PDF), Bulletin of the American Mathematical Society 53 (6): 509
  10. Laczkovich, Miklós (1997), "On Lambert's proof of the irrationality of π", American Mathematical Monthly 104 (5): 439–443, JSTOR 2974737
  11. Gauss, Carl Friedrich (1811–1813), "Disquisitiones generales circa seriem infinitam \scriptstyle1+\frac{\alpha\beta}{1.\gamma}x+\frac{\alpha(\alpha+1)\beta(\beta+1)}{1.2.\gamma(\gamma+1)}xx+\frac{\alpha(\alpha+1)(\alpha+2)\beta(\beta+1)(\beta+2)}{1.2.3.\gamma(\gamma+1)(\gamma+1)}x^3+etc", Commentationes Societatis Regiae Scientiarum Gottingensis recentiores (in Latin) 2