Annales de mathématiques pures et appliquées/Tome 06/Analise transcendante, article 5
# Annales de mathématiques pures et appliquées/Tome 06/Analise transcendante, article 5

### Réflexions sur la méthode qui sert de base au précédent mémoire, et applications diverses de cette méthode ;

La méthode dont M. Kramp vient de faire usage, dans le précédent mémoire, pour résoudre le problème des quadratures, est extrêmement remarquable, et nous paraît tout-à-fait digne de l’attention des géomètres. Elle semble devoir être très-féconde en applications curieuses et utiles ; et nous n’hésitons pas à la regarder comme une des plus belles et des plus ingénieuses inventions d’analise qui aient eu lieu dans ces derniers temps.

L’esprit de cette méthode consiste proprement à chercher, à dessein, des résultats moins approchés que celui dont on est déjà en possession, et à les employer à perfectionner celui-là. C’est exactement prendre de l’élan ; c’est reculer pour mieux sauter. Les détails dans lesquels nous allons entrer pourront faire entrevoir de combien d’applications variées cette méthode peut être susceptible ; ils montreront en même-temps que l’approximation qu’elle est capable de fournir, dans tous les cas, n’a pour ainsi dire d’autre limite que celles de la patience du calculateur. Mais, avant d’entrer en matière, arrêtons-nous encore un moment sur le problème des quadratures.

I. Quelque rapide que puisse être un procédé approximatif, ce procédé doit être jugé imparfait, s’il ne renferme pas en soi quelque moyen d’apprécier l’erreur à laquelle son usage peut exposer. Or, telle serait la méthode des quadratures, développées dans le précédent mémoire, si on ne lui faisait pas subir une légère modification. Cette modification consiste à substituer successivement aux trapèzes des rectangles inscrits et des rectangles circonscrits. Cela conduira à deux résultats, l’un plus grand et l’autre plus petit que le véritable, et dont la différence donnera conséquemment la limite de l’erreur dont chacun d’eux se trouvera affecté. À la vérité, toutes choses égales d’ailleurs, ces résultats seront moins approchés que ceux qu’on déduirait de l’usage des trapèzes ; mais il nous paraît qu’on ne doit pas balancer à sacrifier quelque chose du côté de la précision et de la rapidité, lorsqu’il s’agit de remplir une condition sans laquelle aucun procédé approximatif ne saurait être employé avec quelque sécurité. Nous verrons d’ailleurs bientôt que cet inconvénient disparaît presque totalement, par un emploi convenable de la méthode.

Ceci suppose, au surplus, qu’entre les limites de l’intégrale, les ordonnées de la courbe qu’il s’agit de quarrer sont toujours croissantes ou toujours décroissantes ; mais on sait que, dans le cas contraire, on peut toujours décomposer l’intégrale en plusieurs parties telles que, pour chacune d’elles, cette condition se trouve remplie.

Nous appliquerons uniquement ces réflexions au cas où le diviseur général est 6, et ses aliquotes 1, 2, 3, 6. Soient α, β, γ, δ, ε, ζ, η, les sept ordonnées équidistantes que, pour fixer les idées, nous supposerons perpétuellement croissantes ; prenons de plus pour unité, comme dans le précédent mémoire ; l’intervalle qui sépare les ordonnées extrêmes. En considérant les rectangles inscrits dont les bases sont successivement 1/6, 1/3, 1/2, 1, nous aurons, pour la somme de leurs aires,

=1/6, 6

… 1/6 (α+β+γ+δ+ε+ζ)=1/6 a′,

=1/3, 3

… 2/6 (α+γ+ε)=1/6 b′,

=1/2, 2

… 3/6 (α+δ)=1/6 c′,

=1, 1

… 6/6 α=1/6 d′.

Si nous passons ensuite aux rectangles circonscrits, nous trouverons les sommes d’aires ainsi qu’il suit

=1/6, 6

… 1/6 (β+γ+δ+ε+ζ+η)=3/6 a′,

=1/3, 3

… 2/6 (γ+ε+η)=1/6 b′,

=1/2, 2

… 3/6 (δ+η)=1/6 c′,

=1, 1

… 1/6 η=1/6 d′.

Nous aurons toujours d’ailleurs la formule

A=(1296a′−567b′+112c′−d′)/5040 ;

