from visual import *
from visual.graph import *

pic=gdisplay()
myplot=gdots(gdisplay=pic,radius=0.0001)

def logistic(x,r):
    return 4*r*x*(1-x)


MAX=1000
record=64
y=arange(0,record)

for r in arange(0.0,1.0,0.01):
    print "r is ",r
    x=0.5
    count=0
    i=0;
    for count in range(0,MAX):
  	    x=logistic(x,r)
    for count in range(0,record):
            x=logistic(x,r)
            myplot.plot(pos=(r,x))