# bistabfor.ode
# periodically forced bistable system
x'=x*(1-x)*(x-a)+c*u
u'=u*(1-u^2-v^2)-w*v
v'=v*(1-u^2-v^2)+w*u
init u=1,x=.043
par c=.1,a=.25,w=1
@ dt=.01,total=6.281,xlo=0,xhi=6.5
@ dsmax=.05,parmin=-1,parmax=1
@ autoxmax=1,autoxmin=-1,autoymin=-.5,autoymax=2
done