from pylab import *                                     # Load matplotlib

ion()

N = 100
x = zeros(N)
ii = arange(1,N+1,1)
mu = 2.8 # single attractor
#mu = 3.3 # double attractor
#mu = 3.5  # quadruple attractor
#mu = 3.55  # 8-ple atrtractor
#mu = 3.8   # chaotic regime
#mu = 3.828427 # triple attracrtor

x[0] = 0.8
for i in arange(0,N-1,1):
    x[i+1] = mu*x[i]*(1 - x[i])

plot(ii,x, '-o', lw=1)                                                      

draw()                                                 # Make screen plot