Difference between revisions of "Pi"
(→Definition) |
(→The Beukers integral) |
||
(8 intermediate revisions by 7 users not shown) | |||
Line 1: | Line 1: | ||
+ | '''Pi''' is an [[irrational number]] (in fact, [[transcendental number]], as proved by Lindemann in 1882) denoted by the Greek letter <math>\pi </math>. | ||
+ | |||
== Definition == | == Definition == | ||
− | + | Pi is the [[ratio]] of the [[circumference]] ([[perimeter]]) of a given [[circle]] to its [[diameter]]. It is approximately equal to 3.141592653. The number pi is one of the most important [[constant]]s in all of mathematics and appears in some of the most surprising places, such as in the sum <math>\sum_{n=1}^\infty \frac{1}{n^2}=\frac{\pi^2}{6}</math>. Some common [[fraction]]al approximations for pi are <math>\frac{22}{7} \approx 3.14285</math> and <math>\frac{355}{113} \approx 3.1415929</math>. | |
The number pi often shows up in problems in [[number theory]], particularly [[algebraic number theory]]. For example, many [[class number]] formulae involve pi. | The number pi often shows up in problems in [[number theory]], particularly [[algebraic number theory]]. For example, many [[class number]] formulae involve pi. | ||
Line 12: | Line 14: | ||
== Other interesting properties == | == Other interesting properties == | ||
− | *Letting <math>\theta = \pi</math> in | + | *Letting <math>\theta = \pi</math> in [[Euler's identity]] gives <math>e^{\pi i} + 1 = 0</math>, which is considered to be one of the most beautiful results in mathematics since it involves five of the greatest [[constant]]s: [[e]], pi, [[i]], [[unity | 1]], and [[zero (constant)| 0]]. |
+ | |||
+ | ==Proof that <math>\pi</math> is not Liouvillian== | ||
+ | In this section we'll outline the Beukers proof that <math>\pi</math> is not Liouvillian, i.e., that <math>\left|\pi-\frac pq\right|\ge \frac 1{q^M}</math> for large <math>q</math> with some fixed positive <math>M</math>. The best known value of <math>M</math> is currently (June, 1993) <math>M=8.02</math> and is due to Hata. This article can be read right after its parent article [[rational approximation of famous numbers]], but an unexperienced reader may prefer to read the proof that [[ln 2 is not Liouvillian]] first. | ||
+ | |||
+ | ==Beukers' Proof== | ||
+ | ===Preliminaries=== | ||
+ | Note that if <math>x</math> is Liouvillian, then so is any rational multiple of any positive integer power of <math>x</math> (this is a simple exercise we leave to the reader), so, it will suffice to show that <math>\frac{\pi^2}6</math> | ||
+ | is not Liouvillian. | ||
+ | |||
+ | ===[[Euler's identity]]=== | ||
+ | <math>\frac{\pi^2}6=\sum_{m=1}^{\infty}\frac 1{m^2}</math> | ||
+ | ===The Beukers integral=== | ||
+ | Consider the integrals <math>I(k,\ell):=\iint_{(0,1)\times(0,1)}\frac {x^ky^\ell}{1-xy}\,dx\,dy</math>. Expanding <math>\frac 1{1-xy}=\sum_{j=0}^\infty x^jy^j</math> and integrating term by term, we get <math>I(k,\ell)=\sum_{j=0}^{\infty}\frac1{(k+j+1)(\ell+j+1)}.</math> | ||
+ | |||
+ | If <math>k=\ell</math>, we get the sum of all inverse squares starting with <math>\frac 1{(k+1)^2}</math>, so in this case, | ||
+ | |||
+ | <math> | ||
+ | I(k,\ell)=\frac {\pi^2}{6}-\frac 1{1^2}-\frac 1{2^2}-\dots-\frac 1{k^2}. | ||
+ | </math> | ||
+ | |||
+ | On the other hand, if <math>k>\ell</math>, we get | ||
+ | |||
+ | <math>I(k,\ell)=\frac 1{k-\ell}\sum_{j=0}^\infty\left[\frac 1{\ell+j+1}-\frac {1}{k+j+1}\right]=\frac 1{k-\ell}\left[\frac 1{\ell+1}+\frac 1{\ell+2}+\dots+\frac 1k\right]</math> | ||
+ | |||
+ | Similarly, if <math>\ell>k</math>, we get | ||
+ | |||
+ | <math>I(k,\ell)=\frac 1{\ell-k}\left[\frac 1{k+1}+\frac 1{k+2}+\dots+\frac 1\ell\right]</math> | ||
+ | |||
+ | Let now <math>D_n</math> be the [[least common multiple]] of <math>1,2,\dots,n</math>. It follows immediately from the above formulae for <math>I(k,\ell)</math> that for every polynomial <math>R(x,y)</math> with integer coefficients whose [[polynomial|degree]] in both <math>x</math> and <math>y</math> | ||
+ | does not exceed <math>n</math>, we have | ||
+ | |||
+ | <math>D_n^2\iint_{(0,1)\times(0,1)}\frac {R(x,y)}{1-xy}\,dx\,dy=P_n-Q_n\frac{\pi^2}{6}</math> | ||
+ | where <math>P_n,Q_n\in \mathbb Z</math>. | ||
+ | |||
+ | The trick is to choose a polynomial <math>R(x,y)</math> that makes the integral very small. We put | ||
+ | |||
+ | <math>R(x,y)=\frac 1{n!}(1-y)^n\left(\frac d{dx}\right)^n [x^n(1-x)^n].</math> | ||
+ | |||
+ | Note that <math>R(x,y)</math> has integer coefficients as a product of two polynomials with integer coefficients. | ||
+ | |||
+ | ===Estimate of <math>D_n</math>=== | ||
+ | Since the largest possible power of a given [[prime number|prime]] <math>p\le n</math> that can divide one of the numbers <math>1,2,\dots,n</math> is <math>\left\lfloor\frac{\log n}{\log p}\right\rfloor</math> where <math>\lfloor\cdot\rfloor</math> is the [[floor function]], we have <math>D_n=\prod_{p\mathrm{\ prime,}\,p\le n}p^{\left\lfloor\frac{\log n}{\log p}\right\rfloor}\approx n^{\pi(n)}</math> where <math>\pi(n)</math> is the number of primes not exceeding <math>n</math>. According to the [[Prime Number Theorem|prime number theorem]], we have <math>\pi(n)\approx \frac n{\log n}</math>, so <math>D_n\approx e^n.</math> (The exact meaning of the <math>\approx</math> sign here is that, as <math>n\to\infty</math>, the ratio of the <math>n</math>-th roots of the left and the right sides tends to <math>1</math>. Note also that we need the full strength of the prime number theorem here. The elementary Chebyshev's estimate <math>D_n\le n^{\sqrt n}\cdot 4^n</math> is no longer sufficient). | ||
+ | |||
+ | ===Estimates of the integral <math>\iint_{(0,1)\times(0,1)}\frac {R(x,y)}{1-xy}\,dx\,dy</math>=== | ||
+ | Integrating <math>n</math> times by parts with respect to <math>x</math>, we see that our integral equals <math>\iint_{(0,1)\times(0,1)}\left[\frac {x(1-x)y(1-y)}{1-xy}\right]^n\,\frac {dx\,dy}{1-xy}</math>. We want to estimate <math>\max_{0<x,y<1}\frac {x(1-x)y(1-y)}{1-xy}</math>. The first remark is that it is attained on the diagonal <math>x=y</math>. It becomes clear if we rewrite the fraction as <math>\frac {xy(1-x-y+xy)}{1-xy}</math> and recall that for a fixed product <math>xy</math> the minimum of the sum <math>x+y</math> is achieved when <math>x=y</math>. So, we need to find the maximum of <math>\frac {x^2(1-x)^2}{1-x^2}=\frac {x^2(1-x)}{1+x}</math> over the interval <math>(0,1)</math>. Taking the logarithmic derivative and solving the resulting quadratic equation <math>\frac 2x-\frac 1{1-x}-\frac 1{1+x}=0</math>, we see that the maximum is attained at the point <math>x=\frac{\sqrt 5-1}{2}</math> and equals <math>\frac{5\sqrt 5-11}2</math>. It means that the integral does not exceed | ||
+ | <math>\left(\frac{5\sqrt 5-11}2\right)^n I(0,0)\approx \left(\frac{5\sqrt 5-11}2\right)^n</math>. | ||
+ | |||
+ | On the other hand, to estimate the integral from below, it suffices to notice that the integrand is at least <math>\frac 4{81}</math> for <math>\frac 1 3\le x,y\le\frac 2 3</math>. Thus, the integral is not much less than <math>\left(\frac 4{81}\right)^n</math>. | ||
+ | |||
+ | ===A small miracle=== | ||
+ | The miracle that makes the proof work is that <math>e^2\cdot\frac{5\sqrt 5-11}2</math> is strictly less than <math>1</math> (just check it on your calculator or prove it by hand yourself). Thus, for large <math>n</math>, the whole product <math>D_n^2\iint_{(0,1)\times(0,1)}\frac {R(x,y)}{1-xy}\,dx\,dy</math> | ||
+ | can be estimated from above and from below by the <math>n</math>-th powers of some constants strictly between <math>0</math> and <math>1</math>. Using the main theorem from [[rational approximation of famous numbers|the parent article]], we see that it remains only to show that <math>Q_n</math> grow not faster than some geometric progression. | ||
− | == | + | ===Estimate of <math>Q_n</math> and the end of the proof=== |
− | + | Note that <math>Q_n</math> is just minus the product of <math>D_n^2</math>, which grows like <math>e^{2n}</math> and the sum of the coefficients of the polynomial <math>R(x,y)</math> at the products <math>x^ky^k</math> with <math>k=0,\dots,n</math>. It follows that <math>|Q_n|</math> does not exceed the sum of coefficients of <math>(1-y)^n</math>, which is <math>2^n</math>, times the maximal coefficient of <math>\frac 1{n!}\left(\frac d{dx}\right)^n[x^n(1-x)^n]</math>. But even the sum of the absolute values of the coefficients of this polynomial does not exceed <math>8^n</math>. So <math>|Q_n|</math> grows not much faster than <math>(16e^2)^n</math> and we are done. | |
− | |||
[[Category:Constants]] | [[Category:Constants]] |
Latest revision as of 21:11, 16 September 2022
Pi is an irrational number (in fact, transcendental number, as proved by Lindemann in 1882) denoted by the Greek letter .
Contents
Definition
Pi is the ratio of the circumference (perimeter) of a given circle to its diameter. It is approximately equal to 3.141592653. The number pi is one of the most important constants in all of mathematics and appears in some of the most surprising places, such as in the sum . Some common fractional approximations for pi are
and
.
The number pi often shows up in problems in number theory, particularly algebraic number theory. For example, many class number formulae involve pi.
Approximating pi
can be calculated in several ways, and it can also be approximated. One way to approximate
is to inscribe a unit circle in a square of side length 2. Using a computer, random points are placed inside the square. Because the area of the circle is
, and the area of the square is 4, the ratio of the amount of points inside the circle to the total number of points approximates
. This can simply be multiplied by 4 to approximate
.
One formula for pi is . This can be computed to the desired degree of accuracy.
Other interesting properties
- Letting
in Euler's identity gives
, which is considered to be one of the most beautiful results in mathematics since it involves five of the greatest constants: e, pi, i, 1, and 0.
Proof that
is not Liouvillian
In this section we'll outline the Beukers proof that is not Liouvillian, i.e., that
for large
with some fixed positive
. The best known value of
is currently (June, 1993)
and is due to Hata. This article can be read right after its parent article rational approximation of famous numbers, but an unexperienced reader may prefer to read the proof that ln 2 is not Liouvillian first.
Beukers' Proof
Preliminaries
Note that if is Liouvillian, then so is any rational multiple of any positive integer power of
(this is a simple exercise we leave to the reader), so, it will suffice to show that
is not Liouvillian.
Euler's identity
The Beukers integral
Consider the integrals . Expanding
and integrating term by term, we get
If , we get the sum of all inverse squares starting with
, so in this case,
On the other hand, if , we get
Similarly, if , we get
Let now be the least common multiple of
. It follows immediately from the above formulae for
that for every polynomial
with integer coefficients whose degree in both
and
does not exceed
, we have
where
.
The trick is to choose a polynomial that makes the integral very small. We put
Note that has integer coefficients as a product of two polynomials with integer coefficients.
Estimate of ![$D_n$](//latex.artofproblemsolving.com/2/e/a/2eacbcb991af87d495f8689ac0e952b51317dbaa.png)
Since the largest possible power of a given prime that can divide one of the numbers
is
where
is the floor function, we have
where
is the number of primes not exceeding
. According to the prime number theorem, we have
, so
(The exact meaning of the
sign here is that, as
, the ratio of the
-th roots of the left and the right sides tends to
. Note also that we need the full strength of the prime number theorem here. The elementary Chebyshev's estimate
is no longer sufficient).
Estimates of the integral ![$\iint_{(0,1)\times(0,1)}\frac {R(x,y)}{1-xy}\,dx\,dy$](//latex.artofproblemsolving.com/8/e/3/8e3c0a4faa934ca8328a000ccf677c51ff7f6437.png)
Integrating times by parts with respect to
, we see that our integral equals
. We want to estimate
. The first remark is that it is attained on the diagonal
. It becomes clear if we rewrite the fraction as
and recall that for a fixed product
the minimum of the sum
is achieved when
. So, we need to find the maximum of
over the interval
. Taking the logarithmic derivative and solving the resulting quadratic equation
, we see that the maximum is attained at the point
and equals
. It means that the integral does not exceed
.
On the other hand, to estimate the integral from below, it suffices to notice that the integrand is at least for
. Thus, the integral is not much less than
.
A small miracle
The miracle that makes the proof work is that is strictly less than
(just check it on your calculator or prove it by hand yourself). Thus, for large
, the whole product
can be estimated from above and from below by the
-th powers of some constants strictly between
and
. Using the main theorem from the parent article, we see that it remains only to show that
grow not faster than some geometric progression.
Estimate of
and the end of the proof
Note that is just minus the product of
, which grows like
and the sum of the coefficients of the polynomial
at the products
with
. It follows that
does not exceed the sum of coefficients of
, which is
, times the maximal coefficient of
. But even the sum of the absolute values of the coefficients of this polynomial does not exceed
. So
grows not much faster than
and we are done.