function jac=jac_es_stiff3(t,x)
jac= [ -(1-x(2)), x(1) + 0.8;
    1.e2*( 1-x(2) ), 1.e2*( -x(1) - 0.9)];