function y = ese11odeERE(ts)
% 
% Eulero Rosenbrock esponenziale
%
% y = ese11odeERE(ts)
%
% ts numero di time step
%
A = 16*toeplitz([-2,1,0]);
y0 = [sqrt(2)/2;1;sqrt(2)/2];
tstar = 1;
k = tstar/ts;
y = y0;
for n = 1:ts
  Jn = A+diag(2*y);
  y = y+k*phi1m(k*Jn)*(A*y+y.^2);
end
