40. Nous n’avons considéré, dans la première approximation des mouvements des corps célestes, que les forces principales qui les animent, et nous en avons déduit les lois du mouvement elliptique. Nous aurons égard, dans les recherches suivantes, aux forces perturbatrices de ce mouvement. L’action de ces forces ne fait qu’ajouter des petits termes aux équations différentielles du mouvement elliptique, dont nous avons donné précédemment les intégrales finies : il faut maintenant déterminer par des approximations successives les intégrales des mêmes équations augmentées des termes dus à l’action des forces perturbatrices. Voici, pour cet objet, une méthode générale, quels que soient le nombre et le degré des équations différentielles dont on se propose de trouver des intégrales de plus en plus approchées.
Supposons que l’on ait, entre les variables et la variable dont l’élément est regardé comme constant, les équations différentielles
étant des fonctions de et de leurs différences jusqu’à l’ordre inclusivement, et étant un très-petit
