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

## ANALISE TRANSCENDANTE.

### Intégration par approximation de toute équation différentielle quelconque ;

I. Le mémoire actuel sera destiné à intégrer, par approximation, l’équation différentielle qui suit :

y+A d⁡y/d⁡t=T;

dans laquelle la lettre A désigne un coefficient quelconque constant, et la lettre T une fonction quelconque de t. L’approximation sera semblable à celle que nous avons employée pour intégrer la différentielle Xd⁡x, dans divers mémoires déjà publiés dans ce recueil. Ainsi elle doit, dans les cas ordinaires, savoir, dans ceux qui sont sans asymptotes, sans points d’inflexion ni de rebroussement, faire connaître, dès le premier essai, l’intégrale demandée, jusqu’à cinq, et dès le second jusqu’à dix ou douze décimales. Il sera facile ensuite d’appliquer la méthode à des cas plus compliqués. Ainsi ; l’équation plus générale

y+A d⁡y/d⁡t+B d²⁡y/d⁡t²+C d³⁡y/d⁡t³+…=T,

dans laquelle A, B, C, … sont des constantes, rentrera dans celle de l’équation qui va nous occuper, et s’exécutera par des moyens analogues.

2. Il est bon de remarquer que l’équation

y+A d⁡y/d⁡t=T,

se réduit presque d’elle-même à

y+d⁡y/d⁡x=Q,

en posant simplement, t=Ax. Si, au contraire, la proposée est

y−A d⁡y/d⁡t=T,

en posant également t=Ax, elle deviendra

y−d⁡y/d⁡x=Q;

Q étant, dans l’une et l’autre, une fonction connue de x. Nous nous occuperons donc uniquement, dans tout ce qui va suivre, des deux équations

y+d⁡y/d⁡x=Q, y−d⁡y/d⁡x=Q;

ce qui introduira dans nos calculs des simplifications notables.

3. La valeur rigoureuse de y est

y=e⁻ˣ∫eˣQd⁡x=e⁻ˣ(C+X)=Ce⁻ˣ+Xe⁻ˣ;

en désignant par C une constante arbitraire, et par X une certaine fonction de x, que le calcul nous fera connaître, et dont la détermination est précisément l’objet principal qui doit nous occuper. L’autre partie Ce⁻ˣ, appartient généralement à toutes les équations différentielles de cette classe ; de sorte qu’après avoir trouvé l’intégrale particulière Xe⁻ˣ, il ne s’agira, pour la rendre complète, que de lui ajouter le terme Ae⁻ˣ,. La méthode que nous allons enseigner suppose qu’avant tout on soit instruit des limites entre lesquelles l’intégrale doit être prise. En supposant qu’il faille la prendre entre x=a et x=b, il faudra partager l’intervalle entier b−a en un certain nombre de parties égales. Le nombre en est arbitraire ; mais, plus il sera grand, et plus on approchera de l’intégrale demandée. Il faut cependant bien se garder de croire qu’en prenant les nombres arbitraires en progression arithmétique, et en supposant, par exemple, x successivement égal 0, 1, 2, 3, 4, 5, 6, … la suite des erreurs qui en résultent soit constamment décroissante : la courbe qui aurait ces erreurs pour ordonnées irait en serpentant des deux côtés de l’axe, mais en se rapprochant toujours de cet axe, avec lequel elle coïnciderait à l’infini ; ce qui est conforme à la nature de la chose, et ce que j’ai bien directement prouvé d’ailleurs, par la table des erreurs de ϖ/4 (Annales, tom. VI, pag. 379, 887). Il faut savoir de plus qu’en excluant complètement les asymptotes, les points d’inflexion et les points de rebroussement, les résultats commencent à devenir incertains, lorsque les limites sont trop voisines des points où ces circonstances se rencontrent ; c’est ainsi, par exemple, qu’en appliquant les formules à la logistique ; et en employant les coordonnées ordinaires, on trouve des erreurs assez considérables, mais qui disparaissent pourtant en prenant d’autres coordonnées, plus éloignées de la direction de l’asymptote. Les courbes qui s’intègrent le mieux par cette méthode, sont les courbes rentrantes sur elles-mêmes ; ce sont les plus employées dans la pratique, et ce sont en même temps celles auxquelles la nouvelle méthode est le plus applicable. Il est possible, au reste, que le mémoire actuel paraisse inutile à bien des personnes ; attendu que, puisque nous avons ici

