Applied Mathematics
Vol.09 No.05(2018), Article ID:84855,16 pages
10.4236/am.2018.95036

Global Dynamic Analysis of a Vector-Borne Plant Disease Model with Discontinuous Treatment

Hengmin Lv1*, Lizhi Fei1,2, Zhen Yuan1, Fumin Zhang3

1Department of Basic Course Education, Ji’an College, Ji’an, China

2College of Mathematics, Sichuan University, Chengdu, China

3Key Laboratory of Jiangxi Province for Numerical Simulation and Emulation Techniques, Gannan Normal University, Ganzhou, China

Copyright © 2018 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).

http://creativecommons.org/licenses/by/4.0/

Received: May 2, 2018; Accepted: May 25, 2018; Published: May 28, 2018

ABSTRACT

This paper proposes a vector-borne plant disease model with discontinuous treatment strategies. Constructing Lyapunov function and applying non-smooth theory to analyze discontinuous differential equations, the basic reproductive number R 0 is proved, which determines whether the plant disease will be extinct or not. If R 0 < 1 , the existence and global stability of disease-free equilibrium is established; If R 0 > 1 , there exists a unique endemic equilibrium which is globally stable. The numerical simulations are provided to verify our theoretical results, which indicate that after infective individuals reach some level, strengthening treatment measures is proved to be beneficial in controlling disease transmission.

Keywords:

Vector-Borne, Plant Disease Model, Basic Reproduction Number, Discontinuous Treatment

1. Introduction

The plants play an important role in our lives, as most of our daily food, clothing and building materials come from plants. With the change of environment, there are outbreaks of plant diseases, which seriously affect the health of plants and people’s life, such as huanglongbing [1] , Blackleg [2] . We know that plant diseases have been responsible for the death and suffering of millions of people and countless animals [3] . Controlling the outbreak and spread of plant diseases has become the common goal of scientists. Available control measures include biological, cultural, and chemical methods [4] . Chemical control is a quite effective method, but the residues of chemical drugs have a direct negative impact on environment and thus are not encouraged. To effectively control plant disease and to reduce the harm to environment, it is crucial to understand disease transmission dynamics.

The prevention and control of plant infectious diseases is of vital importance in agricultural production [5] . To work with the plant disease, we first understand how they spread. There are many ways that plant viruses interact with the vectors; this transmission works in the following way. The vectors consume sap from an infected host through their stylet. When the infected vector contacts a healthy plant, some virus particles leave the vector and invade the plant [6] . So the vector-borne is a very important part of the transmission of plant diseases.

Treatment plays a very important role to control the spread of diseases. In recent years, many researchers [7] [8] [9] [10] have studied some mathematical models incorporating treatment. For example, in [11] , Wang and Ruan studied an epidemic model, and provide the limited resources for the treatment of patients. In [12] , Wang proposed constant treatment, which simulates a limited capacity for treatment. In practice, when the number of infectives is large, the constant treatment is suitable for hypothesis of model. Recently, discontinuous treatment strategies are proposed by [13] [14] . The results show that discontinuous treatment strategies would be accord with real condition. Applying this discontinuous treatment strategy makes the mathematical model a discontinuous system. At the same time, some non-smooth analysis techniques [15] are used for this system.

In [16] , Shi and Zhao presented a vector-borne plant disease model, but they do not studied treatment to the infected plant host. Treating infected plant is a quite effective method which to control the outbreak of the plant disease. Although continuous treatment is an effective method, the outbreak of the plant disease is periodic, and continuing treatment can be a huge waste of resources. In order to be realistic, we built a vector-borne plant disease model with discontinuous treatment.

The paper is organized as follows. In the next section, we will construct the model and introduce the rational assumptions for model. In Section 3, positivity of the solution for the model will be clearly discussed. We obtain the existence of possible equilibria, the basic reproductive number, and the stability of equilibria in Section 4. In Section 5 and Section 6, we summarize our main results and main results are numerically simulated.

2. Model and Preliminaries

To construct the model, the following assumptions are being made by Shi et al. in [16] .

