  We present TES, a new n-body integration code for the accurate and rapid
propagation of planetary systems in the presence of close encounters. TES
builds upon the classic Encke method and integrates only the perturbations to
Keplerian trajectories to reduce both the error and runtime of simulations.
Variable step size is used throughout to enable close encounters to be
precisely handled. A suite of numerical improvements are presented that
together makes TES optimal in terms of energy error. Lower runtimes are found
in all test problems considered when compared to direct integration using
ias15. TES is freely available.
