close all;





mu = linspace(0.8,1,500);





N_T = 150;


N_SS = 100;





hold on;





for j = 1:length(mu)


    [i x] = iterate(mu(j),0.2,N_T+N_SS);


    plot(mu(j),x([N_T:N_T+N_SS]),'k.');


end