void rk4(double a, double b, double h, int N,  double y0, double rkout[]);