en y faisant donc successivement les deux substitutions, il viendra

∫Xd⁡x> (82α+216(β+ζ)+27(γ+ε)+272δ)/840.

∫Xd⁡x< (216(β+ζ)+27(γ+ε)+272δ+82η)/840.

La différence

41(η−α)/420,

est la limite de l’erreur que pourra entraîner l’emploi de l’une ou de l’autre de ces deux formules, dont la demi-somme est précisément la formule de M. Kramp, ainsi que ce cela doit être.

Si l’on applique ces formules aux deux exemples de l’auteur, c’est-à-dire, à la recherche du logarithme naturel de 2; et à celle du nombre ϖ/4 ; comme, dans l’un et dans l’autre cas, on a α=1 et η=0, 5, on aura η−α=−1/2 ; de sorte que la limite de l’erreur est 41/840 ; ou environ 1/20. Nous allons voir au surplus que la résolution du problème des quadratures peut encore être présentée sous une autre forme qui, sans exiger un grand nombre de divisions de l’étendue de l’intégrale, est néanmoins susceptible d’une approximation presque illimitée.

Supposons toujours qu’il soit question d’obtenir ϖ/4 ou, ce qui revient au même, l’intégrale de d⁡x/(1+x²), entre 0 et 1 en posant y=1/(1+x²), et divisant d’abord l’intervalle en cinq parties égales seulement ; nous aurons

x=0, 1/5, 2/5, 3/5, 4/5, 1

y=1, 25/26, 25/29, 25/34, 25/41, 1/2.

Comme ici les ordonnées sont continuellement décroissantes, les rectangles inscrits, auxquels nous nous bornerons, et qui, ayant, pour base commune, auront successivement pour hauteur les cinq dernières ordonnées, seront

Premier… 1/5.25/26 =5/26 =0, 1923077, Deuxie^{`}me… 1/5.25/29 =5/29 =0, 1724138, Troisie^{`}me… 1/5.25/34 =5/34 =0, 1470588, Quatrie^{`}me… 1/5.25/41 =5/41 =0, 1219512, Cinquie^{`}me… 1/5.1/2 =1/10 =1, 000000.}somme =0, 6337315.

En multipliant ce résultat par 4, on obtiendra pour première valeur approchée du nombre ϖ

ϖ=2, 9349260=A.

Pour obtenir une valeur plus approchée, cherchons-en une suite d’autres qui le soient moins. Soit d’abord divisée l’étendue de l’intégrale en quatre parties égales ; nos quatre rectangles inscrits seront alors tels qu’il suit :

Premier… 1/4.16/17 =4/17 =0, 2352941, Deuxie^{`}me… 1/4.16/20 =1/5 =0, 2000000, Troisie^{`}me… 1/4.16/25 =4/25 =0, 1600000, Quatrie^{`}me… 1/4.1/2 =1/8 =0, 1250000.}somme =0, 7202941.

Ce résultat, multiplié par 4, donnera pour seconde valeur moins approchée de ϖ

ϖ=2, 8811764=B.

Divisons le même intervalle en trois parties égales seulement, les rectangles inscrits résultans seront

Premier… 1/3.9/10 =3/20 =0, 3000000, Deuxie^{`}me… 1/3.9/13 =3/13 =0, 2307692, Troisie^{`}me… 1/3.1/2 =1/6 =0, 1666667.}somme =0, 6974359.

Ce résultat, multiplié par 4, donnera pour troisième valeur, moins approché que la précédente, du nombre ϖ

ϖ=2, 7897436=C.

Divisons ensuite cet intervalle en deux parties égales seulement, les deux rectangles inscrits correspondans seront

Premier… 1/2.4/5 =2/5 =0, 4000000, Deuxie^{`}me… 1/2.1/2 =1/4 =0, 2500000.}somme =0, 6500000.

résultat qui, multiplié par 4, donne, pour quatrième valeur encore moins approchée du nombre ϖ,

ϖ=2, 6000000=D.

Considérant enfin l’intervalle entier, nous aurons pour le rectangle inscrit 1.1/2=1/2=0, 5000000 qui, multiplié par 4, donnera, pour la dernière valeur, la moins approchée de ϖ,

ϖ=2, 0000000=E.

