% Euler method for y'=-y+t
%%%%%%%%%%%%%%%%%%%%%%%%%%
y0=5; tmax=20; N=10; h=tmax/N; 
T=zeros(1,N+1); Y=T; Ye=T; 
% Fill in the arrays T et Y
T(1)=0;Y(1)=y0; Ye(1)=y0;
for i=2:N+1 
  T(i)=i*h; 
  Y(i)=Y(i-1)*(1-h)+h*T(i-1); 
  Ye(i)=T(i)-1+(y0+1)*exp(-T(i));
end
plot(T,Ye,T,Y)