y=e⁻ˣ∫eˣQd⁡x,

il s’ensuit que l’intégration à effectuer n’est autre que celle de la formule générale ∫Xd⁡x, dans laquelle X a pris la forme particulière eˣQ; et qu’ainsi tout se réduit, pour avoir y, à diviser par eˣ cette même intégrale que nous avons déjà enseignée à déterminer dans nos précédens mémoires. Mais d’abord le produit eˣQ diffère considérablement de la simple fonction Q, et doit, par suite, introduire une différence notable dans l’intégrale. En outre, l’intégration de eˣQd⁡x est un passage nécessaire pour parvenir à l’intégration de l’équation

y+P d⁡y/d⁡x=Q,

ainsi qu’à celles d’autres équations d’une forme plus compliquée.

4. Commençons par supposer le nombre arbitraire égal à cinq unités : ce qui donne

y=A+Bx+Cx²+Dx³+Ex⁴+Fx⁵,

d’où

d⁡y/d⁡x=B+2Cx+3Dx²+4Ex³+5Fx⁴;

et par conséquent

Q=y+d⁡y/d⁡x=(A+B)+(B+2C)x+(C+3D)x²

+(D+4E)x³+(E+5F)x⁴+Fx⁵.

En supposant successivement ici à la variable x les valeurs entières et positives 0, 1, 2, 3, 4, 5; et en désignant par q₀, q₁, q₂, q₃, q₄, q₅, celles qui en résultent pour Q, nous aurons

q₀=A+B, q₁=A+2B+ 3C+4D+ 5E+ 6F, q₂=A+3B+ 8C+ 20D+48E+ 112F, q₃=A+4B+15C+ 54D+ 189E+ 648F, q₄=A+5B+24C+112D+ 512E+2304F, q₅=A+6B+35C+200D+1125E+6250F.

En ôtant chacune de ces équations de celle qui la suit immédiatement et dénotant par Δq₀, Δq₁, Δq₂, Δq₃, Δq₄, les différences qui en résultent, on aura

Δq₀=B+ 3C+ 4D+5E+ 6F, Δq₁=B+ 5C+16D+ 43E+ 106F, Δq₂=B+ 7C+34D+141E+ 536F, Δq₃=B+ 9C+58D+323E+1656F, Δq₄=B+11C+88D+613E+3946F,

Dénotant de même par Δ²q₀, Δ²q₁, Δ²q₂, Δ²q₃, les différence. consécutives de celles-ci, divisées par deux, il viendra

Δ²q₀=C+ 6D+ 19E+50F, Δ²q₁=C+ 9D+ 49E+ 215F, Δ²q₂=C+12D+ 91E+ 560F, Δ²q₃=C+15D+145E+1145F.

Dénotant en outre par Δ³q₀, Δ³q₁, Δ³q₂, les différences consécutives de celles-ci, divisées par trois, nous aurons

Δ³q₀=D+10E+ 55F, Δ³q₁=D+14E+115F, Δ³q₂=D+18E+195F.

Dénotant encore par Δ⁴q₀, Δ⁴q₁, les différences consécutives de ces dernières, divisées par quatre, il viendra

Δ⁴q₀=E+15F, Δ⁴q₁=E+20F,

Dénotant enfin par Δ⁵q₀ la différence de ces deux-ci, divisées par cinq, on aura

Δ⁵q₀=F.

5. En prenant seulement la première équation de chaque série, et supprimant les indices, désormais inutiles, nous aurons, pour le diviseur cinq,

Δ⁵q=F, Δ⁴q=E+15F, Δ³q=D+10E+55F, Δ²q=C+ 6D+19E+50F, Δq=B+ 3C+ 4D+ 5E+6F, q=A+B.

