Module for the Lotka-Volterra Model

Check out the new Numerical Analysis Projects page.

 

The Lotka-Volterra Equations. The "Lotka-Volterra equations" refer to two coupled differential equations

    [Graphics:Images/Lotka-VolterraMod_gr_1.gif]  
    
    [Graphics:Images/Lotka-VolterraMod_gr_2.gif]

There is one critical point which occurs when [Graphics:Images/Lotka-VolterraMod_gr_3.gif] and it is [Graphics:Images/Lotka-VolterraMod_gr_4.gif].  

Background. The Runge-Kutta method is used to numerically solve O.D.E.'s over [Graphics:Images/Lotka-VolterraMod_gr_5.gif].

Program (Runge-Kutta Method)  To compute a numerical approximation for the solution of the initial value problem [Graphics:Images/Lotka-VolterraMod_gr_6.gif] with [Graphics:Images/Lotka-VolterraMod_gr_7.gif] over [Graphics:Images/Lotka-VolterraMod_gr_8.gif] at a discrete set of points using the formula  

    [Graphics:Images/Lotka-VolterraMod_gr_9.gif],  

where [Graphics:Images/Lotka-VolterraMod_gr_10.gif], [Graphics:Images/Lotka-VolterraMod_gr_11.gif], [Graphics:Images/Lotka-VolterraMod_gr_12.gif], and [Graphics:Images/Lotka-VolterraMod_gr_13.gif].

[Graphics:Images/Lotka-VolterraMod_gr_14.gif]

Example 1.  Solve the I.V.P.   [Graphics:Images/Lotka-VolterraMod_gr_15.gif]  over  [Graphics:Images/Lotka-VolterraMod_gr_16.gif].  Use the Runge-Kutta method.
Solution 1.

 

Extension to 2D. The Runge-Kutta method is easily extended to solve a system of D.E.'s over the interval  [Graphics:Images/Lotka-VolterraMod_gr_24.gif].

Program (Runge-Kutta Method in 2D space)  To compute a numerical approximation for the solution of the initial value problem  

    
[Graphics:Images/Lotka-VolterraMod_gr_25.gif]  with  [Graphics:Images/Lotka-VolterraMod_gr_26.gif],  
    
[Graphics:Images/Lotka-VolterraMod_gr_27.gif]  with  [Graphics:Images/Lotka-VolterraMod_gr_28.gif],  

over the interval  
[Graphics:Images/Lotka-VolterraMod_gr_29.gif]  at a discrete set of points.

[Graphics:Images/Lotka-VolterraMod_gr_30.gif]

Note.  The Runge-Kutta method in 2D is a "vector form" of the one-dimensional method, here the function f is replaced with F.

 

Example 2.  Lotka-Volterra Model.  Solve the I.V.P.  

    [Graphics:Images/Lotka-VolterraMod_gr_31.gif]  with  [Graphics:Images/Lotka-VolterraMod_gr_32.gif],  
and
    [Graphics:Images/Lotka-VolterraMod_gr_33.gif]  with  [Graphics:Images/Lotka-VolterraMod_gr_34.gif].    

Use several intervals  
[Graphics:Images/Lotka-VolterraMod_gr_35.gif].
2 (a).  Use the interval   [Graphics:Images/Lotka-VolterraMod_gr_36.gif].
2 (b).  Use the interval   [Graphics:Images/Lotka-VolterraMod_gr_37.gif].
2 (c).  Use the interval   [Graphics:Images/Lotka-VolterraMod_gr_38.gif].
Can you discover if the solution form an "orbit."
Solution 2.

 

Example 3.  Lotka-Volterra Model.  Solve the I.V.P.  

    [Graphics:Images/Lotka-VolterraMod_gr_78.gif]  with  [Graphics:Images/Lotka-VolterraMod_gr_79.gif],  
and
    [Graphics:Images/Lotka-VolterraMod_gr_80.gif]  with  [Graphics:Images/Lotka-VolterraMod_gr_81.gif].    

Combine the system of D. E.'s to form a separable first-order differential equation and solve the D. E..
Solution 3.

 

Example 4.  For the  I.V.P.  

    [Graphics:Images/Lotka-VolterraMod_gr_113.gif]   with   [Graphics:Images/Lotka-VolterraMod_gr_114.gif],  
and
    [Graphics:Images/Lotka-VolterraMod_gr_115.gif]   with   [Graphics:Images/Lotka-VolterraMod_gr_116.gif].    

Show that the numerical solution in Example 3 and the analytic solution in Example 4 are in agreement.
Solution 4.

 

Example 5.  For the  I.V.P.  

    [Graphics:Images/Lotka-VolterraMod_gr_127.gif]   with   [Graphics:Images/Lotka-VolterraMod_gr_128.gif],  
and
    [Graphics:Images/Lotka-VolterraMod_gr_129.gif]   with   [Graphics:Images/Lotka-VolterraMod_gr_130.gif].    

The implicit solution is  
[Graphics:Images/Lotka-VolterraMod_gr_131.gif].
Determine if y can be solved as a function of x.
Solution 5.

 

 

Research Experience for Undergraduates

The Lotka-Volterra Model  The Lotka-Volterra Model  
Internet hyperlinks to web sites and a bibliography of articles.  
   

Downloads (The Lotka-Volterra Model The Lotka-Volterra Model).  
Download this Mathematica notebook.  

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2003