double rk4 ( double t0, double u0, double dt, double f ( double t, double u ) );