Votre navigateur n'accepte pas le Javascript.La navigation sur ce site risque de ne pas fonctionner correctement.

253. La formule de Wallis

Vous vous intéressez dans cet article à obtenir le nombre $\pi$ comme limite d’une suite.

Utilisez les intégrales de Wallis

Pour tout entier naturel $n$, vous posez $I_n = \int_0^{\pi/2} \sin^n t\dt.$

Vous calculez d’abord $I_0.$

\begin{align*}
I_0 &= \int_0^{\pi/2} \sin^0 t\dt\\
&= \int_0^{\pi/2} 1\dt\\
&=\frac{\pi}{2}.
\end{align*}

Vous enchaînez avec $I_1.$

\begin{align*}
I_1 &= \int_0^{\pi/2} \sin^1 t\dt\\
&= \int_0^{\pi/2} \sin t\dt\\
&=\left[-\cos t\right]_0^{\pi/2}\\
&=-\left[\cos t\right]_0^{\pi/2}\\
&=-\left(0-1\right)\\
&=1.
\end{align*}

Soit maintenant $n$ un entier naturel fixé. Vous allez utiliser un outil très important, qui est une intégration par parties.

\begin{align*}
I_{n+2} &= \int_0^{\pi/2} \sin^{n+2} t\dt\\
&=  \int_0^{\pi/2}\sin t \times \sin^{n+1} t\dt\\
&= \left[-\cos t \times \sin^{n+1} t\right]_0^{\pi/2}-  \int_0^{\pi/2} (-\cos t)  (n+1) \sin^{n} t \cos t\dt\\
&=0-  \int_0^{\pi/2} (-\cos t)  (n+1) \sin^{n} t \cos t\dt\\
&=   (n+1) \int_0^{\pi/2} \cos^2 t  \sin^{n} t \dt\\
&=   (n+1) \int_0^{\pi/2} (1-\sin^2 t)  \sin^{n} t \dt\\
&=(n+1) \int_0^{\pi/2} (\sin^n t-\sin^{n+2} t)  \dt\\
&=(n+1)I_n - (n+1)I_{n+2}.
\end{align*}

Cette relation permet d’obtenir :

\begin{align*}
(n+2)I_{n+2} &= (n+1)I_n.
\end{align*}

Ainsi :

\boxed{\forall n\in\N, I_{n+2} =\frac{n+1}{n+2} I_n.}

Le reste dépend de la parité de l’entier $n.$

Pour tout $n\in\N$ calculez $I_{2n}$

Déjà, vous avez $I_0 = \frac{\pi}{2}.$

Puis $I_2 = \frac{1}{2}I_0 = \frac{1}{2}\times \frac{\pi}{2}.$

Ensuite $I_4 =\frac{3}{4}I_2 = \frac{1\times 3}{2\times 4} \times \frac{\pi}{2}.$

Il est d’usage d’écrire ce dernier calcul en utilisant des factorielles :

\begin{align*}
I_4 &= \frac{1\times 3}{2\times 4} \times \frac{\pi}{2}\\
&= \frac{1\times 2 \times 3\times 4}{(2\times 4)^2} \times \frac{\pi}{2}\\
&= \frac{4 !}{(2\times 4)^2} \times \frac{\pi}{2}\\
&= \frac{4 !}{(1\times 2 \times 2^2)^2} \times \frac{\pi}{2}\\
&= \frac{4 !}{(2 ! \times 2^2)^2} \times \frac{\pi}{2}\\
&= \frac{4 !}{(2 !)^2 \times 2^4} \times \frac{\pi}{2}.
\end{align*}

Qu’en est-il pour $I_6$ ? Vous effectuez une démarche similaire :

\begin{align*}
I_6 &= \frac{5}{6}\times I_4 \\
&= \frac{5}{6}\times \frac{1\times 3}{2\times 4} \times \frac{\pi}{2}\\
&= \frac{1\times 3\times 5}{2\times 4 \times 6} \times \frac{\pi}{2}\\
&= \frac{1\times 2\times 3\times 4\times 5\times 6}{(2\times 4 \times 6)^2} \times \frac{\pi}{2}\\
&= \frac{6 !}{(2\times 4 \times 6)^2} \times \frac{\pi}{2}\\
&= \frac{6 !}{(1\times 2 \times 3 \times 2^3)^2} \times \frac{\pi}{2}\\
&= \frac{6 !}{(3 ! \times 2^3)^2} \times \frac{\pi}{2}\\
&= \frac{6 !}{(3 !)^2 \times 2^6} \times \frac{\pi}{2}.
\end{align*}

Vous êtes maintenant prêt à démontrer par récurrence que $\forall n\in\N, I_{2n} = \frac{(2n) !}{(n !)^2\times 2^{2n}}\times \frac{\pi}{2}.$

