| The integration of Newton's equations of motion for self-gravitating systems, particularly in the context of our Solar System's evolution, remains a paradigm for complex dynamics. We implement Stormer's multistep method in backward difference, summed form and perform arithmetic according to what we call ``significance ordered computation''. We achieve results where the truncation error of our thirteenth order integrator resides below machine (double) precision and roundoff error accumulation is strictly random and not systematic. |
Keywords
celestial mechanics, stormer, optimal error growth
Math Review Classification
Primary 65L06
; Secondary 70F05
Last Updated
October 29
Length
13
Availability
This article is available in: