Visualizing a Continuous Predator-Prey Model

On this page, plots of direction fields and solution curves are generated for the phase plane associated with a system of two 1st-order autonomous ordinary differential equations. Plots of the solution curves as functions of time are given below the phase plane, too. This page focuses on a particular form of differential equations which are most often used to model the population numbers for a predator-prey system. The equations are of the form:

dM/dt = r M (1 - M/K) - b M N
dN/dt = c M N - s M

M(0) = Mo
N(0) = No

The populations, M and N, are usually in units relative to a large population, say in thousands or millions.

The user inputs positive values of the parameters r, K, b, and c.  Also, grid bounds, time bounds, and the initial conditions must be specified. NOTE: If you want to ignore the logistic modeling of the the growth of prey, then enter "Infinity" exactly as spelled for the value of K.

r =
K =
b =
c =
s =
Mo =
No =
tmin =
tmax =
Mmin =
Mmax =
Nmin =
Nmax =