//
//Runge-Kutta 2th order functiom to solve ODE                                                             
//                                                                                                            
void rk2(double a, double b, double h, int N,double y0, double rk2out[]);