Il est évident qu’aucune des quantités A, B, C, D, E n’est la valeur de ϖ, et qu’elles sont toutes plus petites que cette valeur ; mais, si on les considère comme répondant respectivement aux indices, 5, 4, 3, 2, 1, il est clair que la valeur de ϖ répondra à l’indice ∞ ; puisque, pour cet indice, on sera dans le même cas que si l’on avait considéré une infinité de rectangles inscrits infiniment petits. Donc, à l’inverse, si l’on considère respectivement A, B, C, D, E comme une suite de termes répondant aux indices 1/5, 1/4, 1/3, 1/2, 1; le terme de cette suite répondant à l’indice 1/∞ ou zéro sera la valeur exacte de ϖ ; et, comme il en sera encore évidemment de même en rendant tous les indices 60 fois plus grands ; il s’ensuit que, si l’on construit une courbe telle qu’aux abscisses 12, 15, 20, 30, 60 répondent respectivement les ordonnées A, B, C, D, E, la valeur de ϖ sera l’ordonnée de cette courbe repondant à l’abscisse zéro.

Or, on a vu, dans le précèdent mémoire, qu’en supposant, pour plus de simplicité, que cette courbe est parabolique, et que son équation ne renferme que des puissances paires de l’abscisse, si a, b, c, d, e représentent les quarrés des abscisses qui répondent respectivement aux ordonnées A, B, C, D, E, on doit avoir sensiblement

ϖ=bcdeA/(b−a)(c−a)(d−a)(e−a)+cdeaB/(c−b)(d−b)(e−b)(a−b)+deabC/(d−c)(e−c)(a−c)(b−c) +eabcD/(e−d)(a−d)(b−d)(c−d)+abcdE/(a−e)(b−e)(c−e)(d−e).(V)

Faisant donc, dans cette formule, A=144, b=225, c=400, d=900, e=3600, elle deviendra, toutes réductions faites

ϖ=(1953125A−2097152B+531441C−24576D+42E)/7.8.8.9.9.10 ;

formule dans laquelle il n’est plus question que de substituer les valeurs ci-dessus. On trouve ainsi

1953125A=5732277, 3437500, 531441C=1482584, 1285276, 42E=84, 0000000;}72144945, 4722776.

2097152B=6042264, 8496128, 24576D=63897, 6000000;} 6106162, 4496128.

7.8.8.9.9.10ϖ=1108783, 0226648.

ϖ=3, 0555088=A′.

Cette valeur est encore peu approchée ; mais on doit en être peu surpris, si l’on considère que d’abord nous avons substitué des rectangles aux trapèzes, et qu’en autre nous n’en avons employé que cinq au plus.

On se tromperait toutefois si l’on se figurait que c’est là tout le degré d’approximation auquel il soit possible de parvenir, avec d’aussi faibles moyens. On peut, en effet, traiter ce nouveau résultat A′ comme nous avons traité le premier A; c’est-à-dire, chercher des résultats moins approchés que lui et les employer à le perfectionner.

Supposons donc que nous n’ayons pas été au-delà de quatre divisions ; c’est-à-dire, faisons abstraction de la valeur A; nous pourrons alors considérer B, C, D, E, comme répondant respectivement aux indices 1/4, 1/3, 1/2, 1, en multipliant par 12, comme répondant aux indices 3, 4, 6, 12; nous aurons alors à employer la formule

ϖ=cdeB/(c−b)(d−b)(e−b)+debC/(d−c)(a−c)(b−c)+ebcD/(e−d)(b−d)(c−d)+bcdE/(b−e)(c−e)(d−e) ;(IV)

dans laquelle il faudra faire b=9, c=16, d=36, e=144; ce qui donnera

ϖ=(8192B−6561C+896D−7E)/5.7.8.9.

En substituant donc nous aurons

8192B=23602, 5970688, 896D=2329, 6000000;}25932, 1970688,

6561C=18303, 5077596, 7E=14, 0000000;}18317, 5077596;

5.7.8.9 ϖ=7614, 6893092;

ϖ=3, 0217025=B′.

En n’allant pas au-delà de trois divisions ; c’est-à-dire, en faisant abstraction des valeurs A, B, nous pourrons considérer C, D, E comme répondant respectivement aux indices 1/3, 1/2, 1, ou, en multipliant par 6, comme repondant aux indices 2, 3, 6 ; nous aurons alors à employer la formule

