Efficient automatic integration of ordinary differential equations with discontinuities

This paper describes a method of automatically detecting and accurately locating discontinuities which occur in many applications of ordinary differential equations. The integration formula is a Runge-Kutta so chosen that accurate values between integration points can be found by Hermite interpolation.