% 用Nicolis方法求最大Lyapunov指數

clear;
close all;
x(1)=0.6;
y(1)=0.4;
dx=[0.002;0.01];
d0=norm(dx);
e=1/d0*dx;
lle=0;
k=300;
for i=1:k
    J=[-2.8*x(i) 1;0.3 0];
    vec=J*e;
    d=norm(vec);
    e=1/d*vec;
    lle=lle+1/k*log(d);
    x(i+1)=1+y(i)-(1.4).*x(i).*x(i);
    y(i+1)=0.3*x(i);
end
lle