# Modified Lotka-Volterra equations ch1-modified-LV.ode 
dx/dt = x*(1-x/K)*phi(x) - x*y
dy/dt = rho*(x*y - y)
phi(x)=(x-eps)/(x+eps)

# Set the default parameters and initial conditions
param K=2.0,eps=.10,rho=1.0
init x=3,  y=0.2

# Set the axes to do a phase portrait.  Change the next line to xplot=t, yplot=x to plot x versus t. 
@ xplot=x, yplot=y

# Set the stopping time, time step, and default viewing window to the first quadrant
@ total=50,dt=.005,xlo=0.0,ylo=0.0,xhi=4.0,yhi=3.0

# Set the maximum number of data points to store and the maximum allowable magnitude of the variables
@ maxstor=100010, bound=10000

# silence the nuisance beeping sound
bell=0

done