Journal of Nonlinear Dynamics
Volume 2015 (2015), Article ID 948728, 13 pages
http://dx.doi.org/10.1155/2015/948728

Research Article

Prey-Predator Model with Two-Stage Infection in Prey: Concerning Pest Control

Swapan Kumar Nandi,¹ Prasanta Kumar Mondal,² Soovoojeet Jana,³ Palash Haldar,⁴ and T. K. Kar¹

¹Department of Mathematics, Indian Institute of Engineering Science and Technology, Shibpur, Howrah, West Bengal 711103, India
²Shibpur Sri Ramkrishna Vidyalaya, 92 Kshetra Banerjee Lane, Shibpur, Howrah 711102, India
³Department of Mathematics, Abhedananda Mahavidyalaya, Sainthia, West Bengal 731234, India
⁴Department of Mathematics, Sundarban Hazi Desarat College, Pathankhali, West Bengal 743611, India

Received 31 May 2015; Accepted 12 November 2015

Academic Editor: Ivo Petras

Copyright © 2015 Swapan Kumar Nandi et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

A prey-predator model system is developed; specifically the disease is considered into the prey population. Here the prey population is taken as pest and the predators consume the selected pest. Moreover, we assume that the prey species is infected with a viral disease forming into susceptible and two-stage infected classes, and the early stage of infected prey is more vulnerable to predation by the predator. Also, it is assumed that the later stage of infected pests is not eaten by the predator. Different equilibria of the system are investigated and their stability analysis and Hopf bifurcation of the system around the interior equilibriums are discussed. A modified model has been constructed by considering some alternative source of food for the predator population and the dynamical behavior of the modified model has been investigated. We have demonstrated the analytical results by numerical analysis by taking some simulated set of parameter values.

1. Introduction

It is of great interest to consider the effects of pesticide on ecological communities from both an environmental and economical point of view with the rapid development of modern technology, industry, and agriculture. The main foods of mankind are different types of crops or vegetables. Due to the limitations of land for planting, our main interest is not to destroy even a single mole of food. But in nature many times the leaves or branches of those crop trees are consumed by so many small insects, animals, weeds, and so forth for their survival indicating the agriculture resource would severely be damaged from the environment. These insects or animals or weeds are known as pests. In a word an organism which harms agriculture by means of feeding on crops or parasitizing livestock is known to us as a pest. The term pest not only is used to refer to harmful animals but also relates to all other harmful organisms, namely, fungi, bacteria, virus parasites, and so forth. It is also possible that an animal may be a pest in one case but beneficial or domesticated in another situation. But in this paper by the word pest we mean those livelihoods which are involved in the cause of the damage to the human populations directly or indirectly.

The control of pest is a worldwide problem. The most common technique of the pest control is use of chemical insecticides to control the pest population. The use of chemical insecticides is an attempt to control pest directly at low cost (see Uboh et al. [1]). But in many research works it is shown that these chemicals have many environmental effects that include chemical residing in crop and in the agricultural ecosystems. In nature, other than pest, there are some other types of population species known as natural enemy, which predates those pests. For an example, there are some birds, frog, and so forth which feed on agricultural pests and thus those species can be treated as predator species to the pests. So these populations can be used as one of the biological controls for the pest. A prey-predator model can be used to study pest control phenomenon (see the works of Tang et al. [2], Ghosh et al. [3], Ghosh and Bhattacharya [4], Pei et al. [5], Liang et al. [6], Apreutesei [7], Kar et al. [8], Jana and Kar [9, 10], and Mondal et al. [11] and references therein). To control the pest species without using any pesticide control, use of natural enemy is one of the important ways. Again the pest population may be affected due to the infection of some bacterial or viral diseases. So another important natural way to control pest population is the consideration of the effect of some infectious diseases among the targeted pest population. Sun and Chen [12] and Tan and Chen [13] used infection to control the pest population. Some viral or bacterial infection caused the reduction of pest species significantly. Wang and Song [14] used mathematical model by infected pest to control the pest population. Hadeler and Freedman [15], Venturino [16, 17], Kar and Mondal [18], Liu and Wang [19], Guo and Chen [20], Y. Zhang and Q. Zhang [21], Wang et al. [22], and Wang and Chen [23] also worked on prey-predator model with the disease in the prey populations. In our present paper we consider two stages of infected pest, namely, infected pest in the first stage and infected pest in the last stage. The two stages of the infection are considered on the assumption that the predator population would consume infected pest of the first stage but the predator would not consume the pest population at the second stage.

