Runge–Kutta–Merson

D02BGF	ODEs, IVP, Runge–Kutta–Merson method, until a component attains given value (simple driver)
D02BHF	ODEs, IVP, Runge–Kutta–Merson method, until function of solution is zero (simple driver)