Si, au lieu du nombre cinq, nous eussions pris tout autre nombre, le nombre douze, par exemple, nous aurions trouvé, par un semblable calcul,

Δ¹²q=N, Δ¹¹q=M+78N, Δ¹⁰q=L+66M+2365N, Δ⁹q=K+55L+1650M+36135N, Δ⁸q=I+45K+1110L+20130M+301587N, Δ⁷q=H+36I+714K+10500L+12866M+1595240N, Δ⁶q=G+28H+434I+5040K+49287L+430584M+3477496N, Δ⁵q=F+21G+245H+2170I+16401K+112035L+764505M+4340160N, Δ⁴q=E+15F+125G+805H+4501I+25079K+111805L+520905M+2360501N, Δ³q=D+10E+55F+240G+931H+3374I+11719K+39580L+131131M+428538N, Δ²q=C+6D+19E+50F+121G+280H+631I+1398K+3061L+6644M+24323N, Δ¹q=B+3C+4D+5E+6F+7G+8H+9I+10K+11L+12M+13N, q=A+B.

Cette table n’est pas seulement applicable au nombre douze ; elle l’est également à tous les diviseurs inférieurs, en s’arrètant dans chacune des formules à la lettre qui marque la limite de la division ; à la lettre G, par exemple, si le diviseur est six ; à la lettre I, si ce diviseur est huit ; et ainsi des autres.

6. Il reste donc à déterminer, à l’aide de ces équations, les valeurs de N, M, L, K, I, …A, pour les substituer dans celle de y, afin de présenter cette intégrale sous la forme d’une série disposée suivant les différences des différens ordres de la quantité q. Comme le calcul est très-facile, il suffira d’en offrir ici les résultats. Ces résultats sont

N=Δ¹²q, M=Δ¹¹q−78Δ¹²q, L=Δ¹⁰q−66Δ¹¹q+2783Δ¹²q, K=Δ⁹ q−55Δ¹⁰q+1980Δ¹¹q−60500Δ¹²q, I=Δ⁸ q−45Δ⁹q+1365Δ¹⁰q−35970Δ¹¹q+901923Δ¹²q, H=Δ⁷ q−36Δ⁸q+906Δ⁹q−20370Δ¹⁰q+445533Δ¹¹q−9852942Δ¹²q, G=Δ⁶ q−28Δ⁷q+574Δ⁸q−10878Δ⁹q+205863Δ¹⁰q−4020786Δ¹¹q+82310129Δ¹²q, F=Δ⁵ q−21Δ⁶q+343Δ⁷q−5404Δ⁸q+87717Δ⁹q−150450Δ¹⁰q+27541646Δ¹¹q−539856504Δ¹²q, E=Δ⁴ q−15Δ⁵q+190Δ⁶q−2450Δ⁷q+33789Δ⁸q−505869Δ⁹q+8246195Δ¹⁰q−146117730Δ¹¹q+2804540596Δ¹²q,

D=Δ³q−10Δ⁴q+95Δ⁵q−985Δ⁶q+11424Δ⁷q−148288Δ⁸q+2141600Δ⁹q−34157480Δ¹⁰q+597224496Δ¹¹q−11369080360Δ¹²q, C=Δ²q−6Δ³q+41Δ⁴q−335Δ⁵q+3229Δ⁶q−36036Δ⁷q+457932Δ⁸q−6534384Δ⁹q+103499016Δ¹⁰q−1802302128Δ¹¹q+34227784920Δ¹²q, B=Δq−3Δ²q+14Δ³q−88Δ⁴q+694Δ⁵q−6578Δ⁶q+72792Δ⁷q−920904Δ⁸q+13109088Δ⁹q−207360912Δ¹⁰q+3608233056Δ¹¹q−68495486640Δ¹²q, A=q−Δq+3Δ²q−14Δ³q+88Δ⁴q−694Δ⁵q+6578Δ⁶q−72792Δ⁷q+920904Δ⁸q−13109088Δ⁹q+207360912Δ¹⁰q−3608233056Δ¹¹q+68495486640Δ¹²q.

