###########################################################
#            Problem definition
###########################################################
#
#
#       BIFURCATION PARAMETER
p  mu=-15
#
#       other parameter values defining model
#
p a=1,eta=0.5,uhat=1.5
#
#       function definitions
#
F(x)=a*((x-uhat)^2-eta^2)
G(x,c)=x^3-3*(x+1)+mu
#
#      differential equations
#
x'=y
y'=-F(x)*y-G(x,mu)
#
#      initial conditions (optional)
#
x(0)=3.0000
y(0)=0.0000
#
@ total=400,xlo=-3,xhi=4,ylo=-4,yhi=4,dt=0.05,maxstor=10000
#
d