function z=f2(t,x,lambda)
% f(t,x) definisce l'eq diff x'=f(t,x), dipendente dal parametro lambda
% 
% La soluzione dell'IVP 
% x'=f(t,x), x(0)=x0
% e' la funzione definita dalla function sol2(t, lambda, x0)
%
% Il problema e' molto "non-stabile" se lambda > 0 grande
% (esempio: lambda=9, x0=0.9 e x0=1)
% stabile se lambda < 0
% (esempio: lambda=-9, x0=0.9 e x0=1, oppure x0=1.125 e x0=1.225)
%
% Anche qui si osserva il problema della A-stabilita'
% con lambda=-10 si deve scegliere h<2/10
% [t,x]=feuler(@f,[0 100],1,500,-10);plot(t,x,t,sol(t,-10,1))
% [t,x]=feuler(@f,[0 100],1,600,-10);plot(t,x,t,sol(t,-10,1))
z=lambda*(x-exp(-t));