Again a prey-predator system plays an important role when an alternative source of food to the predator population is supplied (see Pahari and Kar [24], Srinivasu and Prasad [25],   Srinivasu [26],  Sabelis and van Rijn [27], Harwood and Obryeki [28], etc.). To control the pest, the effect of alternative food takes some attention compared to ordinary predator-prey dynamics in the context of biological conservation. For both the conservation of predator and reduction of the pest, an alternative source of food must be given to the natural enemy of the pest (see Jana and Kar [10]).

The paper is organized in the following schemes: In Section 2, we formulate the pest control model on the basis of some assumptions and in Section 3, we analyze the complex dynamics of that model. In Section 4, we analyze the model if predator populations would be provided with some supplementary food. In Section 5, we present some numerical simulations for validating our model and in the final section we give some outcomes from our whole paper.

2. Model Formulation and Description

The following assumptions are made to build the mathematical model for pest control.

P1. In this paper, we consider mainly two types of species: one is prey, that is, pest, denoted by , and the other is predator, denoted by But in the presence of disease the pest population is divided into two classes, namely, susceptible and infected . Again the infected pest population is classified into two stages: one is the early stage or first stage of infected pest and the other is the later stage or last stage of infected pest . Therefore, at any time the total pest population is .

P2. We consider that the susceptible pest population grows logistically in the absence of the disease (see Das et al. [29]). So if is the intrinsic growth rate and is the carrying capacity, then the growth equation of susceptible pest is

P3. Infected pest populations can be divided into two classes, namely, infected pest in the first stage and infected pest in the last stage. The susceptible pest becomes infected through direct contact with the early stage and later stage of infected pest at rates and , respectively, following the mass action law. After primary infection the pest population becomes more infected by the same disease because the infection takes some time to spread among the population. In this scenario, is assumed as active rate of the disease. The infected pest dies linearly in two stages at rates and , respectively. Therefore the model for the pest population can be written as

P4. The predator population consumes only the susceptible pest and the early stage of infected pest and consciously avoids the latter stage of infected pest for its better survival. We choose the predation function as Holling type II functional response , where is the capturing rate and is the half saturation constant, to predate the susceptible pest as some searching time is needed for capturing and handling the prey. But as the primary infected pest population is somewhat weaker than the susceptible pest, no more searching time is needed by the predators for catching the infected pest and thereby the predation function is assumed as Holling type I functional response in the form , where is the maximum capturing rate.

P5. Conversion factors for the predator are assumed as and from the susceptible and primary infected pest, respectively. Also, to improve this model, we incorporate the density dependent factor for predators, where is the density dependent mortality rate. Let be the death rate of the predator.

Thus, keeping all the above assumptions in mind, we can obtain the following ecoepidemic model as follows:

with the initial conditions ,  ,  , and  .

3. Analysis of the System

In this section we will study the different dynamical behavior like uniform boundedness, existence of different equilibria and their stability, and so forth.

3.1. Boundedness of the System

Lemma 1. System (3) is uniformly bounded in the region .

Proof. See Appendix A.

3.2. Existence of Equilibria

System (3) has the following equilibria:

(i)the trivial equilibrium ;(ii)the axial equilibrium ;(iii)the predator-free equilibrium point , where , , and ;(iv)the infected pest-free equilibrium , where and is a positive root of the following equation:

 where ,  ,  , and   ;(v)the interior equilibrium point ,    where ,  ,  , and is a positive root of the following equation:

 where ,    , and   ,  ,  .

For positivity of , we assume . The infection-free equilibrium will be feasible if and . In fact, (4) has exactly one root under these parametric conditions.

For (5), we can deduce the following cases.

Case 1. Let us assume that . When , then the constant coefficient is negative. So, by Descartes’ rule of sign (5) has exactly one positive root; that is, unique positive interior equilibrium exists.

Case 2. If , then and , where  . Therefore, (5) will have two positive roots provided . In this case, the largest value of will be considered and thus the unique interior equilibrium will be obtained.
Thus, the positive value of is given by

Furthermore, positivity of implies that .

