ABSTRACT We present heyoka, a new, modern and general-purpose implementation of Taylor’s integration method for the numerical solution ordinary differential equations. Detailed tests focused on difficult high-precision gravitational problems in astrodynamics celestial mechanics show how our integrator is competitive with often superior to state-of-the-art specialized symplectic non-symplectic i...