#include "f.hh"

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

  for(int i=0; i<N; i++)
    {
      rk2out[i]=y;
      k1=h*f(x,y);
      y = y + h*f(x+ h/2, y+ k1/2);
      x=x+h;
    }
}