3.3. Stability Analysis

Theorem 2. The trivial equilibrium is always unstable.

Proof. All the eigenvalues at the trivial equilibrium are ,  ,  , and . Therefore, the trivial equilibrium is unstable.

Remark 3. It is concluded from the analysis of Theorem 2 that the trivial equilibrium is a saddle point with one-dimensional unstable and three-dimensional stable manifold. Biologically, it means that the predator population may become extinct from the environment in absence of pest population.

Theorem 4. The axial equilibrium is locally asymptotically stable if .

Proof. The characteristic equation at is

where ,  , and .
If , then two eigenvalues are always real and negative, and another two eigenvalues will be negative or having negative real parts, hence the theorem.

Remark 5. The above stability analysis biologically indicates that the pest population increases rapidly towards its carrying capacity in absence of disease and predators while the environmental carrying capacity for susceptible pest population is below some threshold value . Moreover, the stability of indicates that the predator-free equilibrium does not exist.

Now, we discuss the stability criteria of predator-free equilibrium. The characteristic equation at is

where ,  ,  , and  .

If , then one root of (8) will be real and negative. The other three roots are given by the second part of (8). Clearly, all are positive and function of So, we assume as bifurcation parameter. Note that .

Hence for any , the predator-free equilibrium is locally asymptotically stable if and only if . But the system enters into a Hopf bifurcation around if .

We shall discuss the Hopf bifurcation using the approach of Bhattacharyya and Bhattacharya [30].

As and ,  . Again, ,  , and   indicate that Since is a continuous function in , by intermediate value theorem there exists such that and . As is a monotonic increasing function and is a decreasing function in , then the critical value is unique. Now, we can easily calculate the first- and second-order derivative of with respect to , and then the following two results are obtained:

The above results imply that

Thus, for all . So, by Routh-Hurwitz criteria, the predator-free equilibrium is locally asymptotically stable in and unstable in . According to the definition of Beretta et al. [31], Guckenheimer and Holmos [32], and Bhattacharyya and Bhattacharya [30], a single Hopf bifurcation occurs at and for decreasing , it approaches a periodic solution. Therefore, this discussion provides the following theorem.

Theorem 6. If , then the predator-free equilibrium is locally asymptotically stable for and unstable for . Furthermore, the system enters into a single Hopf bifurcation around when .

Remark 7. Theorem 6 biologically indicates that the infected pest populations coexist with the susceptible pest population in absence of predator and exhibits oscillatory balance behavior, and for large value of the carrying capacity the susceptible pest population may persist for a long time without using pesticide. Also, if is the period of the bifurcating periodic orbits close to , then we get

Theorem 8. The disease-free equilibrium is locally asymptotically stable if , for and , where the symbols are defined in the proof.

Proof. See Appendix B.

Now, the global stability analysis of the system around the disease-free equilibrium will be examined in the next theorem using the approach of reference name.

Remark 9. From Theorem 8, it is evident that the disease may become extinct from pest population and in this situation predator population can not consume pest easily. Therefore, the interaction between susceptible pest and predator will be ensured during long time for their survival.

Theorem 10. The interior equilibrium is locally asymptotically stable if all for and But the system enters into the Hopf bifurcation around when passes through

Proof. The details of the proof are given in Appendix C.

Remark 11. Theorem 10 provides us with a threshold parameter such that, for , the interior equilibrium is stable and the system bifurcates towards a periodic solution of small amplitude when is increased. In this scenario, we can observe an oscillatory nature of the system. The period of the bifurcating periodic orbits close to is given by

4. Modified Model

Here we study the system when some additional food is to be supplied from some alternating source to the predator populations and assume that the maximum growth rate due to this food is . So here predator populations get their food from two sources, namely, from the prey population and from the alternative source of food. Therefore only the final differential equation of the predator population will be changed and other equations will remain in the same form. So combining these with the previous assumptions we construct our final revised model as follows:

The initial conditions are ,  ,  , and .

4.1. Equilibria of Modified System

System (13) has six possible feasible equilibria as follows:

(i)trivial equilibrium ,(ii)the axial equilibrium ,(iii)the predator-free equilibrium ,(iv)the infected pest-free equilibrium ,(v)the pest-free equilibrium , (where ),(vi)the interior equilibrium , where ,     ,  , and is a positive root of the following equation:

 where ,   , and   − ,  ,  .

