nmax = 200;
nout = 100;

A = [0:0.001:4];
X = 0.2 * ones(nmax,length(A));

for n = 1: nmax - 1 do
  X(n+1,:) =  A .* X(n,:) .* (1.0 - X(n,:));
end

plot(A,X(nout:nmax,:),'.r','markersize',1)