An efficient integrator that uses Gauss-Radau spacings

This describes our integrator RADAU, which has been used by several groups in the U.S.A., in Italy, and in the U.S.S.R, over the past 10 years in the numerical integration of orbits and other problems involving numerical solution of systems of ordinary differential equations. First- and second-order equations are solved directly, including the general second-order case. A self-starting integrator, RADAU proceeds by sequences within which the substeps are taken at Gauss-Radau spacings. This allows rather high orders of accuracy with relatively few function evaluations. After the first sequence the information from previous sequences is used to improve the accuracy. The integrator itself chooses the next sequence size. When a 64-bit double word is available in double precision, a 15th-order version is often appropriate, and the FORTRAN code for this case is included here. RADAU is at least comparable with the best of other integrators in speed and accuracy, and it is often superior, particularly at high accuracies.

[1]  T. Levi-Civita.  Traiettorie singolari ed urti nel problema ristretto dei tre corpi , 1904 .

[2]  G. M. Clemence,et al.  Coordinates of the five outer planets, 1653-2060 , 1951 .

[3]  R. Newton.  Periodic Orbits of a Planetoid Passing Close to Two Gravitating Masses , 1959 .

[4]  John C. Butcher,et al.  Integration processes based on Radau quadrature formulas , 1964 .

[5]  A. Stroud,et al.  Gaussian quadrature formulas , 1966 .

[6]  J. Stoer,et al.  Numerical treatment of ordinary differential equations by extrapolation methods , 1966 .

[7]  Victor Szebehely,et al.  Treatment of Close Approaches in the Numerical Integration of the Gravitational Problem of N-Bodies (Papers appear in the Proceedings of IAU Colloquium No. 10 Gravitational N-Body Problem (ed. by Myron Lecar), R. Reidel Publ. Co. , Dordrecht-Holland.) , 1971 .

[8]  E. Fehlberg,et al.  Classical eight- and lower-order Runge-Kutta-Nystroem formulas with stepsize control for special second-order differential equations , 1972 .

[9]  M. Lecar,et al.  Gravitational N-Body Problem , 1972 .

[10]  Fred T. Krogh,et al.  On Testing a Subroutine for the Numerical Integration of Ordinary Differential Equations , 1973, JACM.

[11]  E. Everhart.  Implicit single-sequence methods for integrating orbits , 1974 .

[12]  Z. Sekanina,et al.  New osculating orbits for 110 comets and analysis of original orbits for 200 comets , 1978 .

[13]  G. Weiss,et al.  Numerical integration of stellar orbits , 1978 .

[14]  K. Innanen,et al.  A comparison of five algorithms for numerical orbit computation in galaxy models , 1978 .

[15]  W. Fricke,et al.  Sun and Planetary System , 1982 .

[16]  Methods of Computation of the Perturbed Motion of Small Bodies in the Solar System , 1982 .

[17]  B. Marsden,et al.  New original and future cometary orbits , 1983 .

[18]  Giovanni B. Valsecchi,et al.  The Long-Term Evolution Project , 1985 .

[19]  G. Valsecchi,et al.  Dynamics of comets: Their origin and evolution , 1985 .