Among all these equilibria the coordinates of the first four equilibria are the same as system (3). The main difference between system (3) and system (13) is the fact that in the later system there is one more equilibrium point which does not exist in previous system (3). The equilibrium point is the new equilibrium point and around this equilibrium the system will be totally pest-free. Further this equilibrium point exists if the maximum growth rate of the predator population due to the application of additional food is of higher value than the death rate of the predator population. In the next theorem, we discuss the local asymptotic stability of the pest-free equilibrium .

Theorem 12. System (13) is locally asymptotically stable around the pest-free equilibrium if .

Proof. The characteristic equation of the system around is as follows:

Therefore it is clear that if and , that is, if holds, then all the eigenvalues are negative, hence the theorem.

5. Numerical Simulation

In this section we give some numerical simulations regarding the stability and bifurcation of system (3) around the interior equilibrium . For simulation let us take ,  ,  ,  ,  ,  ,  ,  ,  ,  ,   ,   ,  , and . We see that all the conditions for existence of interior equilibrium hold and the interior equilibrium point is . Now, all the eigenvalues at are , , . Therefore, according to Routh-Hurwitz criterion the interior equilibrium is locally asymptotically stable, depicted in Figure 1. This figure implies that all the trajectories oscillate initially and after some times they all converge their equilibrium points.

Figure 1: Stable solution curves.

Now if the numerical value of the parameter is increased then the oscillation in the solution curve increases. In Figure 2, we see that for oscillation increases but system (3) remains asymptotically stable around Again for and the other parameter value the same as previous we have an unstable interior equilibrium (see Figure 3). Further increment of , say, , and the other parameter value the same as that before make the system extinction of predator population and system is unstable around (see Figure 4).

Figure 2: Stable solution curves and oscillation increases for .

Figure 3: Unstable solution curves for .

Figure 4: Predator extinct for .

Next to compare systems (3) and (13), taking all the other parameters fixed if we take the value of and make the increment in oscillation but system (14) remains asymptotically stable at the interior equilibrium (see Figure 5) but further increment of , say, , and all other parameters having the same values make the system unstable around (see Figure 6). For further increment of , say, for , the new pest-free equilibrium is created and it becomes asymptotically stable and this phenomenon is depicted in Figure 7. Therefore the major important observation of the modified model is that if the additional food is supplied as a sufficient amount then the system may be pest-free.

Figure 5: Stable solution curve when we supplied an alternating source of food (i.e., revised model) with and .

Figure 6: Unstable solution curves of the revised model when and .

Figure 7: Pest-free solution curves of the revised model for and .

6. Conclusions

In the present paper, we have formulated and discussed a prey-predator model by dividing the prey compartment into three classes, where the disease spreads and developed in two stages. Five equilibria are obtained, namely, trivial, axial, predator-free, infection-free, and interior equilibrium, and the dynamical behavior of the system is investigated. The important critical value of the carrying capacity is obtained and it has been found that when the environmental carrying capacity for susceptible pest population is below , the pest population sharply increases in absence of infection and predator. Furthermore, we have seen that the system enters into a single Hopf bifurcation around the predator-free equilibrium and for large value of the carrying capacity the susceptible pest population may persist for a long time without using pesticide.

Also, our global stability analysis of the system around the disease-free equilibrium indicates that the interaction between susceptible pest and predator will be ensured during long time for their survival. We proved that the interior equilibrium is locally asymptotically stable under some parametric conditions and the system enters into the Hopf bifurcation around the interior equilibrium while the density dependent mortality rate passes through the threshold value .

Our main objective is to provide a pest-free system. In this regard, stability of the system at the pest-free equilibrium point is the most targeted job. But in the first model (i.e., system (3)), we observe that there is no pest-free equilibrium point whereas the modified system has a pest-free equilibrium point. In this regard, we may conclude that the additional food on the natural enemy of the pest population has a strong impact to remove the pest population.

Appendices

A. Uniform Boundedness of the System

Let us assume and define a function

Taking the time derivative of (A.1) along the solution of (3), we obtain

Therefore, for each ,

Let . Then from (A.3), we obtain

Now, we assume and in this case, we consider

Then for , we have

The maximum value of is (say).

