Example 6.  Assume that the initial number of foxes and rabbits are  [Graphics:Images/Lotka-VolterraMod_gr_293.gif] and  [Graphics:Images/Lotka-VolterraMod_gr_294.gif], respectively,  and that the coefficients  [Graphics:Images/Lotka-VolterraMod_gr_295.gif],  are used to form the system of D. E.'s
        [Graphics:Images/Lotka-VolterraMod_gr_296.gif]  
and    
        [Graphics:Images/Lotka-VolterraMod_gr_297.gif]
Solve the system of D. E.'s for x(t)and x(t) over the interval  [Graphics:Images/Lotka-VolterraMod_gr_298.gif].  

Solution 6.

Enter the functions  [Graphics:../Images/Lotka-VolterraMod_gr_299.gif]  and  [Graphics:../Images/Lotka-VolterraMod_gr_300.gif]  and form the vector function  [Graphics:../Images/Lotka-VolterraMod_gr_301.gif].  

[Graphics:../Images/Lotka-VolterraMod_gr_302.gif]



[Graphics:../Images/Lotka-VolterraMod_gr_303.gif]
[Graphics:../Images/Lotka-VolterraMod_gr_304.gif]
[Graphics:../Images/Lotka-VolterraMod_gr_305.gif]
[Graphics:../Images/Lotka-VolterraMod_gr_306.gif]

Compute the Runge-Kutta solution.

[Graphics:../Images/Lotka-VolterraMod_gr_307.gif]


[Graphics:../Images/Lotka-VolterraMod_gr_308.gif]

[Graphics:../Images/Lotka-VolterraMod_gr_309.gif]
[Graphics:../Images/Lotka-VolterraMod_gr_310.gif]
[Graphics:../Images/Lotka-VolterraMod_gr_311.gif]


[Graphics:../Images/Lotka-VolterraMod_gr_312.gif]


[Graphics:../Images/Lotka-VolterraMod_gr_313.gif]
[Graphics:../Images/Lotka-VolterraMod_gr_314.gif]
[Graphics:../Images/Lotka-VolterraMod_gr_315.gif]


If we graph the solution over a larger time interval the graph traces over itself, because the solution is periodic.

[Graphics:../Images/Lotka-VolterraMod_gr_316.gif]


[Graphics:../Images/Lotka-VolterraMod_gr_317.gif]

[Graphics:../Images/Lotka-VolterraMod_gr_318.gif]
[Graphics:../Images/Lotka-VolterraMod_gr_319.gif]
[Graphics:../Images/Lotka-VolterraMod_gr_320.gif]

The graph of the number of prey is formed by using the first coordinate [Graphics:../Images/Lotka-VolterraMod_gr_321.gif].  

[Graphics:../Images/Lotka-VolterraMod_gr_322.gif]


[Graphics:../Images/Lotka-VolterraMod_gr_323.gif]

[Graphics:../Images/Lotka-VolterraMod_gr_324.gif]

The graph of the number of predators is formed by using the second coordinate [Graphics:../Images/Lotka-VolterraMod_gr_325.gif].  

[Graphics:../Images/Lotka-VolterraMod_gr_326.gif]

[Graphics:../Images/Lotka-VolterraMod_gr_327.gif]

[Graphics:../Images/Lotka-VolterraMod_gr_328.gif]

The periodic behavior is clearly present in the latter two graphs.  Notice the phase lag for the number of prey.  Did you expect this to happen?

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004