MÉMOIRE
SUR
LA MÉTHODE D’INTERPOLATION.
(Nouveaux Mémoires de l’Académie royale des Sciences et Belles-Lettres
de Berlin, années 1792 et 1793.)
La méthode d’interpolation est, après les logarithmes, la découverte la plus utile qu’on ait faite dans le calcul ; elle est surtout ; comme les logarithmes, d’un usage immense dans l’Astronomie, où elle sert non-seulement pour remplir dans les Tables les lieux intermédiaires entre ceux qu’on a calculés directement, mais encore pour suppléer dans une suite d’observations à celles qui manquent. Lorsque les nombres donnés, entre lesquels il s’agit d’insérer des nombres intermédiaires, sont en progression arithmétique, il est naturel de supposer que les termes intermédiaires forment aussi une même progression arithmétique avec les nombres donnés. Il n’y a donc alors qu’à insérer des moyens arithmétiques entre les nombres donnés c’est en quoi consiste la méthode des parties proportionnelles dont l’usage paraît connu de tout temps.
Mais cette méthode si simple ne peut avoir lieu qu’autant que les nombres donnés croissent ou décroissent également, c’est-à-dire par des différences constantes. Si ces différences ne sont pas constantes, on ne peut pas supposer non plus que celles des termes intermédiaires le soient, et la question se réduit alors à trouver la loi de l’augmentation ou diminution des nombres donnés pour pouvoir y assujettir aussi les nombres intermédiaires c’est l’objet de la méthode d’interpolation.
Ce qui se présente de plus simple dans cette recherche, c’est d’examiner si les différences des nombres donnés forment elles-mêmes une progression arithmétique ; dans ce cas il est visible qu’on peut âppliquer la méthode des parties proportionnelles à la suite des différences, ensuite il n’y aura qu’à remonter de cette suite à celle des nombres cherchés. De même, si les différences des nombres donnés ne formant pas une progression arithmétique, les différences de ces différences, qu’on appelle différences secondes, en forment elles-mêmes une, on pourra trouver les termes intermédiaires de cette dernière suite, et remonter de là successivement à celle des différences premières et enfin à celle des nombres à interpoler. C’est sur ce principe qu’est fondée la Théorie ordinaire de l’interpolation, laquelle se réduit ainsi à la solution de ce Problème :
Étant donnée une suite de termes dont les différences d’un ordre quelconque soient constantes, trouver un nombre quelconque de termes intermédiaires qui suivent la même loi.
Mouton, Astronome de Lyon, est le premier qui ait envisagé l’interpolation sous ce point de vue, dans son Ouvrage intitulé Observationes diametrorum Solis et Lunæ, etc., et imlorimé à Lyon en 1670. Ayant entrepris de calculer une Table des déclinaisons du Soleil pour chaque degré et minute de longitude, il a vu qu’on pouvait se contenter de la calculer directement de degré en degré, et de l’étendre ensuite de minute en minute par la méthode des différences. Il explique par différents exemples l’usage de cette méthode, et il donne d’après François Regnaud, à qui il avait proposé cette question, un procédé général, mais très-long, pour trouver la loi des termes à interpoler dans une Table dont les termes ont des différences constantes d’un ordre donné.
Il paraît néanmoins que la Théorie de l’interpolation est plus ancienne ; et il est bien naturel de penser en effet qu’elle a dû se présenter aux premiers calculateurs des Tables trigonométriques et logarithmiques, vu les secours prodigieux qu’elle offre dans ces sortes de calculs. Aussi je trouve que Henri Briggs, qui a calculé le premier les logarithmes des nombres naturels depuis l’unité jusqu’à
et depuis
jusqu’à
propose, pour remplir la lacune, une méthode d’interpolation fondée sur la considération des différences successives, qu’il dit avoir employée avec succès dans la construction du Canon trigonométrique pour les sinus et tangentes des degrés et centièmes de degré. (Voyez le Chapitre XIII de son Arithmetica logarithmica, et le XIIe de sa Trigonometria Britannica.) Cette méthode, dont Briggs donne les résultats sans démonstration, a été ensuite généralisée par Cotes dans sa Canonotechnia sive constructio tabularum per differentias ; mais ce dernier a également supprimé la démonstration de ses formules, et je ne connais personne qui jusqu’ici ait entrepris d’y suppléer ; ce qui vient peut-être de ce que d’un côté ces formules sont un peu compliquées, et de l’autre, de ce que l’usage en a été abandonné depuis que Newton en a proposé de plus simples, fondées sur la considération des courbes paraboliques.
Tout le monde connait la formule de Newton pour trouver une ordonnée quelconque d’une courbe parabolique, par les différences successives des ordonnées équidistantes ; c’est celle dont on se sert journellement en Astronomie pour interpoler les lieux des Planètes. Cette formule donne en effet tout de suite le terme que l’on cherche, et est par conséquent très-utile pour calculer la valeur de quelques termes ; mais, comme elle demande un calcul particulier pour chaque terme, elle est peu commode pour construire des Tables, où l’on aurait un grand nombre de termes consécutifs à interpoler ; au lieu que la méthode de Mouton, qui consiste à déterminer par les différences de la série donnée celles de la série interpolée, et à remonter de celles-ci aux termes de cette série, a l’avantage de réduire tout le calcul à des additions successives ; ce qui doit la rendre aussi propre pour la construction des Tables que celle de Newton l’est peu.
Mais la difficulté est de trouver ces différences, quel que soit le nombre des termes à interpoler et l’ordre des différences constantes de la série primitive. Mouton a fait le premier cette remarque importante que, dans une série qui conduit à des différences constantes d’un ordre quelconque, si l’on prend les termes de deux en deux, ou de trois en trois, ou, on a une série du même genre, mais dont la différence constante est égale à la différence constante de la première série multipliée par le nombre
ou
ou
élevé à une puissance égale à l’exposant de l’ordre de la différence constante. D’où il a conclu réciproquement, qu’en insérant dans une série à différences constantes des termes intermédiaires qui divisent chaque intervalle en
parties, la nouvelle série aura des différences constantes du même ordre que celles de la série donnée, et qui seront égale à celles-ci divisées par les nombres
élevés à la puissance dont l’exposant indiquera l’ordre des différences constantes. On peut donc trouver ainsi la différence constante de la série interpolée, et il ne s’agit plus que d’avoir les différences des ordres inférieurs, et même il suffit d’avoir les premiers termes des suites de ces différences. Mouton ne donne pour cela que quelques règles particulières sans démonstration, et il renvoie pour la solution générale à la méthode de Regnaud, qui consiste à construire une série dont les différences d’un ordre donné soient constantes, à prendre cette série pour la série interpolée, et à comparer ses termes pris de deux en deux, de trois en trois, ou,
avec ceux de la série donnée à interpoler, suivant qu’on veut partager l’intervalle d’un terme à l’autre en deux, en trois ou en un nombre plus grand de parties ; mais, quoiqu’on puisse toujours de cette manière résoudre la question dans chaque cas, il reste néanmoins à trouver des formules générales. Lalande est, je crois, le seul qui en ait cherché pour les cas de la seconde et de la troisième différence constante, dans un Mémoire sur les interpolations imprimé parmi ceux de l’Académie des Sciences de Paris pour 1761. Sa méthode est à proprement parler celle de Regnaud réduite en Analyse, et, si l’on voulait l’appliquer aux séries dont les différences constantes seraient d’un ordre supérieur au troisième, on tomberait dans des calculs longs et peut-être impraticables ; à plus forte raison serait-il comme impossible de parvenir par cette voie à la solution générale pour des différences constantes d’un ordre quelconque. J’ai donné en 1772 une formule qui renferme cette solution dans un Mémoire qui a pour titre Sur une nouvelle espèce de calcul, et qui est imprimé dans le volume de cette année[1] ; mais, comme je ne pensais pas alors à l’usage dont elle pouvait être pour les interpolations, je ne cherchai pas à développer cette formule.
Je crois donc devoir revenir sur cet objet, et donner ici une solution complète du Problème dont il s’agit, pour servir de Supplément au Mémoire cité. La méthode que j’ai employée dans ce Mémoire est fondée sur une analogie singulière qui a lieu entre les exposants des puissances et les indices du rang des termes dans une série, ou de l’ordre des différences de ces termes, en vertu de laquelle on peut traiter ces indices comme si c’étaient des exposants, et y appliquer les mêmes règles. Je ne m’arrêterai pas ici à démontrer cette analogie, dont on peut aisément vérifier les résultats par les méthodes connues ; je me contenterai de m’en servir, ainsi que je l’ai déjà fait, comme d’un instrument propre à découvrir des formules qu’on ne pourrait trouver par les moyens ordinaires qu’en connaissant leur forme d’avance ; je l’emploierai même ici d’une manière plus simple et plus uniforme, ce qui servira à donner encore à ce nouvel instrument d’Analyse un plus grand degré de perfection.
1. Pour distinguer les indices des exposants, je les placerai au bas des quantités, comme ceux-ci le sont au haut ; on suivra d’ailleurs le même algorithme pour les uns et les autres.
Soient donc
![{\displaystyle \mathrm {T_{0},\ \ T_{1},\ \ T_{2},\ \ T_{3}} ,\ldots ,\ \ \mathrm {T} _{n},\ \ \mathrm {T} _{n+1},\ \ \mathrm {T} _{n+2},\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5353f9f8ac1f3eae04bc9bb289655ee772b379c1)
les termes consécutifs d’une série quelconque, et
![{\displaystyle \mathrm {D_{1},\ \ D_{2},\ \ D_{3}} ,\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/94bda5305fefbca6b37c50df2bf5c1f8f16b83dc)
les différences successives de ces termes, c’est-à-dire
la différence première
la différence seconde
la différence troisième
et ainsi de suite. On aura sur-le-champ l’expression de la différence
ième par l’opération suivante.
Puisque
![{\displaystyle \mathrm {D_{1}=T_{1}-T_{0}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ab13d87280c6fe57238ee25245a2427d8b725ac)
en élevant les deux membres à l’indice
on aura
![{\displaystyle \mathrm {D_{m}=(T_{1}-T_{0})_{m}} \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fb21d006fdd256312470c8128114272dd1bd8ae3)
développant le second membre comme un binôme, en mettant les exposants au bas pour servir d’indices, on aura
(A)
|
|
|
formule connue et qu’on peut trouver par induction.
Si l’on fait
on a
![{\displaystyle \mathrm {D_{0}=T_{0}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3a82c9bdd2c097981233b049737169a528f1568e)
d’où il s’ensuit que le terme
qui doit précéder le terme
dans la série des différences est égal à ![{\displaystyle \mathrm {T} _{0}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0c43c6f66096518ea981f5aca1785d07265864d6)
Si l’on fait
négatif, les différences se changent en sommes ; de sorte que désignant les sommes de différents ordres par
![{\displaystyle \mathrm {S_{1},\ \ S_{2},\ \ S_{3}} ,\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/238a84aae7fc19d66473a47ec4bd8211a7591444)
on aura
![{\displaystyle \mathrm {S_{1}=D_{-1},\quad S_{2}=D_{-2}} ,\ldots ,\quad \mathrm {S} _{m}=\mathrm {D} _{-m}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b03b62c9145ec4634eb3cf73a3b26f2e0a6228d6)
par conséquent, en changeant
en
on aura
![{\displaystyle \mathrm {S} _{m}=\mathrm {T} _{-m}+m\mathrm {T} _{-m-1}+{\frac {m(m+1)}{2}}\mathrm {T} _{-m-2}+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ecd4eab547d8a4e62daac2e9f4695f7edabcf616)
formule connue aussi.
2. Si au contraire on veut avoir l’expression d’un terme quelconque
par le moyen des différences, on reprendra l’équation
![{\displaystyle \mathrm {D_{1}=T_{1}-T_{0}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9ab13d87280c6fe57238ee25245a2427d8b725ac)
laquelle donne
![{\displaystyle \mathrm {T_{1}=T_{0}+D_{1}} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/f9d528b66dda64604b0869ea2abafa388174546c)
ou bien, en mettant pour
![{\displaystyle \mathrm {T} _{0}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e76e38423f8dd877e7088f59f89ca2378acfdc7d)
le terme
![{\displaystyle \mathrm {D} _{0}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dfccac599043e427829c5ccb1042e5e491c2d330)
qui lui est égal,
![{\displaystyle \mathrm {T_{1}=D_{0}+D_{1}} \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dbbda9858220df9b17d8d73495a182ed1ab7a76e)
élevant à l’exposant ou indice
on aura de même
![{\displaystyle \mathrm {T} _{n}=(\mathrm {D_{0}+D_{1}} )_{n}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d2b0c1c7a079201cdd75e31b3ef3aef65249246f)
et développant par la formule du binôme, en ayant soin de mettre les exposants au bas pour servir d’indices,
(B)
|
|
|
où l’on se souviendra que
est la même chose que
C’est la formule différentielle de Newton, dont on se sert communément pour les interpolations, en donnant à l’indice
des valeurs fractionnaires.
Lorsque
est un nombre entier, on peut construire cette formule par des additions successives ; la manière la plus simple est de ranger les différences
![{\displaystyle \mathrm {D_{0},\ \ D_{1},\ \ D_{2},\ \ D_{3}} ,\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/ba8c265a21a1723f5421164508c4726380c93847)
dans une ligne horizontale, et de former au-dessous successivement d’autres lignes correspondantes, en faisant chaque terme de ces lignes égal à la somme de celui qui est au-dessus et de celui qui est à la droite de celui-ci dans la ligne immédiatement supérieure ; la première ligne verticale contiendra par ordre les termes
ou
On peut ainsi, quand la série des différences se termine, ce qui a lieu pour toutes les suites qui conduisent à des différences constantes, construire une Table aussi étendue que l’on voudra, pour avoir successivement tous les termes de la suite proposée.
Mais ce moyen mécanique ne peut plus être employé lorsque le nombre
est fractionnaire. Alors il faut calculer chaque terme séparément par la formule (B), en donnant successivement à
les valeurs convenables ainsi, s’il était question d’interpoler partout un terme entre deux termes consécutifs de la série
![{\displaystyle \mathrm {T_{0},\ \ T_{1},\ \ T_{2},} \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b98de019185791df9497d044eaac40ccdd6071dd)
il faudrait trouver les valeurs des termes
![{\displaystyle \mathrm {T} _{\frac {1}{2}},\ \ \mathrm {T} _{\frac {3}{2}},\ \ \mathrm {T} _{\frac {5}{2}},\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5e7780656ed5a9bd5961d65f14ef7f5e979d6e1c)
en faisant successivement
![{\displaystyle n={\frac {1}{2}},\ \ {\frac {3}{2}},\ \ {\frac {5}{2}},\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/604a57a78908808698a7af2b5b79c279b8f4f608)
et ainsi des autres cas semblables. La méthode de Mouton a l’avantage de réduire tous ces cas au premier, et de rendre par conséquent le procédé des additions successives applicable aux interpolations, comme nous le verrons plus bas.
3. Si l’on voulait que la formule
fût ordonnée suivant les puissances mêmes de
ce qui en rendrait peut-être l’usage plus commode dans le cas de
fractionnaire, rien ne serait plus aisé que de lui donner cette forme par le moyen de notre Analyse.
On n’a qu’à reprendre l’équation
![{\displaystyle \mathrm {T_{1}=D_{0}+D_{1}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c58b49038f39cdecffd918d0536641a597e99da8)
et l’écrire d’abord ainsi
![{\displaystyle \mathrm {T_{1}=D_{0}(1+D_{1})} \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/150108bfd850777ca7b5e5de3af561d654d94d3b)
ce qui est la même chose, puisque
![{\displaystyle \mathrm {D_{0}\times D_{1}=D_{0+1}=D_{1}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc553033a169912d40bb2937c937fce0c27b642b)
en observant pour les indices les lois des exposants ; on transformera ensuite
en
étant le nombre dont le logarithme hyperbolique désigné par
est l’unité ; ce qui donnera
![{\displaystyle \mathrm {T} _{1}=\mathrm {D} _{0}\times e^{\log(1+\mathrm {D} _{1})},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/40a9c378aedc6a2893a09ce0bcb7aef13a80fde8)
et, élevant les deux membres à l’exposant ou indice ![{\displaystyle n,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/397bfafc701afdf14c2743278a097f6f2957eabb)
![{\displaystyle \mathrm {T} _{n}=\mathrm {D} _{0}\times e^{n\log(1+\mathrm {D} _{1})}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5c8191d74c3739415c0dd4a21773ddd12ee798c4)
or on a par les formules connues donc on aura
![{\displaystyle e^{n\log(1+\mathrm {D} _{1})}=1+n\log(1+\mathrm {D} _{1})+{\frac {n^{2}}{2}}\left[\log(1+\mathrm {D} _{1})\right]^{2}+{\frac {n^{3}}{2.3}}\left[\log(1+\mathrm {D} _{1})\right]^{3}+\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0d17d3dc2be59f209c485fd3e77fe58a08069e29)
donc on aura
(C)
|
|
|
j’ai négligé ici de multiplier par
les termes qui contiennent des puissances de
parce que cette multiplication n’y apporterait aucun changement, l’indice
ne pouvant rien changer à ceux auxquels il se trouverait ajouté.
Maintenant il n’y aura qu’à développer
par les formules ordinaires, et l’on aura
![{\displaystyle {\begin{aligned}\mathrm {T} _{n}=&\mathrm {D} _{0}+n\mathrm {\left(D_{1}-{\frac {1}{2}}D_{2}+{\frac {1}{3}}D_{3}-\ldots \right)} \\&+\,\ {\frac {n^{2}}{2}}\mathrm {\left(D_{1}-{\frac {1}{2}}D_{2}+{\frac {1}{3}}D_{3}-\ldots \right)} ^{2}\\&+\,{\frac {n^{3}}{2.3}}\mathrm {\left(D_{1}\,-{\frac {1}{2}}D_{2}+{\frac {1}{3}}D_{3}-\ldots \right)} ^{3}\\&+\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b46175dd83ae7e8393cb07b58422c8db1529b8d9)
où il reste encore à développer les puissances de
![{\displaystyle \mathrm {D_{1}-{\frac {1}{2}}D_{2}+{\frac {1}{3}}D_{3}} -\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/07d92fa32d37ea71fdbd191d6048ce7acc49ae0a)
à la manière ordinaire, en traitant les indices comme des exposants. Cette formule ne sera donc que le développement de la formule (B) en effectuant les multiplications et ordonnant suivant les puissances de
mais, de la manière dont elle est présentée, on y voit la loi des termes qu’il serait très-difficile de découvrir par le développement dont il s’agit.
4. Au reste la formule (C) donne le moyen de trouver d’une manière générale l’expression de la différence d’un ordre quelconque de toute série algébrique dont le terme général
serait de la forme
![{\displaystyle \mathrm {T} _{n}=\mathrm {P} _{0}+n\mathrm {P} _{1}+{\frac {n^{2}}{2}}\mathrm {P} _{2}+{\frac {n^{3}}{2.3}}\mathrm {P} _{3}+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4b42f5f36e934fc13e547c178b276ded08436c5a)
étant des coefficients donnés, et
étant supposé successivement
Car, en comparant cette expression de
à celle de la formule (C), on aura
![{\displaystyle \mathrm {P_{0}=D_{0}=T_{0},\quad P_{1}=\log(1+D_{1}),\quad P_{2}=\left[\log(1+D_{1})\right]^{2}} ,\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/087bb7079bf2a1d9930607a98a0ce8b8c51f8ca0)
et, en général,
![{\displaystyle \mathrm {P} _{s}=\left[\log(1+\mathrm {D} _{1})\right]^{s}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1f290f4f93117038e05cdb9fa478a296b0cf83e9)
Cette équation donne
![{\displaystyle \mathrm {P_{1}=\log(1+D_{1}),\quad 1+D_{1}} =e^{\mathrm {P} _{1}}=\mathrm {1+P_{1}+{\frac {1}{2}}P_{2}+{\frac {1}{2.3}}P_{3}} +\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0f89d72c075fa1d11a553703b14a4dd0d565ab26)
donc
![{\displaystyle \mathrm {D_{1}=P_{1}+{\frac {1}{2}}P_{2}+{\frac {1}{2.3}}P_{3}} +\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/02572a0946c7aaafe8f0b43aa15319fb732b841d)
et de là, pour un indice quelconque ![{\displaystyle s,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e5748cdb81bf00075de8e7e6828c343687513830)
![{\displaystyle \mathrm {D} _{s}=\left(\mathrm {P_{1}+{\frac {1}{2}}P_{2}+{\frac {1}{2.3}}P_{3}+{\frac {1}{2.3.4}}P_{4}} +\ldots \right)^{s},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f3b8a19708ed5de320de5c52c0c685b8d66c5a57)
où il n’y aura qu’à développer la puissance
en convertissant les exposants de
en indices, et observant relativement à ces indices les mêmes lois que pour les exposants. Cette formule peut être utile dans quelques occasions ; elle l’est surtout pour transformer une série de la forme en une série équivalente de la forme
![{\displaystyle \mathrm {P} _{0}+n\mathrm {P} _{1}+{\frac {n^{2}}{2}}\mathrm {P} _{2}+{\frac {n^{3}}{2.3}}\mathrm {P} _{3}+\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/65f1a1fc0ac5b8ae45c2103097ed457aaeda2493)
en une série équivalente de la forme
![{\displaystyle \mathrm {D} _{0}+n\mathrm {D} _{1}+{\frac {n(n-1)}{2}}\mathrm {D} _{2}+{\frac {n(n-1)(n-2)}{2.3}}\mathrm {D} _{3}+\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/77813eafb62fb8f228d03cc23c5868a28a7e6697)
5. Considérons maintenant deux séries correspondantes, l’une représentée par
![{\displaystyle \mathrm {T_{0},\ \ T_{1},\ \ T_{2},\ \ T_{3}} ,\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/b628bd80ac4a3b7040c57a7e8021dca412bb6e4c)
dont les différences successives soient
![{\displaystyle \mathrm {D_{0},\ \ D_{1},\ \ D_{2},\ \ D_{3}} ,\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/ba8c265a21a1723f5421164508c4726380c93847)
comme plus haut, l’autre représentée par
![{\displaystyle t_{0},\ \ t_{1},\ \ t_{2},\ \ t_{3},\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/a67098c6aeea0f0e61ec40cf2e1446ad81813bca)
et dont les différences successives soient
![{\displaystyle d_{0},\ \ d_{1},\ \ d_{2},\ \ d_{3},\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5dd5dd2ff27facd6eabb5215a848f07b7fa95698)
supposons que les termes de la première série soient identiques avec les termes de la seconde pris à des intervalles égaux, de manière que l’on ait, en général,
![{\displaystyle \mathrm {T} _{s}=t_{ms},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6545d14d19bb423e40839699509971b9bd04d7c1)
c’est-à-dire que si par exemple
on ait
![{\displaystyle \mathrm {T} _{0}=t_{0},\ \ \mathrm {T} _{1}=t_{2},\ \ \mathrm {T} _{2}=t_{4},\ \ \mathrm {T} _{3}=t_{6},\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2b5bc07f3693dfbf0d4180ad21fed90e58a4e9d4)
si
on ait
![{\displaystyle \mathrm {T} _{0}=t_{0},\ \ \mathrm {T} _{1}=t_{3},\ \ \mathrm {T} _{2}=t_{6},\ \ \mathrm {T} _{3}=t_{9},\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/222634614307b7afe22f794cf53c3b1987933f18)
et ainsi de suite ; on propose de trouver la relation entre les différences
et les différences ![{\displaystyle d_{0},d_{1},d_{2},\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a87b4fdf8cf0c1577b34c59f4bf2a51e8497183)
L’équation
![{\displaystyle \mathrm {T} _{s}=t_{sm},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fe386fc84aec7b0273cd0ea2a8467f73abf295c0)
se réduit, en extrayant la racine
à
![{\displaystyle \mathrm {T} _{1}=t_{m},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/88658f49ea9eac0951107ebe73662a36274c0a96)
et, mettant pour
sa valeur
(2) de même que pour
sa valeur
elle donnera
![{\displaystyle \mathrm {D_{0}+D_{1}} =(d_{0}+d_{1})_{m}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5beb99ded895e2e2d3beaf08b4a980def453c542)
d’où l’on tirera la valeur de
en
ou de
en
suivant qu’on vou-
dra avoir les différences de la première série exprimées par celles de la seconde, ou réciproquement.
6. Nous aurons d’abord
![{\displaystyle \mathrm {D} _{1}=(d_{0}+d_{1})_{m}-\mathrm {D} _{0}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b798695db0578ad58e745f30b84da14e47f26f83)
mais
![{\displaystyle \mathrm {D} _{0}=\mathrm {T} _{0}=t_{0}=d_{0}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1a057bf02a485a9d9145a831b94a4c18645590f5)
donc
![{\displaystyle \mathrm {D} _{1}=(d_{0}+d_{1})_{m}-d_{0},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a4dbc81fc258121a187221e94871ce3077d5e4eb)
et, élevant à l’indice ![{\displaystyle s,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e5748cdb81bf00075de8e7e6828c343687513830)
![{\displaystyle \mathrm {D} _{s}=\left[(d_{0}+d_{1})_{m}-d_{0}\right]_{s}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ee7960f290188c5375b49a97551f6147c70120d2)
Or
se développe dans la série
![{\displaystyle d_{0}+md_{1}+{\frac {m(m-1)}{2}}d_{2}+{\frac {m(m-1)(m-2)}{2.3}}d_{3}+\ldots \,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e5d0a4079de18f818337d35b512549967d055cca)
donc, faisant cette substitution, on aura
(D)
|
|
|
on développera le second membre de cette équation comme si c’était la puissance
ième d’un polynôme en
en ayant soin de placer toujours au bas de la lettre
les exposants qui devraient affecter la quantité
suivant les règles des exposants.
Ainsi, si l’on suppose, pour abréger,
![{\displaystyle {\begin{aligned}&\left[mx+{\frac {m(m-1)}{2}}x^{2}+{\frac {m(m-1)(m-2)}{2.3}}x^{3}+\ldots \right]^{s}\\&\quad =\mathrm {A} x^{s}+\mathrm {B} x^{s+1}+\mathrm {C} x^{s+2}+\mathrm {E} x^{s+3}+\ldots ,&\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/986c54ac0ee8afc3192346b7022afac8fa48a02c)
on aura sur-le-champ :
![{\displaystyle \mathrm {D} _{s}=\mathrm {A} d_{s}+\mathrm {B} d_{s+1}+\mathrm {C} d_{s+2}+\mathrm {E} d_{s+3}+\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a19acaa6e64c9a8f365538696e6aea63e4bef81f)
Cette formule peut servir, comme l’on voit, à trouver les différences
de tous les ordres d’une série proposée dans laquelle on ne prendrait les termes que de deux en deux, ou de trois en trois, ou, en général, de
![{\displaystyle m}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a07d98bb302f3856cbabc47b2b9016692e3f7bc)
en
À l’égard des coefficients
on peut les trouver par le développement de la puissance
du polynôme, ou bien par la comparaison des termes qui contiendront les mêmes puissances de
après avoir pris les différentielles logarithmiques des deux membres. Ce dernier procédé donne les formules suivantes dont la loi se présente d’elle-même
![{\displaystyle {\begin{aligned}&\mathrm {A} =m^{s},\\&\mathrm {B} =s{\frac {m-1}{2}}\mathrm {A} ,\\&\mathrm {C} ={\frac {2s}{2}}{\frac {(m-1)(m-2)}{2.3}}\mathrm {A} +{\frac {s-1}{2}}{\frac {m-1}{2}}\mathrm {B} ,\\&\mathrm {E} ={\frac {3s}{3}}{\frac {(m-1)(m-2)(m-3)}{2.3.4}}\mathrm {A} +{\frac {2s-1}{3}}{\frac {(m-1)(m-2)}{2.3}}\mathrm {B} +{\frac {s-2}{3}}{\frac {m-1}{2}}\mathrm {C} ,\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9daef5105f99fd1d3b79a2648850e2674c602b13)
7. La même équation
![{\displaystyle \mathrm {D_{0}+D_{1}} =(d_{0}+d_{1})_{m},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a78ad862dfbb4a86edd9d34e10fd96032a1f61ff)
en extrayant pour ainsi dire la racine
donnera celle-ci
![{\displaystyle d_{0}+d_{1}=(\mathrm {D_{0}+D_{1}} )_{\frac {1}{m}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/714bccc70938502ce474c4757c070ae01333ddf8)
d’où l’on tire, à cause de
![{\displaystyle d_{1}=(\mathrm {D_{0}+D_{1}} )_{\frac {1}{m}}-\mathrm {D} _{0},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f8e2a30cf731cf8a3526a1cbda85b6551558e58c)
et de là
![{\displaystyle d_{s}=\left[(\mathrm {D_{0}+D_{1}} )_{\frac {1}{m}}-\mathrm {D} _{0}\right]_{s},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5052eca89676c94d4b35baf5eac2f7de0056a584)
formule qu’on développera comme la précédente ; et il est clair qu’il n’y
aura pour cela qu’à changer, dans les formules du numéro précédent,
![{\displaystyle \mathrm {D} }](https://wikimedia.org/api/rest_v1/media/math/render/svg/3c894024162d89aa8502758b2151cb579806ea9a)
en
![{\displaystyle d}](https://wikimedia.org/api/rest_v1/media/math/render/svg/e85ff03cbe0c7341af6b982e47e9f90d235c66ab)
et
![{\displaystyle m}](https://wikimedia.org/api/rest_v1/media/math/render/svg/0a07d98bb302f3856cbabc47b2b9016692e3f7bc)
en
![{\displaystyle {\frac {1}{m}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aab3d48b046a3fde9bc23b7c2f10d866dc9900bf)
Ainsi l’on aura, en général,
![{\displaystyle d_{s}=a\mathrm {D} _{s}+b\mathrm {D} _{s+1}+c\mathrm {D} _{s+2}+e\mathrm {D} _{s+3}+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f4439077fe1c665cca4fa25d1ac2b63a1390e5c8)
en faisant
![{\displaystyle {\begin{aligned}&a={\frac {1}{m^{s}}},\\&b=s{\frac {1-m}{2m}}a,\\&c={\frac {2s}{2}}{\frac {(1-m)(1-2m)}{2.3m^{2}}}a+{\frac {s-1}{2}}{\frac {1-m}{2m}}b,\\&e={\frac {3s}{3}}{\frac {(1-m)(1-2m)(1-3m)}{2.3.4m^{3}}}a+{\frac {2s-1}{3}}{\frac {(1-m)(1-2m)}{2.3m^{2}}}b+{\frac {s-2}{3}}{\frac {1-m}{2m}}c,\\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f509611f36e6a5967dadb6e02a8e0187bd3e67af)
Ces formules renferment la solution générale du Problème de Mouton. En effet, si l’on a à interpoler la série
![{\displaystyle \mathrm {T_{0},\ \ T_{1},\ \ T_{2}} ,\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/456ea06f2526be8c262a342bda550696bdd94013)
dont les différences
![{\displaystyle \mathrm {D_{0},\ \ D_{1},\ \ D_{2}} ,\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/dd734e575a3558360723c50d69ed65605992f2b1)
sont connues, et qu’on veuille partager en
parties chaque intervalle d’un terme à l’autre en y insérant
termes intermédiaires qui soient liés entre eux et avec ceux de la série donnée par une même loi, en supposant que
![{\displaystyle t_{0},\ \ t_{1},\ \ t_{2},\ \ t_{3},\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/a67098c6aeea0f0e61ec40cf2e1446ad81813bca)
soit la nouvelle série interpolée, il faudra par les conditions du Problème que les termes de cette série, pris de
en
coïncident avec les termes de la série donnée, en sorte qu’on ait
![{\displaystyle t_{0}=\mathrm {T} _{0},\quad t_{m}=\mathrm {T} _{1},\quad t_{2m}=\mathrm {T} _{2},\quad t_{3m}=\mathrm {T} _{3},\ldots ,\quad t_{sm}=\mathrm {T} _{s}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eeb08428a168f6c4a3fb7c19b417fb6e1b25fb24)
ce qui est l’équation fondamentale d’où nous sommes partis pour la so-
lution du Problème que nous venons de traiter (5). Ainsi l’on aura par les formules ci-dessus l’expression de la différence d’un ordre quelconque
![{\displaystyle s}](https://wikimedia.org/api/rest_v1/media/math/render/svg/01d131dfd7673938b947072a13a9744fe997e632)
de la série interpolée, au moyen des différences données de la série proposée. On connaîtra donc par là les différences successives
![{\displaystyle d_{1},\ \ d_{2},\ \ d_{3},\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/59b9acd673a9407d587b820ba2c231ba4f5a8c6a)
de la série
![{\displaystyle t_{0},\ \ t_{1},\ \ t_{2},\ \ t_{3},\ldots }](https://wikimedia.org/api/rest_v1/media/math/render/svg/a67098c6aeea0f0e61ec40cf2e1446ad81813bca)
dont le premier terme
et au moyen de ces différences on pourra trouver successivement tous les termes de cette série par des additions successives, comme nous l’avons vu dans le no 2, parce qu’en employant ces différences, un terme quelconque
de la série dont il s’agit sera exprimé par la formule
![{\displaystyle t_{0}+nd_{1}+{\frac {n(n-1)}{2}}d_{2}+{\frac {n(n-1)(n-2)}{2.3}}d_{3}+\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fc605fef1f7618146e2c4092d3921674160889f6)
dans laquelle
sera un nombre entier comme dans le no 2.
8. Si
est la dernière différence de la série donnée
![{\displaystyle \mathrm {T_{0},\ \ T_{1},\ \ T_{2}} ,\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/456ea06f2526be8c262a342bda550696bdd94013)
en sorte que
![{\displaystyle \mathrm {D} _{r+1}=0,\quad \mathrm {D} _{r+2}=0,\ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/aedefab6d7c7108153f4e1a87f36ef192602adb1)
on aura :
1o En faisant
![{\displaystyle d_{r}=a\mathrm {D} _{r}\quad {\text{et}}\quad a={\frac {1}{m^{r}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/eecf13c26d4d21be38766f6d634cb83993bbe739)
ce qui s’accorde-avec ce que Mouton avait trouvé par induction.
2o En faisant
![{\displaystyle d_{r-1}=a\mathrm {D} _{r-1}+b\mathrm {D} _{r},\quad {\text{et}}\quad a={\frac {1}{m^{r-1}}},\quad b={\frac {(r-1)(1-m)}{2m^{r}}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/361cbe3b29e8f5ffef30a3063e101f9b4a6c20a0)
3o En faisant
![{\displaystyle d_{r-2}=a\mathrm {D} _{r-2}+b\mathrm {D} _{r-1}+c\mathrm {D} _{r},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/14c80e4af01bf7ce029acb02420d2586174a80b0)
et
![{\displaystyle a={\frac {1}{m^{r-2}}},\quad b={\frac {(r-2)(1-m)}{2m^{r-1}}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6d2cb00ed82fa24b83acf93cefd098f7449e48c3)
![{\displaystyle c={\frac {(r-2)(1-m)(1-2m)}{2.3m^{r}}}+{\frac {(r-2)(r-3)(1-m)^{2}}{8m^{r}}}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f108c92447f3401d3d5b5d149491aa475c8cef42)
ainsi de suite. À l’égard des différences plus hautes que la rième, elles seront toutes nulles, comme on le voit par la formule en faisant
Ainsi il suffira de chercher les valeurs de
jusqu’à
et l’on pourra ensuite employer les additions successives pour former et continuer aussi loin qu’on voudra la Table des termes ![{\displaystyle t_{0},t_{1},t_{2},\ldots .}](https://wikimedia.org/api/rest_v1/media/math/render/svg/018e43602f49d279d5bc23f23b593e0fa8ee7ad8)
9. Si dans la formule (D) du no 6 on fait
infiniment grand, et par conséquent les différences représentées par
infiniment petites, on aura la relation entre les différences finies et les différentielles d’une fonction quelconque. Considérant donc une fonction
de
dont les différences finies soient désignées par
et les différentielles par
celles de
l’étant de même par
et
il n’y aura qu’à faire
et écrire
au lieu de
et
en plaçant suivant l’usage les indices de l’ordre des différences et des différentielles au haut des lettres
et
On aura ainsi
(E)
|
|
|
où il n’y aura qu’à développer le second membre comme une puissance
du polynôme
en ayant soin d’écrire
au lieu de ![{\displaystyle \left(d^{m}y\right)^{n}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a2c3d8ba025666c3c69a7c777fffb5913ba153ca)
Par cette formule on peut trouver les différences finies de tous les ordres au moyen des différentielles, et, si l’on fait
négatif, ce qui changera les différences en sommes et les différentielles en intégrales, on aura la somme d’un ordre quelconque des termes d’une série par les intégrales, et les différentielles de la fonction qui exprime le terme général de la série.
10. Réciproquement, si dans la même formule (D) on change
en
et
en
comme nous l’avons dit dans le no 7, et que par conséquent suivant l’hypothèse du numéro précédent, on écrive
au lieu de
au lieu de
et
au lieu de
qui devient alors une quantité infiniment petite, on aura
(F)
|
|
|
où il faudra de même développer le second membre comme une puissance
en ayant soin d’écrire
au lieu de ![{\displaystyle (\mathrm {D} ^{m}y)^{n}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1fa33f1ba0dec43020b8101478808892bdd0faf8)
Cette formule donnera les différentielles de tous les ordres d’une fonction quelconque par le moyen de ses différences finies, et, si l’on fait
négatif, elle donnera les intégrales de la fonction par le moyen des sommes et des différences des valeurs de la même fonction répondantes aux valeurs de
dont
est la différence finie ; ce qui revient à la détermination des aires par les ordonnées équidistantes. [Voyez là-dessus le Mémoire cité dans le volume de 1772[2].]
11. On peut traiter par les mêmes-principes les séries doubles, triples,
c’est-à-dire celles dont les termes varient de deux ou de pluf sieurs manières différentes, et qui forment des Tables à double, à triple entrée,
Soit, par exemple, la série double
![{\displaystyle {\begin{array}{rrrrrr}\mathrm {T} _{0,0}&\mathrm {T} _{1,0}&\mathrm {T} _{2,0}&\mathrm {T} _{3,0}&\mathrm {T} _{4,0},&\ldots \\\mathrm {T} _{0,1}&\mathrm {T} _{1,1}&\mathrm {T} _{2,1}&\mathrm {T} _{3,1}&\mathrm {T} _{4,1},&\ldots \\\mathrm {T} _{0,2}&\mathrm {T} _{1,2}&\mathrm {T} _{2,2}&\mathrm {T} _{3,2}&\mathrm {T} _{4,2},&\ldots \\\mathrm {T} _{0,3}&\mathrm {T} _{1,3}&\mathrm {T} _{2,3}&\mathrm {T} _{3,3}&\mathrm {T} _{4,3},&\ldots \\\mathrm {T} _{0,4}&\mathrm {T} _{1,4}&\mathrm {T} _{2,4}&\mathrm {T} _{3,4}&\mathrm {T} _{4,4},&\ldots \\\ldots &\ldots &\ldots &\ldots &\ldots ,&\ldots \end{array}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2ef58385247f13555e906614c810e92801d6cdb4)
dans laquelle un terme quelconque comme
a deux indices, le premier
pour marquer son rang dans la direction horizontale, et le second
pour marquer son rang dans la direction verticale.
Si l’on désigne de même par
![{\displaystyle \mathrm {D_{0,0},\ \ D_{1,0},\ \ D_{2,0},\ldots ,\quad D_{0,1},\ \ D_{1,1},\ \ D_{2,1},\ldots ,\quad D_{0,2},\ \ D_{1,2},\ \ D_{2,2}} ,\ldots ,\quad \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c40fe0881894817397882f5d52b84e1a85a8ee5b)
et, en général, par
n les différences successives de ces termes, de manière que le premier indice
indique des différences prises dans le sens horizontal, et le second
des différences prises dans le sens vertical, on aura
![{\displaystyle {\begin{aligned}&\mathrm {D_{1,0}=T_{1,0}-T_{0,0}} \\&\mathrm {D_{2,0}=T_{2,0}-2T_{1,0}+T_{0,0}} \\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\&\mathrm {D_{0,1}=T_{0,1}-T_{0,0}} \\&\mathrm {D_{0,2}=T_{0,2}-2T_{0,1}+T_{0,0}} \\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \\&\mathrm {D_{1,1}=T_{1,1}-T_{0,1}-(T_{1,0}-T_{0,0})=T_{1,1}-T_{1,0}-T_{0,1}+T_{0,0}} \\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots ,\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bf3584bba5ef3bfe9d6e3b6e60cc9bca3eecec76)
et ainsi de suite ; et l’on pourra trouver tout de suite l’expression de
par une opération semblable à celle du no 1.
12. Car, puisque
![{\displaystyle \mathrm {D_{1,0}=T_{1,0}-T_{0,0}} \quad {\text{et}}\quad \mathrm {D_{0,1}=T_{0,1}-T_{0,0}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2e936762724b47c53555a5b77db0d255d3d0c1e7)
on aura d’abord, en élevant aux indices
et ![{\displaystyle n,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/397bfafc701afdf14c2743278a097f6f2957eabb)
![{\displaystyle \mathrm {D} _{m,0}=\mathrm {\left(T_{1,0}-T_{0,0}\right)} _{m},\quad \mathrm {D} _{0,n}=\mathrm {\left(T_{0,1}-T_{0,0}\right)} _{n},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/53dee7fed33ce0a3ed91a4f6296102d3f364d626)
et, multipliant ces équations l’une par l’autre, en observant de rapporter toujours les premiers indices aux premiers et les seconds aux seconds, il viendra
![{\displaystyle \mathrm {D} _{m,n}=\mathrm {\left(T_{1,0}-T_{0,0}\right)} _{m}\mathrm {\left(T_{0,1}-T_{0,0}\right)} _{n},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8054930dcc8e012dc7436e7b407b38d034d3b662)
c’est-à-dire, en développant d’abord les deux facteurs,
(G)
|
|
|
Maintenant il n’y aura plus qu’à faire le produit de ces quantités, en observant la règle prescrite, et l’on aura
![{\displaystyle {\begin{aligned}\mathrm {D} _{m,n}=&\mathrm {T} _{m,n}-m\mathrm {T} _{m-1,n}+{\frac {m(m-1)}{2}}\mathrm {T} _{m-2,n}-\ldots \\&-n\mathrm {T} _{m,n-1}-nm\mathrm {T} _{m-1,n-1}-n{\frac {m(m-1)}{2}}\mathrm {T} _{m-2,n-2}+\ldots \\&+{\frac {n(n-1)}{2}}\mathrm {T} _{m,n-2}-{\frac {n(n-1)}{2}}m\mathrm {T} _{m-1,n-2}+\ldots \\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/54fb941bdc4ea5d37d010e8abb1b55955c6bd78e)
13. Réciproquement on aura l’expression d’un terme quelconque
par les différences, au moyen des équations
![{\displaystyle \mathrm {T_{1,0}=D_{0,0}+D_{1,0}} \quad {\text{et}}\quad \mathrm {T_{0,1}=D_{0,0}+D_{0,1}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/987a8a034b2ac6a20231fe68c2f4731d19a3d0bd)
à cause de
![{\displaystyle \mathrm {T_{0,0}=D_{0,0}} ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/bea8c530e106236844c17e0adc1f82d84eb39852)
comme il résulte de la formule précédente ; car, en élevant d’abord la première à l’indice
et la seconde à l’indice
les multipliant ensuite l’une par l’autre, on aura
![{\displaystyle \mathrm {T} _{m,n}=\mathrm {\left(D_{0,0}+D_{1,0}\right)} _{m}\mathrm {\left(D_{0,0}+D_{0,1}\right)} _{n},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1db89ea2db9f58d373718f5e40e387a8e265e11c)
d’où l’on tirera comme dans le numéro précédent
![{\displaystyle {\begin{aligned}\mathrm {T} _{m,n}=&\left[\mathrm {D} _{0,0}+m\mathrm {D} _{1,0}+{\frac {m(m-1)}{2}}\mathrm {D} _{2,0}+\ldots \right]\\\times &\left[\mathrm {D} _{0,0}+n\ \mathrm {D} _{0,1}+{\frac {n\ (n\ -1)}{2}}\mathrm {D} _{0,2}+\ldots \right]\\=&\mathrm {D} _{0,0}+m\mathrm {D} _{1,0}+{\frac {m(m-1)}{2}}\mathrm {D} _{2,0}+\ldots \\&+n\mathrm {D} _{0,1}+nm\mathrm {D} _{1,1}+n{\frac {m(m-1)}{2}}\mathrm {D} _{2,1}+\ldots \\&+{\frac {n(n-1)}{2}}\mathrm {D} _{0,2}+{\frac {n(n-1)}{2}}m\mathrm {D} _{1,2}+{\frac {n(n-1)}{2}}{\frac {m(m-1)}{2}}\mathrm {D} _{2,2}+\ldots \\&\ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots \ldots .\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/7553b2b438ebd779b2673dc1b37f8acbcbb3cefb)
14. On peut aussi étendre par des opérations semblables les formules des nos 6 et 7 pour les interpolations aux séries doubles. En eflet, si l’on imagine la nouvelle série
![{\displaystyle t_{0,0},\ \ t_{1,0},\ \ t_{2,0},\ldots ,\quad t_{0,1},\ \ t_{0,2},\ \ t_{0,3},\ldots ,\quad t_{1,1},\ \ t_{2,1},\ \ t_{3,1},\ldots ,\quad \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/014cebe6b0a0dacb8c7aa04fa23c0e9156f1786c)
dont les différences soient désignées par
![{\displaystyle d_{0,0},\ \ d_{1,0},\ \ d_{2,0},\ldots ,\quad \ldots ,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/8e14655d196d71197a78e34a3fa0ec7ac143faae)
et qui soit telle qu’un terme quelconque
soit égal à
en faisant
et
on aura
![{\displaystyle t_{m,0}=\mathrm {T} _{1,0},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/6b4ac9eeda6eee344af5491a8361da6d0424710d)
et, faisant
on aura
![{\displaystyle t_{0,n}=\mathrm {T} _{0,1}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1b3b1e9d7b56ea32770c5bcc9e5c216b45f4de85)
deux équations qui, étant élevées l’une à l’indice et l’autre à l’indice
et multipliées ensuite l’une par l’autre, redonnent l’équation générale. Or ces deux équations, en y substituant
à la place de
à la place de
et de même
à la place de
à la place de
par conséquent
à la place de
à la place de
deviennent
![{\displaystyle \mathrm {D} _{0,0}+\mathrm {D} _{1,0}=(d_{0,0}+d_{1,0})_{m}\quad {\text{et}}\quad \mathrm {D} _{0,0}+\mathrm {D} _{0,1}=(d_{0,0}+d_{0,1})_{n}\,;}](https://wikimedia.org/api/rest_v1/media/math/render/svg/2357436e29d8951ef48b67af131536c8b2d3c4f6)
d’où, à cause de
![{\displaystyle \mathrm {D} _{0,0}=\mathrm {T} _{0,0}=t_{0,0}=d_{0,0}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/260bd18e299124b76f34acfe24c298dc53fd2cb1)
comme il résulte des formules ci-dessus, on aura
![{\displaystyle {\begin{alignedat}{2}\mathrm {D} _{1,0}=&(d_{0,0}+d_{1,0})_{m}-d_{0,0},\qquad &\mathrm {D} _{s,0}=&\left[(d_{0,0}+d_{1,0})_{m}-d_{0,0}\right]_{s},\\\mathrm {D} _{0,1}=&(d_{0,0}+d_{0,1})_{n}\ -d_{0,0},&\mathrm {D} _{0,r}=&\left[(d_{0,0}+d_{0,1})_{n}\ -d_{0,0}\right]_{r}\,;\end{alignedat}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9a1b3290efdc9e5b2ad59f37e1032190a19d928f)
par conséquent, en multipliant l’une par l’autre,
![{\displaystyle \mathrm {D} _{s,r}=\left[(d_{0,0}+d_{1,0})_{m}-d_{0,0}\right]_{s}\times \left[(d_{0,0}+d_{0,1})_{n}-d_{0,0}\right]_{r},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c6c9418df1dcce3e31233d9bc6973eb82bf90e84)
et, développant les deux facteurs,
(H)
|
|
|
De sorte qu’il n’y aura plus qu’à développer suivant les indices
et
comme on l’a fait dans le no 6, et ensuite effectuer la multiplication comme dans les deux numéros précédents. On aura ainsi pour les séries doubles une formule semblable à celle du no 6 pour les séries simples.
Pareillement, par l’extraction des indices
et
on aura
![{\displaystyle {\begin{aligned}d_{0,0}+d_{1,0}=&\mathrm {\left(D_{0,0}+D_{1,0}\right)} _{\frac {1}{m}},\\d_{0,0}+d_{0,1}=&\mathrm {\left(D_{0,0}+D_{0,1}\right)} _{\frac {1}{n}}\,;\end{aligned}}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/9f4e53025784a1fe3235780df360c28aaf1c7c4b)
et de là on tirera par un procédé semblable au précédent
![{\displaystyle d_{s,r}=\left[\mathrm {\left(D_{0,0}+D_{1,0}\right)} _{\frac {1}{m}}-\mathrm {D} _{0,0}\right]_{s}\times \left[\mathrm {\left(D_{0,0}+D_{0,1}\right)} _{\frac {1}{n}}-\mathrm {D} _{0,0}\right]_{r},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/dcf93bd0c5ffd50a5e841fb0eca17ce1365079f5)
savoir
(I)
|
|
|
où il faudra encore développer le second membre comme dans la formule (H) ci-dessus. Cette expression de
sera analogue à celle de
du no 7, et pourra servir aux mêmes usages pour l’interpolation des Tables à double entrée.
15. Enfin, si l’on suppose
et
infinis, les différences désignées par
deviendront infiniment petites, et l’on aura pour les fonctions à deux variables des formules analogues à celles des nos 9 et 10. Ainsi, en regardant
comme fonction de
et
et désignant par la caractéristique
les différences finies et par la caractéristique
les différentielles, on fera
![{\displaystyle m={\frac {\mathrm {D} x}{dx}},\quad n={\frac {\mathrm {D} y}{dy}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/52c2e399c85a81eea881c7b54bc13a84d61a6e9f)
et la formule (H) se transformera en celle-ci
(K)
|
|
|
dont il faudra développer le second membre par les méthodes ordinaires, en regard’ant et traitant les exposants des différentielles de
comme des exposants de puissances. Cette formule donnera la valeur de la différence finie de
de l’ordre
par rapport à
et de l’ordre
par rapport à
exprimée par les différentielles de
relativement à
et ![{\displaystyle y.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/83f72471aff7c6fbb27df0f971283a068efe091f)
16. Faisant les mêmes substitutions dans la formule (I), et écrivant les indices de l’ordre des différences et des différentielles en forme d’exposants suivant la notation reçue, on aura cette transformée
![{\displaystyle {\frac {d^{s+r}z}{dx^{s}dy^{r}}}=\left({\frac {\mathrm {D} ^{1,0}z-{\frac {1}{2}}\mathrm {D} ^{2,0}z+{\frac {1}{3}}\mathrm {D} ^{3,0}z-\ldots }{\mathrm {D} x}}\right)^{s}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3fa73188d38c7f7f7be467d1e2e01f206917adde)
![{\displaystyle \times \left({\frac {\mathrm {D} ^{0,1}z-{\frac {1}{2}}\mathrm {D} ^{0,2}z+{\frac {1}{3}}\mathrm {D} ^{0,3}z-\ldots }{\mathrm {D} y}}\right)^{r},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5cc4375a17728d4f257c7546036b5ee592b9c0db)
qu’on développera suivant les méthodes ordinaires, en observant relativement aux exposants des différences les mêmes règles que pour ceux des puissances, et ayant soin de calculer séparément les premiers exposants qui sont relatifs à la variable
et les seconds qui se rapportent à la variable ![{\displaystyle y.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/83f72471aff7c6fbb27df0f971283a068efe091f)