%Euler Method

function [y] = Euler(x0,h,tmax)
   x(1)=x0;
   time(1)=0;
   n = tmax/h+1;
  for i=2:n
	  time(i)= (i-1)*h;
      x(i)=x(i-1)+h*f(x(i-1),time(i-1));
  end
  plt = plot(time,x);
  set(plt,'Color','blue')
  y = x(n);  