7. La loi que suivent les coefficiens de Δq, Δ²q, Δ³q, de la table précédente est extrêmement remarquable ; et ils ont avec ceux des facultés numériques, dont j’ai traité dans mon Analise des réfractions (pag. 71) une analogie \operatorname{sin}gulière, dont je n’ai pu encore me rendre compte. Voici en quoi cette analogie consiste. Supposons, par exemple, que l’on demande les coefficiens de Δ⁶q dans les valeurs de A, B, C, D, E, F, G? On rassemblera les coefficiens de la faculté à exposant six, que l’on trouvera être 1, 15, 85, 225, 274, 120; on les multipliera respectivement par lesfacultés 6!, 5!, 4!, 3!, 2!, 1! ou 720, 120, 24, 6, 2, 1; ce qui donnera les produits

720, 1800, 2040, 1350, 548s, 120;

prenant successivement le premier seul, puis la somme deux premiers, puis celle des trois premiers, et ainsi de suite, jusqu’au dernier, on formera la nouvelle suite

720, 2520, 4560, 5910, 6458, 6578,

dont les termes, divisés respectivement par les mêmes facultés 6!, 5!, 4!, 3!, 2!, 1!, donneront les quotiens

1, 21, 190, 985, 3229, 6078,

qui sont précisément les coefficiens de Δ⁶q.

Supposons encore que l’on demande les coeffiens de Δ¹²q ? Il faudra d’abord écrire ceux de la faculté à exposant 12 ; ce sont

1…1, 2…66, 3…1925, 4…32670, 5…357423, 6…2637558, 7…13339535, 8…45995730, 9…105258076, 10…150917976, 11…120543840, 12…39916800;

on les multipliera respectivement par les facultés

12!=479001600, 11!=39916800, 10!=3628800, 9!=362880, 8!=40320, 7!=5040, 6!=720, 5!=120, 4!=24, 3!=6, 2!=2, 1!=1;

ce qui donnera les produits

1.ᵉʳ 479001600, 2…2630508800, 3…6980440000, 4…11855289600, 5…14411295360, 6…13293292320, 7…9604465200, 8…5519487600, 9…252619824, 10…905807856, 11…242087680, 12…39916800,

prenant ensuite le premier, puis la somme des deux premiers puis celle des trois premiers, et ainsi de suite, il viendra

1.ᵉʳ 479001600, 2…3113510400, 3…10098950400, 4…21954240000, 5…36365535360, 6…49658827680, 7…59263292880, 8…64782780480, 9…67308974304, 10…68214482160, 11…68455569840, 12…68495486640;

divisant enfin ces sommes par les mêmes facultés qui avaient d’abord été employées comme multiplicateurs, on obtiendra, pour la série des coefficiens de Δ¹²q, dans nos formules,

N, …1, M, …78, L, …2783, K, …60570, I, …9019239, H, …9852942, G, …82310129,

F, …539856504, E, …2804540596, D, …11369080360, C, …34227784920, B, …68495486640.

8. Un second théorème, qui n’est pas moins digne de remarque, quoiqu’il n’ait point encore pour lui une démonstration rigoureuse, mais qui est fondé sur une induction plus que suffisante, et duquel d’ailleurs je me propose de m’occuper encore, c’est que toutes ces séries de coefficiens qui multiplient les différences d’un même ordre dans nos formules (6) jouissent sensiblement de la propriété de se reproduire eux-mêmes en les divisant respectivement par les nombres naturels 1, 2, 3, 4, …; et se rapprochent en cela des termes de la série hypergéométrique ordinaire 1, 1/1, 1/2, 1/6, 1/24, 1/120, 1/720, … qui multiplient respectivement les termes 1, x, x², x³, x⁴, … dans le développement de eˣ, et qui jouissent rigoureusement de cette propriété.

Prenons, par exemple, du plus grand au plus petit, les coefficiens de Δ⁷q, et divisons-les respectivement par 1, 2, 3, 4, …; nous aurons

1/1.72792=72793, 1/2.72792=36396, 1/3.36036=12012, 1/4.11424=2856, 1/5. 2450=490, …………

qui sont à peu près ces mêmes coefficiens, le premier excepté.

