> restart;
> x:=0.57;
> a:=3.7;
> y:=a*x*(1-x);
> m:=[ ];
> for k to 200 do m:=[op(m),[x,x],[x,y]]: x:=y:y:=a*x*(1-x): od:
> plot(m,scaling=constrained);
>