ϖ=deC/(d−c)(e−c)+ecD/(e−d)(c−d)+cdE/(c−e)(d−e) ;(III)

dans laquelle il faudra faire a=4, b=9, c=16; ce qui donnera

ϖ=(243C−128D+5E)/4.5.6.

En substituant donc, nous aurons

243C=677, 9076948, 5E=10, 0000000;}677, 9076948

128D=332, 8000000

4.5.6 ϖ=355, 1076948

ϖ=2, 9592308=C′.

En ne faisant ensuite que deux divisions, nous aurons

ϖ=eD/(e−d)+dE/(d−e), (III)

ϖ=4D−E/3 ;(II’)

ce qui donnera

4D=10, 4000000E= 2, 0000000

3ϖ=8, 4000000

ϖ=2, 8000000=D′.

En ne considérant enfin qu’un seul rectangle inscrit, nous aurons de nouveau, comme ci-dessus,

Nous pouvons présentement traiter les valeurs de moins en moins approchées A′, B′, C′, D′, E′, comme nous avions traité les valeurs A, B, C, D, E, c’est-à-dire, les substituer dans la formule

ϖ=(1953125A′−2097152B′+531441C′−24576D′+42E′)/7.8.8.9.9.10 ;

ce qui donnera

1953125A′=5967790, 6250000, 531441C′=1572656, 5755828, 42E′=84, 0000000;}7540531, 2005828,

2097152B′=6336969, 4412800, 24576D′=68812, 8000000;} 6405782, 2412800;

7.8.8.9.9.10 ϖ=1134748, 9593028.

ϖ=3, 1270639=A″.

Voilà présentement une valeur un peu plus approchée que les valeurs A′ et A; or, de même que nous avons déduit A′, B′, C′, D′, E′ de A, B, C, D, E nous pourrons déduire A″, B″, C″, D″, E″ de A′, B′, C′, D′, E′; et, en continuant toujours ainsi, nous parviendrons à des valeurs de plus en plus approchées ; à la vérité, le procédé peut paraître un peu long ; mais il l’aurait été beaucoup moins, si nous ne nous étions ; dès l’abord, bornés à dessein à cinq divisions de l’intégrale.

En procédant, comme il vient d’être dit, on aura

A″=3, 1270639, B″=3, 1083625, C″=3, 0891090,

D″=3, 0666667, E″=2, 0000000;

d’où on conclura, par la formule (V’), une nouvelle valeur de ϖ.

Si, ne connaissant pas, à l’avance, la valeur exacte de ϖ, on voulait juger du degré d’approximation obtenu après un certain nombre de pareilles opérations, il ne s’agirait que de faire un semblable calcul sur les rectangles circonscrits. On n’adopterait alors dans la valeur de ϖ que les chiffres décimaux communs aux deux résultats. Nous allons voir, au surplus, qu’en suivant toujours l’esprit de la même méthode on peut se procurer bien plus rapidement une valeur approchée du nombre ϖ, et ce sera là notre première application.

II. Supposons, pour un moment, que la géométrie n’offre absolument aucun moyen de calculer, même par approximation, les périmètres des polygones réguliers au-delà de six côtés. Nous allons voir que, tandis que les procédés ordinaires, étendus jusqu’au polygone de 96 côtés, donnent une valeur qui n’est exacte que dans les deux premiers chiffres décimaux, notre méthode, au contraire, bornée à l’hexagone, donne un résultat qui n’est fautif que dans la sixième décimale seulement.

Observons auparavant que deux diamètres qui se confondent, dans un cercle dont le rayon est un, forment un véritable polygone régulier inscrit de deux côtés, dont le périmètre est quatre. En conséquence, nous aurons les demi-périmètres des polygones réguliers inscrits au cercle dont le rayon est un ainsi qu’il suit :

… 2 =2, 0000000=E,

……… 3/2 √3 =2, 5980762=D,

…… 2√2 =2, 8284272=C,

……… 5/4 √(10−2√5)=2, 9389265=B,

……… 3 =3, 0000000=A.

Si l’on considère les demi-périmètres A, B, C, D, E comme répondant respectivement aux indices 1/6, 1/5, 1/4, 1/3, 1/2, 1, le nombre ϖ devra répondre à l’indice 1/∞=0, et il en sera encore de même, en prenant les indices 60 fois plus grands 10, 12, 15, 20, 30. Faisant donc, dans la formule (V) a=100, b=144, c=225, d=400, e=900, elle deviendra

