clear all

e=0.1;

Nt=400;

tt=linspace(0,2*pi,Nt);
delta=tt(2)-tt(1)

Phi=eye(2);

for it=1:Nt-1

    ti=(tt(it)+tt(it+1))/2;
    c=cos(ti);s=sin(2*ti);
    A=[0,1;-1+4*e*s,-e*(4*c^2-1)];
    eA=expm(A*delta);
    Phi=eA*Phi;

end


lambda=exp(-2*pi*e)

Phi
eig(Phi)




