x'=x*(x*(1-x)-y)
y'=y*(x-a)
par a=1.2
@ xlo=0,ylo=0,xhi=1.3,yhi=.5
@ init x=.5,y=0.5
@ xp=x,yp=y
@ total=100
done

# Choose an IC and reach steady state
# integrate using Integrate - Last
# Xmin=-0.1, Xmax=1.3, Ymin=-.1, Ymax=1.2
# make Ds negative, Dsmax=0.02

# hit Run
# hit Grab, and navigate right until Hopf bifurcation.
# hit Enter to grab, and hit Run again.