En opérant de la même manière sur les coefficiens de Δ¹⁰q nous trouverons

1/1.207360912=207360912, 1/2.207360912=103680456, 1/3.103499016=34499672, 1/4. 34107480=8539370, 1/5. 8246193=1649239, …………

qui sont à peu près ces mêmes coefficiens.

Enfin, en opérant encore ainsi sur les coefficiens de Δ¹²q, nous aurons

1/1.68495486640=68495486640, 1/2.68495486640=34247743320, 1/3.34227784920=11409261640, 1/4.11367080360=2842270090, 1/5. 2804540596=560908119, 1/6. 539806504=89976084, …………;

où la même loi se manifeste également. La démonstration rigoureuse de ce théorème serait sans doute difficile ; mais en attendant, nous l’adopterons, avec d’autant plus de fondement qu’il nous conduira à des résultats exacts et décisifs.

9. L’intégration complète de l’équation y+d⁡y/d⁡x=Q, équation dans laquelle la lettre Q désigne une fonction quelconque de x, suppose deux choses : d’abord la nouvelle fonction de x dont la différentiation nous ramènera à l’équation proposée, et ensuite la constante, multipliée par une certaine autre fonction de x. Nous avons vu que ce dernier produit restait le même quelle que pût être la fonction Q; en conséquence, pour le déterminer, il n’y a qu’à voir ce qu’il deviendra dans la supposition la plus simple qu’on puisse adopter pour Q, qui est celle de Q=0. L’équation sera alors

y+d⁡y/d⁡x=0;

puis donc que nous avons suppose

y=A+Bx+Cx²+Dx³+Ex⁴+Fx⁵+…,

d’où

d⁡y/d⁡x=Bx+2Cx+3Dx²+4Ex³+5Fx⁴+…;

Il faudra résoudre l’équation

0=(A+B)+(B+2C)x+(C+3D)x²+(D+4E)x³+…

Cela conduira aux équations déjà trouvées (6) avec cette seule différence qu’ici q, Δq, Δ²q, Δ³q, … seront zéro.

Il en faudra seulement rejeter la première équation N=Δ¹²q, qui devenant, dans le cas actuel, N=0, donnerait zéro pour valeur de tous les autres coefficiens, tandis qu’il faut nécessairement laisser du jeu à la constante qu’on se propose d’ajouter. Cette légère attention nous met dans la position d’avoir cette constante, qu’il eût été bien difficile de trouver d’une autre manière quelconque.

11. Avec cette attention, les équations trouvées (6) nous donneront

M=−78N, L=+2783N, K=−60500N, I=+901923N, H=−9852942N, G=+82310129N, F=−539856504N, E=+2804540596N, D=−11369080360N, C=+34227784920N, B=−68495486640N, A=+68495486640N.

Il ne restera plus qu’à diviser la dernière de ces équations par chacune de celles qui la précèdent, pour avoir tous ces coefficiens l’un après l’autre. Le premier terme A sera arbitraire ; il formera la constante du problème ; et l’on aura pour les autres

68495486640B=−68495486640A, 68495486640C=+34227784920A, 68495486640D=−11369080360A, 68495486640E=+2804540596A, 68495486640F=−539856504A, 68495486640G=+82310129A, 68495486640H=−9852942A,

68495486640I = +901923A, 68495486640K=−60500A, 68495486640L =+2783A, 68495486640M=−78A, 68495486640N =+A.

En conséquence du théorème énoncé ci-dessus (8), on voit fort bien ce que ces rapports compliqués deviendraient dans l’infini ; on aurait alors

B=−A, 2C=+A, 6D=−A, 24E=+A, 120F=−A, ……

ce qui donnerait

y=A(1−x/1+x²/2−x³/6+x⁴/24−x⁵/120+…)=Ae⁻ˣ,

conformément aux vrais principes du calcul intégral. On voit en même temps l’identité absolue entre le coefficient constant A et le premier terme de la série, qui répond à x=0.

12. Reste donc à trouver l’autre fonction de x, dont la différentiation nous conduit proprement à l’équation proposée

y+d⁡y/d⁡x=Q,

