NOTE VI.
SUR LA MÉTHODE D’APPROXIMATION TIRÉE DES SÉRIES RÉCURRENTES.
1. Reprenons l’équation
![{\displaystyle x^{m}-\mathrm {A} x^{m-1}+\mathrm {B} x^{m-2}-\mathrm {C} x^{m-3}+\ldots =0,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d4b78efd381566ee5786c1c7f1106e4cf2a9a182)
dont on a désigné les racines par
on aura (Note II), par la nature de ces racines, l’équation identique
![{\displaystyle x^{m}-\mathrm {A} x^{m-1}+\mathrm {B} x^{m-2}-\mathrm {C} x^{m-3}+\ldots =(x-\alpha )(x-\beta )(x-\gamma )(x-\delta )\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ba3d7afde4dd810fd0f05013e521eebea46f993c)
laquelle doit avoir lieu, quelle que soit la valeur de
.
L’identité de l’équation subsistera donc encore en mettant
au lieu de
quelles que soient les valeurs de
et
donc aussi, si après la substitution on développe suivant les puissances de
les termes affectés de
fourniront d’autres équations identiques ; ce seront les équations que nous avons appelées dérivées dans la Théorie des fonctions.
La première de ces équations dérivées sera
![{\displaystyle mx^{m-1}-(m-1)\mathrm {A} x^{m-2}+(m-2)\mathrm {B} x^{m-3}-\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ecee501d7f56525f0567051b7cc3f00c177dc4c)
![{\displaystyle =(x-\beta )(x-\gamma )\ldots +(x-\alpha )(x-\gamma )\ldots +(x-\alpha )(x-\beta )\ldots +\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/40d65b094d3aeb17f190d5d2b8f0681536c5d9d2)
Divisons cette équation par l’équation identique ci-dessus ; on aura
![{\displaystyle {\frac {mx^{m-1}-(m-1)\mathrm {A} x^{m-2}+(m-2)\mathrm {B} x^{m-3}-\ldots }{x^{m}-\mathrm {A} x^{m-1}+\mathrm {B} x^{m-2}-\mathrm {C} x^{m-3}+\ldots }}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4dc1aa98c60eec606fd06a2fdcbe8ced33bca52a)
![{\displaystyle ={\frac {1}{x-\alpha }}+{\frac {1}{x-\beta }}+{\frac {1}{x-\gamma }}+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2d2b944fc7e6e961730e8668d3df12624b3a9a0a)
équation qui doit être aussi identique, quelle que soit la valeur de
. Donc elle le sera encore si l’on en développe les deux membres en séries qui procèdent suivant les puissances positives ou négatives de
.
2. Développons d’abord suivant les puissances négatives ; la fraction qui forme le premier membre deviendra
![{\displaystyle {\frac {\mathrm {P} }{x}}+{\frac {\mathrm {Q} }{x^{2}}}+{\frac {\mathrm {R} }{x^{3}}}+{\frac {\mathrm {S} }{x^{4}}}+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/844e4bc27303c402b7f6a5573732ffe8655c42e9)
et, pour trouver les valeurs des coefficients
il n’y a qu’à multiplier par le dénominateur
et comparer ensuite les termes avec ceux du numérateur
on aura ainsi
![{\displaystyle {\begin{aligned}\mathrm {P} =&m,\\\mathrm {Q} =&\mathrm {AP} -(m-1)\mathrm {A} ,\\\mathrm {R} =&\mathrm {AQ-BP} +(m-2)\mathrm {B} ,\\\mathrm {S} =&\mathrm {AR-BQ+CP} -(m-3)\mathrm {C} ,\\.\ldots &\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6db0a09cad1e7679b2e6440f6af819be728b7a0c)
où l’on voit que la suite des quantités
devient après le
ième terme une suite récurrente, dont l’échelle de relation est ![{\displaystyle \mathrm {A,\ -B,\ C,\ -D} ,\ \ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4c34601df4b8026b918fb9ee290a8d2217eef9ae)
Développant de même les fractions qui forment le second membre, il deviendra
![{\displaystyle {\frac {m}{x}}+(\alpha +\beta +\gamma +\ldots ){\frac {1}{x^{2}}}+\left(\alpha ^{2}+\beta ^{2}+\gamma ^{2}+\ldots \right){\frac {1}{x^{3}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/53ca1685488cca5449d93f75ca86e2ad4f2858fb)
![{\displaystyle +\left(\alpha ^{3}+\beta ^{3}+\gamma ^{3}+\ldots \right){\frac {1}{x^{4}}}+\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8f5332291d81890b71baccf79f82fe1f767b5a9f)
Maintenant, la comparaison des termes semblables des deux membres de l’équation donne
![{\displaystyle {\begin{aligned}\mathrm {P} =&m,\\\mathrm {Q} =&\alpha \ \ +\beta \ \,+\gamma \ \ +\ldots ,\\\mathrm {R} =&\alpha ^{2}+\beta ^{2}+\gamma ^{2}+\ldots ,\\\mathrm {S} =&\alpha ^{3}+\beta ^{3}+\gamma ^{3}+\ldots ,\\.\ldots &\ldots \ldots \ldots \ldots \ldots \ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ee8e062f5ab9c7b1d58e63565ad925c2e33925f2)
et, en général, un terme quelconque ; dont le quantième sera
à compter de
sera égal à
![{\displaystyle \alpha ^{\mu }+\beta ^{\mu }+\gamma ^{\mu }+\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ce180d9425adb70e6dcecd43ab1908b1678a8a26)
C’est l’expression du terme général de la série.
On a par là la démonstration la plus simple de la loi donnée par Newton pour la somme des puissances des racines. Mais les formules précédentes sont surtout utiles pour approcher de la valeur de la plus grande des racines
En effet, il est clair que, si toutes ces racines sont réelles et que
soit par exemple la plus grande des racines, soit qu’elle soit positive ou négative, la puissance
surpassera d’autant plus les puissances semblables des autres racines, et même la somme de ces puissances, que l’exposant
sera plus grand ; d’où il s’ensuit que, si
et
sont des termes consécutifs de la série
on aura à très-peu près
et cette valeur de la racine
sera d’autant plus approchée que les termes dont il s’agit seront plus éloignés du commencement de la série.
3. Si parmi les racines
il y en avait d’imaginaires, on aurait, par exemple,
![{\displaystyle \beta =\varpi +\rho {\sqrt {-1}},\quad \gamma =\varpi -\rho {\sqrt {-1}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/939a7a056d39107e30718453e57576bdb3e0ecb7)
alors, faisant
et
on aurait
![{\displaystyle \beta =\Pi \left(\cos \varphi +\sin \varphi {\sqrt {-1}}\right),\quad \gamma =\Pi \left(\cos \varphi -\sin \varphi {\sqrt {-1}}\right)\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a49bfaeb4533c6bc25bfb9cc53b4a038fb5e1dfe)
donc, par le théorème connu,
![{\displaystyle {\begin{aligned}\beta ^{\mu }=&\Pi ^{\mu }\left(\cos \mu \varphi +\sin \mu \varphi {\sqrt {-1}}\right),\\\gamma ^{\mu }=&\Pi ^{\mu }\left(\cos \mu \varphi -\sin \mu \varphi {\sqrt {-1}}\right),\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8f942bc7aba9ea26e252f0a3bc841b708a77fe95)
et par conséquent
![{\displaystyle \beta ^{\mu }+\gamma ^{\mu }=2\Pi ^{\mu }\cos \mu \varphi .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8f4ad7f1c6ba3d6b3c9acbb7fa763a9b30f67272)
Ainsi, pourvu que la racine
soit en même temps plus grande que
ou
c’est-à-dire plus grande que
la puissance
surpassera aussi la somme de pareilles puissances de
et ![{\displaystyle \gamma .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f423d4c0d1a3f651562797e2198c75a3f65e09fe)
Donc la méthode ne sera en défaut, à cause des racines imaginaires, qu’autant qu’il s’en trouvera dans lesquelles le produit réel des deux racines correspondantes sera plus grand que le carré de la plus grande des racines réelles ; et, dans ce cas, la série, au lieu de s’approcher et de se confondre à la fin avec une série géométrique, s’en éloignera continuellement.
4. Cette méthode rentre évidemment dans celle que Daniel Bernoulli a déduite de la considération des suites récurrentes, et qu’Euler a exposée en détail dans son Introduction. Dans celle-ci, on donne à la fraction génératrice de la série, pour numérateur, un polynôme quelconque d’un degré moindre que le dénominateur, ce qui rend les
premiers termes de la série entièrement arbitraires. Cette fraction se décompose dans les fractions simples
![{\displaystyle {\frac {a}{x-\alpha }}+{\frac {b}{x-\beta }}+{\frac {c}{x-\gamma }},\quad \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9db718825f8e44a6c1d6c51a6034d2772a75cab6)
d’où résulte, pour les termes de la série, cette expression générale
![{\displaystyle a\alpha ^{\mu }+b\beta ^{\mu }+c\gamma ^{\mu }+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ab157c769e9a489b11ba9576016398e398b7d99b)
laquelle donne également, lorsque la racine
est beaucoup plus grande que chacune des autres,
pour la valeur approchée de
quelle que soit la valeur du coefficient
Mais l’indétermination des premiers termes de la série, au lieu d’être un avantage de cette méthode, est plutôt un inconvénient ; car s’il arrive que les deux racines
soient égales, alors les deux termes
prennent en général la forme
![{\displaystyle (a'+b'\mu )\alpha ^{\mu },}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b86dcfacce5adf216f660159f398c2ac478a8e71)
et, si les trois racines
sont égales, les trois termes
prennent la forme
![{\displaystyle \left(a'+b'\mu +c'\mu ^{2}\right)\alpha ^{\mu },}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1b579d235eb994183b154ca5815f8d2b16eba9de)
et ainsi de suite ; d’où il est aisé de voir que, lorsque la plus grande racine
est une racine double ou triple, etc., la série converge bien moins rapidement vers une série géométrique. En prenant pour numérateur la fonction prime du dénominateur, ainsi que nous l’avons fait ci-dessus, tous les coefficients
deviennent égaux à l’unité,
et, dans le cas des racines égales
![{\displaystyle \alpha }](https://wikimedia.org/api/rest_v1/media/math/render/svg/b79333175c8b3f0840bfb4ec41b8072c83ea88d3)
et
![{\displaystyle \beta ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/59ab677d974cccb0132cac08bd67fc8ac765627e)
les deux termes
![{\displaystyle \alpha ^{\mu }+\beta ^{\mu }}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9b0080a87ea8242a7149c669897c9ed5770ebdd7)
deviennent simplement
![{\displaystyle 2\alpha ^{\mu },s}](https://wikimedia.org/api/rest_v1/media/math/render/svg/05f49cacc55f01689529f7ccb17b2a372c24a4fb)
et ainsi des autres ; de sorte que les racines égales n’influent en rien sur la convergence de la série.
5. Pour donner un exemple de ce que nous venons de dire, je prendrai celui de l’article 346 de l’Introduction d’Euler. L’équation à résoudre est
![{\displaystyle x^{3}-3x^{2}+4=0\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/264d28ba609044c89b600d5757aed698c427bbd1)
Euler prend
et
pour les trois premiers termes, et il forme par l’échelle de relation,
la série récurrente
![{\displaystyle 1,\ \ 3,\ \ 9,\ \ 23,\ \ 57,\ \ 135,\ \ 313,\ \ 711,\ \ \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/abcda65e1e1e4813f64e5ea4800b74e5f16f91bb)
dans laquelle il observe que le quotient de chaque terme, divisé par le précédent, est toujours plus grand que
racine double, et en même temps la plus grande.
Si l’on emploie les formules données ci-dessus, en faisant
![{\displaystyle m=3,\quad \mathrm {A} =3,\quad \mathrm {B} =0,\quad \mathrm {C} =4,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/18e495de0125e74357515d9c961b4fae9532ea84)
tous les termes
se trouvent multiples de
de sorte que, rejetant ce facteur pour plus de simplicité, on trouve par la même échelle de relation, mais en partant des termes
la série
![{\displaystyle 1,\ \ 1,\ \ 3,\ \ 5,\ \ 11,\ \ 21,\ \ 43,\ \ 85,\ \ 171,\ \ 341,\ \ \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/55a63f364ad49501810a19fdac5db63af2197347)
où l’on voit que le quotient de chaque terme, divisé par celui qui le précède, converge très-rapidement vers la racine double ![{\displaystyle 2.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9f2b3373a07e65d3312989163b5ebd400af86480)
6. Nous avons développé plus haut l’équation identique
![{\displaystyle {\frac {mx^{m-1}-\ldots }{x^{m}-\mathrm {A} x^{m-1}+\ldots }}={\frac {1}{x-\alpha }}+{\frac {1}{x-\beta }}+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/ea23efa65770ccc3e631a72b7f6bef3b780cbe48)
suivant les puissances négatives de
développons-lâ maintenant suivant les puissances positives. Pour cela, soient
![{\displaystyle a-bx+cx^{2}-dx^{3}+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/8ac7165ea9e17da703f42fbbe82fba0701b60d8b)
les derniers termes du polynôme
![{\displaystyle x^{m}-\mathrm {A} x^{m-1}+\mathrm {B} x^{m-2}-\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fece4c491cc4afdba00452a49345bd556b6b86c0)
on mettra le premier membre de l’équation identique sous la forme
![{\displaystyle {\frac {-b+2cx-3dx^{2}+4ex^{3}-\ldots }{a-bx+cx^{2}-dx^{3}+ex^{4}-\ldots }},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/78265aae1327bb6aab1136c91cddafc0263232ea)
et le développement de cette fraction suivant les puissances croissantes de
sera de la forme
![{\displaystyle -\mathrm {P} '-\mathrm {Q} 'x-\mathrm {R} 'x^{2}-\mathrm {S} 'x^{3}-\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c9ccbcb1dd176847b051119bf3c1249a4b3c6f8d)
en multipliant par le dénominateur et comparant les termes, on trouvera
![{\displaystyle {\begin{aligned}a\mathrm {P} '&=&b,\\a\mathrm {Q} '&=b\mathrm {P} '-2c,\\a\mathrm {R} '&=b\mathrm {Q} '-c\mathrm {P} '+3d,\\a\mathrm {S} '&=b\mathrm {R} '-c\mathrm {Q} '+d\mathrm {P} '-4c,\\.\ldots &\ldots \ldots \ldots \ldots \ldots \ldots \ldots ..,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/933eca8e2ba19bb5ace4c1799fd49a3b0a6c296e)
ce qui donne une série récurrente dont l’échelle est
![{\displaystyle {\frac {b}{a}},\quad -{\frac {c}{a}},\quad {\frac {d}{a}},\quad \ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7ed00cf1c66f8393434b667bbffe060d3d5c5ac6)
Le second membre de la même équation, étant développé pareillement suivant les puissances croissantes de
donnera la série
![{\displaystyle -\left({\frac {1}{\alpha }}+{\frac {1}{\beta }}+{\frac {1}{\gamma }}+\ldots \right)-\left({\frac {1}{\alpha ^{2}}}+{\frac {1}{\beta ^{2}}}+{\frac {1}{\gamma ^{2}}}+\ldots \right)x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/527f98bbc2d2336ab14372a8f21cc22875b85b9a)
![{\displaystyle -\left({\frac {1}{\alpha ^{3}}}+{\frac {1}{\beta ^{3}}}+{\frac {1}{\gamma ^{3}}}+\ldots \right)x^{2}-\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d07befb34e87d30dce7598780e044388b91c3ff8)
de sorte qu’on aura par la comparaison
![{\displaystyle {\begin{aligned}{\frac {1}{\ \alpha \ }}+{\frac {1}{\ \beta \ }}+{\frac {1}{\ \gamma \ }}+\ldots &=\mathrm {P} ',\\{\frac {1}{\alpha ^{2}}}+{\frac {1}{\beta ^{2}}}+{\frac {1}{\gamma ^{2}}}+\ldots &=\mathrm {Q} ',\\{\frac {1}{\alpha ^{3}}}+{\frac {1}{\beta ^{3}}}+{\frac {1}{\gamma ^{3}}}+\ldots &=\mathrm {R} ',\\\ldots \ldots \ldots \ldots \ldots \ldots \ldots &\ldots ...\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/268bbe3323367a4f9acf95b9fab26dfef06d4d1f)
Ces formules renferment la loi des sommes des puissances réciproques des racines.
Il est évident que, si
est la plus petite racine, soit positive ou négative, les puissances
surpasseront d’autant plus la somme des pareilles puissances des autres racines, que
sera plus petite que chacune des autres racines
Par conséquent, si
et
sont deux termes consécutifs de la série
le quotient
approchera d’autant plus de la valeur de la plus petite racine réelle de l’équation, que ces termes seront plus éloignés du commencement de la série. Ainsi cette série servira à trouver la plus petite racine, comme la première
sert à trouver la plus grande ; et, à l’égard des racines imaginaires, on prouvera de la même manière qu’elles n’empêcheront pas l’approximation vers la plus petite racine réelle, pourvu que le carré de cette racine soit en même temps plus petit que chacun des produits réels des racines imaginaires correspondantes.
7. On pourrait donc employer cette méthode d’approximation pour chacune des racines réelles d’une équation quelconque si l’on connaissait d’avance une valeur approchée
de cette racine, telle que la différence entre cette valeur et la vraie valeur de la racine fût moindre en quantité, c’est-à-dire abstraction faite des signes, que la différence entre la même valeur et chacune des autres racines réelles, et en même temps moindre que la racine carrée de chacun des produits des racines imaginaires correspondantes, s’il y en
diminuées de la même valeur ; car alors, en nommant
la valeur approchée de la racine cherchée et faisant
on aura une transformée en
dont la plus petite racine pourra se déterminer par la méthode précédente, et cette racine, jointe à la première valeur approchée, donnera la racine cherchée. Mais on ne saurait trouver les premières valeurs qu’en faisant usage des méthodes que nous avons données, et, ces valeurs étant une fois connues, il est bien plus exact d’employer la méthode d’approximation du Chapitre III ; aussi ne suis-je entré dans ce détail sur la méthode d’approximation tirée des séries récurrentes que pour ne rien laisser à désirer sur le sujet dont il s’agit.
8. Si l’on veut appliquer la méthode précédente à l’exemple de Newton, on prendra d’abord la transformée (Note précédente)
![{\displaystyle p^{3}+6p^{2}+10p-1\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/35b691558ccdc240999d48f70fef738d6b9f4652)
et, comme on sait que la racine réelle est moindre que
il s’ensuit que le produit des deux autres racines, qu’on sait être imaginaires, sera
puisque le dernier terme
est le produit des trois racines ; ainsi l’on est assuré que le carré de la racine cherchée est beaucoup moindre que le produit des deux racines imaginaires. On formera donc la série récurrente par le moyen de la fraction
et l’on aura les termes
![{\displaystyle 10,\ \ 112,\ \ 1183,\ \ 12512,\ \ 132330,\ \ \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/faa6c843e1201f9bf24288e595dc18feb13b11f8)
qu’on peut continuer aussi loin qu’on veut par l’échelle de relation
chacun de ces termes, divisé par le suivant, donnera les fractions
![{\displaystyle {\frac {10}{112}},\quad {\frac {112}{1183}},\quad {\frac {1183}{12512}},\quad \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bc878cc117c447f6845d948e8966211f29c2cd11)
qui, étant réduites en décimales, deviennent
![{\displaystyle 0{,}089,\ \ \quad 0{,}09467,\ \ \quad 0{,}094549,\ \ \quad 0{,}0945515,\ \ \ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/04d4e75c61232bc61ef1ddcf43e3b96b4f06824f)
Or, nous avons vu dans la Note précédente que la méthode de Newton donne pour la valeur de
la série convergente
![{\displaystyle 0{,}1,\quad 0{,}946,\quad 0{,}09455147,\quad \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a67f131ddf8fe806f219a0e85fee8eed87dca934)
d’où l’on peut juger de l’accord des deux méthodes. En effet, nous ferons voir plus bas que ces méthodes, quoique fondées sur des principes différents, reviennent à peu près au même dans le fond et donnent des résultats semblables.