function [y,t]=euler(f,a,b,y0,n)
h=(b-a)/n; 
t=a:h:b;
y=zeros(size(t)); 
y(1)=y0;
for k=1:n
y(k+1)=y(k)+h*f(t(k),y(k));
end
t
y