nmax = 200; nout = 100; Mu = [1:0.001:2]; X = 0.5 * ones(nmax,length(Mu)); for m = 1:length(Mu) do for n = 1: nmax - 1 do if X(n,m) < 0.5 then X(n+1,m) = Mu(m) * X(n,m); else X(n+1,m) = Mu(m) * (1 - X(n,m)); end end end plot(Mu,X(nout:nmax,:),'.r','markersize',1) xlabel("$\mu$"); ylabel("$x_n$");