clear all;
clc;
f=2300;
fs=300;
a = linspace(1,4,f);
for j=1:f;
    lyap=0;
       x=0.1;
       for i=1:f
           y=a(j)*x*(1-x);
           if i>fs
               lyap=lyap+log(abs(a(j)*(1-2*y)));
           end
           x=y;
       end
       ly(j)=lyap/(f-fs);
       
       if lyap<-3
           lyap=-3;
       end
end
plot(a,ly)
grid
%AXIS([3 4 -3 1]);