Pour tout entier naturel $n$, vous notez $\mathscr{P}(n)$ la propriété : « $I_{2n} = \frac{(2n) !}{(n !)^2\times 2^{2n}}\times \frac{\pi}{2}.$ »

Initialisation. Pour $n=0$, $I_0 = \frac{\pi}{2}.$

D’autre part, $\frac{(0 !)}{(0 !)^2\times 2^{0}}\times \frac{\pi}{2} = \frac{\pi}{2}$ donc $\mathscr{P}(0)$ est vérifiée.

Hérédité. Soit $n\in\N$, supposez $\mathscr{P}(n).$

\begin{align*}
I_{2n+2} &= \frac{2n+1}{2n+2} I_{2n}\\
&= \frac{2n+1}{2n+2} \times \frac{(2n) !}{(n !)^2\times 2^{2n}}\times \frac{\pi}{2}\\
&= \frac{(2n+1)(2n+2)}{(2n+2)^2} \times \frac{(2n) !}{(n !)^2\times 2^{2n}}\times \frac{\pi}{2}\\
&= \frac{(2n+2) !}{(2n+2)^2 (n !)^2\times 2^{2n}}\times \frac{\pi}{2}\\
&= \frac{(2n+2) !}{(2(n+1))^2 (n !)^2\times 2^{2n}}\times \frac{\pi}{2}\\
&= \frac{(2n+2) !}{2^2 \times (n+1)^2 (n !)^2\times 2^{2n}}\times \frac{\pi}{2}\\
&= \frac{(2n+2) !}{ (n+1)^2 (n !)^2\times 2^{2n+2}}\times \frac{\pi}{2}\\
&= \frac{(2n+2) !}{ ((n+1)\times n !)^2\times 2^{2n+2}}\times \frac{\pi}{2}\\
&= \frac{(2n+2) !}{((n+1) !)^2\times 2^{2n+2}}\times \frac{\pi}{2}.
\end{align*}

Ainsi, la propriété $\mathscr{P}(n+1)$ est bien vérifiée.

Vous avez donc établi le résultat suivant :

\boxed{\forall n\in\N, I_{2n} = \frac{(2n) !}{(n !)^2\times 2^{2n}}\times \frac{\pi}{2}.}

Pour tout $n\in\N$ calculez $I_{2n+1}$

Déjà, vous avez $I_1 = 1.$

Puis $I_3 = \frac{2}{3}I_1 = \frac{2}{3}.$

Ensuite, $I_5 = \frac{4}{5}I_3 = \frac{4}{5}\times \frac{2}{3} = \frac{2\times 4}{3\times 5}.$

Vous écrivez $I_5$ avec des factorielles :

\begin{align*}
I_5 &=  \frac{2\times 4}{3\times 5} \\
&=   \frac{(2\times 4)^2}{2\times 3 \times 4\times 5} \\
&=   \frac{(1\times 2 \times 2^2)^2 \times 2^2}{5 !} \\
&=   \frac{(2 !)^2 \times 2^4}{5 !}.
\end{align*}

Pour tout entier naturel $n$, vous notez $\mathscr{P}(n)$ la propriété : « $I_{2n+1} = \frac{(n !)^2\times 2^{2n}}{(2n+1) !}.$ »

Initialisation. Pour $n=0$, $I_1 = 1.$

$ \frac{(0 !)^2\times 2^0}{(2\times 0+1) !} = 1$ donc $\mathscr{P}(0)$ est vérifiée.

Hérédité. Soit $n\in\N$, supposez $\mathscr{P}(n).$

\begin{align*}
I_{2n+3} &= \frac{2n+2}{2n+3} I_{2n+1}\\
&=\frac{2n+2}{2n+3} \times \frac{(n !)^2\times 2^{2n}}{(2n+1) !}\\
&=\frac{(2n+2)^2}{(2n+3)(2n+2)} \times \frac{(n !)^2\times 2^{2n}}{(2n+1) !}\\
&=\frac{2^2(n+1)^2}{(2n+3)(2n+2)} \times \frac{(n !)^2\times 2^{2n}}{(2n+1) !}\\
&=\frac{2^2(n+1)^2}{(2n+3)(2n+2)} \times \frac{(n !)^2\times 2^{2n}}{(2n+1) !}\\
&=\frac{(n+1)^2(n !)^2\times 2^{2n+2}}{(2n+1) !(2n+2)(2n+3)}\\
&=\frac{(n+1)^2(n !)^2\times 2^{2n+2}}{(2n+3) !}\\
&=\frac{((n+1)n !)^2\times 2^{2n+2}}{(2n+3) !}\\
&=\frac{((n+1) !)^2\times 2^{2n+2}}{(2n+3) !}.
\end{align*}

