void heun(double a, double b, double h, int N, double y0, double heunout[]);