a=0;
b=4;
h=0.5;
N=(b-a)/h;
f=@(x,y)exp(-4*x)-2*y;
x(1)=a;
y(1)=1;
for n=1:N
  z(n)=f(x(n),y(n));
  x(n+1)=x(n)+h;
  y(n+1)=y(n)+h*z(n);
end
plot(x,y)
x
y
