function z=es_stiff3(t,x)
% Modellizza reazione chimica tra 2 sostanze
% L'IVP con condizione iniziale [1; 0 ] e t in [0 50] e' stiff
% 
% Intervallo critico [0 1.e-2]: ode45 e ode78 di Octave falliscono.
%
% Confrontare ode25 passand, oppure no, lo Jacobiano:

% options=odeset('Stats','on')
% S2r=ode2r(@es_stiff3,[0 50], [1; 0],options);
%
% options=odeset('Jacobian',@jac_es_stiff3,'Stats','on')
% S2rj=ode2r(@es_stiff3,[0 50], [1; 0],options);
% 
z=[ -(1-x(2))*x(1) + 0.8*x(2);
    1.e2 *( (1 - x(2)) * x(1) - 0.9 * x(2))];

    
