clear
close all
n=2; % 維數
%e=ones(n,1);
x(1)=0.6;
y(1)=0.4;
lya1=0;lya2=0;
k=300;
for i=1:k
    J=[-2.8*x(i) 1;0.3 0];
    eig_value=eig(J);
    lamd1=max(abs(eig_value));lamd2=min(abs(eig_value));
    lya1=lya1+1/100*log(lamd1);lya2=lya2+1/100*log(lamd2);
    x(i+1)=1+0.3.*y(i)-1.4.*x(i).*x(i);
    y(i+1)=x(i);
end