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