clf;
Nt=800;
thetamin=2; thetamax=4; 
theta=linspace(thetamin,thetamax,Nt);
figure(1);
axis([thetamin thetamax 0 1]);
hold;
x0=0.33;
Niter=600; plotstart=300;
e=ones(1,Niter-plotstart+1);
for i=1:Nt,
      X=simpleiter(x0,theta(i),Niter);
      plot(theta(i)*e,X(plotstart:Niter),'.','MarkerSize',0.5);
end;
hold