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