close all 
clear all
avalues=(1+sqrt(8)-0.0002):0.0001:(1+sqrt(8)+0.0002);
N=1500; a=avalues; x=0.1;
X=zeros(N,length(a));
for n=1:.5*N 
    x=a.*x.*(1-x); 
    X(n,:)=x;
end
figure (1), hold on
for n=.5*N:N 
    x=a.*x.*(1-x); 
    X(n,:)=x;
    plot(a,x,'.','MarkerSize',0.01,'Color',[0,0,0]) 
end
hold off
figure (2), hold on
plot(X(1151:1310,2))
hold off
figure (3), hold on
plot(X(1151:1310,3))