clear all
close all

% Nonlinear system >>> Very sensitive

% a_{n+1} = r(1-a_n) a_n

% a_0 = 0.1

r = 1.546;
%r = 2.75;
%r = 3.25;
%r= 3.525;
%r = 3.555;
%r = 3.75;

a(1) = 0.1;

for k=1:200
    a(k+1) = r*(1-a(k))*a(k);
end


plot(a,'o')
xlim([1 200])
ylim([0 1])