clear all; clc; close all;
% DIAGRAMA DE ORBIT

a = 2.9:.001:3.9; % parámetro de bifurcación
aa=[];yy=[];

for i=1:length(a)
  x(1) = 0.1;
for k = 1:100
  x(k+1) = a(i)*x(k)*(1-x(k));
end

y = [x(91) x(92) x(93) x(94) x(95) x(96) x(97) x(98) x(99) x(100)];
yy=[yy y];
aa=[aa a(i)*ones(1,10)];

end

plot(aa,yy,'k.')
xlabel('a')
ylabel('x')