clear all

e=0.25;
                                                                   
Nt=100;
TT=linspace(0,2*pi,Nt);
dt=TT(2)-TT(1);

OM=linspace(1.5,1.51,201);
EE=linspace(0,0.25,21);
for ie=1:length(EE)
    ie
    e=EE(ie);
for io=1:length(OM)
    om=OM(io);
    Phi=eye(2);


for it=1:Nt-1
tt=TT(it);
c=cos(tt);
om2t=om^2*(1+e*c);
A=[0 1;-om2t 0];
EA=expm(A*dt);
Phi=EA*Phi;
end;

m=abs(eig(Phi));
[sd,ind]=sort(-m);
m=m(ind);
l1(ie,io)=log(m(1))/(2*pi);
l2(ie,io)=log(m(2))/(2*pi);
end
end


figure(10)
%plot(OM,l1,'b',OM,l2,'r')
contourf(OM,EE,l1);colorbar