clear
close all
n=2; % 維數
%e=ones(n,1);
x(1)=0.6;
y(1)=0.4;
V=diag(ones(1,n));
S=V;
b=0; 
k=300;
for i=1:k
    J=[-2.8*x(i) 1;0.3 0];
    B=J*V*S;
    [V,S,U]=svd(B);
    a_max=max(diag(S));
    S=(1/a_max)*S;
    b=b+log(a_max);
    x(i+1)=1+(0.3).*y(i)-(1.4).*x(i).*x(i);
    y(i+1)=x(i);
end
for i=1:n
    Lyapunov(i)=(log(S(i,i))+b)/k;
    fprintf('%f\n',Lyapunov(i));
end