void rk4 ( void dydt ( double t, double u[], double f[] ), double tspan[2], 
  double y0[], int n, int m, double t[], double y[] );
void timestamp ( );