A Mapping for the Gravitational Few-Body Problem with Dissipation

Recently a new class of numerical integration methods – ``mixed variable symplectic integrators’’ – has been introduced for studying long-term evolution in the conservative gravitational few-body problem. These integrators are an order of magnitude faster than conventional ODE integration methods. Here I present a simple modification of this method to include small non-gravitational forces. The new scheme provides a similar advantage of computational speed for a larger class of problems in Solar System dynamics.