// 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
clf
plot2d(aa,yy,style=-0)
xlabel('a')
ylabel('x')
