# logbif.ode
# the logistic map bifurcation diagram
f(x)=a*x*(1-x)
x'=f(x)
a'=a
init x=.1
init a=2
@ maxstor=100000,total=500,trans=350,meth=discrete
@ xlo=2,xhi=4.001,ylo=0,yhi=1.001,xp=a,yp=x
done