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

296. Déterminez sans division si un nombre est premier (2/2)

Soit $n$ un nombre impair supérieur ou égal à $7$ et qui n’est pas un carré.

L’objectif de cet article est de démontrer que la méthode suivie dans le contenu rédigé dans l'article 295 fonctionne bien.

Elle fournit toujours une réponse en un nombre fini d’étapes et permet de conclure que :

  • soit le nombre $n$ est premier :
  • soit il ne l’est pas – on dit qu’il est composé – et il vous permet d’écrire $n$ comme le produit de deux entiers explicites qui sont différents de $1$ et de $n.$

Le point de départ

Vous considérez l’ensemble :

A = \{m\in\N, m^2>n\}.

Comme l’entier $n$ est supérieur ou égal à $7$, l’entier $n-1$ est supérieur ou égal à $6.$

Par produit, vous déduisez :

n(n-1)\geq 42.

Vous développez et isolez le carré de $n$ :

\begin{align*}
n^2-n&\geq 42\\
n^2&\geq n+42\\
n^2&\geq n+1\\
n^2&>n.
\end{align*}

Vous avez obtenu l’appartenance de $n$ à l’ensemble $A$, donc $A$ est non vide.

$A$ étant une partie de $\N$ non vide, vous déduisez que $A$ admet un plus petit élément.

Dans la suite, vous noterez :

\boxed{k=\mathrm{Min}\  \{m\in\N, m^2>n\}.}

Vous aurez aussi besoin du nombre $K=\frac{n+1}{2}.$ Comme $n$ est impair supérieur ou égal à $7$, $K$ est un entier supérieur ou égal à $4$ :

\boxed{K = \frac{n+1}{2}.}

Cas de la factorisation avec $1$ et $n$ et ses conséquences

La factorisation suivante de $n$ semble n’avoir que peu d’intérêt :

n = n\times 1.

Ceci n’est qu’un leurre.

Vous allez maintenant utiliser les identités remarquables bien connues des Babyloniens, qui ont formé une ancienne civilisation de la Mésopotamie. Ils ont contribué au développement des mathématiques et avaient des connaissances sur les identités remarquables.

Comme :

\begin{align*}
(n+1)^2&=n^2+2n+1\\
(n-1)^2&=n^2-2n+1\\
\end{align*}

il vient par soustraction :

(n+1)^2-(n-1)^2 = 4n.

Or, $4$ est lui-même un carré, donc :

\begin{align*}
n &= \frac{(n+1)^2}{4}-\frac{(n-1)^2}{4} \\
&=\left(\frac{n+1}{2}\right)^2-\left(\frac{n-1}{2}\right)^2\\
&=K^2-\left(\frac{n-1}{2}\right)^2.
\end{align*}

De ce résultat, vous déduisez que $K^2 – n$ est un carré et que ce carré vaut $\left(\frac{n-1}{2}\right)^2.$ Comme $n\geq 7$ il vient $n-1\geq 6$ donc $\frac{n-1}{2}\geq 3$ et par suite $\left(\frac{n-1}{2}\right)^2 \geq 9.$

Donc $K^2 – n\geq 9$ et $K^2 > n.$

Par suite, $K\in A$ et par définition du minimum de $A$, il vient :

\boxed{k\leq K.}

Étude du cas d’égalité

Supposez un instant que :

k = K.

Alors $K=\frac{n+1}{2}$ est le minimum de $A.$

Or, comme $n$ est impair supérieur ou égal à $7$, $\frac{n-1}{2}$ est un entier supérieur ou égal à $3$. Comme il est strictement inférieur à $k$, vous déduisez que $\frac{n-1}{2}\notin A.$ Donc :

\begin{align*}
\left(\frac{n-1}{2}\right)^2&\leq n\\
\frac{(n-1)^2}{4}&\leq n\\
(n-1)^2&\leq 4n\\
n^2-2n+1&\leq 4n\\
n^2-6n&\leq -1\\
n^2-6n+9&\leq 8\\
(n-3)^2 &< 9\\
n-3&<3\\
n&<6.
\end{align*}

Ceci est absurde puisque $n\geq 7.$ Il est donc impossible d’avoir $k=K$ et par suite :

\boxed{k < K.}

Le test par la méthode de Fermat

De ce qui précède, l’ensemble $B$ défini par :

B = \{\ell^2-n, \ell\in\llbracket k, K-1\rrbracket\}.

est non vide.

Deux cas se produisent.

Cas 1 : l’ensemble $B$ contient un carré

Sous cette hypothèse, il existe un entier $\ell \in \llbracket k, K-1\rrbracket$ et un entier naturel $m$ tels que :

\ell^2-n  = m^2.

Alors :

\begin{align*}
n &= \ell^2-m^2\\
&=(\ell+m)(\ell-m).
\end{align*}

D’une part, $\ell$ et $m$ sont positifs donc la somme $\ell +m$ est positive.

Comme $n\in\N^{*}$ est un multiple de $\ell+m$, la somme $\ell + m$ est strictement positive donc $\ell+m \geq 1.$

Il en résulte par la règle des signes que la différence $\ell-m$ est strictement positive aussi.

Supposez un instant que $\ell -m = 1.$ Alors : $m = \ell-1.$ Du coup :

\begin{align*}
n &= \ell^2-(\ell-1)^2\\
&=\ell^2-(\ell^2-2\ell+1)\\
&=2\ell-1.
\end{align*}

Ainsi :

\begin{align*}
\ell &= \frac{n+1}{2}\\
&=K.
\end{align*}

Cela est absurde, puisque $\ell\leq K-1.$

Il en résulte que $\ell-m$ est un entier supérieur ou égal à $2.$

Le produit :

\begin{align*}
n =(\ell+m)(\ell-m)
\end{align*}

montre que $\ell -m \leq n.$

En effet si tel n’était pas le cas, vous auriez :

\ell-m > n\\
\ell+m \geq 1

et par produit l’inégalité serait stricte :

(\ell+m)(\ell-m)>n

ce qui serait absurde.

Supposez enfin que $\ell-m = n.$

Alors nécessairement $\ell+m = 1.$

Or $\ell^2\geq k^2>n\geq 7$ donc $\ell^2 > 4$ donc $\ell > 2.$ La positivité de $m$ fournit $\ell + m >2$ ce qui contredit $\ell+m=1.$

Il a été montre que $\ell-m \in\llbracket 2, n-1\rrbracket.$ On dit que $\ell-m$ est un diviseur propre de $n.$ En effet, $n$ admet un troisième diviseur qui est différent de $1$ et de $n$.

Par suite, $n$ n’est pas premier et vous avez une factorisation par ce diviseur propre :

\begin{align*}
n =(\ell+m)(\ell-m).
\end{align*}

Cas 2 : l’ensemble $B$ ne contient pas de carré

Vous supposez, en raisonnant par l’absurde, que $n$ n’est pas un nombre premier.

Il existe un entier $d$ compris entre $2$ et $n-1$ tel que $n$ soit un multiple de $d.$

Du coup il existe aussi un entier naturel $d’$ tel que :

n = dd'.

Supposez $d’\in\{1, n\}.$ Alors $d\in\{1,n\}$ ce qui est absurde, donc $d’\in \llbracket 2, n-1\rrbracket.$

Si $d$ était pair, alors il existerait un entier naturel $d_0$ tel que $d = 2d_0$ et par suite :

n=2d_0d'.

Comme $d_0d’\in\N$ vous déduisez que $n$ est pair, ce qui est absurde.

De même, si $d’$ était pair, alors $n$ serait pair, ce qui est absurde.

Donc $d$ et $d’$ sont deux entiers impairs appartenant à l’ensemble $\llbracket 2, n-1\rrbracket.$

Le reste de la division euclidienne de $d$ par $2$ ne pouvant être nul, c’est donc qu’il est égal à $1.$ Le résultat est aussi identique pour $d’.$

Ainsi, il existe deux entiers naturels $u$ et $u’$ tels que :

\begin{align*}
d&=2u+1\\
d'&=2u'+1.
\end{align*}

Comme $d\geq 2$, vous déduisez $2u+1\geq 2$ donc $2u\geq 1$ et donc $u\neq 0$ donc $u\geq 1.$

De même, vous avez aussi $u’\geq 1.$

Utilisant à nouveau $n = dd’$ il vient :

