Abstract
The numerical integration of systems of differential equations that possess integrals is often approached by using the integrals to reduce the number of degrees of freedom or by using the integrals as a partial check on the resulting solution, retaining the original number of degrees of freedom.
Another use of the integrals is presented here. If the integrals have not been used to reduce the system, the solution of a numerical integration may be constrained to remain on the integral surfaces by a method that applies corrections to the solution at each integration step. The corrections are determined by using linearized forms of the integrals in a least-squares procedure.
The results of an application of the method to numerical integrations of a gravitational system of 25-bodies are given. It is shown that by using the method to satisfy exactly the integrals of energy, angular momentum, and center of mass, a solution is obtained that is more accurate while using less time of calculation than if the integrals are not satisfied exactly. The relative accuracy is ascertained by forward and backward integrations of both the corrected and uncorrected solutions and by comparison with more accurate integrations using reduced step-sizes.
Similar content being viewed by others
References
Aarseth, S. J.: 1966,Monthly Notices Roy. Astron. Soc. 132, 35.
Fehlberg, E.: 1966,Z. Angew. Math. Mech. 46, 1.
Forsyth, A. R.: 1930,Geometry of Four Dimensions,I, Cambridge, p. 80.
Gottlieb, R.: 1970, private communication.
Lanczos, C.: 1949,The Variational Principles of Mechanics, Toronto, p. 43.
Lecar, M.: 1968,Bull. Astron. 3, 91.
Miller, R. H.: 1964,Astrophys. J. 140, 250.
Szebehely, V. G.: 1968,Bull. Astron. 3, 33.
Author information
Authors and Affiliations
Rights and permissions
About this article
Cite this article
Nacozy, P.E. The use of integrals in numerical integrations of theN-body problem. Astrophys Space Sci 14, 40–51 (1971). https://doi.org/10.1007/BF00649193
Issue Date:
DOI: https://doi.org/10.1007/BF00649193