
z=.4999; a=4; x=z;
for i=1:10000;
    z=a*z*(1-z);
    x=[x z];
end
figure
subplot(1,2,1)
plot(0:10000,x,'.')
title('x0=0.4999')
z=.5; a=4; x=z;
for i=1:10000;
    z=a*z*(1-z);
    x=[x z];
end
subplot(1,2,2)
plot(0:10000,x,'.')
title('x0=0.5000')