ϖ=(1469664A−1953125B+720966C−72171D+1056E)/6.6.6.7.10.11 ;

ce qui donnera, en substituant,

1469664A=4408992, 000000000, 720966C=2039001, 8547712, 1056E=2112, 0000000;}6450105, 8547712,

1953125B=5740090, 8203125, 72171D=187505, 7554302;}5927596, 5757427;

6.6.6.7.10.11 ϖ=522509, 2790285.

ϖ=3, 1415902=A′

résultat qui ne commence à être fautif qu’à la sixième décimale.

Si, ne connaissant point à l’avance la valeur exacte du nombre ϖ, on voulait juger du degré d’approximation de ce résultat, il suffirait de faire un semblable calcul relativement aux polygones circonscrits ; et l’on n’admettrait ensuite, dans la valeur approchée du nombre ϖ, que les chiffres décimaux communs aux deux résultats.

On aurait tort de penser au surplus que l’approximation à laquelle nous venons de parvenir est toute celle que peut donner la considération des cinq premiers polygones réguliers ; si, en effet, nous nous arrêtons successivement au 4.ᵉ, au 3.ᵉ, au 2.ᵉ, et au 1.ᵉʳ en désignant respectivement par B′, C′, D′, E′ les valeurs approchées de ϖ résultant de leur considération ; nous pourrons considérer A′, B′, C′, D′, E′, comme des termes répondant respectivement aux indices 1/5, 1/4, 1/3, 1/2, 1; et, en cherchant, comme ci-dessus, le terme qui répond à l’indice 1/∞ ou 0, ce sera une valeur plus approchée de ϖ.

III. Notre deuxième application aura encore pour objet la recherche du nombre ϖ, mais nous y procéderons de manière à faire voir comment la méthode dont nous cherchons ici à étendre l’usage s’applique à la sommation des séries convergentes, dont on connaît seulement un petit nombre des premiers termes, sans que même il soit aucunement besoin d’en connaître la loi.

Prenons la série connue de Leibnitz.

ϖ=4(1−1/3+1/5−1/7+1/9−1/11+…);

en réduisant chaque terme de rang pair avec le terme de rang impair qui le précède immédiatement, elle deviendra

ϖ=8( 1/1.3+1/5.7+1/9.11+1/13.15+…);

Nous allons essayer de la sommer au moyen de ses six premiers termes seulement.

On a

1.ᵉʳ terme 1/1.3=0, 3333333, 1.ᵉʳ terme=0, 3333333=F, 2.ᵉ… 1/5.7=0, 0285714, Som.des deux 1.ᵉʳˢ=0, 3619047=E, 3.ᵉ… 1/9.11=0, 0101010, des trois 1.ᵉʳˢ=0, 3720057=D, 4.ᵉ… 1/13.15=0, 0051282, des quatre 1.ᵉʳˢ=0, 3771339=C, 5.ᵉ… 1/17.19=0, 0030960, des cinq 1.ᵉʳˢ=0, 3802299=B, 6.ᵉ… 1/21.23=0, 0020704, des six 1.ᵉʳˢ=0, 3823003=A.

Si nous considérons ces nombres A, B, C, D, E, F comme répondant respectivement aux indices 1/6, 1/5, 1/4, 1/3, 1/2, 1, il est évident que ϖ/8 devra réppndre à l’indice 1/∞=0 et qu’il en sera encore de même pour les indices 60 fois plus grands 10, 12, 15, 20, 30, 60. Faisant donc, dans la formule du n.ᵒ 13 du précédent mémoire, a=100, b=144, c=225, d=400, e=900, f=3600, on aura

4.7.9.9.10.10.11 ϖ=

181398528A−244140625B+92274688C

−9743085D+168960E−66F;

ce qui donnera en substituant

181398528A=69348711, 6739584, 92274688C=34799912, 9567232, 168960E=61147, 4181120;}104209772, 0487936,

244140625B=92829565, 4296875, 9743085D=3624483, 1555845, 66F=22, 0000000;}96454070, 5852720;

4.7.9.9.10.10.11 ϖ= 7755701, 4635216;