Comme on a supposé

y=A+Bx+Cx²+Dx³+Ex⁴+Fx⁵+…,

on a, en mettant à la place de A, B, C, D, les valeurs trouvées (6),

y=q+(x−1)Δq+(x²−3x+3)Δ²q+(x³−6x²+14x−14)Δ³q+(x⁴−10x³+41x²−88x+88)Δ⁴q+(x⁵−15x⁴+95x³−335x²+694x−694)Δ⁵q+(x⁶−21x⁵+190x⁴−985x³+3229x²−6578x+6578)Δ⁶q+………………………………………………

Mettant ici, à la place de x, les valeurs successives 2, 3, 4, … on aura

Pour

x=2, y=q+ Δq+Δ²q, =3, y=q+2Δq+ 3Δ²q+Δ³q, =4, y=q+3Δq+ 7Δ²q+10Δ³q+8Δ⁴q, =5, y=q+4Δq+13Δ²q+31Δ³q+ 48Δ⁴q+ 26Δ⁵q, =6, y=q+5Δq+21Δ²q+70Δ³q+172Δ⁴q+266Δ⁵q+194Δ⁶q,

7, y=q+6Δq+31Δ²q+133Δ³q+ 452Δ⁴q+1126Δ⁵q+1790Δ⁶q+ 1142Δ⁷q, 8, y=q+7Δq+43Δ²q+226Δ³q+ 984Δ⁴q+3386Δ⁵q+8546Δ⁶q+13672Δ⁷q+9736Δ⁸q, 9, y=q+8Δq+57Δ²q+355Δ³q+1888Δ⁴q+8306Δ⁵q+28862Δ⁶q+73494Δ⁷q+119112Δ⁸q,

13. En mettant ici, à la place des différences successives de q leurs valeurs en q₀, q₁, q₂, q₃, … savoir ;

Δq=q₁− q₀, 2Δ²q=q₂−2q₁+q₀, 6Δ²q=q₃−3q₂+ 3q₁−q₀, 24Δ²q=q₄−4q₃+ 6q₂− 4q₁+ q₀, 120Δ²q=q₅−5q₄+10q₃−10q₂+5q₁−q₀,

on trouvera

Pour

x=2…2y=q₀+q₂, 3, …2y= 2q₀− 3q₁+6q₂+q₃, 4, …6y=q₀− 2q₁+3q₂+ 2q₃+ 2q₄, 5, … 60y= 7q₀−25q₁+ 50q₂−40q₃+65q₄+ 13q₅, 6, …360y=19q₀−72q₁+135q₂−80q₃+45q₄+216q₅+97q₆+…

7… 2520y=108q₀− 623q₁+1764q₂− 2835q₃+ 3220q₄−1953q₅+ 2268q₆+ 5719q₇, 8…20160y=300q₀−1536q₁+3416q₂− 2688q₃− 1680q₄+8960q₅−72249q₆+15744q₇+48689q₈, 9…15120y=260q₀−2115q₁+8218q₂−19278q₃+30744q₄−34030q₅+28560q₆−14778q₇+14148q₈+3391q₉,

Dans le cas parliculier où les quantités q₀, q₁, q₂, q₃, … seraient égales entre elles, on aurait, dans toutes ces équations, y=q; ce qui pourra servir au besoin à vérifier l’exactitude de nos formules.

14. Exemple I. Soit l’équation

y+d⁡y/d⁡x=a,

On aura, dans ce cas, Q=a; ainsi q₀, q₁, q₂, q₃, … sont ici tous égaux à a; en conséquense, toutes les formules donnent pour intégrale complète

y=Ae⁻ˣ+a.

15. Exemple II. Soit l’équation

y=d⁡y/d⁡x=x.

On aura Q=x; donc q₀=0, q₁=1, q₂=2, q₃=3, … et toutes les formules s’accordent également à donner

y=Ae⁻ˣ+x−1;

ce qui est rigoureusement conforme aux principes du calcul ; on tire, en effet,

