ANALISE TRANSCENDANTE.
Essai sur un nouveau mode d’exposition des principes
du calcul différentiel, du calcul aux différences et
de l’interpolation des suites, considérées comme dérivant
d’une source commune ;
Par
M. Ampère, de l’Académie royale des sciences
de Paris, de celles d’Édimbourg, de Cambridge, de Genève, etc.,
Professeur de physique au Collége de France.
≈≈≈≈≈≈≈≈≈
I. Le calcul différentiel, le calcul aux différences et les diverses méthodes d’interpolation reposent également sur un petit nombre de formules générales, applicables à toutes les fonctions, et qu’on n’a démontrées jusqu’ici que par des considérations souvent compliquées, presque toujours déduites de principes éloignée, ou d’inductions de nature à laisser des doutes sur leur généralité. On remarque surtout, dans l’exposition de ces diverses branches d’analise, une variété de procédé et de raisonnemens qui ne laisse que difficilement apercevoir leur liaison réciproque, et l’identité des principes dont elles, ne sont pourtant, en quelque sorte, que des traductions variées.
Frappés de ces considérations, nous avons pensé faire une chose qui pourrait intéresser les géomètres, en déduisant toutes ces diverses formules de quelques théorèmes nouveaux ou peu connus, que nous démontrerons d’abord, et dont il nous suffira ensuite de traduire les énoncés dans l’algorithme reçu, relatif aux divers cas particuliers, pour en voir éclore, d’une manière tout-à-fait naturelle, les formules connues qui répondent à chacun d’eux.
2. Soit
une fonction de ce de forme quelconque, de manière que
soient respectivement ce que devient cette fonction, lorsqu’on y fait, tour-à-tour,
Soient posés successivement
![{\displaystyle {\frac {\operatorname {f} b-\operatorname {f} a}{b-a}}=\operatorname {f} _{1}(a,b),\ {\frac {\operatorname {f} c-\operatorname {f} a}{c-a}}=\operatorname {f} _{1}(a,c),\ {\frac {\operatorname {f} d-\operatorname {f} a}{d-a}}=\operatorname {f} _{1}(a,d),\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/6589fcd14a24cb38fa10ef6caf3c0d53d0719a1d)
![{\displaystyle {\frac {\operatorname {f} _{1}(a,c)-\operatorname {f} _{1}(a,b)}{c-b}}=\operatorname {f} _{2}(a,b,c),\qquad {\frac {\operatorname {f} _{1}(a,d)-\operatorname {f} _{1}(a,b)}{d-b}}=\operatorname {f} _{2}(a,b,d),\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/cd6d6c237eab8a7353865178281f69a2554efdc8)
![{\displaystyle {\frac {\operatorname {f} _{2}(a,b,d)-\operatorname {f} _{2}(a,b,c)}{d-c}}=\operatorname {f} _{3}(a,b,c,d),\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/248ef50ff6e34fdc57632ed7e69ddb7e25159a5e)
Les fonctions
sont ce que nous appellerons à l’avenir les fonctions interpolaires des différens ordres des quantités ![{\displaystyle a,b,c,d,\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/ddc8541a9372c87dce345ae3ca7b3f7019fa7647)
3. La première remarque que nous ferons au sujet de ces sortes de fonctions, c’est qu’elles sont toujours symétriques, de telle sorte qu’on y peut intervertir, comme on voudra, l’ordre des élémens
qui concourent à leur formation, sans qu’elles eu éprouvent aucun changement. En effet, on a, en premier lieu,
![{\displaystyle \operatorname {f} _{1}(a,b)={\frac {\operatorname {f} b-\operatorname {f} a}{b-a}}={\frac {\operatorname {f} a}{a-b}}+{\frac {\operatorname {f} b}{b-a}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/cd5b8487ee555e0967770dcce24789024f21f7cc)
fonction évidemment symétrique en
et
et on pourrait en dire autant de toutes les fonctions interpolaires à deux lettres ou du premier ordre.
On aura donc, semblablement,
![{\displaystyle \operatorname {f} _{1}(a,c)={\frac {\operatorname {f} a}{a-c}}+{\frac {\operatorname {f} c}{c-a}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9eec2ef678107c5a5e23f7f6e7667cfc15ae3cfd)
mais,
![{\displaystyle \operatorname {f} _{2}(a,b,c)={\frac {\operatorname {f} _{1}(a,c)-\operatorname {f} _{1}(a,b)}{c-b}}={\frac {\operatorname {f} _{1}(a,b)}{b-c}}+{\frac {\operatorname {f} _{1}(a,c)}{c-b}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/783326eaa7f33b21472916a778eb78b4da0fc96e)
mettant donc pour
et
les valeurs ci-dessus, il viendra
![{\displaystyle \operatorname {f} _{2}(a,b,c)={\frac {\operatorname {f} a}{(a-b)(b-c)}}+{\frac {\operatorname {f} b}{(b-a)(b-c)}}+{\frac {\operatorname {f} a}{(a-c)(c-b)}}+{\frac {\operatorname {f} c}{(c-a)(c-b)}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2756e1cef20c46fd9f3d1ea195843774ca30ae3d)
ou, en réduisant à une seule les deux fractions de même numérateur,
![{\displaystyle \operatorname {f} _{2}(a,b,c)={\frac {\operatorname {f} a}{(a-b)(a-c)}}+{\frac {\operatorname {f} b}{(b-a)(b-c)}}+{\frac {\operatorname {f} c}{(c-a)(c-b)}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b71e5172e184077783ea47cc1e47174d2c6e1d10)
fonction également symétrique, et on démontrerait la même chose de toutes les fonctions interpolaires à trois lettres, ou du second ordre.
On pourrait, par des transformations analogues, étendre progressivement le même principe aux fonctions interpolaires des ordres supérieurs ; mais, afin de ne pas nous borner à une simple induction, à l’égard d’un principe qui est fondamental, dans la théorie qui nous occupe, prouvons que, si la symétrique se soutient jusqu’aux fonctions du
ième ordre, inclusivement, elle aura lieu également pour celles du
ième.
Supposons donc que, pour les
lettres
on ait trouvé
![{\displaystyle \operatorname {f} _{n-1}(a,b,c,\ldots p,q,r)=\left\{{\begin{aligned}&{\frac {\operatorname {f} a}{(a-b)(a-c)\ldots (a-q)(a-r)}}\\\\+&{\frac {\operatorname {f} b}{(b-a)(a-c)\ldots (b-q)(b-r)}}\\&\ldots \ldots \ldots \ldots \ldots \ldots \\+&{\frac {\operatorname {f} q}{(q-a)(q-b)\ldots (q-p)(q-r)}}\\\\+&{\frac {\operatorname {f} r}{(r-a)(r-b)\ldots (r-p)(r-q)}}\end{aligned}}\right\}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0875c50c4142bbe6f243f7e66387e1add187b540)
nous aurons semblablement, pour les
lettres ![{\displaystyle a,b,c,\ldots p,q,s,\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/9b848bd1c9319916c3f1ff6551aa84b46e690b21)
![{\displaystyle \operatorname {f} _{n-1}(a,b,c,\ldots p,q,s)=\left\{{\begin{aligned}&{\frac {\operatorname {f} a}{(a-b)(a-c)\ldots (a-q)(a-s)}}\\\\+&{\frac {\operatorname {f} b}{(b-a)(b-c)\ldots (b-q)(b-s)}}\\&\ldots \ldots \ldots \ldots \ldots \ldots \\+&{\frac {\operatorname {f} q}{(q-a)(q-b)\ldots (q-p)(q-s)}}\\\\+&{\frac {\operatorname {f} s}{(s-a)(s-b)\ldots (s-p)(s-q)}}\end{aligned}}\right\}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8d1a16482a73f87322eee012c6eb66ee34ce603d)
fonctions qui sont symétriques l’une et l’autre. Mais, en vertu de la définition des fonctions interpolaires, on a, pour les
lettres ![{\displaystyle a,b,c,\ldots p,q,s,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/67e2cf8ae9ebf0aa630b60ae4a4350d7d375e571)
![{\displaystyle {\begin{aligned}\operatorname {f} _{n}(a,b,c,\ldots p,r,s)&={\frac {\operatorname {f} _{n-1}(a,b,c,\ldots p,q,s)-\operatorname {f} _{n-1}(a,b,c,\ldots p,q,r)}{s-r}},\\\\&={\frac {\operatorname {f} _{n-1}(a,b,c,\ldots p,q,r)}{r-s}}+{\frac {\operatorname {f} _{n-1}(a,b,c,\ldots p,q,s)}{s-r}}\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a50f6faffabe238b1556ecd8c69675688baec089)
il viendra donc, en substituant,
![{\displaystyle \operatorname {f} _{n}(a,b,c,\ldots q,r,s)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ea2ac82844287929e2f2701a49ad4f43ee17c1e1)
![{\displaystyle =\left\{{\begin{aligned}&{\frac {\operatorname {f} a}{(a-b)(a-c)\ldots (a-q)(a-r)}}+{\frac {\operatorname {f} a}{(a-b)(a-c)\ldots (a-q)(a-s)}}\\\\+&{\frac {\operatorname {f} b}{(b-a)(b-c)\ldots (b-q)(b-r)}}+{\frac {\operatorname {f} b}{(b-a)(b-c)\ldots (b-q)(b-s)}}\\+&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots +\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\+&{\frac {\operatorname {f} q}{(q-a)(q-b)\ldots (q-p)(q-r)}}+{\frac {\operatorname {f} q}{(q-a)(q-b)\ldots (q-p)(q-s)}}\\\\+&{\frac {\operatorname {f} r}{(r-a)(r-b)\ldots (r-p)(r-q)}}+{\frac {\operatorname {f} s}{(s-a)(s-b)\ldots (s-p)(s-q)}}\end{aligned}}\right\}\,:}](https://wikimedia.org/api/rest_v1/media/math/render/svg/522f969a65f05542daab966498d41df755aec8a6)
où, en réduisant à une seule les fractions de même numérateur,
![{\displaystyle \operatorname {f} _{n}(a,b,c,\ldots p,q,s)=\left\{{\begin{aligned}&{\frac {\operatorname {f} a}{(a-b)(a-c)\ldots (a-r)(a-s)}}\\\\+&{\frac {\operatorname {f} b}{(b-a)(b-c)\ldots (b-r)(b-s)}}\\&+\ldots \ldots \ldots \ldots \ldots \ldots \\+&{\frac {\operatorname {f} r}{(r-a)(r-b)\ldots (r-q)(r-s)}}\\\\+&{\frac {\operatorname {f} r}{(s-a)(s-b)\ldots (s-q)(s-r)}}\end{aligned}}\right\}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a60d9f8f65920445bfe744262b8e2b778abb76c)
fonction également symétrique. Il demeure donc établi par là que, si la symétrie se soutient jusqu’à l’emploi de la
.ième lettre, inclusivement, elle aura lieu encore après l’introduction de la
ième, quel que soit
puis donc que cette symétrique a lieu en effet pour des fonctions interpolaires formées de deux et de trois lettres, il s’ensuit qu’elle doit être regardée comme un fait analitique généralement démontré.
4. D’après le mode de génération des fonctions interpolaires, il existe une relation fort simple entre deux d’entre elles de même ordre, ne différant l’une de l’autre que par un seul des élémens qui les composent. Puisqu’en effet on a
![{\displaystyle \operatorname {f} _{n+1}(a,b,c,d,\ldots )={\frac {\operatorname {f} _{n}(b,c,d,\ldots )-\operatorname {f} _{n}(a,c,d,\ldots )}{b-a}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/afecf3f0b997a4db51e9a0e0fb3c52fc1cbdd83c)
on en conclura
![{\displaystyle \operatorname {f} _{n}(a,c,d,\ldots )-\operatorname {f} _{n}(b,c,d,\ldots )=(a-b)\operatorname {f} _{n+1}(a,b,c,d,\ldots ).\qquad }](https://wikimedia.org/api/rest_v1/media/math/render/svg/33d2f21b0a63a02efd56ce15e50e6d3d8310aa85)
(1)
On peut également obtenir une relation très-remarquable entre trois de ces fonctions interpolaires de même ordre, ne différant que par l’exclusion donnée tour-à-tour à une lettre, sur trois d’entre elles. On a en effet, par la formule (1)
![{\displaystyle {\begin{aligned}&\operatorname {f} _{n}(a,c,d,\ldots )-\operatorname {f} _{n}(b,c,d,\ldots )=(a-b)\operatorname {f} _{n+1}(a,b,c,d,\ldots ),\\&\operatorname {f} _{n}(a,b,d,\ldots )-\operatorname {f} _{n}(a,c,d,\ldots )=(b-c)\operatorname {f} _{n+1}(a,b,c,d,\ldots ).\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/13a68ab8fcf6801e63d4e0f2bbd6c1e7887128a7)
d’où, en retranchant du produit de la première par
le produit de la seconde par
et réduisant,
![{\displaystyle (a-b)\operatorname {f} _{n}(a,b,d,\ldots )+(b-c)\operatorname {f} _{n}(b,c,d,\ldots )+(c-a)\operatorname {f} _{n}(a,c,d,\ldots )=0\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bdd8fb869da09e3a74cbfcf1d88077098e67db51)
(2)
formule que l’on retiendra facilement dans sa mémoire, en remarquant 1.o que la somme des trois coefficiens binômes est nulle ; 2.o que la lettre qui manque dans chacun de ces coefficiens est aussi celle qui manque dans la fonction qu’il multiplie. Remarquons, au surplus, que, bien que nous ayons supposé qu’il s’agissait des trois premières lettres, le théorème n’en demeure pas moins établi généralement pour trois lettres quelconques ; puisqu’à cause de la symétrie des fonctions interpolaires, on peut toujours amener ces trois lettres à être les trois premières.
5. Les fonctions interpolaires des différens ordres sont, comme l’on voit, complètement déterminées, tant que les élémens
dont elles se composent sent tous différens les uns des autres ; mais, si tous ou partie d’entre eux sont égaux, il arrive que toutes ou partie des fonctions
qu’il faut successivement calculer, pour parvenir à la valeur de
se présentent sous la forme
Or, on sait que, lorsqu’une fonction quelconque, qui n’est susceptible que d’une valeur unique, se présente sous cette forme, en vertu de l’égalité de deux des élémens
et
dont elle se compose, trois cas seulement peuvent se présenter. 1.o il peut arriver que la différence
décroissant indéfiniment, la fonction décroisse aussi indéfiniment, de manière à pouvoir devenir moindre que toute grandeur donnée ; ce qu’on exprime en disant qu’elle devient nulle quand
2.o il peut arriver, au contraire, que la différence
décroissant indéfiniment, la fonction croisse indéfiniment, de manière à pouvoir surpasser toute grandeur donnée ; ce qu’on exprime en disant qu’elle devient infinie, quand
3.o enfin il peut se faire que
décroissant indéfiniment, la fonction ne décroisse ni ne croisse indéfiniment, mais seulement de manière à tendre sans cesse vers une grandeur constante, vers une limite finie, dont elle pourra différer de moins de toute grandeur donnée, en prenant
d’une petitesse suffisante ; et alors cette limite fixe sera dite la valeur de la fonction qui répond
Or, nous allons voir qu’à l’égard des fonctions interpolaires, ce troisième cas est, généralement parlant, le seul qui puisse avoir lieu, ou, en d’autres termes, qu’en général ces sortes de fonctions ne sauraient jamais devenir nulles, par l’effet de l’égalité de tous ou partie des élémens dont elles se composent.
Pour le prouver, imaginons qu’on partage l’intervalle
en un nombre arbitraire
de parties égales ; posons, pour abréger,
et soient
![{\displaystyle a_{1}=a+{\frac {k}{m}},\ a_{2}=a_{1}+{\frac {k}{m}},\ a_{3}=a_{2}+{\frac {k}{m}},\ldots b=a_{m-1}+{\frac {k}{m}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7e41de2532ee6171cae2b79df217cee82b197bf2)
nous aurons (1)
![{\displaystyle {\begin{aligned}\operatorname {f} _{n-1}(a_{1},c,d,\ldots )-\operatorname {f} _{n-1}(a,c,d,\ldots )&={\frac {k}{m}}\operatorname {f} _{n}(a_{1},a,c,d,\ldots ),\\\\\operatorname {f} _{n-1}(a_{2},c,d,\ldots )-\operatorname {f} _{n-1}(a_{1},c,d,\ldots )&={\frac {k}{m}}\operatorname {f} _{n}(a_{1},a_{2},c,d,\ldots ),\\\\\operatorname {f} _{n-1}(a_{3},c,d,\ldots )-\operatorname {f} _{n-1}(a_{2},c,d,\ldots )&={\frac {k}{m}}\operatorname {f} _{n}(a_{2},a_{3},c,d,\ldots ),\\\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots &\ldots \ldots \ldots \ldots \ldots \ldots \ldots \\\operatorname {f} _{n-1}(b,c,d,\ldots )-\operatorname {f} _{n-1}(a_{m-1},c,d,\ldots )&={\frac {k}{m}}\operatorname {f} _{n}(a_{m-1},b,c,d,\ldots )\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7f1e66dc3346aa7d8d235814bd92489d27e7cd91)
en ajoutant, réduisant et se rappelant que (1)
![{\displaystyle \operatorname {f} _{n-1}(b,c,d,\ldots )-\operatorname {f} _{n-1}(a,c,d,\ldots )=k\operatorname {f} _{n}(a,b,c,d,\ldots ),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f257dd6e5995e40c0d21c1b88d031624e1ebf342)
il viendra, en divisant par ![{\displaystyle k,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/44e185ab9c990830d5055fa3ae698a4225ce67e0)
![{\displaystyle \operatorname {f} _{n}(a,b,c,d,\ldots )={\frac {1}{m}}\left\{{\begin{aligned}&\operatorname {f} _{n}(a,a_{1},c,d,\ldots )\\\\+&\operatorname {f} _{n}(a_{1},a_{2},c,d,\ldots )\\\\+&\operatorname {f} _{n}(a_{2},a_{3},c,d,\ldots )\\+&\ldots \ldots \ldots \ldots \ldots \\+&\operatorname {f} _{n}(a_{m-1},b,c,d,\ldots )\end{aligned}}\right\}\,;\quad (3)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/feedb9d1d434256335e462cdae6627a31ba33340)
c’est-à-dire que
est la moyenne arithmétique entre toutes les autres fonctions interpolaires du même ordre déduites de celle-là, en y mettant successivement
et
et
et
et
en place de
et ![{\displaystyle b.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aef051eb30c89e5493d672f6479566c673b0890a)
Il suit de là qu’excepté le cas où toutes les fonctions dont la somme divisée par
compose le second membre de l’équation (3) seraient égales entre elles, auquel cas chacune d’elles serait égale à celle qui forme le premier membre, il y en aura toujours de plus grandes et de plus petites que celle-là. Or, on peut toujours prendre le nombre
assez grand pour rendre et conséquemment les différences consécutives
d’une petitesse indéfinie ; donc, si une fonction interpolaire telle que
pouvait devenir nulle ou infinie, lorsqu’on y fait
on pourrait toujours prendre pour
un assez grand nombre pour rendre toutes les fonctions qui composent le second membre de l’équation (3) moindres ou plus grandes qu’une grandeur donnée quelconque, et, en particulier, moindres ou plus grandes que celle qui forme son premier membre, ce qui rendrait cette équation absurde.
On voit donc que la fonction
continue à avoir une valeur déterminée, qui n’est ni nulle ni infinie, lorsqu’on suppose les deux premiers élémens
et
égaux entre eux ; et il en sera de même encore, dans le cas de l’égalité entre deux autres quelconques de ses élémens ; puisqu’en vertu de la symétrie de la fonction, on pourra toujours amener ces deux élémens à être les premiers. Il en serait encore évidemment de même dans le cas de l’égalité de plus de deux élémens entre eux, ou de l’égalité entre plusieurs groupes d’élémens ; et, bien qu’alors le procédé que nous avons indiqué pour parvenir à ces sortes de fonctions n’en assigne plus la valeur, on pourra néanmoins les représenter et les employer dans les calculs, comme si cette valeur était connue. Observons d’ailleurs que rien n’empêche que les fonctions interpolaires, comme toutes les autres fonctions, ne deviennent nulles ou infinies, pour certaines valeurs particulières des élémens qui les composent ; mais cela n’arrivera jamais, tant que ces élémens conserveront leur indétermination.
6. De la formule (1) on peut conclure
![{\displaystyle {\begin{aligned}&\operatorname {f} x=\operatorname {f} a+(x-a)\operatorname {f} _{1}(a,x),\\&\operatorname {f} _{1}(a,x)=\operatorname {f} _{1}(a,b)+(x-b)\operatorname {f} _{2}(a,b,x),\\&\operatorname {f} _{2}(a,b,x)=\operatorname {f} _{2}(a,b,c)+(x-c)\operatorname {f} _{3}(a,b,c,x),\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\&\operatorname {f} _{n-1}(a,b,c,\ldots p,x)=\operatorname {f} _{n-1}(a,b,c,\ldots p,q)+(x-q)\operatorname {f} _{n}(a,b,c,\ldots p,q,x)\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3e71c33348dbb2f97298f2345e2a2712acfc1ae3)
ce qui donne, par des substitutions successives
![{\displaystyle \operatorname {f} x=\left\{{\begin{aligned}&\operatorname {f} a\\+&(x-a)\operatorname {f} _{1}(a,b)\\+&(x-a)(x-b)\operatorname {f} _{2}(a,b,c)\\+&(x-a)(x-b)(x-c)\operatorname {f} _{3}(a,b,c,d)\\+&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\+&(x-a)(x-b)(x-c)\ldots (x-p)\operatorname {f} _{n-1}(a,b,c,\ldots p,q)\\+&(x-a)(x-b)(x-c)\ldots (x-p)(x-q)\operatorname {f} _{n}(a,b,c,\ldots p,q,r).\end{aligned}}\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7134eb4d0edfd7eaf0a66a35bf19609d33c92194)
Cette valeur de
contient une fonction interpolaire de
, à son dernier terme ; mais, si les élémens
ont des valeurs comprises entre des limites peu étendues, et que
soit lui-même compris entre eux, la série sera assez convergente pour qu’on puisse se permettre d’en négliger le dernier terme, et on aura alors sensiblement
![{\displaystyle \operatorname {f} x=\left\{{\begin{aligned}&\operatorname {f} a\\+&a(x-a)\operatorname {f} _{1}(a,b)\\+&a(x-a)(x-b)\operatorname {f} _{2}(a,b,c)\\+&a(x-a)(x-b)(x-b)\operatorname {f} _{3}(a,b,c,d)\\+&a(x-a)(x-b)(x-c)\operatorname {f} _{3}(a,b,c,d)\\+&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\+&a(x-a)(x-b)(x-c)\ldots (x-p)\operatorname {f} _{n-1}(a,b,c,\ldots p,q)\end{aligned}}\right\}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7ab216611e63aa72d91f44ca39e07344ffc03ac0)
(4)
On reconnaît ici la formule ordinaire d’interpolation ; ce qui justifie la dénomination d’interpolaires que nous avons donnée aux fonctions dont elle se compose.
7. Posons
et désignons
par
![{\displaystyle \Delta y\ }](https://wikimedia.org/api/rest_v1/media/math/render/svg/261921101a00fb25e849f6d64d58221666753c06)
la quantité dont
![{\displaystyle \ \ y\quad }](https://wikimedia.org/api/rest_v1/media/math/render/svg/4374433723c2710527aaa84a6a6e10a0a5ea9161)
augmente, lorsque
![{\displaystyle x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/87f9e315fd7e2ba406057a97300593c4802b53e4)
devient
![{\displaystyle x+\Delta x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/79b92751e02d2719e1c5b73a57784c25d0a434d2)
;
par
![{\displaystyle \Delta ^{2}y}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95f6d3d67ede37da51f969d627d48b8a5d02cdee)
la quantité dont
![{\displaystyle \Delta y\ \,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4507680cf323a8a9b1301fa5e77d9cb988a60938)
augmente, lorsque
![{\displaystyle x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/87f9e315fd7e2ba406057a97300593c4802b53e4)
devient
![{\displaystyle x+\Delta x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/79b92751e02d2719e1c5b73a57784c25d0a434d2)
;
par
![{\displaystyle \Delta ^{3}y}](https://wikimedia.org/api/rest_v1/media/math/render/svg/474392a7ebada0cdaf5a1155313ca9c75fe7d752)
la quantité dont
![{\displaystyle \Delta ^{2}y}](https://wikimedia.org/api/rest_v1/media/math/render/svg/95f6d3d67ede37da51f969d627d48b8a5d02cdee)
augmente, lorsque
![{\displaystyle x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/87f9e315fd7e2ba406057a97300593c4802b53e4)
devient
![{\displaystyle x+\Delta x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/79b92751e02d2719e1c5b73a57784c25d0a434d2)
;
par
![{\displaystyle \Delta ^{n}y}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4f5f1e8a0979ad396cd84f5aa379a3f2c1f5a982)
la quantité dont
![{\displaystyle \Delta ^{n-1}y}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e29463e9a18c63936a3ac34d10f5104fd3391421)
augmente, lorsque
![{\displaystyle x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/87f9e315fd7e2ba406057a97300593c4802b53e4)
devient
![{\displaystyle x+\Delta x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/79b92751e02d2719e1c5b73a57784c25d0a434d2)
;
nous aurons, par la définition des fonctions interpolaire,
![{\displaystyle \quad y\ =\operatorname {f} x,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/04d769a9a8f082ec37d6810773786721574a0669)
![{\displaystyle \Delta y\ =\operatorname {f} (x+\Delta x)-\operatorname {f} x=\Delta x.\operatorname {f} _{1}(x,x+\Delta x),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3b0b3a6451bdb811fac2d0dfedd92bf9394bed99)
![{\displaystyle \Delta ^{2}y=\Delta x\left\{\operatorname {f} _{1}(x+\Delta x,x+2\Delta x)-\operatorname {f} _{1}(x,x+\Delta x)\right\}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c782e121547a5455d29d1536d1b377d2431f3240)
L’analogie conduit à soupçonner qu’en général
![{\displaystyle \Delta ^{n}y=n!\Delta x^{n}\{\operatorname {f} _{n}(x,x+\Delta x,x+2\Delta x,\ldots x+n\Delta x).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2d8e45fa8d96e6635ca027ff95053d3bbffe93e5)
Pour changer ce soupçon en certitude, il suffit de prouver qu’il devra en être ainsi ; si l’on a
![{\displaystyle \Delta ^{n-1}y=(n-1)!\Delta x^{n-1}\operatorname {f} _{n-1}(x,x+\Delta x,x+2\Delta x,\ldots x+(n-1)\Delta x)\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/adb708c2ea600ebd35c6a291f8760dd4d544bc4d)
or, en adoptant cette hypothèse, on aura
![{\displaystyle \Delta ^{n}y=(n-1)!\Delta x^{n-1}\left\{\operatorname {f} _{n-1}(x,x+\Delta x,x+2\Delta x,\ldots x+n\Delta x)\right.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5eb8c944c20f8a9d549ef8c1f201773c850b5cff)
![{\displaystyle \left.-\operatorname {f} _{n-1}(x,x+\Delta x,\ldots x+(n-1)\Delta x)\right\}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/05681ef3fd050defcd80537d68b6134a5776f327)
c’est-à-dire,
![{\displaystyle \Delta ^{n}y=n!\Delta x^{n}\{\operatorname {f} _{n}(x,x+\Delta x,x+2\Delta x,\ldots x+n\Delta x)\,;\qquad }](https://wikimedia.org/api/rest_v1/media/math/render/svg/ad1fa2c73b84437a53b0e1c7f7394f9cd295779a)
(5)
comme nous l’avions annoncé.
8. Reprenons la formule (4)
![{\displaystyle \operatorname {f} x=\operatorname {f} a+(x-a)\operatorname {f} _{1}(a,b)+(x-a)(x-b)\operatorname {f} _{2}(a,b,c)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ddbeea2b6d52b8c989e4584d5b59cd220a771a14)
![{\displaystyle +(x-a)(x-b)(x-c)\operatorname {f} _{3}(a,b,c,d)+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/bd6f6d076b2c2a2eb5b0d2dbeddb521d43a68570)
supposons les élémens
équidifférens, et soient posés
![{\displaystyle b-a=c-b=d-c=\ldots =k\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/809e68006a9a3d8ac83d046f7b2bc548ac81cbe6)
soient faits, de plus
![{\displaystyle y_{p}=\operatorname {f} a,\quad y_{p+1}=\operatorname {f} b,\quad y_{p+2}=\operatorname {f} c,\ldots ,y_{p+n}=\operatorname {f} x\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ee223605f63c28c28f93c3e2fb7519afa9ea74b1)
nous aurons
![{\displaystyle {\frac {\Delta y_{p}}{k}}={\frac {\operatorname {f} b-\operatorname {f} a}{b-a}}=\operatorname {f} _{1}(a,b),\ {\frac {\Delta ^{2}y_{p}}{k^{2}}}=2{\frac {\operatorname {f} _{1}(b,c)-\operatorname {f} _{1}(a,b)}{c-a}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0c3f88f9a3d4cadca0e900f45cc4a4cf7a888739)
![{\displaystyle =2\operatorname {f} _{2}(a,b,c),\ {\frac {\Delta ^{3}y_{p}}{k^{3}}}=\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/2719324bbbb77fe53da7cf71d5b54d921c81b749)
![{\displaystyle {\frac {\Delta y_{p+1}}{k}}={\frac {\operatorname {f} c-\operatorname {f} b}{c-b}}=\operatorname {f} _{1}(b,c),\ {\frac {\Delta ^{2}y_{p+1}}{k^{2}}}=2{\frac {\operatorname {f} _{1}(c,d)-\operatorname {f} _{1}(b,c)}{d-b}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5055fd9f315765666113382f42b4e4d0c20da4ac)
![{\displaystyle =2\operatorname {f} _{2}(b,c,d),\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/b3bebfd970452302f2f2f69d57b47cea46b3d808)
![{\displaystyle {\frac {\Delta y_{p+3}}{k}}={\frac {\operatorname {f} d-\operatorname {f} c}{d-c}}=\operatorname {f} _{1}(c,d),\ldots \ldots \ldots \ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/cae568a1afb738ba5f412a7f46a4accbf3fe61ad)
. . . . . . . . . . . . . . . . . . . .
d’où, en substituant dans la formule (4),
![{\displaystyle y_{p+n}=y_{p}+{\frac {x-a}{1}}.{\frac {\Delta y_{p}}{k}}+{\frac {x-a}{1}}.{\frac {x-b}{2}}.{\frac {\Delta ^{2}y_{p}}{k^{2}}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/11ea81e3be31604047dbf011ab8cfcaa36afa421)
![{\displaystyle +{\frac {x-a}{1}}.{\frac {x-b}{2}}.{\frac {x-c}{3}}.{\frac {\Delta ^{3}y_{p}}{k^{3}}}+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/2a138cb105530ca7e62435f71b88844f9d7a57b1)
(6)
Si nous posons
nous aurons
![{\displaystyle {\begin{aligned}&x-b=(x-a)-(b-a)=nk-k=(n-1)k,\\&x-c=(x-a)-(c-a)=nk-2k=(n-2)k,\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c89e2a1ca811ce13ed726fa2be876d8fe7729df3)
de sorte qu’alors cette formule deviendra
![{\displaystyle y_{p+n}=y_{p}+{\frac {n}{1}}.\Delta y_{p}+{\frac {n}{1}}.{\frac {n-1}{2}}.\Delta ^{2}y_{p}+{\frac {n}{1}}.{\frac {n-1}{2}}.{\frac {n-2}{3}}.\Delta ^{3}y_{p}+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/82d388b39a682bd2f955960329e257b6906403d7)
(7)
et, par suite, en supposant
nul,
![{\displaystyle y_{n}=y+{\frac {n}{1}}.\Delta y+{\frac {n}{1}}.{\frac {n-1}{2}}.\Delta ^{2}y+{\frac {n}{1}}.{\frac {n-1}{2}}.{\frac {n-2}{3}}.\Delta ^{3}y+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/693f910d5f578a883333607c3a25f3559867b0b7)
(8)
ce sont les formules connues d’interpolation, pour le cas de l’équidifférence des valeurs de la variable indépendante.
9. D’après ce qui a été démontré ci-dessus (5), on a
![{\displaystyle \operatorname {f} _{n}(a,b,c,\ldots r,s)=\left\{{\begin{aligned}&{\frac {\operatorname {f} a}{(a-b)(a-c)\ldots (a-r)(a-s)}}\\\\+&{\frac {\operatorname {f} b}{(b-a)(b-c)\ldots (b-r)(b-s)}}\\\\+&{\frac {\operatorname {f} c}{(c-a)(c-b)\ldots (c-r)(c-s)}}\\+&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\+&{\frac {\operatorname {f} s}{(s-a)(s-b)\ldots (s-q)(s-r)}}\end{aligned}}\right\}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b6112fb5ea18cc3eadbe9381dc9eaa8f2b333140)
supposant alors, comme ci-dessus,
![{\displaystyle b-a=c-b=d-c=\ldots =s-r=k,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/68cd65b7ca52acc6a2d4e07c67b2b8ee4b44176d)
et posant en outre
![{\displaystyle y_{p}=\operatorname {f} a,\ y_{p+1}=\operatorname {f} b,\ y_{p+2}=\operatorname {f} c,\ldots y_{p+n-1}=\operatorname {f} r,\ y_{p+n}=\operatorname {f} s,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f4af9eec915a74af351166642fa4352ce0e5ef93)
on trouvera, en substituant et renversant l’ordre des termes du second membre,
![{\displaystyle \operatorname {f} _{n}(a,a+k,a+2k,a+3k,\ldots a+nk)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8e1520f4c28722ee0d2fdf9fde6363ab369fa414)
![{\displaystyle ={\frac {y_{p+n}}{n!k^{n}}}-{\frac {y_{p+n-1}}{(n-1)!1!k^{n}}}+{\frac {y_{p+n-2}}{(n-2)!2!k^{n}}}-{\frac {y_{p+n-3}}{(n-3)!3!k^{n}}}+}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b5e57990612e3b11a18ae540e9cecfba19241802)
![{\displaystyle \ldots \pm {\frac {y_{p+1}}{1!(n-1)!k^{n}}}\mp {\frac {y_{p}}{n!k^{n}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f209baf1873ecb90254e0eba824e2a0b71e6fcf2)
mais si, dans la formule (5), on change
en
en
et qu’on y fasse, en outre,
elle donnera
![{\displaystyle \operatorname {f} _{n}(a,a+k,a+2k,a+3k,\ldots a+nk)={\frac {\Delta ^{n}y_{p}}{n!k^{n}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3033cac26c3ef394cb2ff7b14bc2dfaf77e355ef)
donc, en égalant ces deux valeurs et multipliant par ![{\displaystyle n!,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fd61e8437239aa98a3c40ddf512d8819b0e0ff6d)
![{\displaystyle \Delta ^{n}y_{p}=y_{p+n}-{\frac {n}{1}}y_{p+n-1}+{\frac {n}{1}}.{\frac {n-1}{2}}y_{p+n-2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c618a14e0a6bcbed8b07484b1997dc8d08229a63)
![{\displaystyle -{\frac {n}{1}}.{\frac {n-1}{2}}.{\frac {n-2}{3}}y_{p+n-3}+\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/08daf9e9e757ad25e80db48d54a5b3d156ad80ca)
(9)
d’où, en supposant
nul
![{\displaystyle \Delta ^{n}y=y_{n}-{\frac {n}{1}}y_{n-1}+{\frac {n}{1}}.{\frac {n-1}{2}}y_{n-2}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1f89b983f6095a98627b81be335375c550216af2)
![{\displaystyle -{\frac {n}{1}}.{\frac {n-1}{2}}.{\frac {n-2}{3}}y_{n-3}+\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/49ece5635bc86dafc1fd857e5d1748d7d6691d84)
(10)
formules inverses des formules (7) et (8).
10. Revenons au cas ou plusieurs des élémens de la fonction
sont égaux entre eux. Nous avons observé (5) que, bien qu’alors cette fonction se présentât sous la forme
elle n’en avait pas moins une valeur déterminée qui, en général, n’était ni nulle ni infinie, et que, bien qu’alors notre procédé ne fut plus propre à en assigner la valeur, les mêmes notations n’en étaient pas moins propres à la représenter ; en continuant donc à les employer, nous aurons d’abord (1)
![{\displaystyle \operatorname {f} _{1}(b,b)-\operatorname {f} _{1}(a,b)=(b-a)\operatorname {f} _{2}(a,b,b),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9c47b8371b5957ae8bf949753825ecca44c3f4ce)
![{\displaystyle \operatorname {f} _{1}(a,b)-\operatorname {f} _{1}(a,a)=(b-a)\operatorname {f} _{2}(a,a,b),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c8d420132f98810df87124579de59017dfb84c59)
d’où, en ajoutant et réduisant,
![{\displaystyle \operatorname {f} _{1}(b,b)-\operatorname {f} _{1}(a,a)=(b-a)\left\{\operatorname {f} _{2}(a,b,b)+\operatorname {f} _{2}(a,a,b)\right\}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/92b0fd35d3e7998c88d2c5cdfef35c911054172d)
Nous aurons aussi
![{\displaystyle {\begin{aligned}\operatorname {f} _{2}(b,b,b)-\operatorname {f} _{2}(a,b,b)&=(b-a)\operatorname {f} _{3}(a,b,b,b),\\\operatorname {f} _{2}(a,b,b)-\operatorname {f} _{2}(a,a,b)&=(b-a)\operatorname {f} _{3}(a,a,b,b),\\\operatorname {f} _{2}(a,a,b)-\operatorname {f} _{2}(a,a,a)&=(b-a)\operatorname {f} _{3}(a,a,a,b),\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fd4479ee0953a037b45901cf20e8c3d595bb0952)
d’où en ajoutant encore et réduisant
![{\displaystyle \operatorname {f} _{2}(b,b,b)-\operatorname {f} _{2}(a,a,a)=(b-a)\left\{\operatorname {f} _{3}(a,b,b,b)+\operatorname {f} _{3}(a,a,b,b)+\operatorname {f} _{3}(a,a,a,b)\right\}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e791015b5161fb2010c3c7820b2e32e4ca2249e3)
En continuant de la même manière, on aura en général
![{\displaystyle \operatorname {f} _{n-1}(b,b,b,\ldots b)-\operatorname {f} _{n-1}(a,a,a,\ldots a)=(b-a)\left\{{\begin{aligned}&\operatorname {f} _{n}(a,b,b,\ldots b,b,b)\\+&\operatorname {f} _{n}(a,a,b,\ldots b,b,b)\\+&\operatorname {f} _{n}(a,a,a,\ldots b,b,b)\\+&\ldots \ldots \ldots \ldots \ldots \\+&\operatorname {f} _{n}(a,a,a,\ldots a,b,b)\\+&\operatorname {f} _{n}(a,a,a,\ldots a,a,b)\end{aligned}}\right\}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9b037f959015a83c11d46168a82445a35688db85)
si, dans cette formule, on fait
d’où
et qu’on change ensuite
en
, elle deviendra
![{\displaystyle \operatorname {f} _{n-1}(x+k,x+k,x+k,\ldots x+k,x+k)-\operatorname {f} _{n-1}(x,x,x,\ldots x,x)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4cadee8ae9c101571eccc11a1441c73254175a57)
![{\displaystyle =k\left\{{\begin{aligned}&\operatorname {f} _{n}(x,x+k,x+k,\ldots x+k,x+k,x+k)\\+&\operatorname {f} _{n}(x,x,x+k,\ldots x+k,x+k,x+k)\\+&\operatorname {f} _{n}(x,x,x,\ldots x+k,x+k,x+k)\\+&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \\+&\operatorname {f} _{n}(x,x,x,\ldots x,x+k,x+k)\\+&\operatorname {f} _{n}(x,x,x,\ldots x,x+k)\end{aligned}}\right\}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc3430d3ab5b7edbdb3664dbe62d7d7e2c9baeca)
formule qui va nous servir tout-à-l’heure.
11. En conservant toujours les mêmes notations, on a, par la formule fondamentale,
![{\displaystyle {\begin{aligned}&\operatorname {f} (x+k)=\operatorname {f} x+k\operatorname {f} _{1}(x,x+k),\\&\operatorname {f} _{1}(x,x+k)=\operatorname {f} _{1}(x,x)+k\operatorname {f} _{2}(x,x,x+k),\\&\operatorname {f} _{2}(x,x,x+k)=\operatorname {f} _{2}(x,x,x)+k\operatorname {f} _{3}(x,x,x,x+k),\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\&\operatorname {f} _{n-1}(x,x,\ldots x,x+k)=\operatorname {f} _{n-1}(x,x,x,\ldots x,x)+k\operatorname {f} _{n}(x,x,x,\ldots x,x,x+k)\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/043c230ab9d1bded205a5e5ba3ab542ed882b8a4)
prenant la somme des produits respectifs de ces équations par
et réduisant, on trouvera
![{\displaystyle \operatorname {f} (x+k)=\left\{{\begin{aligned}&\operatorname {f} x\\+&k\operatorname {f} _{1}(x,x)\\+&k^{2}\operatorname {f} _{2}(x,x,x)\\+&k^{3}\operatorname {f} _{3}(x,x,x,x)\\+&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\+&k^{n-1}\operatorname {f} _{n-1}(x,x,x,\ldots x,x)\\+&k^{n}\operatorname {f} _{n}(x,x,x,\ldots x,x,x+k)\end{aligned}}\right\}\,;\quad }](https://wikimedia.org/api/rest_v1/media/math/render/svg/152b49528ae1ee11710bfafd032d9c3b5bfe30be)
(12)
développement de
en série où on pourra toujours prendre
assez petit, sans être nul, pour que la série soit convergente, puisque les fonctions qui multiplient les diverses puissances de
ont toujours des valeurs finies. En outre, si l’on parvient à assigner deux limites entre lesquelles se trouve compris le dernier terme du développement, le seul qui renferme à sous le signe de fonction, on connaîtra aussi par là les limites de l’erreur commise en bornant la série aux seuls termes qui précèdent celui-là.
Examinons présentement, d’une manière plus particulière, la nature des multiplicateurs des diverses puissances de
dans la formule (12). Soit, en général, désigné par
la limite vers laquelle tend sans cesse le rapport
![{\displaystyle {\frac {\operatorname {f} ^{(p-1)}(x+k)-\operatorname {f} ^{(p-1)}x}{k}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9453b3cc3e90d415a36d902bdf0c3a35d10648b6)
lorsque
tend vers zéro ; où, en d’autres termes, ce que devient ce rapport, lorsque
devient infiniment petit ou nul. Alors
seront ce qu’on appelle les dérivées successives de la fonction
On aura d’abord
![{\displaystyle \operatorname {f} _{1}(x,x)=\operatorname {f} 'x\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/26c6cd7e08eec24adf1c51dffd530bef56631b2f)
car on a
![{\displaystyle \operatorname {f} _{1}(x,x+k)={\frac {\operatorname {f} (x+k)-\operatorname {f} x}{k}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a724db43b4b1f920a1459d77470fe074c6d44c74)
et les deux quantités
et
sont respectivement ce que deviennent les deux membres de cette équation, lorsqu’on suppose ![{\displaystyle k=0.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/274c79b3ac027e0d0cd04cf86ae43c15567ba0bc)
On aura en conséquence,
![{\displaystyle {\frac {\operatorname {f} '(x+k)-\operatorname {f} 'x}{k}}={\frac {\operatorname {f} _{1}(x+k,x+k)-\operatorname {f} _{1}(x,x)}{k}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ae0a40c0f4057067764eba6e83ea8a541ff72bc3)
mais, en vertu de la formule (11)
![{\displaystyle \operatorname {f} _{1}(x+k,x+k)-\operatorname {f} _{1}(x,x)=k\left\{\operatorname {f} _{2}(x,x+k,x+k)+\operatorname {f} _{2}(x,x,x+k)\right\}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/de0d1f9da1bcc05ab150ae51df4d41324fbf3de3)
donc
![{\displaystyle {\frac {\operatorname {f} '(x+k)-\operatorname {f} 'x}{k}}=\operatorname {f} _{2}(x,x+k,x+k)+\operatorname {f} _{2}(x,x,x+k)\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4b3e1d40e44441d5733274849fd8cb0f778de49e)
donc aussi
car ce sont là les limites respectives vers lesquelles tendent les deux membres de cette équation, à mesure que
tend vers zéro.
On aura, en conséquence,
![{\displaystyle {\frac {\operatorname {f} ''(x+k)-\operatorname {f} ''x}{k}}=2.{\frac {\operatorname {f} _{2}(x+k,x+k,x+k)-\operatorname {f} _{2}(x,x,x)}{k}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/376b8a3acca3ed8655cf4c3fb91c8aac41237015)
mais, en vertu de la formule (11)
![{\displaystyle \operatorname {f} _{2}(x+k,x+k,x+k)-\operatorname {f} _{2}(x,x,x)=k\left\{{\begin{aligned}&\operatorname {f} _{3}(x,x+k,x+k,x+k)\\+&\operatorname {f} _{3}(x,x,x+k,x+k)\\+&\operatorname {f} _{3}(x,x,x,x+k)\end{aligned}}\right\}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/264d79175d11d4c88baf73eef8084fde43aea827)
donc
![{\displaystyle {\frac {\operatorname {f} ''(x+k)-\operatorname {f} ''x}{k}}=}](https://wikimedia.org/api/rest_v1/media/math/render/svg/08ea6318c44582396978b2b5baa85f8899095624)
![{\displaystyle 2.\left\{\operatorname {f} _{3}(x,x+k,x+k,x+k)+\operatorname {f} _{3}(x,x,x+k,x+k)+\operatorname {f} _{3}(x,x,x,x+k)\right\}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4221cc6d5aebc8894d00bbd829b7629075268685)
donc aussi
car ce sont là les limites respectives vers lesquelles tendent les deux membres de cette équation à mesure que
tend vers zéro.
Comme rien n’empêche de poursuivre ce raisonnement aussi loin qu’on voudra, il s’ensuit qu’on a
![{\displaystyle \operatorname {f} _{1}(x,x)={\frac {\operatorname {f} 'x}{1}},\quad \operatorname {f} _{2}(x,x,x)={\frac {\operatorname {f} ''x}{1.2}},\quad \operatorname {f} _{3}(x,x,x,x)={\frac {\operatorname {f} '''x}{1.2.3}},\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/6065abbdaf1b4b82f544c0f5f1e6a0106d9d90af)
et, en général,
![{\displaystyle \operatorname {f} _{n}(x,x,x,\ldots x,x)={\frac {\operatorname {f} ^{(n)}x}{n!}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/016e760b3716cd5ce6b6b7cd8449315a1f1811b1)
valeurs qui, substituées dans la formule (12), donneront
![{\displaystyle \operatorname {f} (x+k)=\operatorname {f} x+{\frac {k}{1}}\operatorname {f} 'x+{\frac {k^{2}}{1.2}}\operatorname {f} ''x+{\frac {k^{3}}{1.2.3}}\operatorname {f} '''x+\ldots +{\frac {k^{n-1}}{(n-1)!}}\operatorname {f} ^{(n-1)}x}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ab1709f4932ebd6f1d5918a4bd7c363f2cb35e9d)
![{\displaystyle +{\frac {k^{n}}{\operatorname {f} }}_{n}(x,x,x,\ldots x,x,x+k).\qquad }](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc6871519a37de36d0531c9c06ac335508556b4f)
(13)
On reconnaît ici la série de Taylor, qui se trouve ainsi démontrée.
12. On a cette suite d’équations (1)
![{\displaystyle {\begin{aligned}&\operatorname {f} _{1}(x,x)=\operatorname {f} _{1}(x,x+k)-k\operatorname {f} _{2}(x,x,x+k),\\&\operatorname {f} _{2}(x,x,x+k)=\operatorname {f} _{2}(x,x+k,x+2k)-2k\operatorname {f} _{3}(x,x,x+k,x+2k),\\&\operatorname {f} _{3}(x,x,x+k,x+2k)=\operatorname {f} _{3}(x,x+k,x+2k,x+3k)\\&\qquad \qquad \qquad \qquad \qquad \qquad -3k\operatorname {f} _{4}(x,x,x+k,x+2k,x+3k),\\\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf155006f235ca63dd74183a81188adc56266e79)
En prenant la somme des produits respectifs de ces équations par
et réduisant, il viendra
![{\displaystyle \operatorname {f} _{1}(x,x)=\operatorname {f} _{1}(x,x+k)-k\operatorname {f} _{2}(x,x+k,x+2k)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/349d9b9961b342aaa676090eca9c501aed089e52)
![{\displaystyle +1.2k^{2}\operatorname {f} _{3}(x,x+k,x+2k,x+3k)-\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/38b7bfb5d749a8d7c924d77eaf874423bcbbe794)
![{\displaystyle \pm (n-1)!k^{n-1}\operatorname {f} _{n}(x,x+k,x+2k,\ldots x+nk)}](https://wikimedia.org/api/rest_v1/media/math/render/svg/67691004aede089da5e903d5b02c8a017213f578)
![{\displaystyle \pm n!k^{n}\operatorname {f} _{n+1}(x,x+k,x+2k,\ldots x+nk)\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ac27b293f26ecb1315bded9b2ee13c74748322b8)
mais on a par la formule (5), en posant
d’où ![{\displaystyle \operatorname {f} _{1}(x,x)=\operatorname {f} 'x={\frac {\operatorname {d} y}{\operatorname {d} x}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/04ff32916fc79ce46749aa79313e879ddf6366f8)
![{\displaystyle {\begin{array}{rc}\operatorname {f} _{1}(x,x+k)&={\frac {\Delta y}{k}},\\\\\operatorname {f} _{2}(x,x+k,x+2k)&={\frac {\Delta ^{2}y}{1.2k^{2}}},\\\\\operatorname {f} _{3}(x,x+k,x+2k,x+3k)&={\frac {\Delta ^{3}y}{1.2.3k^{3}}},\end{array}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b82f1d9fdb640b6e8e40f2ac02f2d218b2ac539b)
il viendra donc, en substituant et multipliant par ![{\displaystyle k}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c3c9a2c7b599b37105512c5d570edc034056dd40)
![{\displaystyle k{\frac {\operatorname {d} y}{\operatorname {d} x}}={\frac {\Delta y}{1}}-{\frac {\Delta ^{2}y}{2}}+{\frac {\Delta ^{3}y}{3}}-\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/dd1b74821687def31084347e8f49cd91b17174f4)
![{\displaystyle \pm {\frac {\Delta ^{n}y}{n}}\mp n!k^{n+1}\operatorname {f} _{n+1}(x,x,x+k,x+2k,\ldots x+nk).}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ed77d005e0f8764a8d326b3aaee244257578717d)
On pourrait étendre indéfiniment ces recherches, mais ce qui précède suffit pour atteindre le but que nous nous étions proposé[1].