  The method of choice for integrating the equations of motion of the general
N-body problem has been to use an individual time step scheme. For the sake of
efficiency, block time steps have been the most popular, where all time step
sizes are smaller than a maximum time step size by an integer power of two. We
present the first successful attempt to construct a time-symmetric integration
scheme, based on block time steps. We demonstrate how our scheme shows a vastly
better long-time behavior of energy errors, in the form of a random walk rather
than a linear drift. Increasing the number of particles makes the improvement
even more pronounced.
