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).

The two real roots of the characteristic equation  [Graphics:../Images/HarvestingModelProof_gr_96.gif],  are  [Graphics:../Images/HarvestingModelProof_gr_97.gif].  

[Graphics:../Images/HarvestingModelProof_gr_98.gif]

[Graphics:../Images/HarvestingModelProof_gr_99.gif]

Thus there are two stationary solutions  [Graphics:../Images/HarvestingModelProof_gr_100.gif]  and  [Graphics:../Images/HarvestingModelProof_gr_101.gif].

We will verify that  [Graphics:../Images/HarvestingModelProof_gr_102.gif]  satisfies the D. E.  [Graphics:../Images/HarvestingModelProof_gr_103.gif].   

Aside.  There are three equivalent forms for the solution for case (ii).

 

[Graphics:../Images/HarvestingModelProof_gr_104.gif]

This third form is the one we will promote.

[Graphics:../Images/HarvestingModelProof_gr_105.gif]



[Graphics:../Images/HarvestingModelProof_gr_106.gif]
[Graphics:../Images/HarvestingModelProof_gr_107.gif]
[Graphics:../Images/HarvestingModelProof_gr_108.gif]
[Graphics:../Images/HarvestingModelProof_gr_109.gif]

[Graphics:../Images/HarvestingModelProof_gr_110.gif]

[Graphics:../Images/HarvestingModelProof_gr_111.gif]
[Graphics:../Images/HarvestingModelProof_gr_112.gif]
[Graphics:../Images/HarvestingModelProof_gr_113.gif]

[Graphics:../Images/HarvestingModelProof_gr_114.gif]
[Graphics:../Images/HarvestingModelProof_gr_115.gif]
[Graphics:../Images/HarvestingModelProof_gr_116.gif]

[Graphics:../Images/HarvestingModelProof_gr_117.gif]
[Graphics:../Images/HarvestingModelProof_gr_118.gif]

It is difficult for Mathematica to find the limit of  x[t]  because of the arbitrary constants (which are assumed to be complex symbols in Mathematica).

So we will need to be clever in taking the limit.

[Graphics:../Images/HarvestingModelProof_gr_119.gif]



[Graphics:../Images/HarvestingModelProof_gr_120.gif]

[Graphics:../Images/HarvestingModelProof_gr_121.gif]
[Graphics:../Images/HarvestingModelProof_gr_122.gif]


[Graphics:../Images/HarvestingModelProof_gr_123.gif]
[Graphics:../Images/HarvestingModelProof_gr_124.gif]

[Graphics:../Images/HarvestingModelProof_gr_125.gif]
[Graphics:../Images/HarvestingModelProof_gr_126.gif]

[Graphics:../Images/HarvestingModelProof_gr_127.gif]

[Graphics:../Images/HarvestingModelProof_gr_128.gif]

[Graphics:../Images/HarvestingModelProof_gr_129.gif]

[Graphics:../Images/HarvestingModelProof_gr_130.gif]
[Graphics:../Images/HarvestingModelProof_gr_131.gif]

[Graphics:../Images/HarvestingModelProof_gr_132.gif]

[Graphics:../Images/HarvestingModelProof_gr_133.gif]

[Graphics:../Images/HarvestingModelProof_gr_134.gif]

[Graphics:../Images/HarvestingModelProof_gr_135.gif]

[Graphics:../Images/HarvestingModelProof_gr_136.gif]

[Graphics:../Images/HarvestingModelProof_gr_137.gif]

This seems to verifies the statement

        If   [Graphics:../Images/HarvestingModelProof_gr_138.gif]  then  [Graphics:../Images/HarvestingModelProof_gr_139.gif].  

To find the solution with the initial condition  [Graphics:../Images/HarvestingModelProof_gr_140.gif],  proceed as follows.

[Graphics:../Images/HarvestingModelProof_gr_141.gif]



[Graphics:../Images/HarvestingModelProof_gr_142.gif]
[Graphics:../Images/HarvestingModelProof_gr_143.gif]

[Graphics:../Images/HarvestingModelProof_gr_144.gif]
[Graphics:../Images/HarvestingModelProof_gr_145.gif]

[Graphics:../Images/HarvestingModelProof_gr_146.gif]
[Graphics:../Images/HarvestingModelProof_gr_147.gif]

[Graphics:../Images/HarvestingModelProof_gr_148.gif]
[Graphics:../Images/HarvestingModelProof_gr_149.gif]

[Graphics:../Images/HarvestingModelProof_gr_150.gif]
[Graphics:../Images/HarvestingModelProof_gr_151.gif]

[Graphics:../Images/HarvestingModelProof_gr_152.gif]

[Graphics:../Images/HarvestingModelProof_gr_153.gif]

[Graphics:../Images/HarvestingModelProof_gr_154.gif]
[Graphics:../Images/HarvestingModelProof_gr_155.gif]

Example (ii).  Find the solution to the D. E.  [Graphics:../Images/HarvestingModelProof_gr_156.gif] using the general solution and compare it with the one found with Mathematica.

Solution (ii).

[Graphics:../Images/HarvestingModelProof_gr_157.gif]

[Graphics:../Images/HarvestingModelProof_gr_158.gif]

[Graphics:../Images/HarvestingModelProof_gr_159.gif]

[Graphics:../Images/HarvestingModelProof_gr_160.gif]

Now find the solution with the initial condition  [Graphics:../Images/HarvestingModelProof_gr_161.gif].

[Graphics:../Images/HarvestingModelProof_gr_162.gif]

[Graphics:../Images/HarvestingModelProof_gr_163.gif]

[Graphics:../Images/HarvestingModelProof_gr_164.gif]

[Graphics:../Images/HarvestingModelProof_gr_165.gif]

[Graphics:../Images/HarvestingModelProof_gr_166.gif]

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews 2004