Differential Equations Project

Computer Lab Modules

 

The Population Model with Harvesting

 

 

Preliminaries. We shall study the population model with harvesting [Graphics:e5.txtgr1.gif].

 

Computer Lab Work.

 

Exercise 1. First, enter the D.E. into Mathematica and solve it.

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr2.gif]

Exercise 2. Second, enter the characteristic equation and find the roots.

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr4.gif]

There are three possibilities, equal real roots, distinct real roots, and complex roots.
We desire that the D. E. has one or two constant solutions, with a real constant.

Exercise 3. Case (i) One critical point. Suppose that [Graphics:e5.txtgr5.gif].

Show that there is one root of the characteristic equation [Graphics:e5.txtgr6.gif], which is [Graphics:e5.txtgr7.gif]

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr8.gif]

Show that there is one constant solution to the D. E. which is [Graphics:e5.txtgr9.gif],

and that [Graphics:e5.txtgr10.gif].

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr11.gif]

Exercise 4. Solve the population model with harvesting [Graphics:e5.txtgr12.gif]
using the constants a = 2, b = 1, k = 1, and explore this situation.

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr13.gif]

The constant solution is:

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr14.gif]

The general solution is:

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr15.gif]

Exercise 5. Plot some solutions to this D. E.
The constants for solutions with the initial condition x[0] = 2, 3, 4, 5, 6 are:

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr16.gif]

The constants for solutions with the initial condition [Graphics:e5.txtgr17.gif] are:

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr18.gif]

Exercise 6. Discuss the graphs in the above plot.
What are the vertical lines ?
What are the curves that lie below x = 1.
What are the curves that lie above x = 1 ? What use are they ?

In order to clear things up, it is necessary to specify the individual domain,
for each of the solutions, then plot all the curves on the same graph.

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr19.gif]

Delete the output graphs given above and report the following composite graph for your report.

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr20.gif]

Exercise 7. Case (ii) Two critical points. Suppose that [Graphics:e5.txtgr21.gif].

Then there are two real roots of the characteristic equation [Graphics:e5.txtgr22.gif],
they are [Graphics:e5.txtgr23.gif].

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr24.gif]

Exercise 8. Solve the population model with harvesting [Graphics:e5.txtgr25.gif]
using the constants a = 4, b = 1, k = 3, and explore this situation.

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr26.gif]

The constant solutions are:

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr27.gif]

The general solution is:

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr28.gif]

Exercise 9. Plot some solutions to this D. E.
The constants for solutions with the initial condition x[0] = 4, 5, 6, 7, 8 are:

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr29.gif]

The constants for solutions with the initial condition [Graphics:e5.txtgr30.gif] are:

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr31.gif]

You need to get the above list of functions. If you can't seem to do it, then type them in !

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr32.gif]

The solutions with the initial condition [Graphics:e5.txtgr33.gif] are:

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr34.gif]

You need to get the above list of functions. If you can't seem to do it, then type them in !

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr35.gif]

Discuss the graphs in the above plot.
What are the vertical lines ?
What are the curves that lie below x = 1.
What are the curves that lie above x = 3 ? What use are they ?

In order to clear things up, it is necessary to specify the individual domain for
some of the solutions, then plot all the curves on the same graph.

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr36.gif]

Delete the output graphs given above and report the following composite graph for your report.

[Graphics:e5.txtgr3.gif][Graphics:e5.txtgr37.gif]

Exercise 10. Look at the above graph and summarize what happens for the various initial conditions x[0]>0.

 

Solutions.

 



Return to the Differential Equations Project

Return to the Numerical Analysis Project

Return to the Complex Analysis Project

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) John H. Mathews, 1998