function y = ese11odeEE(ts)
%
% Eulero Esponenziale
%
% y = ese11odeEE(ts)
%
% ts numero passi temporali
%
A = 16*toeplitz([-2,1,0]);
y0 = [sqrt(2)/2;1;sqrt(2)/2];
tstar = 1;
k = tstar/ts;
y = y0;
P = phi1m(k*A);
for n = 1:ts
  y = y+k*P*(A*y+y.^2);
end
