tocke(R):=block(kill(a),a[0]:0.01, 
a[n]:=R*a[n-1]*(1-a[n-1]),
niz: makelist(float(a[n]),n,0,1000),
niz2:lastn(niz,100),
niz3: makelist([R,niz2[i]],i,1,100));

toc: apply(append,makelist(tocke(i),i,1,4,0.01))$
wxplot2d([discrete,toc],[style,[points,0.2]], [color,black],[x,1,4],[y,0,1]);

kill(a); 
R:3.82;
a[0]:0.01;
a[n]:=R*a[n-1]*(1-a[n-1]);
niz: makelist(float(a[n]),n,0,1000)$

niz2: lastn(niz,100);

niz3: makelist([R,niz2[i]],i,1,100);

wxplot2d([discrete,niz3],[style,[points,1]], [color,red],[x,1,4],[y,0,1]);

kill(a); 
R:3.83;
a[0]:0.01;
a[n]:=R*a[n-1]*(1-a[n-1]);
niz: makelist(float(a[n]),n,0,1000)$

niz2: lastn(niz,100);

niz3: makelist([R,niz2[i]],i,1,100);

wxplot2d([discrete,niz3],[style,[points,1]], [color,red],[x,1,4],[y,0,1]);