\begin{align*}
n&=(2u+1)(2u'+1)\\
&=4uu'+2u+2u'+1\\
&=u^2+u'^2+2uu'+2u+2u'+1-(u^2+u'^2)+2uu'\\
&=(u+u'+1)^2-(u^2+u'^2-2uu')\\
&=(u+u'+1)^2-(u-u')^2.
\end{align*}

Posez $\ell = u+u’+1.$ Vous avez $\ell \in\N$ et $\ell \geq 3.$

Si vous aviez $u=u’$ alors $n$ serait le carré de $\ell$ ce qui est absurde.

Comme $u\neq u’$ vous déduisez $(u-u’)^2\geq 1$ et donc $(u+u’+1)^2-n \geq 1$ donc $(u+u’+1)^2>n$ et par suite $u+u’+1\in A.$ Cela prouve que :

\begin{align*}
k&\leq u+u'+1\\
k&\leq \ell.
\end{align*}

D’autre part, de l’égalité $n = 4uu’+2u+2u’+1$ et en minorant le produit $uu’$ par $1$, vous déduisez :

\begin{align*}
n&\geq 4uu'+2u+2u'+1 
\\
&\geq 4+2u+2u'+1
\\
&\geq 5+2u+2u'.
\end{align*}

Or $n=2K-1$ donc :

\begin{align*}
2K-1&\geq 5+2u+2u'\\
2K&\geq 6+2u+2u'\\
K&\geq 3+u+u'\\
K-1&\geq 2+u+u'\\
K-1&\geq 1+\ell\\
K-1&\geq \ell.
\end{align*}

Vous avez $\ell\in \llbracket k, K-1\rrbracket$ et d’autre part :

\ell^2-n = (u-u')^2.

Donc $B$ contient un carré, ce qui est absurde.

Ainsi, $n$ est premier.

267. Déterminez les racines carrées d’un nombre complexe

Sur un exemple concret, vous allez voir comment il est possible de trouver tous les nombres complexes $z$ tels que :

z^2 = -3-4i.

Justifiez de l’existence d’un tel nombre

L’utilisation de l’artillerie lourde permet de conclure : d’après le théorème de d’Alembert, appelé aussi théorème fondamental de l’algèbre, tout polynôme non constant de $\C[X]$ admet une racine complexe.

Il en est donc ainsi pour le polynôme du second degré $X^2+3+4i.$

Déterminez la partie réelle et la partie imaginaire

Notez $z$ un nombre complexe tel que $z^2 = -3-4i.$ Il existe un couple $(a,b)\in\R^2$ tel que :

z=a+ib.

Vous élevez au carré :

\begin{align*}
z^2 &= (a+ib)^2\\
&= a^2+(ib)^2+2a\times ib\\
&=a^2+i^2b^2+2iab\\
&=a^2-b^2+i(2ab).
\end{align*}

En identifiant la partie réelle et la partie imaginaire de $z^2$ il vient :

\left\{\begin{align*}
a^2-b^2 &= -3\\
2ab &= -4.
\end{align*}\right.

Remarquez que la seconde équation est simplifiable après division par $2$ :

\left\{\begin{align*}
a^2-b^2 &= -3\\
ab &= -2.
\end{align*}\right.

Eliminez une inconnue

Partez de la relation :

a^2-b^2=-3.

Souhaitant utiliser $ab=-2$ vous multipliez par $a$ :

\begin{align*}
a^3-ab^2&=-3a\\
a^3-(ab)b &=-3a\\
a^3-(-2)b&=-3a\\
a^3+2b&=-3a.
\end{align*}

Vous avez abaissé le degré de $b$ qui est passé de $2$ à $1.$ Poursuivez en multipliant encore par $a$ :

\begin{align*}
a^3+2b&=-3a\\
a^4+2ab&=-3a^2\\
a^4-4&=-3a^2\\
a^4+3a^2-4 &= 0.
\end{align*}

Résolvez l’équation bicarrée

L’équation de degré $4$ se ramène à une équation de degré $2$ :

\begin{align*}
a^4+3a^2-4 &= 0\\
(a^2)^2+3a^2-4&=0.
\end{align*}

Posez $A = a^2.$ Le réel $A$ est solution de l’équation :

A^2+3A-4 = 0.

Vous constatez que $A=1$ est racine évidente. L’expression du membre de gauche est factorisable par $A-1.$ Cela conduit à écrire ce qui suit :

\begin{align*}
(A-1)(A+4) &=0.
\end{align*}

Du coup :

A\in\{1, -4\}.

Comme $A$ est un carré, $A$ est nécessairement positif, donc $A = 1.$

Ainsi, $a^2 = 1$ et donc $a\in\{-1,1\}.$

Trouvez les valeurs possibles

Si $a=1$

Comme $ab = -2$ vous déduisez $b=-2$ et donc $z = 1-2i.$

Si $a=-1$

Comme $ab = -2$ vous déduisez $b=2$ et donc $z = -1+2i.$

Si $z$ est une solution de l’équation $z^2 = -3-4i$, alors $z\in\{1-2i, -1+2i\}.$

Effectuez la réciproque

Si $z = 1-2i$, alors :

\begin{align*}
z^2  &= (1-2i)^2\\
&= 1^2+(2i)^2-2\times 1\times 2i\\
&= 1-4-4i\\
&=-3-4i.
\end{align*}

De même, si $z = -1+2i$, alors :

\begin{align*}
z^2 &= (-1+2i)^2\\
&= ((-1)\times (1-2i))^2\\
&=(-1)^2(1-2i)^2\\
&=1\times(-3-4i)\\
&=-3-4i.
\end{align*}

Concluez

Pour tout nombre complexe $z$, l’équivalence suivante est vérifiée :

\boxed{z^2 = -3-4i \Longleftrightarrow z\in\{1-2i, -1+2i\}.}

264. Déterminez la valeur du discriminant universel d’ordre 3

Dans cet article vous travaillez dans l’anneau des polynômes à $3$ indéterminées à coefficients dans $\Z$ qui sera noté $\Z[X_1,X_2,X_3].$

Par définition, le discriminant universel d’ordre $3$ est défini par :

\boxed{\Delta = (X_1-X_2)^2(X_1-X_3)^2(X_2-X_3)^2.}

Les carrés sont pris pour obtenir une invariance de $\Delta$ par toute permutation des indéterminées.

Ainsi, vous savez d’après le théorème sur les polynômes symétriques que $\Delta$ s’écrit comme un polynôme en fonction des polynômes symétriques élémentaires $\sigma_1$, $\sigma_2$ et $\sigma_3$ définis par :

\boxed{
\begin{align*}
\sigma_1 &= X_1+X_2+X_3 \\
\sigma_2 &= X_1X_2+X_1X_3+X_2X_3 \\
\sigma_3 &= X_1X_2X_3.
\end{align*}
}

L’objectif est ici de déterminer un tel polynôme.

Développez le discriminant $\Delta$

Vous développez d’abord le polynôme $R$ défini par $R = (X_1-X_2)(X_1-X_3)(X_2-X_3)$ en prenant soin d’utiliser l’ordre lexicographique sur les puissances respectives des $3$ indéterminées :

\begin{align*}
R &=   \left[(X_1-X_2)(X_1-X_3)\right](X_2-X_3) \\
 &= (X_1^2-X_1X_3-X_1X_2+X_2X_3)(X_2-X_3)\\
 &= (X_1^2-X_1X_2-X_1X_3+X_2X_3)(X_2-X_3)\\
 &= X_1^2X_2-X_1X_2^2-X_1X_2X_3+X_2^2X_3\\
&\quad -X_1^2X_3+X_1X_2X_3+X_1X_3^2-X_2X_3^2\\
 &= X_1^2X_2-X_1^2X_3-X_1X_2^2+X_1X_3^2+X_2^2X_3 -X_2X_3^2.
\end{align*}

En mettant le tout au carré, vous obtenez, en ordonnant les exposants $(4,2,0)$, $(4,1,1)$, $(3,3,0)$, $(3,2,1)$, $(2,2,2)$ avec leurs permutés :

\begin{align*}
\Delta &= R^2 \\
 &=  ( X_1^2X_2-X_1^2X_3-X_1X_2^2+X_1X_3^2+X_2^2X_3 -X_2X_3^2) ^2 \\
&=X_1^4X_2^2+X_1^4X_3^2+X_1^2X_2^4+X_1^2X_3^4+X_2^4X_3^2 +X_2^2X_3^4\\
&\quad -2X_1^4X_2X_3-2X_1^3X_2^3+2X_1^3X_2X_3^2+2X_1^2X_2^3X_3-2X_1^2X_2^2X_3^2\\
&\quad +2X_1^3X_2^2X_3-2X_1^3X_3^3-2X_1^2X_2^2X_3^2+2X_1^2X_2X_3^3\\
&\quad -2X_1^2X_2^2X_3^2-2X_1X_2^4X_3+2X_1X_2^3X_3^2\\
&\quad +2X_1X_2^2X_3^3-2X_1X_2X_3^4\\
&\quad -2X_2^3X_3^3\\
&=X_1^4X_2^2+X_1^4X_3^2+X_1^2X_2^4+X_1^2X_3^4+X_2^4X_3^2 +X_2^2X_3^4\\
&\quad -2X_1^4X_2X_3-2X_1X_2^4X_3-2X_1X_2X_3^4\\
&\quad -2X_1^3X_2^3-2X_1^3X_3^3-2X_2^3X_3^3\\
&\quad +2X_1^3X_2^2X_3+2X_1^3X_2X_3^2+2X_1^2X_2^3X_3+2X_1^2X_2X_3^3+2X_1X_2^3X_3^2+2X_1X_2^2X_3^3\\
&\quad -6X_1^2X_2^2X_3^2.
\end{align*}

Traitez l’exposant $(4,2,0)$ et ses permutés

Pour éliminer la somme :

X_1^4X_2^2+X_1^4X_3^2+X_1^2X_2^4+X_1^2X_3^4+X_2^4X_3^2 +X_2^2X_3^4

vous remarquez que $X_1^4X_2^2 = X_1^2(X_1X_2)^2.$ Il convient de développer $\sigma_1^2\sigma_2^2.$ Vous commencez par le développement du produit suivant :

\begin{align*}
\sigma_1\sigma_2 &= (X_1+X_2+X_3)(X_1X_2+X_1X_3+X_2X_3)\\
&=X_1^2X_2+X_1^2X_3+X_1X_2X_3\\
&\quad +X_1X_2^2+X_1X_2X_3+X_2^2X_3\\
&\quad +X_1X_2X_3+X_1X_3^2+X_2X_3^2\\
&=X_1^2X_2+X_1^2X_3+X_1X_2^2+X_1X_3^2+X_2^2X_3+X_2X_3^2\\
&\quad + 3X_1X_2X_3.
\end{align*}

Puis vous mettez le tout au carré :

\begin{align*}
\sigma_1^2\sigma_2^2 &= (\sigma_1\sigma_2)^2 \\
&= (X_1^2X_2+X_1^2X_3+X_1X_2^2+X_1X_3^2+X_2^2X_3+X_2X_3^2 + 3X_1X_2X_3)^2\\
&=X_1^4X_2^2+X_1^4X_3^2+X_1^2X_2^4+X_1^2X_3^4+X_2^4X_3^2+X_2^2X_3^4+ 9X_1^2X_2^2X_3^2\\
&\quad  +2X_1^4X_2X_3+2X_1^3X_2^3+2X_1^3X_2X_3^2+2X_1^2X_2^3X_3+2X_1^2X_2^2X_3^2+6X_1^3X_2^2X_3\\
&\quad +2X_1^3X_2^2X_3+2X_1^3X_3^3+2X_1^2X_2^2X_3^2+2X_1^2X_2X_3^3+6X_1^3X_2X_3^2\\
&\quad +2X_1^2X_2^2X_3^2+2X_1X_2^4X_3+2X_1X_2^3X_3^2+6X_1^2X_2^3X_3\\
&\quad +2X_1X_2^2X_3^3+2X_1X_2X_3^4+6X_1^2X_2X_3^3\\
&\quad +2X_2^3X_3^3+6X_1X_2^3X_3^2\\
&\quad +6X_1X_2^2X_3^3\\
&=X_1^4X_2^2+X_1^4X_3^2+X_1^2X_2^4+X_1^2X_3^4+X_2^4X_3^2+X_2^2X_3^4\\
&\quad +2X_1^4X_2X_3+2X_1X_2^4X_3+2X_1X_2X_3^4\\
&\quad +2X_1^3X_2^3+2X_1^3X_3^3+2X_2^3X_3^3\\
&\quad +8X_1^3X_2^2X_3+8X_1^3X_2X_3^2 +8X_1^2X_2^3X_3+8X_1^2X_2X_3^3+8X_1X_2^3X_3^2+8X_1X_2^2X_3^3\\
&\quad+ 15X_1^2X_2^2X_3^2.
\end{align*}

Par soustraction avec l’expression développée de $\Delta$ il vient :

\begin{align*}
\Delta - \sigma_1^2\sigma_2^2&= -4X_1^4X_2X_3-4X_1X_2^4X_3-4X_1X_2X_3^4\\
&\quad -4X_1^3X_2^3-4X_1^3X_3^3-4X_2^3X_3^3\\
&\quad -6X_1^3X_2^2X_3-6X_1^3X_2X_3^2-6X_1^2X_2^3X_3-6X_1^2X_2X_3^3-6X_1X_2^3X_3^2-6X_1X_2^2X_3^3\\
&\quad -21X_1^2X_2^2X_3^2.
\end{align*}

Traitez l’exposant $(4,1,1)$ et ses permutés

Pour éliminer les trois permutés :

 -4X_1^4X_2X_3-4X_1X_2^4X_3-4X_1X_2X_3^4

il convient de remarquer que $X_1^4X_2X_3 = X_1^3(X_1X_2X_3).$ Vous effectuez le développement de $\sigma_1^3\sigma_3.$

Tout d’abord :

\begin{align*}
\sigma_1^2 &= (X_1+X_2+X_3)^2\\
&=X_1^2+X_2^2+X_3^2+2X_1X_2+2X_1X_3+2X_2X_3.
\end{align*}

Ensuite :

\begin{align*}
\sigma_1^3 &= \sigma_1^2\times \sigma_1\\
&=(X_1^2+X_2^2+X_3^2+2X_1X_2+2X_1X_3+2X_2X_3)(X_1+X_2+X_3)\\
&=X_1^3+X_1X_2^2+X_1X_3^2+2X_1^2X_2+2X_1^2X_3+2X_1X_2X_3\\
&\quad + X_1^2X_2+X_2^3+X_2X_3^2+2X_1X_2^2+2X_1X_2X_3+2X_2^2X_3\\
&\quad + X_1^2X_3+X_2^2X_3+X_3^3+2X_1X_2X_3+2X_1X_3^2+2X_2X_3^2 \\
&= X_1^3+X_2^3+X_3^3\\
&\quad +3X_1^2X_2+3X_1^2X_3+3X_1X_2^2+3X_1X_3^2+3X_2^2X_3+3X_2X_3^2\\
&\quad +6X_1X_2X_3.
\end{align*}

Du coup :

\begin{align*}
\sigma_1^3\sigma_3 &= (X_1^3+X_2^3+X_3^3)(X_1X_2X_3)\\
&\quad +(3X_1^2X_2+3X_1^2X_3+3X_1X_2^2+3X_1X_3^2+3X_2^2X_3+3X_2X_3^2)(X_1X_2X_3)\\
&\quad +6X_1^2X_2^2X_3^2\\
&= X_1^4X_2X_3+X_1X_2^4X_3+X_1X_2X_3^4 \\
&\quad +3X_1^3X_2^2X_3+3X_1^3X_2X_3^2+3X_1^2X_2^3X_3+3X_1^2X_2X_3^3+3X_1X_2^3X_3^2+3X_1X_2^2X_3^3\\
&\quad +6X_1^2X_2^2X_3^2.
\end{align*}

Vous multipliez par $4$ :

\begin{align*}
4\sigma_1^3\sigma_3 &= 4X_1^4X_2X_3+4X_1X_2^4X_3+4X_1X_2X_3^4 \\
&\quad +12X_1^3X_2^2X_3+12X_1^3X_2X_3^2+12X_1^2X_2^3X_3+12X_1^2X_2X_3^3+12X_1X_2^3X_3^2+12X_1X_2^2X_3^3\\
&\quad +24X_1^2X_2^2X_3^2.
\end{align*}

Vous déduisez ainsi :

\begin{align*}
\Delta - \sigma_1^2\sigma_2^2+4\sigma_1^3\sigma_3 &= -4X_1^3X_2^3-4X_1^3X_3^3-4X_2^3X_3^3\\
&\quad +6X_1^3X_2^2X_3+6X_1^3X_2X_3^2+6X_1^2X_2^3X_3+6X_1^2X_2X_3^3+6X_1X_2^3X_3^2+6X_1X_2^2X_3^3\\
&\quad +3X_1^2X_2^2X_3^2.
\end{align*}

Traitez l’exposant $(3,3,0)$ et ses permutés

Pour éliminer $-4X_1^3X_2^3-4X_1^3X_3^3-4X_2^3X_3^3$ il convient de calculer $\sigma_2^3.$

Tout d’abord :

\begin{align*}
\sigma_2^2&=(X_1X_2+X_1X_3+X_2X_3)^2\\
&=X_1^2X_2^2+X_1^2X_3^2+X_2^2X_3^2\\
&\quad + 2X_1^2X_2X_3+2X_1X_2^2X_3+2X_1X_2X_3^2.
\end{align*}

Puis :

\begin{align*}
\sigma_2^3 &= \sigma_2^2\times \sigma_2\\
&= (X_1^2X_2^2+X_1^2X_3^2+X_2^2X_3^2)(X_1X_2+X_1X_3+X_2X_3)\\
&\quad + (2X_1^2X_2X_3+2X_1X_2^2X_3+2X_1X_2X_3^2)(X_1X_2+X_1X_3+X_2X_3)\\
&= X_1^3X_2^3+X_1^3X_2X_3^2+X_1X_2^3X_3^2\\
&\quad + X_1^3X_2^2X_3+X_1^3X_3^3+X_1X_2^2X_3^3\\
&\quad +X_1^2X_2^3X_3+X_1^2X_2X_3^3+X_2^3X_3^3\\
&\quad +2X_1^3X_2^2X_3+2X_1^2X_2^3X_3+2X_1^2X_2^2X_3^2\\
&\quad +2X_1^3X_2X_3^2+2X_1^2X_2^2X_3^2+2X_1^2X_2X_3^3\\
&\quad +2X_1^2X_2^2X_3^2+2X_1X_2^3X_3^2+2X_1X_2^2X_3^3\\
&=X_1^3X_2^3+X_1^3X_3^3+X_2^3X_3^3\\
&\quad +3X_1^3X_2^2X_3+3X_1^3X_2X_3^2+3X_1^2X_2^3X_3+3X_1^2X_2X_3^3+3X_1X_2^3X_3^2+3X_1X_2^2X_3^3\\
&\quad +6X_1^2X_2^2X_3^2.
\end{align*}

Vous multipliez par $4$ :

\begin{align*}
4\sigma_2^3 &=4X_1^3X_2^3+4X_1^3X_3^3+4X_2^3X_3^3\\
&\quad +12X_1^3X_2^2X_3+12X_1^3X_2X_3^2+12X_1^2X_2^3X_3+12X_1^2X_2X_3^3+12X_1X_2^3X_3^2+12X_1X_2^2X_3^3\\
&\quad +24X_1^2X_2^2X_3^2.
\end{align*}

Et vous déduisez :

\begin{align*}
\Delta - \sigma_1^2\sigma_2^2+4\sigma_1^3\sigma_3 +4\sigma_2^3&= 18X_1^3X_2^2X_3+18X_1^3X_2X_3^2+18X_1^2X_2^3X_3+18X_1^2X_2X_3^3+18X_1X_2^3X_3^2+18X_1X_2^2X_3^3\\
&\quad +27X_1^2X_2^2X_3^2.
\end{align*}

Traitez l’exposant $(3,2,1)$ et ses permutés

Pour éliminer $18X_1^3X_2^2X_3+18X_1^3X_2X_3^2+18X_1^2X_2^3X_3+18X_1^2X_2X_3^3+18X_1X_2^3X_3^2+18X_1X_2^2X_3^3$ il convient de calculer $\sigma_1\sigma_2\sigma_3$ puisque $X_1^3X_2^2X_3 = X_1(X_1X_2)(X_1X_2X_3).$

Le produit $\sigma_1\sigma_2$ a déjà été calculé plus haut :

\begin{align*}
\sigma_1\sigma_2 &=X_1^2X_2+X_1^2X_3+X_1X_2^2+X_1X_3^2+X_2^2X_3+X_2X_3^2\\
&\quad + 3X_1X_2X_3.
\end{align*}

Du coup :

\begin{align*}
\sigma_1\sigma_2\sigma_3 &=(X_1^2X_2+X_1^2X_3+X_1X_2^2+X_1X_3^2+X_2^2X_3+X_2X_3^2)(X_1X_2X_3)\\
&\quad + 3X_1^2X_2^2X_3^2 \\
&= X_1^3X_2^2X_3+X_1^3X_2X_3^2+X_1^2X_2^3X_3+X_1^2X_2X_3^3+X_1X_2^3X_3^2+X_1X_2^2X_3^3\\
&\quad + 3X_1^2X_2^2X_3^2.
\end{align*}

Vous multipliez par $18$ :

\begin{align*}
18\sigma_1\sigma_2\sigma_3 &= 18X_1^3X_2^2X_3+18X_1^3X_2X_3^2+18X_1^2X_2^3X_3+18X_1^2X_2X_3^3+18X_1X_2^3X_3^2+18X_1X_2^2X_3^3\\
&\quad + 54X_1^2X_2^2X_3^2.
\end{align*}

Il vient :

\begin{align*}
\Delta - \sigma_1^2\sigma_2^2+4\sigma_1^3\sigma_3 +4\sigma_2^3 - 18\sigma_1\sigma_2\sigma_3&= -27X_1^2X_2^2X_3^2.
\end{align*}

Concluez

D’après le dernier résultat obtenu, vous avez :

\begin{align*}
\Delta - \sigma_1^2\sigma_2^2+4\sigma_1^3\sigma_3 +4\sigma_2^3 - 18\sigma_1\sigma_2\sigma_3&= -27\sigma_3^2.
\end{align*}

Vous avez ainsi obtenu l’expression du discriminant universel d’ordre $3$ en fonction des polynômes symétriques élémentaires :

\boxed{\Delta =   \sigma_1^2\sigma_2^2-4\sigma_1^3\sigma_3 -4\sigma_2^3 + 18\sigma_1\sigma_2\sigma_3 -27\sigma_3^2.}

263. Calculez une intégrale qui comporte une racine carrée au dénominateur

Cet article va vous permettre de calculer l’intégrale suivante :

I=\displaystyle\int_{-1}^1 \frac{x^3+2x^2+1}{\sqrt{x^2+1}}\,\mathrm{d}x.

Gérez le dénominateur

Le dénominateur qui est $\sqrt{x^2+1}$ pousse à utiliser le changement de variable $x = \mathrm{sh}\,t$, de sorte que $1+x^2 = 1+(\mathrm{sh}\,t)^2 = (\mathrm{ch}\,t)^2.$ La positivité de $1+x^2$ et de $\mathrm{ch}\,t$ aboutit à $\sqrt{x^2+1}=\mathrm{ch}\,t.$

Du coup :

\begin{align*}
 \frac{x^3+2x^2+1}{\sqrt{x^2+1}} &= \frac{(\mathrm{sh} \,t)^3+2(\mathrm{sh}\,t)^2+1}{\mathrm{ch}\,t}.
\end{align*}

Or :

\begin{align*} \mathrm{d}x &= \mathrm{ch}\,t\,\mathrm{d}t. \end{align*}

Ainsi :

\begin{align*}
\frac{x^3+2x^2+1}{\sqrt{x^2+1}}\,\mathrm{d}x &= \frac{(\mathrm{sh} \,t)^3+2(\mathrm{sh}\,t)^2+1}{\mathrm{ch}\,t}\mathrm{ch}\,t\,\mathrm{d}t\\
&=\left[(\mathrm{sh}\,t)^3+2(\mathrm{sh}\,t)^2+1\right]\,\mathrm{d}t.
\end{align*}

Pour le calcul des nouvelles bornes de l’intégrale obtenue par changement de variable, la forme explicite de la fonction réciproque de la fonction $\mathrm{sh}$ va s’avérer précieuse. En effet, quels que soient les réels $a$ et $b$, vous avez $b=\mathrm{sh}\,a \Longleftrightarrow a=\mathrm{Argsh}\,b = \ln\left(b+\sqrt{1+b^2}\right).$

Quand $x=1$ vous avez $t=\mathrm{Argsh}\, 1=\ln\left(1+\sqrt{2}\right).$

Par imparité de la fonction $\mathrm{Argsh}$ il vient, quand $x=-1$, $t=\mathrm{Argsh}(-1) = -\mathrm{Argsh}\, 1 = -\ln\left(1+\sqrt{2}\right).$

Ainsi :

I = \displaystyle\int_{-\ln\left(1+\sqrt{2}\right)}^{\ln\left(1+\sqrt{2}\right)} \left[(\mathrm{sh}\,t)^3+2(\mathrm{sh}\,t)^2+1\right]\,\mathrm{d}t.

Utilisant l’imparité de la fonction $t\mapsto (\mathrm{sh}\,t)^3$ il vient :

\int_{-\ln\left(1+\sqrt{2}\right)}^{\ln\left(1+\sqrt{2}\right)} (\mathrm{sh}\,t)^3\,\mathrm{d}t = 0.

Or :

\begin{align*}
(\mathrm{sh}\,t)^2&=\left(\displaystyle\frac{\mathrm{e}^t-\mathrm{e}^{-t}}{2}\right)^2 \\
 &= \frac{\mathrm{e}^{2t} + \mathrm{e}^{-2t} -2}{4} \\ 
&= \frac{2\mathrm{ch}(2t)-2}{4} \\
 &= \frac{\mathrm{ch}(2t)-1}{2}.
\end{align*}

Si bien que $2(\mathrm{sh}\,t)^2+1=\mathrm{ch}(2t).$

Donc :

I =\displaystyle\int_{-\ln\left(1+\sqrt{2}\right)}^{\ln\left(1+\sqrt{2}\right)} \mathrm{ch}(2t)\,\mathrm{d}t.

La fonction $t\mapsto \mathrm{ch}(2t)$ étant paire, vous avez :

I = 2\int_{0}^{\ln\left(1+\sqrt{2}\right)} \mathrm{ch}(2t)\,\mathrm{d}t.

Vous effectuez le changement de variable $u=2t$ qui donne $\mathrm{d}u=2\mathrm{d}t$ d’où :

\begin{align*}
I &= \int_{0}^{\ln\left(1+\sqrt{2}\right)} \mathrm{ch}(2t)\,2\mathrm{d}t\\
 &=\int_{0}^{2\ln\left(1+\sqrt{2}\right)} \mathrm{ch}\,u\,\mathrm{d}u \\
 &= \left[\mathrm{sh}\,u\right]_0^{2\ln\left(1+\sqrt{2}\right) } \\
 &= \mathrm{sh}\left(2\ln\left(1+\sqrt{2}\right)\right)\\
 &= \mathrm{sh}\left(\ln\left[\left(1+\sqrt{2}\right)^2\right]\right)\\
 &= \mathrm{sh}\left(\ln\left(3+2\sqrt{2}\right)\right) \\
 &= \frac{3+2\sqrt{2} - \frac{1}{3+2\sqrt{2}}}{2} \\
 &= \frac{3+2\sqrt{2} - \frac{3-2\sqrt{2}}{(3+2\sqrt{2})(3-2\sqrt{2})}}{2}\\
 &= \frac{3+2\sqrt{2} - \frac{3-2\sqrt{2}}{1}}{2} \\
 &= \frac{3+2\sqrt{2} - 3+2\sqrt{2}}{2} \\
 &= \frac{4\sqrt{2}}{2} \\
 &= 2\sqrt{2}.
\end{align*}

En définitive :

\boxed{\displaystyle\int_{-1}^1 \frac{x^3+2x^2+1}{\sqrt{x^2+1}}\,\mathrm{d}x = 2\sqrt{2}.}

262. Calculez la valeur d’une intégrale impropre

Dans cet article, vous souhaitez justifier la convergence de l’intégrale impropre $\int_0^{+\infty} \frac{\mathrm{d}x}{x^2+\sqrt{x}}$ et déterminer sa valeur.

Isolez le problème avec $+\infty$

Soit $M$ un nombre réel strictement positif fixé.
L’idée principale est de comprendre pourquoi l’intégrale $I=\displaystyle\int_0^M \frac{\mathrm{d}x}{x^2+\sqrt{x}}$ est bien définie. En effet, il y a un problème en $0$ où la fonction $x\mapsto \frac{1}{x^2+\sqrt{x}}$ n’est pas définie en $0.$

Effectuez un changement de variable

Soit maintenant $\varepsilon$ un réel strictement positif tel que $\varepsilon < M.$

Afin d’éliminer la racine carrée, vous posez $y=\sqrt{x}$ de sorte que $y^2=x$ et $2y\mathrm{d}y=\mathrm{d}x.$
Du coup :

\begin{align*}
\int_{\varepsilon}^M \frac{\mathrm{d}x}{x^2+\sqrt{x}}&=\displaystyle\int_{\sqrt{\varepsilon}}^{\sqrt{M}} \frac{2y\mathrm{d}y}{y^4+y}\ &= 2\int_{\sqrt{\varepsilon}}^{\sqrt{M}} \frac{\mathrm{d}y}{y^3+1}. \end{align*}

La dernière expression fait apparaître l’intégrale de la fonction $x\mapsto \frac{1}{x^3+1}$ qui est une fonction continue sur le segment $[0,\sqrt{M}].$

Ainsi, la limite de $\int_{\sqrt{\varepsilon}}^{\sqrt{M}} \frac{\mathrm{d}y}{y^3+1}$ quand $\varepsilon\to 0$ existe et :

\lim_{\varepsilon\to 0} \int_{\sqrt{\varepsilon}}^{\sqrt{M}} \frac{\mathrm{d}y}{y^3+1} = \int_{0}^{\sqrt{M}} \frac{\mathrm{d}y}{y^3+1}.

Vous en déduisez la convergence de l’intégrale $I = \int_{0}^M \frac{\mathrm{d}x}{x^2+\sqrt{x}}.$

D’après ce qui précède, vous avez $I = \int_{0}^{\sqrt{M}} \frac{\mathrm{d}y}{y^3+1}.$

Effectuez une décomposition en éléments simples

Dans le corps des fractions rationnelles $\R(X)$, la fraction $\displaystyle\frac{1}{X^3+1}$ n’est pas décomposée en éléments simples, puisque le dénominateur $X^3+1$ est un polynôme de degré $3.$

L’équation $X^3+1=0$ admettant une unique solution réelle qui est $X=-1$, vous savez que le polynôme $X^3+1$ est factorisable par $X+1.$

La division euclidienne s’effectue de la façon suivante :

\begin{align*}
 X^3+1 &= (X+1)X^2-X^2+1\\
 &= (X+1)X^2-(X+1)X+X+1\\ 
&= (X+1)X^2-(X+1)X+(X+1)1\\
 &= (X+1)(X^2-X+1). 
\end{align*}

D’après le théorème de décomposition en éléments simples, il existe un triplet de réels $(a,b,c)\in\R^3$ tel que :

\displaystyle\frac{1}{X^3+1} = \frac{a}{X+1}+\frac{bX+c}{X^2-X+1}.

Le plus rapide est de calculer $a$ en multipliant par $X+1$ puis en évaluant en $X=-1$. Vous obtenez :

\begin{align*} a &= \left.\frac{1}{X^2-X+1}\right\vert_{X=-1}\\
 &= \frac{1}{1+1+1}\\
 &=\frac{1}{3}. 
\end{align*}

Il est rapide de trouver $b$ à partir du nombre $a$, vous multipliez le tout par $X$ :

\displaystyle\frac{X}{X^3+1} = \frac{aX}{X+1}+\frac{bX^2+cX}{X^2-X+1}.

Puis vous passez à la limite lorsque $X\to +\infty$. Vous obtenez $0 = a+b$ donc $b=\displaystyle-\frac{1}{3}.$
Comme il ne reste plus qu’un coefficient à trouver, vous pouvez effectuer $X=0$ dans l’égalité :

\displaystyle\frac{1}{X^3+1} = \frac{a}{X+1}+\frac{bX+c}{X^2-X+1}.

D’où $1 = a+c$ et par suite $c=1-a = \frac{2}{3}.$ Ainsi :

\boxed{\displaystyle\frac{1}{X^3+1} = \frac{1}{3(X+1)}+\frac{-X+2}{3(X^2-X+1)}.}

Intégrez les éléments simples

Notez $J = \displaystyle\int_{0}^{\sqrt{M}} \frac{\mathrm{d}x}{x+1}$ et $ K= \displaystyle\int_{0}^{\sqrt{M}} \frac{-x+2}{x^2-x+1}\,\mathrm{d}x.$

Calculez $J$

\begin{align*} J &= \displaystyle \left[\ln(x+1)\right]{0}^{\sqrt{M}}\\
 &=\ln\left(\sqrt{M}+1\right). 
\end{align*}

Calculez $K$

L’idée est d’utiliser la forme canonique du trinôme $X^2-X+1 = \displaystyle\left(X-\frac{1}{2}\right)^2+\frac{3}{4}$ pour préparer un changement de variable :

\begin{align*} K&= \displaystyle\int_{0}^{\sqrt{M}} \frac{-x+2}{x^2-x+1}\,\mathrm{d}x\\ &=\int_{0}^{\sqrt{M}} \frac{-x+2}{\left(x-\frac{1}{2}\right)^2+\frac{3}{4}}\,\mathrm{d}x.
\end{align*}

Vous posez alors $y=\displaystyle x-\frac{1}{2}.$ Alors :

\begin{align*} K &=\displaystyle\int_{-1/2}^{\sqrt{M}-1/2} \frac{-(y+\frac{1}{2})+2}{y^2+\frac{3}{4}}\,\mathrm{d}y \\
 &=\int_{-1/2}^{\sqrt{M}-1/2} \frac{-y+\frac{3}{2}}{y^2+\frac{3}{4}}\,\mathrm{d}y\\
 &=\int_{-1/2}^{\sqrt{M}-1/2} \frac{-y}{y^2+\frac{3}{4}}\,\mathrm{d}y + \int_{-1/2}^{\sqrt{M}-1/2} \frac{\frac{3}{2}}{y^2+\frac{3}{4}}\,\mathrm{d}y\\
 &=-\frac{1}{2}\int_{-1/2}^{\sqrt{M}-1/2} \frac{2y}{y^2+\frac{3}{4}}\,\mathrm{d}y + \frac{3}{2}\int_{-1/2}^{\sqrt{M}-1/2} \frac{1}{y^2+\left(\frac{\sqrt{3}}{2}\right)^2}\,\mathrm{d}y \\
&= -\frac{1}{2}\left[ \ln\left(y^2+\frac{3}{4}\right)\right]_{-1/2}^{\sqrt{M}-1/2} + \frac{\sqrt{3}\times \sqrt{3}}{2} \frac{1}{\frac{\sqrt{3}}{2}}\left[\arctan\left(\frac{y}{\frac{\sqrt{3}}{2}}\right)\right]_{-1/2}^{\sqrt{M}-1/2}\\
 &=-\frac{1}{2}\ln (M-\sqrt{M}+1) + \sqrt{3} \left[\arctan\left(\frac{y}{\frac{\sqrt{3}}{2}}\right)\right]_{-1/2}^{\sqrt{M}-1/2}\\
 &=-\frac{1}{2}\ln (M-\sqrt{M}+1) + \sqrt{3} \arctan\left(\frac{\sqrt{M}-1/2}{\frac{\sqrt{3}}{2}}\right) - \sqrt{3} \arctan\left(\frac{-1/2}{\frac{\sqrt{3}}{2}}\right) \\
 &=-\frac{1}{2}\ln (M-\sqrt{M}+1) + \sqrt{3} \arctan\left(\frac{\sqrt{M}-1/2}{\frac{\sqrt{3}}{2}}\right) + \sqrt{3} \arctan\left(\frac{1}{\sqrt{3}}\right) \\ &=-\frac{1}{2}\ln (M-\sqrt{M}+1) + \sqrt{3} \arctan\left(\frac{\sqrt{M}-1/2}{\frac{\sqrt{3}}{2}}\right) + \frac{\pi\sqrt{3}}{6}.
\end{align*}

Passez à la limite

\begin{align*} I&=\displaystyle\int_0^M \frac{\mathrm{d}x}{x^2+\sqrt{x}}\\
 &= 2\int_0^{\sqrt{M}} \frac{\mathrm{d}x}{x^3+1}\\
 &= \frac{2}{3} J + \frac{2}{3}K \\
 &=\frac{2}{3}\ln\left(\sqrt{M}+1\right) -\frac{1}{3}\ln (M-\sqrt{M}+1) + \frac{2\sqrt{3}}{3} \arctan\left(\frac{\sqrt{M}-1/2}{\frac{\sqrt{3}}{2}}\right) + \frac{\pi\sqrt{3}}{9}\\
 &=\frac{1}{3}\ln\left(M+2\sqrt{M}+1\right) -\frac{1}{3}\ln (M-\sqrt{M}+1)+ \frac{2\sqrt{3}}{3} \arctan\left(\frac{\sqrt{M}-1/2}{\frac{\sqrt{3}}{2}}\right) + \frac{\pi\sqrt{3}}{9} \\
 &=\frac{1}{3}\ln \frac{M+2\sqrt{M}+1}{M-\sqrt{M}+1} + \frac{2\sqrt{3}}{3} \arctan\left(\frac{\sqrt{M}-1/2}{\frac{\sqrt{3}}{2}}\right) + \frac{\pi\sqrt{3}}{9}\\&= \frac{1}{3}\ln \frac{1+\frac{2}{\sqrt{M}}+\frac{1}{M}}{1-\frac{1}{\sqrt{M}}+\frac{1}{M}} + \frac{2\sqrt{3}}{3} \arctan\left(\frac{\sqrt{M}-1/2}{\frac{\sqrt{3}}{2}}\right) + \frac{\pi\sqrt{3}}{9}.
 \end{align*}

Quand $M$ tend vers $+\infty$ vous avez :

\frac{1+\frac{2}{\sqrt{M}}+\frac{1}{M}}{1-\frac{1}{\sqrt{M}}+\frac{1}{M}} \to 1.

Vous déduisez :

\ln \frac{1+\frac{2}{\sqrt{M}}+\frac{1}{M}}{1-\frac{1}{\sqrt{M}}+\frac{1}{M}} \to 0.

Quand $M$ tend vers $+\infty$ vous avez :

\frac{\sqrt{M}-1/2}{\frac{\sqrt{3}}{2}} \to +\infty

Du coup :

\arctan\left(\frac{\sqrt{M}-1/2}{\frac{\sqrt{3}}{2}}\right) \to \frac{\pi}{2}.

Concluez

De ce qui précède, la limite de $\displaystyle\int_0^M \frac{\mathrm{d}x}{x^2+\sqrt{x}}$ quand $M\to +\infty$ existe. Elle est égale à :

\begin{align*}\displaystyle\int_0^{+\infty} \frac{\mathrm{d}x}{x^2+\sqrt{x}} 
&= \frac{2\sqrt{3}}{3}\times \frac{\pi}{2} + \frac{\pi\sqrt{3}}{9}\\ 
&= \frac{\pi\sqrt{3}}{3} + \frac{\pi\sqrt{3}}{9} \\ 
&= \frac{3\pi\sqrt{3}}{9} + \frac{\pi\sqrt{3}}{9} \\ 
&= \frac{4\pi\sqrt{3}}{9}. 
\end{align*}

En définitive :

\boxed{\displaystyle\int_0^{+\infty} \frac{\mathrm{d}x}{x^2+\sqrt{x}} = \frac{4\pi\sqrt{3}}{9}.}

261. Calculez l’intégrale d’une fraction rationnelle avec des termes de deuxième espèce

Pour calculer l’intégrale $\int_0^1 \frac{x^6}{(x^2+1)(x^4+1)}\mathrm{d}x$, vous pouvez utiliser :

  • Une décomposition en éléments simples dans $\R(x)$ avant d’intégrer. Cette décomposition fera apparaître des éléments de deuxième espèce ;
  • Les calculs d’intégrale de la forme $\int_{a}^b \frac{1}{x^2+1}\mathrm{d}x$ qui font intervenir des valeurs de la fonction arctangente. L’égalité $\tan \frac{\pi}{4}=1$ est fondamentale pour mener à bien les calculs, mais elle ne suffira pas, d’autres angles feront leur apparition.

Utilisez un changement de variable

Une première étape consiste à décomposer la fraction rationnelle $\frac{x^6}{(x^2+1)(x^4+1)} = \frac{(x^2)^3}{(x^2+1)((x^2)^2+1)}$ en éléments simples. Vous posez : $X=x^2.$

L’égalité suivante :

\displaystyle\frac{x^6}{(x^2+1)(x^4+1)}=\frac{X^3}{(X+1)(X^2+1)}

va permettre de travailler dans le corps des fractions rationnelles $\R(X).$

Comme le degré du numérateur est égal à celui du dénominateur, il va falloir calculer la partie entière de cette fraction.
Vous effectuez la division euclidienne de $X^3$ par$(X+1)(X^2+1) = X^3+X^2+X+1.$

Comme $X^3 = (X^3+X^2+X+1)\times 1 + (-X^2-X-1)$, vous obtenez après division :

\displaystyle\frac{X^3}{(X+1)(X^2+1)} = 1 + \frac{-X^2-X-1}{(X+1)(X^2+1)}.

Il reste maintenant à séparer les parties polaires de la fraction $\frac{-X^2-X-1}{(X+1)(X^2+1)}.$
Le théorème de décomposition en éléments simples fournit l’existence d’un triplet $(a,b,c)\in\R^3$ tel que :

\displaystyle\frac{-X^2-X-1}{(X+1)(X^2+1)}=\frac{a}{X+1}+\frac{bX+c}{X^2+1}.

Vous multipliez le tout par $X+1$ et évaluez avec $X=-1$ ce qui fournit :

a=\displaystyle\frac{-(-1)^2-(-1)-1}{(-1)^2+1}=\frac{-1+1-1}{1+1}=-\frac{1}{2}.

Vous multipliez le tout par $X^2+1$ et évaluez avec $X=i$ ce qui fournit :

bi+c=\displaystyle\frac{-i^2-i-1}{i+1}=\frac{-i}{1+i}=\frac{-i(1-i)}{2}=\frac{-1-i}{2}.

Du coup :

\begin{align*}
b&=\displaystyle-\frac{1}{2} \\
c&=\displaystyle-\frac{1}{2}.
\end{align*}

Il a été obtenu que :

\displaystyle\frac{-X^2-X-1}{(X+1)(X^2+1)}=-\frac{1}{2(X+1)}-\frac{X+1}{2(X^2+1)}.

Du coup :

\boxed{\displaystyle\frac{X^3}{(X+1)(X^2+1)} = 1 -\frac{1}{2(X+1)}-\frac{X+1}{2(X^2+1)}.}

Traitez les termes de seconde espèce

En revenant dans le corps des fractions rationnelles $\R(x)$,il vient :

\boxed{\displaystyle\frac{x^6}{(x^2+1)(x^4+1)} = 1 -\frac{1}{2(x^2+1)}-\frac{x^2+1}{2(x^4+1)}.}

La fraction $\frac{x^2+1}{x^4+1}$ n’est ni un terme de première espèce, ni un terme de seconde espèce.
En effet, le dénominateur n’est pas irréductible et le numérateur a un degré trop élevé.

Décomposez le polynôme de degré $4$ en produit de deux polynômes

Vous faites apparaître une identité remarquable :

x^4+1=x^4+2x^2+1-2x^2=(x^2+1)^2-2x^2.

Vous préparez ensuite la seconde identité remarquable :

x^4+1=x^4+2x^2+1-2x^2=(x^2+1)^2-(\sqrt{2}x)^2.

Vous déduisez la factorisation recherchée :

\boxed{x^4+1=(x^2+\sqrt{2}x+1)(x^2-\sqrt{2}x+1).}

Séparez les parties polaires

Un réel élevé à la puissance $4$ étant toujours positif, le polynôme $x^4+1$ n’a pas de racine réelle.
Du coup, les deux polynômes $x^2+\sqrt{2}+1$ et $x^2-\sqrt{2}+1$ n’ont pas de racine réelle.
Le théorème de décomposition des fractions rationnelles fournit l’existence et l’unicité d’un quadruplet $(a,b,c,d)\in\R^4$ tel que :

\displaystyle\frac{x^2+1}{x^4+1}=\frac{ax+b}{x^2+\sqrt{2}x+1}+\frac{cx+d}{x^2-\sqrt{2}x+1}.

En substituant $-x$ à $x$ il vient :

\begin{align*}
\displaystyle\frac{(-x)^2+1}{(-x)^4+1}&=\frac{a(-x)+b}{(-x)^2+\sqrt{2}(-x)+1}+\frac{c(-x)+d}{(-x)^2-\sqrt{2}(-x)+1}\\
\displaystyle\frac{x^2+1}{x^4+1}&=\frac{-ax+b}{x^2-\sqrt{2}x+1}+\frac{-cx+d}{x^2+\sqrt{2}x+1}\\
\displaystyle\frac{x^2+1}{x^4+1}&=\frac{-cx+d}{x^2+\sqrt{2}x+1}+\frac{-ax+b}{x^2-\sqrt{2}x+1}.
\end{align*}

Par unicité des coefficients, vous déduisez :

\displaystyle\frac{x^2+1}{x^4+1}=\frac{ax+b}{x^2+\sqrt{2}x+1}+\frac{-ax+b}{x^2-\sqrt{2}x+1}.

Notez maintenant $u$ une racine complexe et non réelle du polynôme $x^2+\sqrt{2}x+1.$ Alors $u^2=-1-\sqrt{2}u.$ En multipliant par $x^2+\sqrt{2}x+1$ et en substituant par $u$ il vient :

\begin{align*}
au+b&=\displaystyle\frac{u^2+1}{u^2-\sqrt{2}u+1} \\
&=\displaystyle\frac{-\sqrt{2}u}{-2\sqrt{2}u}\\&=\frac{-1}{2}.
\end{align*}

Si $a$ n’était pas nul, vous auriez $u=\frac{-\frac{1}{2}-b}{a}$ et $u$ serait réel, ce qui est absurde. Ainsi : $a=0.$
Il vient alors $b=\frac{1}{2}.$ D’où :

\boxed{\displaystyle\frac{x^2+1}{2(x^4+1)}=\frac{1}{4(x^2+\sqrt{2}x+1)}+\frac{1}{4(x^2-\sqrt{2}x+1)}.}

Déterminez la décomposition de la fraction rationnelle de départ

Vous avez obtenu ce qui suit :

\frac{x^6}{(x^2+1)(x^4+1)} = 1 -\frac{1}{2(x^2+1)}-\frac{x^2+1}{2(x^4+1)}

Finalement :

\boxed{\frac{x^6}{(x^2+1)(x^4+1)} = 1 -\frac{1}{2(x^2+1)}-\frac{1}{4(x^2+\sqrt{2}x+1)}-\frac{1}{4(x^2-\sqrt{2}x+1)}.}

Effectuez des calculs intermédiaires

Pour calculer $\int_0^1 \frac{x^6}{(x^2+1)(x^4+1)}\mathrm{d}x$ il sera commode de poser pour la suite :

\begin{align*}
I&=\displaystyle\int_0^1 \frac{1}{x^2+1}\mathrm{d}x\\
J&=\displaystyle\int_0^1 \frac{1}{x^2+\sqrt{2}x+1}\mathrm{d}x\\
K&=\displaystyle\int_0^1 \frac{1}{x^2-\sqrt{2}x+1}\mathrm{d}x.
\end{align*}

Calculez l’intégrale $I$

\begin{align*}
 I&=\int_0^1 \frac{1}{x^2+1}\mathrm{d}x\\
&=\left[\arctan x\right]_0^1 \\
&=\arctan 1 - \arctan 0 \\
&=\frac{\pi}{4}.
\end{align*}

Calculez l’intégrale $J$

L’idée est de se ramener à obtenir $z^2+1$ par changements de variables.

\begin{align*}
 J &= \int_0^1 \frac{1}{x^2+\sqrt{2}x+1}\mathrm{d}x \\
 &= \int_0^1 \frac{1}{\left(x+\frac{\sqrt{2}}{2}\right)^2-\frac{2}{4}+1}\mathrm{d}x \\
 &= \int_0^1 \frac{1}{\left(x+\frac{\sqrt{2}}{2}\right)^2+\frac{1}{2}}\mathrm{d}x 
\end{align*}

Cela conduit à poser $y=x+\frac{\sqrt{2}}{2}$ ce qui donne, avec $\mathrm{d}y=\mathrm{d}x$ :

\begin{align*} J &= \int_{\frac{\sqrt{2}}{2}}^{1+\frac{\sqrt{2}}{2}} \frac{1}{y^2+\frac{1}{2}}\mathrm{d}y \\
 &= \int_{\frac{\sqrt{2}}{2}}^{1+\frac{\sqrt{2}}{2}} \frac{2}{2y^2+1}\\
&= \int_{\frac{\sqrt{2}}{2}}^{1+\frac{\sqrt{2}}{2}} \frac{2}{\left(\sqrt{2}y\right)^2+1}\mathrm{d}y.
\end{align*}

Vous posez alors $z=\sqrt{2}y.$ D’où : $\mathrm{d}z=\sqrt{2}\mathrm{d}y$ et $2\mathrm{d}y=\sqrt{2}\mathrm{d}z$ et par suite :

\begin{align*}
 J &= \int_{1}^{1+\sqrt{2}} \frac{\sqrt{2}}{z^2+1}\mathrm{d}z \\
&= \sqrt{2}\int_{1}^{1+\sqrt{2}} \frac{1}{z^2+1}\mathrm{d}z \\
&= \sqrt{2} \left[\arctan z\right]_{1}^{1+\sqrt{2}} \\
&= \sqrt{2}(\arctan(1+\sqrt{2})-\arctan 1) \\
&= \sqrt{2}\left(\arctan(1+\sqrt{2})-\frac{\pi}{4}\right).
\end{align*}

Calculez l’intégrale $K$

\begin{align*}
 K &= \int_0^1 \frac{1}{x^2-\sqrt{2}x+1}\mathrm{d}x \\
 &= \int_0^1 \frac{1}{\left(x-\frac{\sqrt{2}}{2}\right)^2-\frac{2}{4}+1}\mathrm{d}x \\
 &= \int_0^1 \frac{1}{\left(x-\frac{\sqrt{2}}{2}\right)^2+\frac{1}{2}}\mathrm{d}x.
\end{align*}

Vous posez $y=x-\frac{\sqrt{2}}{2}$ ce qui donne :

\begin{align*}
 K &= \int_{-\frac{\sqrt{2}}{2}}^{1-\frac{\sqrt{2}}{2}} \frac{1}{y^2+\frac{1}{2}}\mathrm{d}y \\
 &= \int_{-\frac{\sqrt{2}}{2}}^{1-\frac{\sqrt{2}}{2}} \frac{2}{2y^2+1}\\
&= \int_{-\frac{\sqrt{2}}{2}}^{1-\frac{\sqrt{2}}{2}} \frac{2}{\left(\sqrt{2}y\right)^2+1}\mathrm{d}y.
\end{align*}

Vous posez encore $z=\sqrt{2}y.$ D’où :

\begin{align*} 
K &= \int_{-1}^{\sqrt{2}-1} \frac{\sqrt{2}}{z^2+1}\mathrm{d}z \\
&= \sqrt{2}\int_{-1}^{\sqrt{2}-1} \frac{1}{z^2+1}\mathrm{d}z \\
&= \sqrt{2} \left[\arctan z\right]_{-1}^{\sqrt{2}-1} \\
 &= \sqrt{2}(\arctan(\sqrt{2}-1)-\arctan (-1)) \\
 &= \sqrt{2}\left(\arctan(\sqrt{2}-1)-\frac{-\pi}{4}\right) \\
&= \sqrt{2}\left(\arctan(\sqrt{2}-1)+\frac{\pi}{4}\right).
\end{align*}

Calculez deux angles

Notez $\theta = \arctan(\sqrt{2}-1)$. Comme $\sqrt{2} > 1$ vous déduisez que $\sqrt{2}-1$ est strictement positif et par suite, $\theta$ est strictement positif. Du coup, $\theta\in\left]0,\frac{\pi}{2}\right[$.
De $\tan \theta = \sqrt{2}-1$ vous déduisez :
$\tan^2 \theta = \left(\sqrt{2}-1\right)^2 = 3-2\sqrt{2}.$

\begin{align*}
\cos^2 \theta &= \frac{1}{1+\tan^2 \theta} \\
 &= \frac{1}{4-2\sqrt{2}} \\
&=\frac{1}{2(2-\sqrt{2})} \\
&=\frac{2+\sqrt{2}}{4}.
\end{align*}

Du coup:

\begin{align*}
\cos 2\theta &= 2\cos^2 \theta - 1\\
 &=\frac{2+\sqrt{2}}{2}-1\\
 &=\frac{\sqrt{2}}{2}\\
 &=\cos \frac{\pi}{4}. 
\end{align*}

Or, les réels $2\theta$ et $\frac{\pi}{4}$ appartiennent à l’intervalle $[0,\pi].$ La fonction cosinus étant strictement décroissante sur cet intervalle, vous déduisez $2\theta = \frac{\pi}{4}$ donc $\theta = \frac{\pi}{8}.$

\boxed{\arctan\left(\sqrt{2}-1\right) = \frac{\pi}{8}.}

Notez maintenant $\alpha = \arctan(\sqrt{2}+1)$. Comme $\sqrt{2}+1>0$, vous avez $\alpha > 0$ donc $\alpha\in\left]0,\frac{\pi}{2}\right[$.
De $\tan \alpha = \sqrt{2}+1$ vous déduisez:

\tan^2 \alpha = \left(\sqrt{2}+1\right)^2 = 3+2\sqrt{2}.
\begin{align*}
\cos^2 \alpha &= \frac{1}{1+\tan^2 \alpha} \\
 &= \frac{1}{4+2\sqrt{2}} \\
&=\frac{1}{2(2+\sqrt{2})} \\
&=\frac{2-\sqrt{2}}{4}.
\end{align*}

Du coup:

\begin{align*}
\cos 2\alpha &= 2\cos^2 \alpha - 1\\
 &=\frac{2-\sqrt{2}}{2}-1\\
 &=\frac{-\sqrt{2}}{2}\\
 &=\cos \frac{3\pi}{4}.
 \end{align*}

Or, les réels $2\alpha$ et $\frac{3\pi}{4}$ appartiennent à l’intervalle $[0,\pi].$ La fonction cosinus étant strictement décroissante sur cet intervalle, vous déduisez $2\alpha = \frac{3\pi}{4}$ donc $\alpha = \frac{3\pi}{8}.$

\boxed{\arctan\left(\sqrt{2}+1\right) = \frac{3\pi}{8}}.

Concluez

De tout ce qui précède :

\begin{align*}
 \int_0^1 \frac{x^6}{(x^2+1)(x^4+1)}\mathrm{d}x \\
&= 1 - \frac{I}{2} - \frac{J}{4}-\frac{K}{4} \\
 &= 1-\frac{\frac{\pi}{4}}{2}-\frac{\sqrt{2}\left(\frac{3\pi}{8}-\frac{\pi}{4}\right)}{4} - \frac{\sqrt{2}\left(\frac{\pi}{8}+\frac{\pi}{4}\right)}{4}\\
 &= 1-\frac{\pi}{8}-\frac{\sqrt{2}\left(\frac{3\pi}{8}-\frac{2\pi}{8}\right)}{4} - \frac{\sqrt{2}\left(\frac{\pi}{8}+\frac{2\pi}{8}\right)}{4}\\
 &= 1-\frac{\pi}{8}-\frac{\pi\sqrt{2}}{32} - \frac{3\pi\sqrt{2}}{32} \\
 &= 1-\frac{\pi}{8}-\frac{4\pi\sqrt{2}}{32} \\ 
&= \frac{8}{8}-\frac{\pi}{8}-\frac{\pi\sqrt{2}}{8}.
\end{align*}

Finalement :

\boxed{\int_0^1 \frac{x^6}{(x^2+1)(x^4+1)}\mathrm{d}x = \frac{8 - \pi(1+\sqrt{2})}{8}.}

257. L’intégrale de Gauss (2/2)

Cet article constitue le prolongement du contenu écrit dans l'article 256.

Pour rappel, il a été établi dans le document précité les résultats suivants.

L’intégrale de Gauss $I = \int_{0}^{+\infty} \e^{-x^2}\dx$ est un réel strictement positif.

De plus :

\begin{align*}
\forall n\in\N, I_{2n} &= \int_{0}^{+\infty} x^{2n}\e^{-x^2}\dx = \frac{(2n)!}{n! \times 2^{2n}}\times I\\
\forall n\in\N, I_{2n+1} &= \int_{0}^{+\infty} x^{2n+1}\e^{-x^2}\dx = \frac{n!}{2}.
\end{align*}

Trouvez une inégalité

Soient $n\in\N$ et $t\in\R$ fixés.

Le contenu écrit dans l'article 256 a montré que les intégrales $\int_{0}^{+\infty} x^{n+2}\e^{-x^2}\dx$, $\int_{0}^{+\infty} x^{n+1}\e^{-x^2}\dx$ et $\int_{0}^{+\infty} x^{n}\e^{-x^2}\dx$ sont convergentes.

Par suite :

\begin{align*}
\lim_{M\to +\infty} \int_{0}^{M} x^{n+2}\e^{-x^2}\dx &= \int_{0}^{+\infty} x^{n+2}\e^{-x^2}\dx  = I_{n+2}\\
\lim_{M\to +\infty} \int_{0}^{M} x^{n+1}\e^{-x^2}\dx &= \int_{0}^{+\infty} x^{n+1}\e^{-x^2}\dx = I_{n+1}\\
\lim_{M\to +\infty} \int_{0}^{M} x^{n}\e^{-x^2}\dx &= \int_{0}^{+\infty} x^{n}\e^{-x^2}\dx = I_{n}.
\end{align*}

D’autre part, pour tout réel $M>0$ :

\begin{align*}
\int_{0}^{M} x^n(x+t)^2\e^{-x^2}\dx &= \int_{0}^{M} x^n(x^2+2tx+t^2)\e^{-x^2}\dx \\
&= \int_{0}^{M} (x^{n+2}+2tx^{n+1}+t^2x^n)\e^{-x^2}\dx \\
&= \int_{0}^{M} x^{n+2}\e^{-x^2}\dx + 2t  \int_{0}^{M} x^{n+1}\e^{-x^2}\dx + t^2  \int_{0}^{M} x^{n}\e^{-x^2}\dx.
\end{align*}

En faisant tendre $M$ vers $+\infty$, il vient, par addition de limites :

\begin{align*}
\lim_{M \to +\infty} \int_{0}^{M} x^n(x+t)^2\e^{-x^2}\dx &= \int_{0}^{+\infty} x^{n+2}\e^{-x^2}\dx + 2t  \int_{0}^{+\infty} x^{n+1}\e^{-x^2}\dx + t^2  \int_{0}^{+\infty} x^{n}\e^{-x^2}\dx \\
&=I_{n+2}+2tI_{n+1}+t^2I_n\\
&= I_n\ t^2+2I_{n+1}\ t+I_{n+2.}
\end{align*}

La fonction $x\mapsto x^n(x+t)^2\e^{-x^2}$ est positive, donc :

\forall M>0,  \int_{0}^{M} x^n(x+t)^2\e^{-x^2}\dx \geq 0.

En faisant tendre $M$ vers $+\infty$ une nouvelle fois, vous obtenez $I_n\ t^2+2I_{n+1}\ t+I_{n+2} \geq 0.$

En définitive :

\boxed{\forall n\in\N, \forall t\in\R, I_n\ t^2+2I_{n+1}\ t+I_{n+2} \geq 0.}

Trouvez une autre inégalité

Soit $n\in\N.$

La fonction $t\mapsto I_n\ t^2+2I_{n+1}\ t+I_{n+2}$ est bien un trinôme du second degré ($I_n$ est non nul) qui est toujours positif.

Par conséquent, le discriminant de ce trinôme est négatif ou nul. Cela donne :

\begin{align*}
(2I_{n+1})^2-4I_n\ I_{n+2} &\leq 0 \\
4I_{n+1}^2-4I_n\ I_{n+2} &\leq 0 \\
I_{n+1}^2-I_n\ I_{n+2} &\leq 0 \\
I_{n+1}^2\leq I_n\ I_{n+2}.
\end{align*}

En définitive :

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

Déterminez la valeur de l’intégrale de Gauss

Première partie

Soit maintenant $n\in\NN.$

D’après le paragraphe précédent :

\begin{align*}
I_{2n}^2&\leq I_{2n+1}\ I_{2n-1} \\
\left( \frac{(2n)!}{n! \times 2^{2n}}\right)^2\times I^2 &\leq  \frac{n!}{2} \times  \frac{(n-1)!}{2}\\
\frac{[(2n)!]^2}{[n!]^2\times 2^{4n}}\times I^2 &\leq \frac{[n!]^2}{4n} \\
 I^2 &\leq \frac{[n!]^4\times 2^{4n}}{[(2n)!]^2\times4n}\\
I^2 &\leq \frac{[n!]^4\times 2^{4n-2}}{[(2n)!]^2\times n}.
\end{align*}

Par conséquent :

\forall n\in\NN, I \leq \frac{[n!]^2\times 2^{2n-1}}{(2n)!\times \sqrt{n}}.

Vous allez maintenant déterminer la limite de $\frac{[n!]^2\times 2^{2n-1}}{(2n)!\times \sqrt{n}}$ quand $n$ tend vers $+\infty.$

Grâce à la formule de Stirling, dont une démonstration se trouve dans les contenus écrits dans l'article 255 et dans l'article 254, vous déduisez :

\begin{align*}
 n! &\underset{n \to+\infty}{\sim} \sqrt{2 \pi n}\left(\frac{n}{\e}\right)^n \\
[n!]^2 &\underset{n \to+\infty}{\sim} 2 \pi n\left(\frac{n}{\e}\right)^{2n}\\
 &\underset{n \to+\infty}{\sim} 2 \pi n^{2n+1}\e^{-2n}\\
(2n)! &\underset{n \to+\infty}{\sim} \sqrt{4 \pi n}\left(\frac{2n}{\e}\right)^{2n} \\
&\underset{n \to+\infty}{\sim} 2\sqrt{ \pi}\sqrt{ n}\ 2^{2n}n^{2n}\e^{-2n} \\
&\underset{n \to+\infty}{\sim} \sqrt{ \pi}\sqrt{ n}\ 2^{2n+1}n^{2n}\e^{-2n} \\
\frac{[n!]^2\times 2^{2n-1}}{(2n)!\times \sqrt{n}} &\underset{n \to+\infty}{\sim}  \frac{2 \pi n^{2n+1}\e^{-2n} \times 2^{2n-1}}{\sqrt{ \pi}\sqrt{ n}\ 2^{2n+1}n^{2n}\e^{-2n} \times \sqrt{n}} \\
&\underset{n \to+\infty}{\sim}  \frac{2 \pi n^{2n+1} \times 2^{2n-1}}{\sqrt{ \pi}n\ 2^{2n+1}n^{2n}} \\
&\underset{n \to+\infty}{\sim}  \frac{2 \sqrt{\pi} \sqrt{\pi} \times 2^{2n-1}}{\sqrt{ \pi}\ 2^{2n+1}} \\
&\underset{n \to+\infty}{\sim}  \frac{ \sqrt{\pi}  \times 2^{2n}}{2^{2n+1}} \\
&\underset{n \to+\infty}{\sim}  \frac{ \sqrt{\pi}  }{2}.
\end{align*}

Ainsi :

\lim_{n\to +\infty} \frac{[n!]^2\times 2^{2n-1}}{(2n)!\times \sqrt{n}} = \frac{\sqrt{\pi}}{2}.

En faisant tendre $n$ vers $+\infty$ dans l’inégalité $I\leq \frac{[n!]^2\times 2^{2n-1}}{(2n)!\times \sqrt{n}}$ vous déduisez $\boxed{I \leq \frac{\sqrt{\pi}}{2}.}$

Seconde partie

Soit maintenant $n\in\NN.$

\begin{align*}
I_{2n+1}^2 &\leq I_{2n}\ I_{2n+2} \\
\left(\frac{n!}{2}\right)^2 &\leq  \left(\frac{(2n)!}{n! \times 2^{2n}}\times I\right)\left( \frac{(2n+2)!}{(n+1)! \times 2^{2n+2}}\times I\right) \\
\frac{[n!]^2}{4} &\leq \frac{[(2n)!]^2 (2n+1)(2n+2)}{[n!]^2(n+1) \times 2^{4n+2}}\times I^2 \\
\frac{[n!]^4 (n+1) \times 2^{4n+2}}{4 [(2n)!]^2 (2n+1)(2n+2)}&\leq I^2\\
\frac{[n!]^2 \sqrt{n+1} \times 2^{2n+1}}{2\times (2n)!  \sqrt{(2n+1)(2n+2)}}&\leq I \\
\frac{[n!]^2 \sqrt{n+1} \times 2^{2n}}{ (2n)!  \sqrt{(2n+1)(2n+2)}}&\leq I \\
\frac{[n!]^2  \times 2^{2n-1}}{ (2n)!  \sqrt{n} } \times \frac{2\sqrt{n}\sqrt{n+1}}{\sqrt{(2n+1)(2n+2)}} &\leq I.
\end{align*}

Par conséquent :

\forall n\in\NN, \frac{[n!]^2  \times 2^{2n-1}}{ (2n)!  \sqrt{n} } \times \frac{2\sqrt{n}\sqrt{n+1}}{\sqrt{(2n+1)(2n+2)}} \leq I.

Pour tout $n\in\NN$ :

\begin{align*}
 \frac{2\sqrt{n}\sqrt{n+1}}{\sqrt{(2n+1)(2n+2)}} &=  \sqrt{\frac{4n(n+1)}{(2n+1)(2n+2)}} \\
&=  \sqrt{\frac{4n^2+4n}{4n^2+6n+2}} \\
&=  \sqrt{\frac{4+\frac{4}{n}}{4+\frac{6}{n}+\frac{2}{n^2}}} \\
\end{align*}

Comme :

\begin{align*}
\lim_{n\to +\infty} 4+\frac{4}{n} &= 4 \\
\lim_{n\to +\infty} 4+\frac{6}{n}+\frac{2}{n^2} &= 4 \\
\end{align*}

Par quotient, il vient :

\lim_{n\to +\infty} \frac{4+\frac{4}{n}}{4+\frac{6}{n}+\frac{2}{n^2}} = 1.

En prenant la racine carrée, vous déduisez :

\begin{align*}
\lim_{n\to +\infty} \sqrt{\frac{4+\frac{4}{n}}{4+\frac{6}{n}+\frac{2}{n^2}}} &= \sqrt{1} =  1\\
\lim_{n\to +\infty} \frac{2\sqrt{n}\sqrt{n+1}}{\sqrt{(2n+1)(2n+2)}}  &= 1.
\end{align*}

Il a été vu dans la première partie que :

\lim_{n\to +\infty} \frac{[n!]^2\times 2^{2n-1}}{(2n)!\times \sqrt{n}} = \frac{\sqrt{\pi}}{2}.

Par produit, vous déduisez :

\lim_{n\to +\infty }\frac{[n!]^2  \times 2^{2n-1}}{ (2n)!  \sqrt{n} } \times \frac{2\sqrt{n}\sqrt{n+1}}{\sqrt{(2n+1)(2n+2)}} = \frac{\sqrt{\pi}}{2}.

En faisant tendre $n$ vers $+\infty$ dans l’inégalité $\frac{[n!]^2 \times 2^{2n-1}}{ (2n)! \sqrt{n} } \times \frac{2\sqrt{n}\sqrt{n+1}}{\sqrt{(2n+1)(2n+2)}} \leq I$ vous déduisez $\boxed{\frac{\sqrt{\pi}}{2} \leq I.}$

Concluez

Les inégalités $\frac{\sqrt{\pi}}{2} \leq I$ et $I \leq\frac{\sqrt{\pi}}{2}$ fournissent $I =\frac{\sqrt{\pi}}{2}.$

L’intégrale de Gauss $I = \int_{0}^{+\infty} \e^{-x^2}\dx$ est égale à $\frac{\sqrt{\pi }}{2}$ :

\boxed{ \int_{0}^{+\infty} \e^{-x^2}\dx = \frac{\sqrt{\pi }}{2}.}

256. L’intégrale de Gauss (1/2)

L’objectif de cet article et du contenu que vous trouverez dans l'article 257 est de déterminer la valeur de l’intégrale de Gauss définie par :

I = \int_{0}^{+\infty} \e^{-x^2}\dx.

Pour y parvenir, vous définissez la suite d’intégrales suivante en posant :

\forall n\in\N, I_n =  \int_{0}^{+\infty} x^n\e^{-x^2}\dx.

Remarquez que pour tout entier naturel $n$, la fonction $x\mapsto x^n\e^{-x^2}$ est positive sur l’intervalle $[0,+\infty[.$

Par conséquent, l’intégrale $I_n$ est bien définie. Elle est égale à ce stade à un réel positif, ou bien à $+\infty.$

Etablissez la convergence de l’intégrale $I_n$

Soit $n$ un entier naturel.

Pour justifier que l’intégrale $I_n$ n’est pas égale à $+\infty$ vous pouvez utiliser des majorations.

Partez du fait que l’exponentielle domine tous les polynômes de n’importe quel degré, en particulier le degré $n+2$, ce qui s’écrit ainsi $\lim_{x\to +\infty} \frac{\e^x}{x^{n+2}} = +\infty.$ Vous déduisez que $\lim_{x\to +\infty} \frac{x^{n+2}}{\e^x} = 0$ ce qui donne $\lim_{x\to +\infty} x^{n+2}\e^{-x} = 0.$

Par conséquent, il existe un réel $A>0$ tel que:

\begin{align*}
\forall x\geq A &, x^{n+2}\e^{-x}\leq 1 \\
\forall x\geq A &, x^{n}\e^{-x}\leq \frac{1}{x^2}.
\end{align*}

Posez $B = A+1.$

Soit maintenant $x$ un réel tel que $x\geq B.$ Comme $x\geq 1$ il vient $x^2\geq x$ donc $-x^2 \leq -x$ et $\e^{-x^2} \leq \e^{-x}.$ Vous déduisez que:

\forall x\geq B, x^n\e^{-x^2}\leq x^n\e^{-x} \leq \frac{1}{x^2}.

Il en résulte que:

\begin{align*}
\int_B^{+\infty} x^n\e^{-x^2}\dx &\leq \int_B^{+\infty} \frac{1}{x^2}\dx\\
&\leq  \left[\frac{-1}{x}\right]_B^{+\infty}\\
&\leq \frac{1}{B}.
\end{align*}

Du coup, pour l’intégrale $I_n$ vous déduisez que :

\begin{align*}
I_n &\leq \int_0^{B} x^n\e^{-x^2}\dx + \int_B^{+\infty} x^n\e^{-x^2}\dx \\
&\leq  \int_0^{B} x^n\e^{-x^2}\dx + \frac{1}{B}.
\end{align*}

L’intégration de la fonction continue $x\mapsto x^n\e^{-x^2}$ sur le segment $[0,B]$ fournit $\int_0^{B} x^n\e^{-x^2}\dx\in \R.$

Il en résulte que $\int_0^{B} x^n\e^{-x^2}\dx + \frac{1}{B} \in \R.$

On ne peut donc avoir $I_n = +\infty$, sinon l’inégalité $I_n \leq \int_0^{B} x^n\e^{-x^2}\dx + \frac{1}{B}$ ne serait pas vérifiée.

Vous en tirez ceci :

\int_{0}^{+\infty} x^n\e^{-x^2}\dx \in \R_{+}.

D’autre part, la fonction $x\mapsto x^n\e^{-x^2}$ est positive, continue et non identiquement nulle sur l’intervalle $[0,1].$ Vous déduisez que :

\int_{0}^{+\infty} x^n\e^{-x^2}\dx \in \R_{+} \geq \int_{0}^{1} x^n\e^{-x^2}\dx > 0.

En définitive, vous avez montré que l’intégrale $I_n$ est un réel strictement positif.

\boxed{\forall n\in\N, \int_{0}^{+\infty} x^n\e^{-x^2}\dx \in \R_{+}^{*}.}

Remarque. En particulier, l’intégrale de Gauss notée $I$, est égale à $I_0.$ C’est donc bien un réel strictement positif aussi.

Trouvez une relation de récurrence

En vous inspirant des intégrations par parties effectuées dans les intégrales de Wallis, vous effectuez le même raisonnement.

Soit $n$ un entier naturel et $M$ un réel supérieur ou égal à $1.$

\begin{align*}
 \int_{0}^{M} x^{n}\e^{-x^2}\dx &= \left[\frac{x^{n+1}}{n+1}\e^{-x^2}\right]_0^M - \int_0^M \frac{x^{n+1}}{n+1}(-2x)\e^{-x^2}\dx \\
&= \frac{M^{n+1}}{n+1}\e^{-M^2} -0 +\frac{2}{n+1} \int_0^M x^{n+2}\e^{-x^2}\dx\\
&= M\times\frac{M^{n}\e^{-M^2}}{n+1} +\frac{2}{n+1} \int_0^M x^{n+2}\e^{-x^2}\dx\\
\end{align*}

Vous observez que :

0\leq M^{n}\e^{-M^2} \leq M^n\e^{-M}.

Via $\lim_{M\to +\infty} M^n\e^{-M} = 0$ et le théorème des gendarmes, vous déduisez $\lim_{M\to +\infty} M^n\e^{-M^2} = 0.$

D’autre part, $\lim_{M\to + \infty} \int_{0}^{M} x^{n}\e^{-x^2}\dx = I_n$ et $\lim_{M\to + \infty} \int_{0}^{M} x^{n+2}\e^{-x^2}\dx = I_{n+2}$ du coup en passant à la limite quand $M\to +\infty$ il vient $I_n = \frac{2}{n+1}I_{n+2}.$

Vous aboutissez à la relation de récurrence suivante :

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

Pour tout $n\in\N$, calculez l’intégrale $I_{2n+1}$

Vous avez :

I_1 =  \int_{0}^{+\infty} x\e^{-x^2}\dx.

Soit $M$ un réel strictement positif. Effectuez le changement de variable suivant :

\begin{align*}
y&=x^2\\
 \dy &=2x\dx\\
\int_0^M x\e^{-x^2}\dx &= \int_0^{{M^2}}\frac{1}{2}\e^{-y}\dy\\
&= \frac{1}{2} \int_0^{{M^2}}\e^{-y}\dy\\
&=\frac{1}{2} \left[-\e^{-y}\right]_0^{{M^2}} \\
&=-\frac{1}{2}\e^{-M^2} + \frac{1}{2}.
\end{align*}

Vous trouvez :

\begin{align*}
\lim_{M\to +\infty} -M^2 &= -\infty \\
\lim_{M\to +\infty} \e^{-{M^2}} &= 0.
\end{align*}

Il en résulte que :

I_1 = \int_{0}^{+\infty} x\e^{-x^2}\dx = \frac{1}{2}.

Utilisant la relation de récurrence sur la suite $(I_n)_{n\geq 0}$ vous déduisez :

\begin{align*}
I_3 &= \frac{2}{2}\times I_1 = \frac{2}{2}\times \frac{1}{2} =  \frac{1!}{2} \\
I_5 &= \frac{4}{2}\times I_3 = \frac{4}{2}\times  \frac{2}{2}\times \frac{1}{2} =  \frac{2!}{2} \\
I_7 &= \frac{6}{2}\times I_5 =  \frac{6}{2}\times \frac{4}{2}\times  \frac{2}{2}\times \frac{1}{2} =  \frac{3!}{2}.
\end{align*}

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

Initialisation. Pour $n=0$, vous trouvez $\frac{0!}{2} = \frac{1}{2} = I_1$ donc $\mathscr{P}(0)$ est vérifiée.

Hérédité. Soit $n\in\N.$ Supposez que $\mathscr{P}(n)$ est vérifiée.

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

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

Vous venez de montrer par récurrence sur $n$ que :

\boxed{\forall n\in\N, I_{2n+1} = \int_{0}^{+\infty} x^{2n+1}\e^{-x^2}\dx = \frac{n!}{2}.}

Pour tout $n\in\N$, calculez l’intégrale $I_{2n}$

Vous allez effectuer les calculs en fonction de l’intégrale de Gauss : $I = \int_{0}^{+\infty} \e^{-x^2}\dx.$ La valeur définitive de cette intégrale sera déterminée dans l'article 257.

Partez de la relation de récurrence $\forall n\in\N, I_{n+2} = \frac{n+1}{2}\times I_n.$

Vous obtenez successivement :

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

Ces calculs préliminaires étant effectués, vous êtes prêts à généraliser le tout en lançant une récurrence.

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

Initialisation. Pour $n=0$, $I_0 = I.$ D’autre part $ \frac{(2\times 0)!}{0! \times 2^{2\times 0}}\times I = 1\times I = I.$

Donc $\mathscr{P}(0)$ est vérifiée.

Hérédité. Soit $n$ un entier naturel. Supposez $\mathscr{P}(n)$.

Vous avez alors :

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

Ainsi avec $\boxed{I = \int_{0}^{+\infty} \e^{-x^2}\dx}$ vous avez :

\boxed{\forall n\in\N, I_{2n} = \int_{0}^{+\infty} x^{2n}\e^{-x^2}\dx = \frac{(2n)!}{n! \times 2^{2n}}\times I.}

255. La formule de Stirling (2/2)

Dans le prolongement du contenu rédigé dans l'article 254, vous allez commencer par une majoration en suivant une démarche similaire.

Etablissez une majoration d’une intégrale

Soit $n$ un entier naturel non nul. Vous souhaitez majorer l’intégrale $\int_{n}^{n+1} \frac{\dt}{t}$ qui est égale à $\ln(n+1)-\ln n.$

Considérez la figure suivante, dans laquelle la courbe $\mathscr{C}$ représente la fonction $x\mapsto \frac{1}{x}.$

04/08/2022 - Capture decran 2022 08 04 a 15.12.03

Dans cette figure, le point $A$ admet pour coordonnées $(n,0).$

Le point $B$ admet pour coordonnées $(n+1,0).$

Le point $C$ situé sur la courbe $\mathscr{C}$ admet pour coordonnées $\left(n+1,\frac{1}{n+1}\right).$

Le point $D$ situé sur la courbe $\mathscr{C}$ admet pour coordonnées $\left(n,\frac{1}{n}\right).$

La convexité de la fonction $x\mapsto \frac{1}{x}$ sur l’intervalle $]0,+\infty[$ permet de déduire que l’intégrale $\int_{n}^{n+1} \frac{\dt}{t}$ est majorée par l’aire du trapèze rectangle $ABCD.$

Des égalités : $AD = \frac{1}{n}$ et $BC = \frac{1}{n+1}$ vous déduisez:

\begin{align*}
\int_{n}^{n+1} \frac{\dt}{t} &\leq \frac{(AD+BC)\times 1}{2} \\
\ln(n+1)-\ln n &\leq  \frac{\frac{1}{n}+\frac{1}{n+1}}{2}\\
\ln\left(\frac{n+1}{n}\right) &\leq  \frac{\frac{1}{n}+\frac{1}{n+1}}{2}\\
\ln\left(\frac{n+1}{n}\right) &\leq  \frac{\frac{2n+1}{n(n+1)}}{2}\\
\ln\left(\frac{n+1}{n}\right) &\leq\frac{2n+1}{2n(n+1)}\\
\frac{n+1}{n} &\leq \e^{\frac{2n+1}{2n(n+1)}}.
\end{align*}

Déduisez-en un encadrement de $n!$

En reprenant les mêmes notations que celles se trouvant dans l'article 254, vous posez :

\forall n\in\NN, u_n = n! \e^n n^{-n-\frac{1}{2}}.

Il a été vu dans l'article 254 que la suite $(u_n)_{n\geq 1}$ est convergente.

Soit $n\in\NN.$ Il vient :

\begin{align*}
\frac{u_n}{u_{n+1}} &= \frac{ n! \e^n n^{-n-\frac{1}{2}}}{ (n+1)! \e^{n+1} (n+1)^{-n-1-\frac{1}{2}}}\\
&= \frac{n^{-n-\frac{1}{2}}}{(n+1)\e (n+1)^{-n-1-\frac{1}{2}}}\\
&= \frac{n^{-n-\frac{1}{2}}}{\e (n+1)^{-n-\frac{1}{2}}}\\
&= \frac{1}{\e}\times \left(\frac{ n }{ n+1 }\right)^{-n-\frac{1}{2}}\\
&= \frac{1}{\e}\times \left(\frac{ n+1 }{ n }\right)^{n+\frac{1}{2}}.
\end{align*}

Comme $\frac{n+1}{n} \geq 1$ et comme $n+\frac{1}{2}\geq 0$, vous déduisez :

\begin{align*}
\frac{n+1}{n} &\leq \e^{\frac{2n+1}{2n(n+1)}} \\
\left(\frac{n+1}{n}\right)^{n+\frac{1}{2}} &\leq \left(\e^{\frac{2n+1}{2n(n+1)}}\right)^{n+\frac{1}{2}}\\
\left(\frac{n+1}{n}\right)^{n+\frac{1}{2}} &\leq  \left(\e^{\frac{2n+1}{2n(n+1)}}\right)^{\frac{2n+1}{2}}\\
\left(\frac{n+1}{n}\right)^{n+\frac{1}{2}} &\leq \e^{\frac{(2n+1)^2}{4n(n+1)}}\\
\frac{1}{\e}\times \left(\frac{n+1}{n}\right)^{n+\frac{1}{2}} &\leq \e^{\frac{(2n+1)^2}{4n(n+1)}}\times \e^{-1}\\
\frac{u_n}{u_{n+1}} &\leq \e^{\frac{(2n+1)^2}{4n(n+1)}-1}\\
\frac{u_n}{u_{n+1}} &\leq \e^{\frac{(2n+1)^2-4n(n+1)}{4n(n+1)}}\\
\frac{u_n}{u_{n+1}} &\leq \e^{\frac{4n^2+4n+1-4n^2-4n}{4n(n+1)}}\\
\frac{u_n}{u_{n+1}} &\leq \e^{\frac{1}{4n(n+1)}}\\
\frac{u_n}{u_{n+1}} &\leq \e^{\frac{1}{4n}-\frac{1}{4(n+1)}}\\
\frac{u_n}{u_{n+1}} \times \e^{\frac{1}{4(n+1)}-\frac{1}{4n}} &\leq 1\\
\frac{u_n\times \e^{\frac{-1}{4n}}}{u_{n+1}\times \e^{\frac{-1}{4(n+1)}}} &\leq 1.
\end{align*}

Par conséquent la suite $\left(u_n\times \e^{\frac{-1}{4n}} \right)_{n\geq 1}$ est croissante.

Or, la suite $(u_n)_{n\geq 1}$ est convergente vers un réel $\ell.$ Comme $\lim_{n\to +\infty} \e^{\frac{-1}{4n}} = 1$ vous déduisez que la suite $\left(u_n\times \e^{\frac{-1}{4n}} \right)_{n\geq 1}$ converge aussi vers $\ell.$

La suite $(u_n)_{n\geq 1}$ étant décroissante et convergente vers $\ell$ vous déduisez $\forall n\in\NN, u_n \geq \ell.$

La suite $\left(u_n\times \e^{\frac{-1}{4n}} \right)_{n\geq 1}$ étant croissante et convergente vers $\ell$ vous déduisez $\forall n\in\NN, u_n\times \e^{\frac{-1}{4n}} \leq \ell.$

Vous avez montré qu’il existe un réel $\ell$ tel que, pour tout $n\in\NN$ :

\begin{align*}
n! \e^n n^{-n-\frac{1}{2}} \times  \e^{\frac{-1}{4n}} &\leq \ell\leq n! \e^n n^{-n-\frac{1}{2}}\\
n! \e^n n^{-n} \times  \e^{\frac{-1}{4n}} &\leq \ell \sqrt{n }\leq n! \e^n n^{-n}\\
n! \e^n  \times  \e^{\frac{-1}{4n}} &\leq \ell \sqrt{n}\ n^n\leq n! \e^n \\
n!  \times  \e^{\frac{-1}{4n}} &\leq \ell \sqrt{n}\ \left(\frac{n}{\e}\right)^n\leq n! \\
 \ell \sqrt{n}\ \left(\frac{n}{\e}\right)^n&\leq n! \leq \ell \sqrt{n}\ \left(\frac{n}{\e}\right)^n \times \e^{\frac{1}{4n}}.
\end{align*}

En prenant $n=1$, vous déduisez l’encadrement suivant :

\begin{align*}
\forall n\in\NN, \e\times  \e^{\frac{-1}{4}} \leq \ell\leq \e \\
\forall n\in\NN, \e^{\frac{3}{4}} \leq \ell\leq \e.
\end{align*}

En définitive de cette section, vous avez montré qu’il existe un nombre $\ell\in\left[\e^{\frac{3}{4}}, \e\right]$ ($\ell$ est ainsi strictement positif) tel que :

\boxed{\forall n\in\NN, \ell \sqrt{n}\ \left(\frac{n}{\e}\right)^n \leq n! \leq \ell \sqrt{n}\ \left(\frac{n}{\e}\right)^n \times \e^{\frac{1}{4n}}.}

Calculez le nombre $\ell$

De l’encadrement précédent, vous déduisez que :

\forall n\in\NN, \ell \leq \frac{n!}{\sqrt{n}\ \left(\frac{n}{\e}\right)^n } \leq \ell \times \e^{\frac{1}{4n}}.

Comme $\lim_{n\to +\infty} \e^{\frac{1}{4n}} = 1$, il vient, d’après le théorème des gendarmes :

\lim_{n\to +\infty} \frac{n!}{\sqrt{n}\left(\frac{n}{\e}\right)^n} = \ell.

D’après le contenu rédigé dans l'article 253 le nombre $\pi$ est obtenu par la limite de la suite suivante définie avec des factorielles (appelée formule de Wallis) :

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

Soit maintenant $n$ un entier naturel non nul. Vous souhaitez forcer l’apparition du nombre $\ell$ donc vous préparez le tout en rajoutant des éléments inspirés de la limite donnant $\ell$ :

\begin{align*}
\frac{(n!)^4\times 2^{4n+1}}{(2n)!(2n+1)!} &= \left(\frac{n!}{\sqrt{n}\left(\frac{n}{\e}\right)^n}\right)^4 \left(\sqrt{n}\left(\frac{n}{\e}\right)^n\right)^4\times 2^{4n+1}\times \frac{\sqrt{2n} \left(\frac{2n}{\e}\right)^{2n}}{(2n)!}\times \frac{1}{\sqrt{2n} \left(\frac{2n}{\e}\right)^{2n}} \times \frac{\sqrt{2n+1} \left(\frac{2n+1}{\e}\right)^{2n+1}}{(2n+1)!}\times \frac{1}{\sqrt{2n+1} \left(\frac{2n+1}{\e}\right)^{2n+1}}. 
\end{align*}

Vous avez :

\begin{align*}
\lim_{n\to +\infty} \left(\frac{n!}{\sqrt{n}\left(\frac{n}{\e}\right)^n}\right)^4 = \ell^4.
\end{align*}

Comme $\ell$ est non nul, vous déduisez :

\begin{align*}
\lim_{n\to +\infty }\frac{\sqrt{2n} \left(\frac{2n}{\e}\right)^{2n}}{(2n)!} = \frac{1}{\ell}\\
\lim_{n\to +\infty } \frac{\sqrt{2n+1} \left(\frac{2n+1}{\e}\right)^{2n+1}}{(2n+1)!} = \frac{1}{\ell}.
\end{align*}

Par produit de limites, vous en tirez que :

\lim_{n\to +\infty} \left(\frac{n!}{\sqrt{n}\left(\frac{n}{\e}\right)^n}\right)^4 \times \frac{\sqrt{2n} \left(\frac{2n}{\e}\right)^{2n}}{(2n)!}\times   \frac{\sqrt{2n+1} \left(\frac{2n+1}{\e}\right)^{2n+1}}{(2n+1)!}  = \ell^2.

Soit $n\in\NN$, vous cherchez à simplifier ce qui suit :

\begin{align*}
\left(\sqrt{n}\left(\frac{n}{\e}\right)^n\right)^4\times 2^{4n+1}\times  \frac{1}{\sqrt{2n} \left(\frac{2n}{\e}\right)^{2n}} \times \frac{1}{\sqrt{2n+1} \left(\frac{2n+1}{\e}\right)^{2n+1}} &=\frac{n^2 n^{4n}\e^{-4n}2^{4n+1}}{\sqrt{2n(2n+1)}\ 2^{2n}n^{2n}\e^{-2n}(2n+1)^{2n+1}\e^{-2n-1}} \\
 &=\frac{\e\ n^{2n+2}2^{2n+1}}{\sqrt{4n^2+2n}\ (2n+1)^{2n+1}}\\
&=\frac{\e\ n^{2n+2}2^{2n+1}}{2n\sqrt{\frac{4n^2+2n}{4n^2}}\ (2n+1)^{2n+1}}\\
&=\frac{\e\ n^{2n+1}2^{2n}}{\sqrt{1+\frac{1}{2n}}\ (2n+1)^{2n+1}}\\
&=\frac{\e\ n^{2n+1}2^{2n}}{(2n)^{2n+1}\sqrt{1+\frac{1}{2n}}\left(\frac{2n+1}{2n}\right)^{2n+1}}\\
&=\frac{\e\ n^{2n+1}2^{2n}}{(2n)^{2n+1}\sqrt{1+\frac{1}{2n}}\left(1+\frac{1}{2n}\right)^{2n}\left(1+\frac{1}{2n}\right)}\\
&=\frac{\e\ n^{2n+1}2^{2n}}{2^{2n+1}n^{2n+1}\sqrt{1+\frac{1}{2n}}\left(1+\frac{1}{2n}\right)^{2n}\left(1+\frac{1}{2n}\right)}\\
&=\frac{\e}{2\sqrt{1+\frac{1}{2n}}\left(1+\frac{1}{2n}\right)^{2n}\left(1+\frac{1}{2n}\right)}.
\end{align*}

 

D’après le contenu rédigé dans l'article 251 le nombre $\e$ vérifie :

\lim_{n\to +\infty} \left(1+\frac{1}{2n}\right)^{2n} = \e.

Comme $\lim_{n\to +\infty} 1+\frac{1}{2n} = 1$, vous déduisez :

\lim_{n\to +\infty}\left(\sqrt{n}\left(\frac{n}{\e}\right)^n\right)^4\times 2^{4n+1}\times  \frac{1}{\sqrt{2n} \left(\frac{2n}{\e}\right)^{2n}} \times \frac{1}{\sqrt{2n+1} \left(\frac{2n+1}{\e}\right)^{2n+1}} = \frac{1}{2}.

Du coup :

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

Par unicité de la limite d’une suite, vous déduisez que :

\begin{align*}
\frac{\ell^2}{2} &=\pi\\
\ell^2 &= 2\pi.
\end{align*}

Comme $\ell$ est positif, il vient $\boxed{\ell = \sqrt{2\pi}}.$

Concluez

D’après l’étude de cet article, vous avez établi l’encadrement suivant :

\boxed{\forall n\in\NN,  \sqrt{2\pi n}\ \left(\frac{n}{\e}\right)^n \leq n! \leq  \sqrt{2\pi n}\ \left(\frac{n}{\e}\right)^n \times \e^{\frac{1}{4n}}.}

De cet encadrement, vous déduisez la formule de Stirling :

\boxed{\lim_{n\to +\infty} \frac{n!}{\sqrt{2 \pi n}\left(\frac{n}{\e}\right)^n} = 1.}

254. La formule de Stirling (1/2)

Le but de cette série d’articles est de présenter une approche permettant de faire apparaître la formule de Stirling.

Vous allez partir de la fonction inverse $x\mapsto \frac{1}{x}$ qui est la dérivée de la fonction logarithme népérien $x\mapsto \ln x$ et utiliser ses propriétés pour en déduire une inégalité faisant apparaître une factorielle.

Etablissez une minoration d’une intégrale

Soit $n$ un entier naturel non nul. Vous souhaitez minorer l’intégrale $\int_{n}^{n+1} \frac{\dt}{t}$ qui est égale à $\ln(n+1)-\ln n.$

Considérez la figure suivante, dans laquelle la courbe $\mathscr{C}$ représente la fonction $x\mapsto \frac{1}{x}.$

03/08/2022 - Capture decran 2022 08 03 a 15.33.59

Le point $A$ admet pour coordonnées $(n,0).$

Le point $B$ admet pour coordonnées $(n+1,0).$

Le point $I$, milieu du segment $[AB]$ admet pour coordonnées $\left(n+\frac{1}{2},0\right).$

Le point $E$ situé sur la courbe $\mathscr{C}$ admet pour coordonnées $\left(n,\frac{1}{n}\right).$

Le point $F$ situé sur la courbe $\mathscr{C}$ admet pour coordonnées $\left(n+1,\frac{1}{n+1}\right).$

Le point $J$ situé sur la courbe $\mathscr{C}$ admet pour coordonnées $\left(n+\frac{1}{2},\frac{1}{n+\frac{1}{2}}\right)$ c’est-à-dire $\left(n+\frac{1}{2}, \frac{2}{2n+1}\right).$

La droite $(CD)$ est la tangente à la courbe $\mathscr{C}$ au point $J$ qui a pour abscisse $n+\frac{1}{2}.$

Le coefficient directeur de cette tangente est $m = \frac{-1}{\left(n+\frac{1}{2}\right)^2}.$

L’aire du trapèze $ABCD$ est inférieure ou égale à $\ln(n+1)-\ln n$ étant donné que la fonction $x\mapsto \frac{1}{x^2}$ est convexe sur l’intervalle $]0,+\infty[.$

Vous simplifiez cette écriture en supprimant les fractions empilées:

\begin{align*}
 m &= \frac{-1}{\left(n+\frac{1}{2}\right)^2}\\
&= \frac{-1}{\left(\frac{2n+1}{2}\right)^2}\\
&= \frac{-1}{\frac{(2n+1)^2}{4}}\\
&= \frac{-4}{(2n+1)^2}.
\end{align*}

L’équation réduite de la droite $(CD)$ est $y=mx+p.$ Cette droite passant par $J$ il vient:

\begin{align*}
y &=mx+p \\
\frac{2}{2n+1} &= m\times \left(n+\frac{1}{2}\right)+p\\
\frac{2}{2n+1} &=  \frac{-4}{(2n+1)^2}\times \left(n+\frac{1}{2}\right)+p\\
2(2n+1)&=-4\left(n+\frac{1}{2}\right)+(2n+1)^2\times p\\
4n+2&=-4n-2+(2n+1)^2\times p\\
8n+4&=(2n+1)^2\times p\\
4(2n+1)&=(2n+1)^2\times p\\
4 &=(2n+1)\times p\\
\frac{4}{2n+1}&=p.
\end{align*}

L’ordonnée du point $C$ est égale à:

\begin{align*}
y_C &= m\times (n+1)+p\\
&= \frac{-4}{(2n+1)^2}\times (n+1)+\frac{4}{2n+1}\\
&=\frac{-4n-4}{(2n+1)^2}+\frac{4(2n+1)}{(2n+1)^2}\\
&=\frac{-4n-4+8n+4}{(2n+1)^2}\\
&=\frac{4n}{(2n+1)^2}.
\end{align*}

La distance verticale entre les points $B$ et $C$ est égale à $\boxed{BC = \frac{4n}{(2n+1)^2}.}$

L’ordonnée du point $D$ est égale à:

\begin{align*}
y_D &= m\times n+p\\
&= \frac{-4}{(2n+1)^2}\times n+\frac{4}{2n+1}\\
&=\frac{-4n}{(2n+1)^2}+\frac{4(2n+1)}{(2n+1)^2}\\
&=\frac{-4n+8n+4}{(2n+1)^2}\\
&=\frac{4n+4}{(2n+1)^2}.
\end{align*}

La distance verticale entre les points $A$ et $D$ est égale à $\boxed{AD = \frac{4n+4}{(2n+1)^2}.}$

Vous déduisez de ces calculs les inégalités:

\begin{align*}
\frac{(AD+BC)\times 1}{2} &\leq \ln(n+1)-\ln n \\
AD+BC &\leq 2\left(\ln(n+1)-\ln n\right)\\
\frac{8n+4}{(2n+1)^2} &\leq 2\left(\ln(n+1)-\ln n\right)\\
\frac{4n+2}{(2n+1)^2} &\leq \ln(n+1)-\ln n\\
\frac{2(2n+1)}{(2n+1)^2} &\leq \ln(n+1)-\ln n\\
\frac{2}{2n+1} &\leq \ln(n+1)-\ln n\\
1&\leq \frac{2n+1}{2} \left( \ln(n+1)-\ln n \right)\\
1&\leq \left(n+\frac{1}{2}\right) \ln\frac{n+1}{n}\\
1&\leq \ln\left(\left(\frac{n+1}{n}\right)^{n+\frac{1}{2}}\right)\\
\e &\leq \left(\frac{n+1}{n}\right)^{n+\frac{1}{2}}.
\end{align*}

Ainsi:

\boxed{\forall n\in\NN, \e \leq \left(\frac{n+1}{n}\right)^{n+\frac{1}{2}}.}

Vers la formule de Stirling

Soit $n$ un entier naturel non nul.

Vous déduisez de ce qui précède que:

\begin{align*}
\e &\leq \left(\frac{n+1}{n}\right)^{n+\frac{1}{2}}\\
1 &\leq \frac{1}{\e} \times   \left(\frac{n+1}{n}\right)^{n+\frac{1}{2}}\\
1 &\leq \frac{\e^n}{\e^{n+1}} \times   \left(\frac{n+1}{n}\right)^{n+\frac{1}{2}}\\
1 &\leq \frac{\e^n}{\e^{n+1}} \times   \left(\frac{n}{n+1}\right)^{-n-\frac{1}{2}}\\
1 &\leq \frac{\e^n}{\e^{n+1}} \times  \frac{n^{-n-\frac{1}{2}}}{(n+1)^{-n-\frac{1}{2}}}\\
1 &\leq \frac{\e^n}{\e^{n+1}} \times  \frac{n^{-n-\frac{1}{2}}}{(n+1)^{-n-1-\frac{1}{2}}}\times \frac{1}{n+1}\\
1 &\leq \frac{\e^n}{\e^{n+1}} \times  \frac{n^{-n-\frac{1}{2}}}{(n+1)^{-n-1-\frac{1}{2}}}\times \frac{n!}{(n+1)!}.
\end{align*}

Cela vous invite à poser $\boxed{\forall n\in\NN, u_n = n! \e^n n^{-n-\frac{1}{2}}.}$ La suite $(u_n)_{n\geq 1}$ est strictement positive, elle est minorée par $0.$

D’après le résultat établi ci-dessus, $\forall n\in\NN, 1\leq \frac{u_n}{u_{n+1}}$ donc la suite $(u_n)_{n\geq 1}$ est décroissante.

La suite $(u_n)_{n\geq 1}$ étant décroissante et minorée il en résulte qu’elle converge.

Prolongement

Pour accéder à la formule de Stirling, allez jeter un oeil dans le contenu rédigé dans l'article 255.