Theory and Proof

for

The Harvesting Model

 

Background

    The exponential model  [Graphics:Images/HarvestingModelProof_gr_1.gif]  is used to study uninhibited population growth and solution is the exponential function  [Graphics:Images/HarvestingModelProof_gr_2.gif].  When the term  [Graphics:Images/HarvestingModelProof_gr_3.gif] is added we obtain the logistic differential equation which is used to model inhibited population growth or bounded population growth.  The logistic differential equation is

            [Graphics:Images/HarvestingModelProof_gr_4.gif].  

One form of the solution is  

            [Graphics:Images/HarvestingModelProof_gr_5.gif].  
            
The terms have been carefully determined so that the initial condition is  

            [Graphics:Images/HarvestingModelProof_gr_6.gif].

The limiting value  L  of  y(t)  is given by  
    
            [Graphics:Images/HarvestingModelProof_gr_7.gif].

The graph is the so called "S-shaped" curve. The choice of parameters  [Graphics:Images/HarvestingModelProof_gr_8.gif]  creates the curve shown below.   

[Graphics:Images/HarvestingModelProof_gr_9.gif]

              [Graphics:Images/HarvestingModelProof_gr_10.gif]

Proof  Logistic Differential Equation  Logistic Differential Equation  

 

 

Harvesting a Logistic Population

    When the harvesting term  -k  is incorporated into into bounded population model we have
    
            [Graphics:Images/HarvestingModelProof_gr_11.gif].  

There are three solution forms for this differential equation, and they correspond to the nature of the stationary solutions  ( x(t) = c).  

 

Definition(Stationary Points)  The stationary points of the D. E.  [Graphics:Images/HarvestingModelProof_gr_12.gif]  are solutions where [Graphics:Images/HarvestingModelProof_gr_13.gif] and are the roots of the characteristic equation  

            [Graphics:Images/HarvestingModelProof_gr_14.gif].

The roots are known to be  [Graphics:Images/HarvestingModelProof_gr_15.gif],  and the stationary solutions are  [Graphics:Images/HarvestingModelProof_gr_16.gif].  

Remark.  Since  x(t)  is a real function, there are no stationary solutions when  [Graphics:Images/HarvestingModelProof_gr_17.gif].  

 

Case (i) If   [Graphics:Images/HarvestingModelProof_gr_18.gif]  there is one stationary solution [Graphics:Images/HarvestingModelProof_gr_19.gif].

    When  [Graphics:Images/HarvestingModelProof_gr_20.gif],  the differential equation has the form  [Graphics:Images/HarvestingModelProof_gr_21.gif]  and the solution is  

            [Graphics:Images/HarvestingModelProof_gr_22.gif].      
[Graphics:Images/HarvestingModelProof_gr_23.gif]

  The solution with the initial condition  [Graphics:Images/HarvestingModelProof_gr_24.gif]  is

            [Graphics:Images/HarvestingModelProof_gr_25.gif].
        
If  [Graphics:Images/HarvestingModelProof_gr_26.gif]  then   [Graphics:Images/HarvestingModelProof_gr_27.gif].

If  [Graphics:Images/HarvestingModelProof_gr_28.gif]  then function  x(t)  has a vertical asymptote at  [Graphics:Images/HarvestingModelProof_gr_29.gif]
and the population   x(t)  becomes extinct at some time  [Graphics:Images/HarvestingModelProof_gr_30.gif]  (where [Graphics:Images/HarvestingModelProof_gr_31.gif]),  i. e.  

            [Graphics:Images/HarvestingModelProof_gr_32.gif].

Proof (i).

 

 

Case (ii)  If  [Graphics:Images/HarvestingModelProof_gr_82.gif] there are two stationary solutions  [Graphics:Images/HarvestingModelProof_gr_83.gif]  and  [Graphics:Images/HarvestingModelProof_gr_84.gif].

    When  [Graphics:Images/HarvestingModelProof_gr_85.gif],  the differential equation has the form  [Graphics:Images/HarvestingModelProof_gr_86.gif]  and the solution is  

            [Graphics:Images/HarvestingModelProof_gr_87.gif]
[Graphics:Images/HarvestingModelProof_gr_88.gif]

  The solution with the initial condition  [Graphics:Images/HarvestingModelProof_gr_89.gif]  is

    [Graphics:Images/HarvestingModelProof_gr_90.gif].

If   [Graphics:Images/HarvestingModelProof_gr_91.gif]   then   [Graphics:Images/HarvestingModelProof_gr_92.gif].  
    
If   [Graphics:Images/HarvestingModelProof_gr_93.gif]   then the population   x(t)  becomes extinct at some time  [Graphics:Images/HarvestingModelProof_gr_94.gif],  i. e.   [Graphics:Images/HarvestingModelProof_gr_95.gif].  

Proof (ii).

 

 

Case (iii)  If  [Graphics:Images/HarvestingModelProof_gr_167.gif] there are no stationary solutions.

When  [Graphics:Images/HarvestingModelProof_gr_168.gif],  the differential equation has the form  [Graphics:Images/HarvestingModelProof_gr_169.gif]  and the solution is  

            [Graphics:Images/HarvestingModelProof_gr_170.gif].   
[Graphics:Images/HarvestingModelProof_gr_171.gif]

  The solution with the initial condition  [Graphics:Images/HarvestingModelProof_gr_172.gif]  is

            [Graphics:Images/HarvestingModelProof_gr_173.gif].

The function  x(t)  has a vertical asymptote at  [Graphics:Images/HarvestingModelProof_gr_174.gif]

so the population   x(t)  becomes extinct at some time  [Graphics:Images/HarvestingModelProof_gr_175.gif]  (where [Graphics:Images/HarvestingModelProof_gr_176.gif].), i.e.  

            [Graphics:Images/HarvestingModelProof_gr_177.gif].

Proof (iii).

 

 

 Return to Numerical Methods - Numerical Analysis

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004