% programa 4.3, W.L.Luyben, pag. 108
% "first-order explicit euler integration" 
% modif. E.Quintero.M 18nov'98
%-------------------------------------------------
% ode is xdot=1-x
% initial conditions: x=0 at t=0
clear
format short e
x=0; t=0; xx=[];
delta=1.6;
while(t<=7)
   xdot=1-x;
   xx=[xx;t,x,xdot,1-exp(-t)];
   x=x+xdot*delta;
   t=t+delta;
end
tt=0:.1:7; 
plot(tt,1-exp(-tt),'-'),grid
hold on
plot(xx(:,1),xx(:,2),'--')
xlabel('tiempo'),ylabel('x')
title('Efecto del paso de integración con el método de Euler')
