void midpoint_explicit ( double *dydt ( double x, double y[] ), 
  double tspan[2], double y0[], int n, int m, double t[], double y[] );
double *r8vec_linspace_new ( int n, double a, double b );