function x = sol2(t, lambda, x0)
% E' la soluzione del problema ai valori iniziali
% x'=lambda(x-exp(-t)), x(0)= x0
% definito dalla function f(t,x,lambda) 

% Il problema e' "molto" nonstabile 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]=euler_forw(@f2,[0 10],1,0.15,-10);plot(t,x,t,sol2(t,-10,1))
% [t,x]=euler_forw(@f2,[0 10],1,0.2,-10);plot(t,x,t,sol2(t,-10,1))

x = (x0 - lambda/(lambda + 1)) * exp(lambda*t);
x = x + lambda/(lambda+1) * exp(-t);

