Considerthepredator–preyproblemdescribedinSection5.10.2.Thisproblemmaybeextended to consider the effect of culling on the interacting populations by subtracting a term from both equations in (5.28) as follows:
dP/dt=K1P −CPQ−S1P
dQ/dt=K2Q−DPQ−S2Q
Here S1 and S2 are constants which provide the culling level for the populations. Use ode45 to solve this problem with K1 =2, K2 =10, C =0.001, and D=0.002 and initial values of the population P =5000 and Q=100. Assuming that S1 and S2 are equal, experiment with values in the range 1 to 2. There is a wealth of experimental opportunity in this problem and the reader is encouraged to investigate different values of S1 and S2.