Then (A.6) can be written as

Applying the theorem of differential inequality [33], we obtain

Thus, for ,

Hence all the solutions of (3) initiating in are confined in the region

B. Local Stability of the Disease-Free Equilibrium

The characteristic equation of system (3) around the disease-free equilibrium is

where

Hence for local stability of the disease-free equilibrium , the following Routh-Hurwitz criterion must be satisfied:

(i),(ii).

C. Local Stability and Hopf Bifurcation of the Interior Equilibrium

The characteristic equation of system (3) around the interior equilibrium is

where

By Routh-Hurwitz criterion, the interior equilibrium is locally asymptotically stable if and .

Now, we shall find out the conditions for which the system enters into Hopf bifurcation around the interior equilibrium .

Let us define a function of as follows:

The necessary and sufficient conditions for Hopf bifurcation to occur are that there exists such that

(i),(ii).

Note that the critical value is obtained from the equation . For , characteristic equation (C.1) becomes

Suppose (C.4) has four roots, say, , with , being the pair of purely imaginary roots. Then we have

where

From the above relations, it is obtained that

If and are complex conjugate, then from above we have , but if they are real roots, then we get and . Further, since is a continuous function of all its roots, then there exists an open interval containing , where and are complex conjugate for

Suppose

Now, we shall verify the transversality condition

Substituting in (C.1) and calculating the derivative, we get

where

Since at , we have

Thus, the Hopf bifurcation conditions are verified.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

References

F. E. Uboh, E. N. Asuquo, M. U. Eteng, and E. O. Akpanyung, “Endosulfan-induces renal toxicity independent of the route of exposure in rats,” American Journal of Biochemistry and Molecular Biology, vol. 1, no. 4, pp. 359–367, 2011.

S. Tang, Y. Xiao, L. Chen, and R. A. Cheke, “Integrated pest management models and their dynamical behaviour,” Bulletin of Mathematical Biology, vol. 67, no. 1, pp. 115–135, 2005.

S. Ghosh, S. Bhattacharyya, and D. K. Bhattacharya, “The role of viral infection in pest control: a mathematical study,” Bulletin of Mathematical Biology, vol. 69, no. 8, pp. 2649–2691, 2007.

S. Ghosh and D. K. Bhattacharya, “Optimization in microbial pest control: an integrated approach,” Applied Mathematical Modelling, vol. 34, no. 5, pp. 1382–1395, 2010.

Y. Pei, X. Ji, and C. Li, “Pest regulation by means of continuous and impulsive nonlinear controls,” Mathematical and Computer Modelling, vol. 51, no. 5-6, pp. 810–822, 2010.

J. Liang, S. Tang, and R. A. Cheke, “An integrated pest management model with delayed responses to pesticide applications and its threshold dynamics,” Nonlinear Analysis: Real World Applications, vol. 13, no. 5, pp. 2352–2374, 2012.

N. C. Apreutesei, “An optimal control problem for a pest, predator, and plant system,” Nonlinear Analysis. Real World Applications, vol. 13, no. 3, pp. 1391–1400, 2012.

T. K. Kar, A. Ghorai, and S. Jana, “Dynamics of pest and its predator model with disease in the pest and optimal use of pesticide,” Journal of Theoretical Biology, vol. 310, pp. 187–198, 2012.

S. Jana and T. K. Kar, “Modeling and analysis of a prey-predator system with disease in the prey,” Chaos, Solitons & Fractals, vol. 47, no. 1, pp. 42–53, 2013.

S. Jana and T. K. Kar, “A mathematical study of a prey-predator model in relevance to pest control,” Nonlinear Dynamics, vol. 74, no. 3, pp. 667–683, 2013.

P. K. Mondal, S. Jana, and T. K. Kar, “A theoretical approach on controlling agricultural pest by biological controls,” Acta Biotheoretica, vol. 62, no. 1, pp. 47–67, 2014.

S. Sun and L. Chen, “Mathematical modelling to control a pest population by infected pests,” Applied Mathematical Modelling, vol. 33, no. 6, pp. 2864–2873, 2009.

Y. Tan and L. Chen, “Modelling approach for biological control of insect pest by releasing infected pest,” Chaos, Solitons & Fractals, vol. 39, no. 1, pp. 304–315, 2009.