(A1) The total of the insect vector population is divided into X and Y, which denotes the densities of the susceptible vector and infective vector at time t, respectively. The total of the plant host population is divided into S, I, and R, which represents the numbers of the susceptible, infective, and recovered host plant population at time t, respectively. At the same time, we assume that the number of plants in one area is fixed. The total number of plants K = S + I + R is a positive constant. In fact, when a plant has died, it would be replaced by a new plant to keep the total number of plants. Further, we assume that those new plants are susceptible, i.e., we chose the birth rate of susceptible plant host as f ( S , I ) = μ K + d I .

(A2) The susceptible plants can be infected not only by the infected insect vectors but also by the infected plants.

(A3) A susceptible vector can be infected only by an infected plant host, and after it is infected, it will hold the virus for the rest of its life. Further, there is no vertical infection being considered.

(A4) The replenishment rate of insect vectors is a positive constant, and all of the new born vectors are susceptible.

According to the principle of the compartmental model, consider the following model with discontinuous treatment:

{ S ˙ = f ( S , I ) μ S ( β P Y + β s I ) S , I ˙ = ( β P Y + β s I ) S ( d + μ + γ ) I h ( I ) , R ˙ = γ I + h ( I ) μ R , X ˙ = Λ β 1 I X m X , Y ˙ = β 1 I X m Y . (2.1)

Here the dimensionless variables and parameters (with parameter values) are given in Table 1.

The function h ( I ) = φ ( I ) I represents the treatment rate. φ ( I ) satisfies the following assumptions. Obviously, the treatment rate should be nondecreasing as the number of infectious individuals is increasing. The following assumption will be needed throughout the paper.

(H1) φ : [ 0, ) [ 0, ) is nondecreasing and has at most a finite number of jump discontinuities in every compact interval. No loss of generality, we always assume that φ is continuous at I = 0 , otherwise we define φ ( 0 ) to be φ ( 0 + ) . Here φ ( 0 + ) denotes the right limit of φ ( I ) as I 0 + .

By adding the fourth and fifth equations of system (2.1), we get

N ˙ = Λ m N (2.2)

where N = X + Y . From Equation (2.2), we easily get N Λ m as t .

Note that S + I + R = K . Since the variable R and X does not appear in the first two equations of model (2.1), meanwhile, let’s substitute X for ( Λ m Y ) in

the fifth equation. We only need to study the first two equations and the fifth

Table 1. Dimensionless variables and parameters (with illustrative parameter values) in system (2.1).

equation of model (2.1), thereby lowering the order of the system to be studied, i.e.

{ S ˙ = μ ( K S ) ( β P Y + β s I ) S + d I , I ˙ = ( β P Y + β s I ) S ω I h ( I ) , Y ˙ = Λ β 1 I m ( β 1 I + m ) Y . (2.3)

where ω = d + μ + γ . Obviously,

Ω = { ( S , I , Y ) R + 3 : 0 S + I K ,0 Y Λ m } (2.4)

is the positively invariant set for system (2.3).

According to the definition of solutions for differential equations with discontinuous right-hand sides in [15] [17] , ( S ( t ) , I ( t ) , Y ( t ) ) is called a solution with initial condition

( S ( 0 ) , I ( 0 ) , Y ( 0 ) ) = ( S 0 , I 0 , Y 0 ) , S 0 , I 0 , Y 0 0 (2.5)

of model (2.3) on [ 0 , T ) , 0 < T , if it is absolutely continuous on any compact subinterval of [ 0, T ) , and almost everywhere on [ 0, T ) (abbreviated to a.e. on [ 0, T ) ) satisfies the following differential inclusion:

{ S ˙ = μ ( K S ) ( β P Y + β s I ) S + d I , I ˙ ( β P Y + β s I ) S ω I c o ¯ [ h ( I ) ] , Y ˙ = Λ β 1 I m ( β 1 I + m ) Y . (2.6)

where c o ¯ [ h ( I ) ] = [ h ( I 0 ) , h ( I + 0 ) ] . Here, h ( I 0 ) and h ( I + 0 ) denote the left limit and the right limit of the function h ( I ) at I, respectively.

From (H1), it is clear that the set map

( S , I , Y ) ( μ ( K S ) ( β P Y + β s I ) S + d I , ( β P Y + β s I ) S ω I c o ¯ [ h ( I ) ] , Λ β 1 I m ( β 1 I + m ) Y ) (2.7)

is an upper semi-continuous set-valued map with non-empty compact convex values. By the measurable selection theorem [15] , if ( S ( t ) , I ( t ) , Y ( t ) ) is a solution of model (2.3) on [ 0, T ) , then there is a measurable function m ( t ) c o ¯ [ h ( I ( t ) ) ] such that

{ S ˙ = μ ( K S ) ( β P Y + β s I ) S + d I , I ˙ ( β P Y + β s I ) S ω I m ( t ) , a . e . on [ 0 , T ) . Y ˙ = Λ β 1 I m ( β 1 I + m ) Y . (2.8)

3. Positivity

In this section, we will prove the positive of the solution to the initial condition of the model (2.3) with positive initial value. First, we will prove the following theorem.

Theorem 3.1. Suppose that assumption (H1) holds and let ( S ( t ) , I ( t ) , Y ( t ) ) be the solution with initial condition (2.5) of model (2.3) on [ 0, T ) . Then ( S ( t ) , I ( t ) , Y ( t ) ) is nonnegative on [ 0, T ) .

Proof: By the definition of a solution of (2.3) in the sense of Filippov, ( S ( t ) , I ( t ) , Y ( t ) ) must be a solution to differential inclusion (2.6). From the first equation of (2.6), we have

[ S 0 + 0 t ( μ k + d I ( u ) ) exp ( 0 u ( μ + β P Y ( ρ ) + β s I ( ρ ) ) d ρ ) d u ] exp ( ( 0 t ( μ + β P Y ( ρ ) + β s I ( ρ ) ) d ρ ) ) > 0 (3.1)

for all t ( 0, T ) .

Based on the previous hypothesis of (H1), we have c o ¯ [ h ( 0 ) ] = 0 and h ( I ) is continuous at I = 0 . Combining the continuity of φ at I = 0 , it may be concluded that there exists a positive constant δ such that φ ( I ) is continuous as | I | < δ . On this account, when | I | < δ the differential inclusion (2.6) becomes the following system of differential equations:

{ I ˙ = ( β P Y + β s I ) S ( ω + φ ( I ) ) I , Y ˙ = Λ β 1 I m ( β 1 I + m ) Y . (3.2)

We divide this into four cases to discuss the positivity of the solutions for (2.6).

1) I 0 = Y 0 = 0 .

From (3.2), we see that I ( t ) = Y ( t ) = 0 for all t [ 0, T ) .

2) I 0 > 0 , Y 0 = 0 .

By the continuity of I ( t ) at t = 0 and d Y d t | t = 0 = Λ β 1 I 0 m > 0 , we conclude I ( t ) > 0 and Y ( t ) > 0 for all t ( 0, T ) . If it is not true, then we can set

t 1 = inf { t : I ( t ) = 0 or Y ( t ) = 0 } ( 0 , T ) . (3.3)

If I ( t 1 ) = 0 , then from d I d t ( ω + φ ( I ) ) I for 0 t t 1 , we have I ( t 1 ) I 0 exp ( ( ω + φ ( I ) ) t 1 ) > 0 . This is a contradiction.

If I ( t 1 ) = 0 , then there is a θ such that t 1 θ > 0 and 0 < I ( t ) < δ on [ t θ , t 1 ) . Therefore, the second equation of (3.2) implies

d Y d t ( β 1 I + m ) Y (3.4)

We have

Y ( t 1 ) Y ( t 1 θ ) exp ( t 1 θ t 1 ( β 1 I ( ξ ) + m ) d ξ ) > 0 (3.5)

This is also a contradiction. Hence, I ( t ) and Y ( t ) are positive for all t ( 0, T ) . The same conclusion can be reached for the following two cases.

3) I 0 = 0 , Y 0 > 0 .

4) I 0 > 0 , Y 0 > 0 . This completes the proof.

4. The Equilibria and Their Stability

In this section, we will discuss the existence of equilibria of system (2.3). First, we prove the existence of endemic equilibrium.

Let ( S ( t ) , I ( t ) , Y ( t ) ) = ( S * , I * , Y * ) is a constant solution of (2.3), where ( S * , I * , Y * ) satisfies the following system:

{ 0 = μ ( K S * ) ( β P Y * + β s I * ) S * + d I * , 0 ( β P Y * + β s I * ) S * ω I * c o ¯ [ h ( I * ) ] , 0 = β 1 I * Λ m ( β 1 I * + m ) Y * . (4.1)

Since h ( 0 ) = 0 , there always exists a disease-free equilibrium P 0 of the model (2.3), where P 0 = ( K , 0 , 0 ) . Next, we consider that the existence of an endemic equilibrium of the model (2.3).

It follows from the first and third equations of (4.1), we conclude that

S * = d I * + μ K μ + β P Y * + β s I * , Y * = Λ β 1 I * m ( β 1 I * + m ) . (4.2)

Substituting (4.2) into the second inclusion of (4.1), we have the follows

A 1 I * 2 + B 1 I * + C 1 A 2 I * 2 + B 2 I * + C 2 ω c o ¯ [ φ ( I ) ] = [ φ ( I * 0 ) , φ ( I * + 0 ) ] (4.3)

where

A 1 = m d β 1 β s

B 1 = d Λ β p β 1 + m 2 d β s + μ m K β 1 β s ,

C 1 = μ K ( Λ β p β 1 + m 2 β s ) ,

A 2 = m β 1 β s ,

B 2 = μ m β 1 + Λ β 1 β p + m 2 β s ,

C 2 = μ m 2 .

Denote

g ( I * ) = A 1 I * 2 + B 1 I * + C 1 A 2 I * 2 + B 2 I * + C 2 ω (4.4)

and let

R 0 = K ( Λ β p β 1 + β s m 2 ) m 2 ( ω + φ ( 0 ) ) (4.5)

We next claim that R 0 is the basic reproductive number for the model (2.3) which will determine the existence of an endemic equilibrium.

Theorem 4.1. Suppose that assumption (H1) holds. If R 0 1 , then there only exists a disease-free equilibrium P 0 ( K ,0,0 ) . If R 0 > 1 , then there exists a unique positive endemic equilibrium P * ( S * , E * , I * ) except P 0 .

Proof: By R 0 1 , we get g ( 0 ) φ ( 0 ) . Since g ( I ) is nonincreasing on I and φ ( I ) is nondecreasing on I. For this reason, the inclusion (4.3) is only valid at I = 0 . Hence, the model (2.3) has a unique disease-free equilibrium as long as R 0 1 .

From (4.4), we have the following

A I 2 + B I + C = 0 (4.6)

where

A = m β 1 β s ( m d ω ) < 0 ,

B = Λ d β 1 β p + m 2 d β s + μ m K β 1 β s ω m μ β 1 Λ ω β 1 β p ω m 2 β s ,

C = μ ω m 2 ( R 0 1 ) + μ m 2 φ ( 0 ) .

If R 0 1 , then C > 0 , and the Equation (4.6) has a unique positive root I, where

I = B + Δ 1 2 2 A , Δ = B 2 4 A C . (4.7)

If R 0 > 1 , then. Meanwhile, the inequality, it implies. Therefore, the set

(4.8)

is bounded and non-empty. We can write

It follows easily that

(4.9)

We claim. Assumption, contrary to our claim, that

From (H1), there exists a such that

(4.10)

This contradicts the definition of. Thus, we have . That is to say, is a positive solution of the inclusion (4.3). We proceed to show that is the only one positive solution of the inclusion (4.3). If the inclusion (4.3) has another positive solution, then there must exist two numbers

(4.11)

which satisfy

(4.12)

and

(4.13)

Subtracting (4.13) from (4.12) gives

(4.14)

which implies

(4.15)

This is a contradiction. Hence, is the unique positive solution of the inclusion (4.3). Combining it with (4.2), we conclude that the is the unique endemic equilibrium of (2.3). The proof is completed.

Next, we prove the global stability of the disease-free equilibrium and the endemic equilibrium. We do this in several steps. We first investigate the local properties of the equilibria of system (2.3).

Theorem 4.2. Assume (H1) holds. The disease-free equilibrium is locally asymptotically stable if, and is unstable if.

Proof: We analyze the stability of the disease-free equilibrium by investigating the eigenvalues of the Jacobian matrix of model (2.3) at. The matrix is

(4.16)

Thus, the characteristic equation at the disease-free equilibrium is

(4.17)

It is easy to see that one of the roots with respect to of (4.14) is. the other two roots are determined by the following characteristic equation of:

(4.18)

From (4.18) and Routh-Hurwitz criteria [18] , it is easily seen that both the real parts of and of are negative when. When, one of and is a number with a positive real part. Thus the disease-free equilibrium is locally asymptotically stable if and unstable if.,

We have shown that there exists a positive endemic equilibrium if and only if in Theorem 4.1. Here, we will establish its local stability.

Theorem 4.3. Suppose that assumption (H1) holds. If, the endemic equilibrium of the system (2.3) is locally asymptotically stable.

Proof: The Jacobian matrix of (2.3) at the endemic equilibrium is

(4.19)

Replacing, by, , respectively. So we have

(4.20)

The characteristic equation of is

(4.21)

where

Since is nondecreasing,. This implies.

Then

Hence, all of the Routh-Hurwitz criteria are satisfied. Thus it follows that the endemic equilibrium of (2.3), which exists if, is always locally asymptotically stable. The proof is completed.,

Next, we will prove global stability of the disease-free equilibrium and endemic equilibrium of (2.3). We need to use the LaSalle-type invariance principle for the differential inclusion (Theorem 3 in [19] ) to prove their global stability.

Let. We obtain the following system analogous to (2.5)

(4.22)

Set

(4.23)

and

For any, there exists an such that

(4.24)

Hence,

(4.25)

When, the nondecreasing of implies

(4.26)

It shows that is a Lyapunov function of (4.23).

Furthermore, when, we have

(4.27)

When, we have, which implies. For any, we set

(4.28)

Hence, the largest weakly invariant subset of is the singleton.

When, we have

(4.29)

From the first equation of (4.23) and x = 0, it may be concluded that. Therefore, we see that the largest weakly invariant subset of is also the singleton. By the LaSalle-type invariance principle, the equilibrium of (4.23) is globally asymptotically stable as. Summarizing the above analysis, we obtain the following theorem.

Next, we demonstrate the global stability of the endemic equilibrium of (2.3). So, we have the following theorem.

Theorem 4.4. Suppose that assumption (H1) holds. If, the disease-free equilibrium of the system (2.3) is globally asymptotically stable.

Proof: Let

(4.30)

Write

and

(4.31)

For any, there exists an such that

(4.32)

Hence

(4.33)

The monotonicity of implies. Thus. This shows that V is a Lyapunov function of (2.24). Define

(4.34)

If, then the first equation of (2.6) implies. Consequently, for any, the largest weakly invariant subset of of (2.6) is the singleton overline. Here

Therefore, is globally asymptotically stable if. This completes the proof.,

5. Numerical Simulation

To make our analysis more intuitive, some numerical simulations of solutions of the model (2.6) is provided which to illustrate the influence of insect vector and discontinuous treatment on the spread of plant disease. We apparent a treatment function satisfying (H1) as follows:

(5.1)

where. The treatment function is applied at the following case: when the infective individuals I attain some threshold, the treatment rate should be strengthened to control the spread of the plant disease.

To better illustrate the effects of non-continuous healing on the spread of plant disease, the following parameters are derived from [16] . Let, , , , , , , , , , then we easily calculate by using (4.5). Figure 1 shows that the infective individuals I and the infective insect vectors Y tend to 0, and it means that the disease goes to extinction. In addition, we find that the peak values of the infective is affected by the different values of. Figure 1 shows that larger values of can reduce the peak values of the infective. Therefore, we can increase the treatment rate to prevent the spread of disease after the number infective individuals reaching some high level. From Figure 1, the infective individuals reach some level and strengthening the treatment rate is also effective for disease control, even though we do not take any treatment measures at the initial time of the diseases outbreak.

If we fixed all parameter values as follows:, , , , , , , , one could easily see that, by using (4.5). Figure 2 shows that the endemic equilibrium is globally asymptotically stable. By (4.5), it’s obvious find that the basic reproductive number is independent of, and the different values of can affect the stability level of the infective. That is to say, larger values of

Figure 1. Stability of disease-free equilibrium. The parameters are fixed as follows:, , , , , , , , , , and the initial values are.The time series charts for and the phase diagram are given in (a), (b), (c), and (d), respectively.

can impact level of the infective. It implies that the strengthening of the treatment rate can effectively control the spread of plant disease after the number of infective individuals has increased to some high level.

6. Discussion

As for the plant infectious disease model, our main object is to investigate the effect of the insect vector and discontinuous treatment function on the dynamics of spreading the plant disease. We calculated the basic reproduction number, which is derived under some reasonable assumptions on the discontinuous treatment function. It is an important threshold parameter which plays an important role in determining the global dynamics of the model (2.6) and whether it persists or dies out of the disease. When, the disease-free equilibrium is globally stable, which means that the disease always dies out, and when, the plant disease will be permanent which means that after some period of time the plant disease will become endemic and it is global stable.

In this paper, we studied the existence, local stability and global stability of the disease-free equilibrium and endemic equilibrium of the system (2.3) in detail. By building a suitable Lyapunov function, and the Jacobian matrix method,

Figure 2. Stability of endemic equilibrium. The parameters are fixed as follows:, , , , , , , ,. The time series charts for S(t), I(t), Y(t) and the phase diagram are given in (a), (b), (c), and (d), respectively.

employing Routh-Hurwitz criteria and LaSalle-type invariance principle, the main results as shown in Theorems 4.2, 4.3 and 4.4 have been derived. Our main results indicate that if, then the disease-free equilibrium is globally asymptotically stable, if, the unique endemic equilibrium is globally asymptotically stable of the system (2.3). From above results, it is easy to find that the basic reproduction number plays an important role in determining the persistence or dying out of the disease.

Acknowledgements

The research have been supported by The Natural Science Foundation of China (11561004), the Science and Technology research project of Jiangxi Provincial Education Department (171373, 171374, GJJ170815), The bidding project of Gannan Normal University (16zb02).

Cite this paper

Lv, H.M., Fei, L.Z., Yuan, Z. and Zhang, F.M. (2018) Global Dynamic Analysis of a Vector-Borne Plant Disease Model with Discontinuous Treatment. Applied Mathematics, 9, 496-511. https://doi.org/10.4236/am.2018.95036

References

  1. 1. Lee, J.A., Halbert, S.E., Dawson, W.O., et al. (2015) Asymptomatic Spread of Huanglongbing and Implications for Disease Control. Proceedings of the National Academy of Sciences, 112, 7605-7610. https://doi.org/10.1073/pnas.1508253112

  2. 2. Smith, A.B., Beeck, C.P., Cowling, W.A., et al. (2013) A Bivariate Mixed Model Approach for the Analysis of Plant Survival Data. Euphytica, 190, 371-383. https://doi.org/10.1007/s10681-012-0791-0

  3. 3. Xia, L., Gao, S., Zou, Q., et al. (2013) Analysis of a Nonautonomous Plant Disease Model with Latent Period. Applied Mathematics and Computation, 223, 147-159. https://doi.org/10.1016/j.amc.2013.08.011

  4. 4. Wang, J., Zhang, F., Wang, L., et al. (2016) Equilibrium, Pseudoequilibrium and Sliding-Mode Heteroclinic Orbit in a Filippov-Type Plant Disease Model. Nonlinear Analysis: Real World Applications, 31, 308-324. https://doi.org/10.1016/j.nonrwa.2016.01.017

  5. 5. Zhao, W., Li, J., Zhang, T., et al. (2017) Persistence and Ergodicity of Plant Disease Model with Markov Conversion and Impulsive Toxicant Input. Communications in Nonlinear Science and Numerical Simulation, 48, 70-84. https://doi.org/10.1016/j.cnsns.2016.12.020

  6. 6. Jackson, M. and Chen-Charpentier, B.M. (2016) Modeling Plant Virus Propagation with Delays. Journal of Computational and Applied Mathematics, 2016, 1-14.

  7. 7. Xue, Y. and Wang, J. (2012) Backward Bifurcation of an Epidemic Model with Infectious Force in Infected and Immune Period and Treatment. Abstract and Applied Analysis, 2012, Article ID: 647853. https://doi.org/10.1155/2012/647853

  8. 8. Hussaini, N. and Winter, M. (2010) Travelling Waves for an Epidemic Model with Non-Smooth Treatment Rates. Journal of Statistical Mechanics, 11, Article ID: 11019. https://doi.org/10.1088/1742-5468/2010/11/P11019

  9. 9. Brauer, F. (2008) Epidemic Models with Heterogeneous Mixing and Treatment. Bulletin of Mathematical Biology, 70, 1869-1885. https://doi.org/10.1007/s11538-008-9326-1

  10. 10. Hu, Z., Liu, S. and Wang, H. (2008) Backward Bifurcation of an Epidemic Model with Standard Incidence Rate and Treatment Rate. Nonlinear Analysis: Real World Applications, 9, 2302-2312. https://doi.org/10.1016/j.nonrwa.2007.08.009

  11. 11. Wang, W. and Ruan, S. (2004) Bifurcations in an Epidemic Model with Constant Removal Rate of the Infectives. Journal of Mathematical Analysis and Applications, 291, 775-793. https://doi.org/10.1016/j.jmaa.2003.11.043

  12. 12. Wang, W. (2006) Backward Bifurcation of an Epidemic Model with Treatment. Mathematical Biosciences, 201, 58-71. https://doi.org/10.1016/j.mbs.2005.12.022

  13. 13. Guo, Z. Huang, L. and Zou, X. (2012) Impact of Discontinuous Treatments on Disease Dynamics in an SIR Epidemic Model. Mathematical Biosciences and Engineering, 9, 97-110. https://doi.org/10.3934/mbe.2012.9.97

  14. 14. Zhang, T., Kang, R., Wang, K. and Liu, J. (2015) Global Dynamics of an SEIR Epidemic Model with Discontinuous Treatment. Advances in Difference Equations, 361, 1-16.

  15. 15. Huang, L., Guo, Z. and Wang, J. (2011) Theory and Applications of Differential Equations with Discontinuous Righthand Sides. Science Press, Beijing. (In Chinese)

  16. 16. Shi, R., Zhao, H. and Tang, S. (2014) Global Dynamic Analysis of a Vector-Borne Plant Disease Model. Advances in Difference Equations, 59, 1-16. https://doi.org/10.1186/1687-1847-2014-59

  17. 17. Filippov, A.F. (1988) Differential Equations with Discontinuous Righthand Sides. Mathematics and Its Applications (Soviet Series). Kluwer Academic, Boston. https://doi.org/10.1007/978-94-015-7793-9

  18. 18. Ma, Z. and Zhou, C. (2016) Methods of Qualitative and Stability of Ordinary Differential Equations. Science Press, Beijing. (In Chinese)

  19. 19. Baciotti, A. and Ceragioli, F. (1999) Stability and Stabilization of Discontinuous Systems and Non-Smooth Lyapunov Function. ESAIM: Control, Optimisation and Calculus of Variations, 4, 361-376. https://doi.org/10.1051/cocv:1999113