x=0.5;
nmax=1000;
a=3.567
y=[x];
for i=1:nmax

x=a*x*(1-x);
y=[y x];

end

plot(y,'+')
axis([0 nmax .3 1]);
