BDF

D02EJF	ODEs, stiff IVP, BDF method, until function of solution is zero, intermediate output (simple driver)
D02NVF	ODEs, IVP, BDF method, setup for D02M–N routine s