function z=es_stiff1(t,x,lambda)
% Definisce l'eq differenziale
% x'=lambda*(x -  cos(t));
% Se lambda e' negativo e grande in modulo, l'IVP
% con x(0)=0, t in [0 5] e' stiff
%
% Pet t "non vicino" a zero, la soluzione e' ben approssimata da cos(t).
% 
% Scegliere ad esempio lambda=-1.e3
% e confrontare ode45 con i risolutori per problemi stiff: ode23s in
% Matlab, oppure odebwe, ode2r, ode5r in Octave
% 
% Ad esempio:
% lambda=-1.e3;
% ops=odeset('Stats','on');
% S45=ode45(@es_stiff1,[0 5], 0, ops, lambda);
% Sbe=odebwe(@es_stiff1,[0 5], 0, ops, lambda);
% S2r=ode2r(@es_stiff1,[0 5], 0, ops, lambda);
% 
% Confrontare le statistiche e disegnare le approssimazioni ottenute.
% Osservare, facendo degli zoom, che ode45 produce una approssimazione 
% che oscilla
z = lambda*(x -  cos(t));
