figure(1)
hold on
%This routine does a bifurcation plot for the discrete logistic
%equation.  The first 201 iterates are discarded, to remove transients.
for jlam=1:601
    lambda=1+(jlam-1)/200;
    xn=0.2;
    for n=1:501
        xnew=lambda*xn*(1-xn);
        if n>201
            plot(lambda,xnew,'.')
        end
        xn=xnew;
    end
end
hold off