Defining the Model
Next we study solutions of the Lotka-Volterra system (similar to the system discussed at the bottom of page 141 in the text).. In this model the prey is assumed to grow logistically in the absence of any predators. Can you see how the rate equations have been changed from the original L-V model to incorporate this assumption? This time represents the hare (rabbit) population and represents 100 times the lynx (fox) population. We are going to work with three different initial conditions.
> rate_eqn1 := diff(h(t),t) = (0.1)*h(t)-.00001*h(t)^2-(0.005)*h(t)*(1/100)*u(t) ;
> rate_eqn2 := diff(u(t),t) = (0.00004)*h(t)*u(t)-(0.04)*u(t) ;
> MODEL := { rate_eqn1, rate_eqn2 };
> vars := { h(t), u(t) } ;
> init1 := [h(0)=2000,u(0)=500];
> init2 := [h(0)=2000,u(0)=1000];
> init3 := [h(0)=2000,u(0)=5000];
> domain := t = 0 .. 380;