x=0.5;
nmax=2000;
% cascade
a=3.567
%window
a=3.843
% intermittency
a=3.82837
y=[x];
for i=1:nmax

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

end

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