# Michaelis-Menten kinetics   ch4-mm.ode 
# See Chapter 4.4.5 (Nullclines and Trapping Regions) in our textbook for details.

dx/dt = -x*(1-y)+y
dy/dt = (1/epsilon)*(x*(1-y)-(1+kappa)*y)

# Set the default values of the parameters:
# Don't put spaces between the parameter name, the equal sign, and the parameter value!
par kappa=4,epsilon=0.1

# Set the default initial conditions:
init x=15  y=0.02

# 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=20,ylo=0.0,yhi=1.0
done