# A reduced activator-inhibitor system   ch6-ai.ode 
# See Chapter 5.6 or 6.3 in our textbook for background on these equations.

dx/dt = (sigma/(1+y))*x^2-x
dy/dt = rho*(x^2-y)

# Set the default values of the parameters:
# Don't put spaces between the parameter name, the equal sign, and the parameter value!
par sigma=2.1,rho=1.1

# Set the default initial conditions:
init x=1.0  y=0.5

# Set the default axes for a phase plane plot of y versus x: 
@ xplot=x, yplot=y

# Set the default stopping time and time step
@ total=20, dt=0.01

# Set the default viewing window
@ xlo=0.0,xhi=2.5,ylo=0.0,yhi=4.5
done