De x=2, 2y=q₀+q₂=2;donc y=1=2−1=x−1;De x=3, 6y=2q₀−3q₁+6q₂+q₃=12;donc y=2=3−1=x−1;De x=4, 6y=q₀−2q₁+3q₂+2q₃+2q₄=18;donc y=3=4−1=x−1;De x=5, 60y=7q₀−25q₁+50q₂−40q₃+55q₄+13q₅=240;donc y=4=5−1=x−1.

il en sera de même des suivans, de sorte qu’on aura généralement et rigoureusement

y=Ae⁻ˣ+x−1.

16. Exemple III. Soit l’équation

y+d⁡y/d⁡x=x².

On aura Q=x², donc

q₀=0, q₁=1, q₂=4, q₃=9, q₄=16, …

cela donne

Pour x=2, 2y= 4, donc y=2=x²−2x+2, x=3, 6y=30, 5=x²−2x+2, x=4, 6y=60, 10=x²−2x+2, x=5, 60y=1020, 17=x²−2x+2, x=6, 360y=9360, 26=x²−2x+2, ……………………

on aura donc généralement

y=Ae⁻ˣ+x²−2x+2;

valeur qui en effet est rigoureuse,

17. Exemple IV. Soit l’équation

y+d⁡y/d⁡x=x³.

On aura ici Q=x³, d’où

q₀=0, q₁=1, q₂=8, q₃=27, q₄=64, …

donc à commencer par la valeur 3

Pour x=3, y=12, 4, 34, 5, 74, 6, 138, 7, 362, ………

Ces nombres étant tous compris sous la formule x³−3x²+6x−6; on aura généralement

y=Ae⁻ˣ+x³−3x²+6x−6,

intégrale qui, en effet, est rigoureusement exacte.

18. Exemple V. Soit l’équation

y+d⁡y/d⁡x=x⁴.

Ayant ici Q=x⁴, on aura

q₀=0, q₁=1, q₂=16, q₃=81, q₄=256, …

donc, en commençant par la valeur 4,

Pour x=4, y=120, 5, 329, 6, 744,

valeurs comprises dans la formule x⁴−4x³+12x²−24x+24, en sorte qu’on aura

y=Ae⁻ˣ+x⁴−4x³+12x²−24x+24;

ce qui est rigoureusement exact.

19. Soit plus généralement l’équation

y+d⁡y/d⁡x=a+bx+cx²+dx³+ex⁴;

on aura

Q=a+bx+cx²+dx³+ex⁴;

d’où

q₀=a, q₁=a+ b+c+ d+ e, q₂=a+2b+ 4c+8d+ 16e, q₃=a+3b+ 9c+ 27d+ 81e, q₄=a+4b+16c+ 64d+256e, q₅=a+5b+25c+125d+625e,

donc, en partant de la valeur 5,

Pour x=5, y=a+4b+17c+ 74d+ 329e, 6, a+5b+26c+138d+ 744e, 7, a+6b+37c+232d+1473e,

résultats qui sont tous compris dans la formule générale

a+b(x−1)+c(x²−2x+2)+d(x³−3x²+6x−6)

+e(x⁴−4x³+12x²−24x+24)

de sorte qu’on doit avoir

y=Ae⁻ˣ+a+b(x−1)+c(x²−2x+2)+d(x³−3x²+6x−6)+e(x⁴−4x³+12x²−24x+24).

20. Les résultats obtenus dans ce mémoire ont tous été exacts et rigoureux ; et ils ont dû l’être à raison de ce que les exposans de toutes les puissances dont se composait la quantité Q étaient entiers et positifs. Dans le mémoire qui suivra celui-ci, nous prendrons pour Q des fonctions quelconques de x; nous leur appliquerons la même méthode ; nous serons conduits à des résultats absolument neufs, et nous aurons lieu d’être satisfaits de leur exactitude ; conséquence nécessaire de l’approximation que nous avons employée. La seule difficulté qui reste sera la détermination de la constante. Dans le cas que nous venons d’exposer, savoir y+d⁡y/d⁡x=Q, cette constante était Ae⁻ˣ; encore ne sommes-nous parvenus à ceci que par une induction très permise : mais qui eut été difficile dans d’autres cas quelconques.

