These IDL routines perform a Runga Kutta integration of a system of first order differential equations.