Hilbert's sixteenth problem: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
en>Clarkcj12
m Reverted edits by 203.17.203.1 (talk): nonconstructive edits (HG)
en>Materialscientist
m Reverted 1 edit by 119.93.109.57 identified as test/vandalism using STiki
 
Line 1: Line 1:
{{DISPLAYTITLE:Proof that 22/7 exceeds {{pi}}}}
Dalton is what's written by his birth certificate remember, though , he never really preferred that name. His friends say it's unhealthy for him but alternatives he loves doing [http://Www.Wikipedia.org/wiki/typically typically] is acting and he's just lately been doing it for a very long time. In my professional life he is without question a filing assistant and his salary has been very really fulfilling. His wife and your pet live in Idaho having said that he needs to shift because of his spouse. Go to his website identify out more: http://[http://www.adobe.com/cfusion/search/index.cfm?term=&circuspartypanama&loc=en_us&siteSection=home circuspartypanama].com<br><br>my webpage :: [http://circuspartypanama.com clash of clans hack]
{{Pi box}}
[[mathematical proof|Proofs]] of the famous mathematical result that the [[rational number]] 22/7 is greater than [[pi|<span class=texhtml>π</span>]] (pi) date back to antiquity. One of these proofs, more recently developed but requiring only elementary techniques from calculus, has attracted attention in modern mathematics due to its [[mathematical elegance]] and its connections to the theory of [[diophantine approximations]]. Stephen Lucas calls this proof, “One of the more beautiful results related to approximating <span class=texhtml>π</span>”.<ref name="Lucas2005">{{Citation
  | last = Lucas
  | first = Stephen
  | title = Integral proofs that&nbsp;355/113&nbsp;>&nbsp;{{pi}}
  | journal = Australian Mathematical Society Gazette
  | volume = 32
  | issue = 4
  | pages = 263–266
  | year = 2005
  | url = http://www.austms.org.au/Publ/Gazette/2005/Sep05/Lucas.pdf
  | mr = 2176249
  | zbl = 1181.11077}}</ref>
Julian Havil ends a discussion of [[continued fraction]] approximations of <span class=texhtml>π</span> with the result, describing it as “impossible to resist mentioning” in that context.<ref>{{Citation
  | last = Havil
  | first = Julian
  | title = Gamma. Exploring Euler's Constant
  | place = Princeton, NJ
  | publisher = Princeton University Press
  | year = 2003
  | page = 96
  | mr = 1968276
  | zbl = 1023.11001
  | isbn = 0-691-09983-9}}</ref>
 
The purpose of the proof is not primarily to convince its readers that 22/7 {{nowrap|(or 3<sup>1</sup>⁄<sub>7</sub>)}} is indeed bigger than&nbsp;{{pi}}; systematic methods of computing the value of {{pi}} exist.
If one knows that {{pi}} is approximately 3.14159, then it trivially follows that {{pi}}&nbsp;<&nbsp;22/7, which is approximately 3.142857. But it takes much less work to show that {{pi}}&nbsp;<&nbsp;22/7 by the method used in this proof than to show that {{pi}} is approximately&nbsp;3.14159.
 
==Background==
22/7 is a widely used [[Diophantine approximation]] of {{pi}}. It is a convergent in the simple [[continued fraction]] expansion of {{pi}}. It is greater than {{pi}}, as can be readily seen in the [[decimal expansion]]s of these values:
 
:<math>\begin{align}
  \frac{22}{7} & = 3. \overline{142\,857}, \\
  \pi\,        & = 3.141 \,592\,65\ldots
  \end{align}</math>
 
The approximation has been known since antiquity. [[Archimedes]] wrote the first known proof that 22/7 is an overestimate in the 3rd century BCE, although he may not have been the first to use that approximation. His proof proceeds by showing that 22/7 is greater than the ratio of the [[perimeter]] of a [[circumscribe]]d [[regular polygon]] with 96 sides to the diameter of the circle. Another rational approximation of {{pi}} that is far more accurate is [[355/113]].
 
== The proof ==
The proof can be expressed very succinctly:
 
:<math> 0 < \int_0^1 \frac{x^4(1-x)^4}{1+x^2} \, dx = \frac{22}{7} - \pi. </math>
 
Therefore 22/7&nbsp;>&nbsp;{{pi}}.
 
The evaluation of this integral was the first problem in the 1968 [[William Lowell Putnam Mathematical Competition|Putnam Competition]].<ref>{{Citation
  | last = Alexanderson
  | first = Gerald L.
  | last2 = Klosinski
  | first2 = Leonard F.
  | last3 = Larson
  | first3 = Loren C. (editors)
  | title = The William Lowell Putnam Mathematical Competition: Problems and Solutions: 1965–1984
  | place = Washington, D.C.
  | publisher = The Mathematical Association of America
  | year = 1985
  | url = http://books.google.com/books?id=HNLRgSGZrWMC&pg=PA9&dq=December-7-1968+Putnam+Mathematical-Competition&ei=DZCfR4iRJJu4sgPRu-CwCg&sig=u4-SIYyVtV3rbwN-p56c42BGUKw
  | isbn = 0-88385-463-5
  | zbl = 0584.00003}}</ref>
It is easier than most Putnam Competition problems, but the competition often features seemingly obscure problems that turn out to refer to something very familiar. This integral also has been used in the entrance examinations for the [[Indian Institutes of Technology]].<ref>[http://entrance.icbse.com/iit-jee/solutions/2010/fiitjee-1.pdf 2010 IIT Joint Entrance Exam], question 38 on page 15 of the mathematics section.</ref>
 
==Details of evaluation of the integral==
That the [[integral]] is positive follows from the fact that the [[integrand]] is a quotient whose numerator and denominator are both [[non-negative]], being sums or products of powers of non-negative [[real numbers]]. Since the integrand is positive, the integral from 0 to 1 is positive because the lower limit of integration is less than the upper limit of integration.
 
It remains to show that the integral in fact evaluates to the desired quantity:
 
: <math>
\begin{align}
0 & < \int_0^1\frac{x^4(1-x)^4}{1+x^2}\,dx \\[8pt]
& = \int_0^1\frac{x^4-4x^5+6x^6-4x^7+x^8}{1+x^2}\,dx\quad\text{(expansion of terms in the numerator)} \\[8pt]
& = \int_0^1 \left(x^6-4x^5+5x^4-4x^2+4-\frac{4}{1+x^2}\right) \,dx \\
& {} \qquad \text{(polynomial long division)} \\[8pt]
& = \left.\left(\frac{x^7}{7}-\frac{2x^6}{3}+ x^5- \frac{4x^3}{3}+4x-4\arctan{x}\right)\,\right|_0^1 \quad \text{(definite integration)} \\[6pt]
& = \frac{1}{7}-\frac{2}{3}+1-\frac{4}{3}+4-\pi\quad (\text{since }\arctan(1) = \pi/4 \text{ and } \arctan(0) = 0) \\[8pt]
& = \frac{22}{7}-\pi. \quad \text{(addition)}
\end{align}
</math>
 
(See [[polynomial long division]].)
 
==Quick upper and lower bounds==
In {{harvtxt|Dalzell|1944}}, it is pointed out that if 1 is substituted for {{math|''x''}} in the denominator, one gets a lower bound on the integral, and if 0 is substituted for {{math|''x''}} in the denominator, one gets an upper bound:<ref>{{citation|last=Dalzell|first=D. P.|year=1944|title=On 22/7|journal=Journal of the London Mathematical Society|volume=19|pages=133–134|mr=0013425|zbl=0060.15306}}.</ref>
 
:<math>\frac{1}{1260} = \int_0^1\frac{x^4 (1-x)^4}{2}\,dx < \int_0^1\frac{x^4 (1-x)^4}{1+x^2}\,dx < \int_0^1\frac{x^4 (1-x)^4}{1}\,dx = {1 \over 630}.</math>
 
Thus we have
:<math>{22 \over 7} - {1 \over 630} < \pi < {22 \over 7} - {1 \over 1260},</math>
 
hence 3.1412 < {{pi}} < 3.1421 in decimal expansion. The bounds deviate by less than 0.015% from&nbsp;{{pi}}. See also {{harvtxt|Dalzell|1971}}.<ref>{{citation|last=Dalzell|first=D. P.|year=1971|title=On 22/7 and 355/113|journal=Eureka; the Archimedeans' Journal|volume=34|pages=10–13|issn=0071-2248}}.</ref>
 
==Proof that 355/113 exceeds {{pi}}==
As discussed in {{harvtxt|Lucas|2005}}, the well-known Diophantine approximation and far better upper estimate [[355/113]] for {{pi}} follows from the relation
 
:<math>0<\int_0^1\frac{x^8(1-x)^8(25+816x^2)}{3164(1+x^2)}\,dx=\frac{355}{113}-\pi.</math>
 
Note that
 
:<math>\frac{355}{113}= 3.141\,592\,92\ldots,</math>
 
where the first six digits after the period agree with those of {{pi}}. Substituting 1 for {{math|''x''}} in the denominator, we get the lower bound
 
:<math>\int_0^1\frac{x^8(1-x)^8(25+816x^2)}{6328}\,dx =\frac{911}{5\,261\,111\,856} = 0.000\,000\,173\ldots,</math>
 
substituting 0 for {{math|''x''}} in the denominator, we get twice this value as an upper bound, hence
 
:<math>\frac{355}{113}-\frac{911}{2\,630\,555\,928}<\pi<\frac{355}{113}-\frac{911}{5\,261\,111\,856}\,.</math>
 
In decimal expansion, this means
'''3.141'''&thinsp;'''592'''&thinsp;57 < {{pi}} < '''3.141'''&thinsp;'''592'''&thinsp;74, where the bold digits of the lower and upper bound are those of&nbsp;{{pi}}.
 
==Extensions==
The above ideas can be generalized to get better approximations of&nbsp;{{pi}}, see also {{harvtxt|Backhouse|1995}}<ref>{{citation|last=Backhouse|first=Nigel|date=July 1995|title=Note 79.36, Pancake functions and approximations to {{pi}}|journal=The Mathematical Gazette|volume=79|issue=485|pages=371–374|jstor=3618318}}</ref> and {{harvtxt|Lucas|2005}} (in both references, however, no calculations are given). For explicit calculations, consider, for every integer {{math|''n'' ≥ 1}},
 
:<math>
\frac1{2^{2n-1}}\int_0^1 x^{4n}(1-x)^{4n}\,dx
<\frac1{2^{2n-2}}\int_0^1\frac{x^{4n}(1-x)^{4n}}{1+x^2}\,dx
<\frac1{2^{2n-2}}\int_0^1 x^{4n}(1-x)^{4n}\,dx,
</math>
 
where the middle integral evaluates to
 
:<math>\begin{align}
&\frac1{2^{2n-2}}\int_0^1\frac{x^{4n}(1-x)^{4n}}{1+x^2}\,dx\\
&\qquad=\sum_{j=0}^{2n-1}\frac{(-1)^j}{2^{2n-j-2}(8n-j-1)\binom{8n-j-2}{4n+j}}
+(-1)^n\biggl(\pi-4\sum_{j=0}^{3n-1}\frac{(-1)^j}{2j+1}\biggr)
\end{align}</math>
 
involving&nbsp;{{pi}}. The last sum also appears in [[Leibniz formula for pi|Leibniz' formula for {{pi}}]]. The correction term and [[error bound]] is given by
 
:<math>\begin{align}\frac1{2^{2n-1}}\int_0^1 x^{4n}(1-x)^{4n}\,dx
&=\frac{1}{2^{2n-1}(8n+1)\binom{8n}{4n}}\\
&\sim\frac{\sqrt{\pi n}}{2^{10n-2}(8n+1)},
\end{align}</math>
 
where the approximation (the tilde means that the quotient of both sides tends to one for large {{math|''n''}}) of the [[central binomial coefficient]] follows from [[Stirling's formula]] and shows the fast convergence of the integrals to&nbsp;{{pi}}.
 
<div style="clear:both;width:95%;" class="NavFrame">
<div class="NavHead" style="background-color:#FFFAF0; text-align:left; font-size:larger;">Calculation of these integrals</div>
<div class="NavContent" style="text-align:left;display:none;">
For all integers {{math|''k'' ≥ 0}} and {{math|''l'' ≥ 2}} we have
 
:<math>\begin{align}
x^k(1-x)^l&=(1-2x+x^2)x^k(1-x)^{l-2}\\
&=(1+x^2)\,x^k(1-x)^{l-2}-2x^{k+1}(1-x)^{l-2}.
\end{align}</math>
 
Applying this formula recursively {{math|2''n''}} times yields
 
:<math>x^{4n}(1-x)^{4n}
=(1+x^2)\sum_{j=0}^{2n-1}(-2)^jx^{4n+j}(1-x)^{4n-2(j+1)}+(-2)^{2n}x^{6n}.</math>
 
Furthermore,
 
:<math>\begin{align}
x^{6n}-(-1)^{3n}
&=\sum_{j=1}^{3n}(-1)^{3n-j}x^{2j}-\sum_{j=0}^{3n-1}(-1)^{3n-j}x^{2j}\\
&=\sum_{j=0}^{3n-1}\bigl((-1)^{3n-(j+1)}x^{2(j+1)}-(-1)^{3n-j}x^{2j}\bigr)\\
&=-(1+x^2)\sum_{j=0}^{3n-1}(-1)^{3n-j}x^{2j},\\
\end{align}</math>
 
where the first equality holds, because the terms for {{math|1 ≤ ''j'' ≤ 3''n'' &minus; 1}} cancel, and the second equality arises from the index shift {{math|''j'' → ''j'' + 1}} in the first sum.
 
Application of these two results gives
 
:<math>\begin{align}\frac{x^{4n}(1-x)^{4n}}{2^{2n-2}(1+x^2)}
&=\sum_{j=0}^{2n-1}\frac{(-1)^j}{2^{2n-j-2}}x^{4n+j}(1-x)^{4n-2j-2}\\
&\qquad{}-4\sum_{j=0}^{3n-1}(-1)^{3n-j}x^{2j}+(-1)^{3n}\frac4{1+x^2}.\qquad(*)
\end{align}</math>
 
For integers {{math|''k'', ''l'' ≥ 0}}, using [[integration by parts]] {{math|''l''}} times, we obtain
 
:<math>\begin{align}
\int_0^1x^k(1-x)^l\,dx
&=\frac l{k+1}\int_0^1x^{k+1}(1-x)^{l-1}\,dx\\
&=\cdots\\
&=\frac l{k+1} \frac{l-1}{k+2}\cdots\frac1{k+l}\int_0^1x^{k+l}\,dx\\
&=\frac{1}{(k+l+1)\binom{k+l}{k}}.\qquad(**)
\end{align}</math>
 
Setting {{math|''k'' {{=}} ''l'' {{=}} 4''n''}}, we obtain
 
:<math>\int_0^1 x^{4n}(1-x)^{4n}\,dx
=\frac{1}{(8n+1)\binom{8n}{4n}}.</math>
 
Integrating (*) from 0 to 1 using (**) and {{math|arctan(1) {{=}} π/4}}, we get the claimed equation involving&nbsp;{{pi}}.
</div>
</div>
 
The results for {{math|''n'' {{=}} 1}} are given above. For {{math|''n'' {{=}} 2}} we get
 
:<math>\frac14\int_0^1\frac{x^8(1-x)^8}{1+x^2}\,dx=\pi -\frac{47\,171}{15\,015}</math>
 
and
 
:<math>\frac18\int_0^1 x^8(1-x)^8\,dx=\frac1{1\,750\,320},</math>
 
hence  '''3.141'''&thinsp;'''592'''&thinsp;31 < {{pi}} < '''3.141'''&thinsp;'''592'''&thinsp;89, where the bold digits of the lower and upper bound are those of&nbsp;{{pi}}. Similarly for {{math|''n'' {{=}} 3}},
 
:<math>\frac1{16}\int_0^1\frac{x^{12}(1-x)^{12}}{1+x^2}\,dx= \frac{431\,302\,721}{137\,287\,920}-\pi</math>
 
with correction term and error bound
 
:<math>\frac1{32}\int_0^1 x^{12}(1-x)^{12}\,dx=\frac1{2\,163\,324\,800},</math>
 
hence '''3.141'''&thinsp;'''592'''&thinsp;'''653'''&thinsp;40 < {{pi}} < '''3.141'''&thinsp;'''592'''&thinsp;'''653'''&thinsp;87. The next step for {{math|''n'' {{=}} 4}} is
 
:<math>\frac1{64}\int_0^1\frac{x^{16}(1-x)^{16}}{1+x^2}\,dx= \pi-\frac{741\,269\,838\,109}{235\,953\,517\,800}</math>
 
with
 
:<math>\frac1{128}\int_0^1 x^{16}(1-x)^{16}\,dx=\frac1{2\,538\,963\,567\,360},</math>
 
which gives '''3.141'''&thinsp;'''592'''&thinsp;'''653'''&thinsp;'''589'''&thinsp;55 < {{pi}} < '''3.141'''&thinsp;'''592'''&thinsp;'''653'''&thinsp;'''589'''&thinsp;96.
 
==See also==
*[[Approximations of π|Approximations of {{pi}}]]
*[[Chronology of computation of π|Chronology of computation of {{pi}}]]
*[[Proof that π is irrational|Proof that {{pi}} is irrational]]
*[[Lindemann–Weierstrass theorem]] (proof that {{pi}} is transcendental)
*[[List of topics related to π|List of topics related to {{pi}}]]
 
==References==
{{reflist|2}}
 
==External links==
* [http://web.archive.org/web/20070730150749/http://www.kalva.demon.co.uk/putnam/putn68.html The problems of the 1968 Putnam competition], with this proof listed as question A1.
* [http://www.cecm.sfu.ca/%7Ejborwein/pi-slides.pdf The Life of Pi] by Jonathan Borwein&mdash;see page 5 for this integral.
 
{{DEFAULTSORT:Proof That 22 7 Exceeds}}
[[Category:Pi]]
[[Category:Article proofs]]

Latest revision as of 12:00, 29 August 2014

Dalton is what's written by his birth certificate remember, though , he never really preferred that name. His friends say it's unhealthy for him but alternatives he loves doing typically is acting and he's just lately been doing it for a very long time. In my professional life he is without question a filing assistant and his salary has been very really fulfilling. His wife and your pet live in Idaho having said that he needs to shift because of his spouse. Go to his website identify out more: http://circuspartypanama.com

my webpage :: clash of clans hack