X. Wang and X. Song, “Mathematical models for the control of a pest population by infected pest,” Computers & Mathematics with Applications, vol. 56, no. 1, pp. 266–278, 2008.

K. P. Hadeler and H. I. Freedman, “Predator-prey populations with parasitic infection,” Journal of Mathematical Biology, vol. 27, no. 6, pp. 609–631, 1989.

E. Venturino, “Epidemics in predator models: disease among the prey,” in Mathematical Population Dynamics: Analysis of Heterogenity, Theory of Epidemics, O. Arino, D. Axelord, M. KImmel, and M. Langlais, Eds., vol. 1, pp. 381–393, Wuertz Publishing, Winnipeg, Canada, 1995.

E. Venturino, “Epidemics in predator-prey models: disease in the predators,” Journal of Mathematics Applied in Medicine and Biology, vol. 19, no. 3, pp. 185–205, 2002.

T. K. Kar and P. K. Mondal, “A mathematical study on the dynamics of an eco-epidemiological model in the presence of delay,” Applications and Applied Mathematics, vol. 7, no. 1, pp. 300–333, 2012.

X. Liu and C. Wang, “Bifurcation of a predator-prey model with disease in the prey,” Nonlinear DynamicsSystems, vol. 62, no. 4, pp. 841–850, 2010.

H. Guo and L. Chen, “Time-limited pest control of a Lotka-Volterra model with impulsive harvest,” Nonlinear Analysis. Real World Applications, vol. 10, no. 2, pp. 840–848, 2009.

Y. Zhang and Q. Zhang, “Dynamical behavior in a delayed stage-structure population model with stochastic fluctuation and harvesting,” Nonlinear Dynamics, vol. 66, no. 1-2, pp. 231–245, 2011.

X. Wang, Y. Tao, and X. Song, “Analysis of pest-epidemic model by releasing diseased pest with impulsive transmission,” Nonlinear Dynamics, vol. 65, no. 1-2, pp. 175–185, 2011.

T. Wang and L. Chen, “Nonlinear analysis of a microbial pesticide model with impulsive state feedback control,” Nonlinear Dynamics, vol. 65, no. 1, pp. 1–10, 2011.

U. K. Pahari and T. K. Kar, “Conservation of a resource based fishery through optimal taxation,” Nonlinear Dynamics, vol. 72, no. 3, pp. 591–603, 2013.

P. D. N. Srinivasu and B. S. R. V. Prasad, “Role of quantity of additional food to predators as a control in predator–prey systems with relevance to pest management and biological conservation,” Bulletin of Mathematical Biology, vol. 73, no. 10, pp. 2249–2276, 2011.

P. D. N. Srinivasu, B. S. R. V. Prasad, and M. Venkatesulu, “Biological control through provision of additional food to predators: a theoretical study,” Theoretical Population Biology, vol. 72, no. 1, pp. 111–120, 2007.

M. W. Sabelis and P. C. J. van Rijn, “When does alternative food promote biological pest control?” in Proceedings Second International Symposium on Biological Control of Arthropods, M. S. Hoddle, Ed., vol. 1, pp. 428–437, USDA Forest Service, 2005.

J. D. Harwood and J. J. Obryeki, “The role of alternating prey in sustaining predator populations,” in Proceedings Second International Symposium on Biological Control of Arthropods, II, M. S. Hoddle, Ed., pp. 4253–4262, 2005.

K. P. Das, S. Roy, and J. Chattopadhyay, “Effect of disease-selective predation on prey infected by contact and external sources,” BioSystems, vol. 95, no. 3, pp. 188–199, 2009.

S. Bhattacharyya and D. K. Bhattacharya, “Pest control through viral disease: mathematical modeling and analysis,” Journal of Theoretical Biology, vol. 238, no. 1, pp. 177–196, 2006.

E. Beretta, G. I. Bischi, and F. Solimano, “Oscillations in a system with material cycling,” Journal of Mathematical Biology, vol. 26, no. 2, pp. 143–167, 1988.

J. Guckenheimer and P. J. Holmes, Nonlinear Oscillations, Dynamical Systems, and Bifurcations of Vector Fields, Springer, New York, NY, USA, 1983.

G. Birkoff and G. C. Rota, Ordinary Differential Equations, Ginn, Boston, Mass, USA, 1982.