Composition methods in the presence of small parameters

We derive numerical methods for arbitrary small perturbations of exactly solvable differential equations. The methods, based in one instance on Gaussian quadrature, are symplectic if the system is Hamiltonian and are asymptotically more accurate than previously known methods.