clear
close all
n=2; % 維數(shù)
%e=ones(n,1);
x(1)=0.6;
y(1)=0.4;
lya1=0;lya2=0;
k=1000;
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/k)*log(lamd1);lya2=lya2+(1/k)*log(lamd2);
    x(i+1)=1+y(i)-1.4*x(i)*x(i);
    y(i+1)=0.3*x(i);
end
lya1,lya2