%%
%close all 
clear all
avalues=2.8:0.002:4;
N=800; a=avalues; x=0.1;
X=zeros(N,length(a));
for n=1:.3*N 
    x=a.*x.*(1-x); 
    X(n,:)=x;
end
figure (2), hold on
for n=.3*N:N 
    x=a.*x.*(1-x); 
    X(n,:)=x;
    plot(a,x,'.','MarkerSize',0.05) 
    axis ([2.8 4 0 1]) 
end
hold off