#
#  Generic Pitchfork Bifurcation
#
#  Parameter default values
#
p mu=-1
#
#  Vector field definition
#
f(x,y)=x*(mu-x^2)
g(x,y)=-y
#
#  ODE definition
#
x'=f(x,y)
y'=g(x,y)
#
#  Program parameters (optional)
#
@ total=20, xlo=-3,xhi=3,ylo=-3,yhi=3
@ maxstor=2000000
@ nplot=1,xpt=x,yp1=y
@ BUT=DynamicsRocks:fq