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