ϖ=3, 1087468=A′.

On ne sera pas surpris du peu d’exactitude de cette valeur, si l’on fait attention à l’extrême lenteur de la série, qui tend sans cesse à n'être plus convergente.

Il est d’ailleurs aisé ici, comme dans les précédens exemples, de se procurer une valeur plus approchée, en en cherchant d’autres qui le soient moins ; si, en effet, on désigne respectivement par B′, C′, D′, E′, F′ les valeurs qu’on obtient, en se bornant successivement à cinq, quatre, trois, deux et un termes, on aura

A′=3, 1087468, B′=3, 0985080, C′=3, 0815390, D′=3, 0493500, E′=2, 9714280, F′=2, 6666666.

En substituant ces valeurs dans la précédente formule à la place de A, B, C, D, E, F et divisant le résultat par 8, attendu qu’elles n’expriment plus ici ϖ/8 mais ϖ lui-même, on obtiendra

ϖ=3, 1360361=A″,

valeur plus approchée que la précédente, et de laquelle il serait facile, par les mêmes moyens, d’en déduire d’autres qui le soient davantage encore. Ainsi, malgré le peu de convergence de la série, il ne faudra qu’un peu de patience pour obtenir, à l’aide de ses six premiers termes seulement, des valeurs de plus en plus approchées de la somme de tous ses termes.

Si, ne connaissant pas à l’avance la valeur rigoureuse du nombre ϖ, on voulait juger de la précision des résultats successivement obtenus, on remarquerait que la série de Leibnitz peut aussi être mise sous cette autre forme

ϖ/8=1/2−( 1/3.5+1/7.9+1/11.13+1/15.17+…);

faisant donc le calcul du nombre ϖ par cette nouvelle série, comme par la première, on n’admettrait, dans sa valeur définitive, que les chiffres décimaux communs aux deux résultats.

IV. Pour troisième application, nous choisirons le problème important et délicat de l’interpolation des suites ; mais ici les formes de l’application de la méthode pouvant être variées d’une multitude de manières différentes ; nous insisterons principalement sur quelques procédés, en nous bornant par rapport aux autres à une briève indication.

On sait que le problème qui nous occupe se réduit à former l’équation d’une courbe parabolique passant par un certain nombre de points dont on connaît les coordonnées, ou du moins s’écartant le moins possible de ces points, que l’on peut supposer n’être qu’à peu près sur la courbe qu’on cherche. Supposons donc, en premier lieu, pour suivre exactement l’esprit du procédé de M. Kramp, que l’on ait sept ordonnées équidistantes α, β, γ, δ, ε, ζ, η; on pourra chetcher successivement l’expression générale de l’ordonnée de la parabole 1.ᵒ du sixième degré, passant par les extrémités de ces ordonnées ; 2.ᵒ du troisième degré, passant par les extrémités des mêmes ordonnées prises de deux en deux seulement ; 3.ᵒ du deuxième degré, passant par leurs extrémités, de trois en trois ; 4.ᵒ enfin du premier degré, passant uniquement par les extrémités des deux ordonnées extrêmes.

Désignant alors respectivement ces expressions par A, B, C, D et les considérant comme répondant aux indices respectifs 1/6, 1/3, 1/2, 1, ou, ce qui revient au même, aux indices 1, 2, 3, 6; le terme qui répondra à l’indice zéro équivaudra sensiblement à la valeur générale de l’ordonnée qui répondrait au cas où on aurait fait entrer en considération une infinité d’ordonnées intermédiaires entre les ordonnées extrêmes ; et sera conséquemment une expression plus exacte de l’ordonnée générale que celle qu’avait fourni la considération des sept points donnés.

L’application de ce procédé exige que le nombre des points donnés diminué d’une unité ; ait le plus grand nombre possible de diviseurs. Voici un autre procédé qui n’est point sujet à cette limitation.

Lorsque les ordonnées données sont équidistantes et en nombre impair, on peut toujours, pour plus de simplicité, prendre leur commune distance pour unité, et supposer en outre que l’ordonnée moyenne répond à l’origine. Supposons donc qu’on ait les cinq ordonnées consécutives βₗₗ, βₗ, β, β′, β″ répondant respectivement aux abscisses −2, −1, ∓0, +1, +2 ; et proposons-nous de trouver l’ordonnée y qui doit répondre à l’abscisse quelconque x.

