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