Shampine and Gordon ODE solver (by John Burkardt)