#include "f.hh"
#include "math.h"

//
// Euler function
//
void euler(double a, double b, double h, int N,double y0, double eulerout[])
{
  double y=y0;
  double x=a;
  double k1;
  for(int i=0; i<N; i++)
    {
      eulerout[i]=y;
      k1=h*f(x,y);
      y = y + k1;
      x=x+h;
    }
}