The equations of motion for N-body systems are usually integrated by means of a central-difference algorithm. An alternative average force algorithm is described in this paper. The new method is shown to have both theoretical and practical advantages when compared to the central-difference method.