%
% risolve y' = l u con dato iniziale y(0) nell'intervallo [0 1]
% usando Eulero esplicito
%
clear all
figure(1);
clf;
set(1,'position',[407 12 612 684]);
set(1,'defaultlinelinewidth',3)
hold on;
%
% dati iniziali
%
global u0 l;
l=1; u0=1;
fplot('esatta(x)',[0 1])
%
n=input('Numero suddivisioni? ');
h=1/n;
%
x0=0;
for i=1:n
  u1=u0+h*l*u0;
  plot([x0 x0+h],[u0 u1],'r*-')
  u0=u1; x0=x0+h;
  pause
end
