function z=f4(t,x,lambda)
% Example 4.6 Ascher Petzold
% IVP x'=lambda*(x-sin(t))+cos(t), x(0)=0
% La soluzione esatta e' x(t)=sin(t), very smooth
% La stima dell'errore locale per la determinazione del passo fallisce,
% in quanto per questo esempio e_i cresce come exp(lambda t_i)*delta,
% dove delta e' l'errore locale di discretizzazione
%
% Con lambra=50, prova a confrontare la soluzione esatta con quella ottenuta con
% [t,x]=ode45(@f4, [0 10], 0,[],50);
% ode45 fallisce
%
% Modificare le tolleranze:
% ops=odeset('AbsTol',1.e-8,'RelTol',1.e-7,'Stats','on');
% [t,x]=ode45(@f4, [0 10], 0,ops,50); 
% 
% Con lambda=1 i risultati sono ragionevoli
z=lambda*(x-sin(t))+cos(t);
