# Chemostat model   ch4-chemostat.ode 
# See Chapter 4.4.1 (Nullclines and Trapping Regions) in our textbook for details.

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

# 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,rho=0.5

# Set the default initial conditions:
init x=0.1  y=0.2

# 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=40, dt=0.01

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