1.ᵒ Ne considérons d’abord que l’ordonnée βₗₗ, et posons

E=βₗₗ.

2.ᵒ Considérons en second lieu les deux ordonnées βₗₗ, βₗ; en désignant par D l’ordonnée générale de la droite qui joint leurs extrémités supérieures, nous aurons

D=(2βₗ−βₗₗ)+(βₗ−βₗₗ)x.

3.ᵒ Considérons ensuite les trois ordonnées βₗₗ, βₗ; en désignant par C l’ordonnée générale de la parabole ordinaire qui joint leurs extrémités supérieures, nous aurons

C=β+1/2 (3β−4βₗ+βₗₗ)x+1/2 (β−2βₗ+βₗₗ)x².

4.ᵒ Appelant de même B l’ordonnée générale de la parabole du troisième degré qui joint les extrémités supérieures des quatre ordonnées βₗₗ, βₗ, β, β′, nons aurons

B=β+1/6 (2β′+3β−6βₗ+βₗₗ)x+1/6 (β′−2β+βₗ)x²

+1/6 (β′−3β+3βₗ−βₗₗ)x³.

5.ᵒ Appelant enfin A l’ordonnée générale de la parabole du quatrième degré qui résulte de l’emploi total de cinq données βₗₗ, βₗ, β, β′, β″ on aura

A=β−1/12 (β″−8β′+8βₗ−βₗₗ)x−1/24 (β″−8β′+14β−8βₗ+βₗₗ)x²

+1/2 (β″−2β′+2βₗ−βₗₗ)x³+1/24 (β″−2β′+2β−2βₗ+βₗₗ)x⁴.

Si présentement on considère les valeurs successives A, B, C, D, E comme répondant successivement aux indices 1/5, 1/4, 1/3, 1/2, 1, ou, ce qui revient au même, aux indices 60 fois plus grands 12, 15, 20, 30, 60, le terme répondant à l’indice 0 sera sensibletnent ce qu’on obtiendrait pour y, en ayant égard à une infinité d’autres ordonnées qui seraient censées suivre β″ suivant la loi qui régit les premières ; on trouve dans ce cas, comme nous l’avons déjà vu,

y=(1953125A−2097152B+531441C−24576D+42E)/7.8.8.9.9.10 ;

et il ne restera plus que les substitutions à exécuter.

On pourra ensuite procéder d’une manière inverse, c’est-à-dire, prendre successivement une, deux, trois, quatre et cinq ordonnées en allant de β″ vers βₗₗ; on obtiendra ainsi une nouvelle expression de y, qui ne différera au surplus de la précédente qu’en ce que ϰ″ et β′ y seront respectivement changés en βₗₗ et βₗ, et réciproquement. Cette nouvelle sera relative à l’hypothèse où l’on aurait eu égard à une infinité d’ordonnées précédant βₗₗ. La demi-somme de ces expressions donnera l’expression la plus convenable à employer. Leur différence qui sera nécessairement très-petite fera connaître sensiblement l’erreur dans laquelle l’emploi de chacune d’elles peut entraîner.

Mais de toutes les manières d’appliquer la nouvelle méthode à l’interpolation des suites la plus exacte paraît devoir être la suivante. Soient n+1 le nombre des valeurs données et correspondantes de x et de y. Soient A, B, C, D, … les fonctions des degrés n, n−1, n−2, n−3… représentant le plus exactement possible les valeurs données ; ces fonctions étant obtenues par la méthode des moindres quarrés, ainsi qu’il a été expliqué dans ce volume (pag. 242 et suiv.). On considérera A, B, C, D, … comme répondant respectivement aux indices 1/n, 1/n−1, 1/n−2, 1/n−3, …;

et cherchant, comme ci-dessus, le terme qui doit répondre à l’indice zéro, on prendra ce terme pour la valeur de γ.

Nous n’entrerons point actuellement dans plus de détails à ce sujet, sur lequel nous pourrons peut-être revenir une autre fois. Il nous suffit pour le présent d’avoir montré que l’analise possède, dans la méthode développée par M. Kramp, un nouvel instrument, susceptible sans doute de perfectionnement ; mais qui, tel qu’il est, peut déjà, dans un grand nombre de circonstances, devenir d’un usage très-précieux.

