void midpoint ( void dydt ( double t, double y[], double dy[] ), 
  double tspan[2], double y0[], int n, int m, double t[], double y[] );