Numerical Method in Classical Dynamics

Abstract

A set of algorithms is presented for efficient numerical calculation of the time evolution of classical dynamical systems. Starting with a first approximation for solving the differential equations that has a "reversible" character, we show how to bootstrap easily to higher order accuracy.The method, first shown for a single particle in one dimension, is then neatly extended to many dimensions and many particles.