Ainsi, la propriété $\mathscr{P}(n+1)$ est bien vérifiée.

Vous avez donc établi le résultat suivant :

\boxed{\forall n\in\N, I_{2n+1} = \frac{(n !)^2\times 2^{2n}}{(2n+1) !}.}

La suite $(I_n)$ est décroissante

Soit $n$ un entier naturel.

Sur l’intervalle $[0,\pi/2]$ la fonction $t\mapsto \sin t$ est positive.

Il en résulte que :

\forall t\in\left[0,\pi/2\right], 0\leq (\sin t)^n.

D’autre part :

\forall t\in \left[0,\pi/2\right], \sin t \leq 1.

En multipliant cette inégalité par $(\sin t)^n$ qui est positif, vous obtenez :

\forall t\in \left[0,\pi/2\right], 0\leq (\sin t)^{n+1} \leq (\sin t)^n.

En intégrant sur l’intervalle $[0,\pi/2]$ il vient :

\begin{align*}
\int_0^{\pi/2} \sin^{n+1} t\dt \leq \int_0^{\pi/2} \sin^{n} t\dt.
\end{align*}

En définitive vous avez montré que :

\boxed{\forall n\in\N, I_{n+1}\leq I_n.}

Déduisez-en la formule de Wallis

Soit $n$ un entier naturel non nul.

Vous avez successivement :

\begin{align*}
I_{2n} &= \frac{(2n) !}{(n !)^2\times 2^{2n}}\times \frac{\pi}{2}\\
I_{2n+1} &= \frac{(n !)^2\times 2^{2n}}{(2n+1) !}\\
\frac{I_{2n}}{I_{2n+1}} &=  \frac{(2n) !}{(n !)^2\times 2^{2n}}\times \frac{\pi}{2} \times \frac{(2n+1) !}{(n !)^2\times 2^{2n}}\\
&=  \frac{(2n) !(2n+1) !}{(n !)^4\times 2^{4n+1}} \times \pi.
\end{align*}

D’autre part, par décroissance de la suite $(I_p)_{p\geq 0}$ vous avez :

\begin{align*}
I_{2n+1}\leq I_{2n}\leq I_{2n-1}.
\end{align*}

Des expressions explicites de la suite susmentionnée vous avez $\forall p\in\N, I_p > 0.$

En divisant par $I_{2n+1}$ il vient :

\begin{align*}
1\leq \frac{I_{2n}}{I_{2n+1}}\leq \frac{I_{2n-1}}{I_{2n+1}}.
\end{align*}

Or :

\begin{align*}
I_{2n+1} &=\frac{2n}{2n+1} I_{2n-1} \\
\frac{2n+1}{2n} &=\frac{ I_{2n-1}}{I_{2n+1}} \\
1+\frac{1}{2n} &=\frac{ I_{2n-1}}{I_{2n+1}}.
\end{align*}

Donc :

1\leq  \frac{(2n) !(2n+1) !}{(n !)^4\times 2^{4n+1}} \times \pi \leq 1+\frac{1}{2n}.

Comme $\lim_{n\to +\infty} 1+\frac{1}{2n} =1 $ l’application du théorème des gendarmes fournit :

$\lim_{n\to +\infty} \frac{(2n) !(2n+1) !}{(n !)^4\times 2^{4n+1}} \times \pi = 1.$

En divisant par $\pi$ vous déduisez alors que :

$\lim_{n\to +\infty} \frac{(2n) !(2n+1) !}{(n !)^4\times 2^{4n+1}} = \frac{1}{\pi}.$

En prenant l’inverse vous obtenez la formule de Wallis :

\boxed{\lim_{n\to +\infty} \frac{(n !)^4\times 2^{4n+1}}{(2n) !(2n+1) !} =\pi.}

Application avec un ordinateur. Pour $n=493$ vous trouverez l’encadrement suivant $3,140<\pi<3,144.$ La formule de Wallis n’est en pratique pas adaptée pour calculer efficacement les décimales de $\pi.$

Prolongement

En utilisant la formule de Wallis et en admettant qu’il existe une constante strictement positive $C$ telle que :

n !   \underset{+\infty}{\sim} C\sqrt{n}\left(\frac{n}{\e}\right)^n,

pourriez-vous calculer la constante $C$ ?

Partagez!

Diffusez cet article auprès de vos connaissances susceptibles d'être concernées en utilisant les boutons de partage ci-dessous.

Aidez-moi sur Facebook!

Vous appréciez cet article et souhaitez témoigner du temps que j'y ai passé pour le mettre en œuvre. C'est rapide à faire pour vous et c'est important pour moi, déposez un j'aime sur ma page Facebook. Je vous en remercie par avance.

Lisez d'autres articles!

Parcourez tous les articles qui ont été rédigés. Vous en trouverez sûrement un qui vous plaira!