DIM prob(1:100)
LET nbin = 100                    ! number of bins
LET delta_x = 1/nbin
LET n = 0                         ! number of iterations
LET x = 0.6
LET nmax = 10000
DO while n <= nmax
   LET x = 4*x*(1-x)
   LET n = n + 1
   LET ibin = x/delta_x
   LET ibin = truncate(ibin,0) + 1
   LET prob(ibin) = prob(ibin) + 1
LOOP
FOR ibin = 1 to nbin
    IF prob(ibin) > 0 then PRINT ibin,prob(ibin)/nmax
NEXT ibin
END