function [x,y]=euler(x0,y0,xn,f,N)
  h=(xn-x0)/N;
  xn=x0;
  yn=y0;
  for n=1:N   
    yn=yn+h*f(xn,yn);
    xn=x0+h*n;
    x(n)=xn;
    y(n)=yn;
  end;
