SUMMARYIn a companion paper, Thomas and Gladwell: the authors designed a class of multistage methods for secondorder systems of ordinary differential equations (ODES), together with local error estimators based on embedding techniques. Here we provide a set of subroutines implementing the algorithms and discuss our numerical experience with the resulting codes.