function z = es_stiff2(t, x)
% Definisce l'eq differenziale
% x' = x^2 - x^3
% L'IVP con
% x(0) = delta,  0<= t <= 2/delta
% e' stiff se delta e' piccolo
% 
% ESEMPI:
% opts=odeset('RelTol',1.e-4,'Stats','on');
% 1:  delta =  0.01;
%     ode45(@es_stiff2, [0 2/delta], delta, opts)
%     ode2r(@es_stiff2, [0 2/delta], delta, opts)
% 2:  delta = 0.00001;
%     ode45(@es_stiff2, [0 2/delta], delta, opts)
% 3:  ode2r(@es_stiff2, [0 2/delta], delta, opts)
%
% Prova anche ode23 con delta = 0.01 (l'accuratezza e' inferiore rispetto a
% ode45) 
% 
z = x^2 - x^3;