Research Article  Open Access
Ruiqing Shi, Ting Lu, Cuihong Wang, "Dynamic Analysis of a FractionalOrder Model for Hepatitis B Virus with Holling II Functional Response", Complexity, vol. 2019, Article ID 1097201, 13 pages, 2019. https://doi.org/10.1155/2019/1097201
Dynamic Analysis of a FractionalOrder Model for Hepatitis B Virus with Holling II Functional Response
Abstract
In this paper, a fractionalorder model is constructed to describe the transmission of Hepatitis B Virus (HBV). Firstly, the existence and uniqueness of positive solutions are proved. Secondly, the basic reproduction number and the sufficient conditions for the existence of two equilibriums are obtained. Thirdly, the stability of equilibriums are analyzed. After that, some numerical simulations are performed to verify the theoretical prediction. Finally, a brief discussion is presented.
1. Introduction
In recent years, more and more attention has been paid to the research of epidemic disease, and it is a great challenge to control the spread of epidemic disease among people. Mathematical models of infectious disease are important tools, which provided theoretical basis for the prevention and control of disease[1–6]. Therefore, many scholars have established mathematical models to simulate the epidemic dynamics of virus in host cells. Nowak and Bangham presented a basic model of the dynamics of HIV infection[7, 8]. This model have since been used for other viral infections, like HBV [9–11] and HCV [12, 13]. The basic model can be formulated as a system of three differential equations:where , , and are the densities of uninfected hepatocytes, infected hepatocytes and free virions, respectively. Uninfected hepatocytes are assumed to be produced by a constant rate , and die at the per capita rate of . It becomes to infected hepatocytes at the rate of , and is the infection rate constant which characterize the infection efficiency. The death rate of infected hepatocytes is . Free virion are assumed to be produced from infected hepatocytes at the rate of , and is the clearance rate of viral particles. All parameters , , , , , and are positive constants.
Mathematical analysis is an effective method to study the state of virus in vivo. According to the World Health Organization (WHO), over two billion people has been or is actively infected by HBV infections [14]. Chronic HBV infection often occurs in our early life and the virus persists in the body due to the destruction of strong antibody or cellular immune responses [15]. Recently, in an HBV model of viral infection, a reversion rate constant has been introduced into the uninfected state to reflect a noncytolytic mechanism [16]. That is to say, the infected hepatocytes may be reverted to be uninfected state by losing all cccDNA from their nucleus [10]. Some research indicates that there is indeed a cytokinemediated cure of infected cells during HBV infection [10, 17–19]. Specifically, in [20] the authors considered the following model:where the term represents the rate at which infected hepatocytes are recovered to uninfected hepatocytes through cure.
As is known to all, differential equations provide an important tool for mathematical models, which are very useful in understanding the dynamic behavior of biological systems. However, many engineering, physical and biological systems have temporal memory [21–23]. Classical integerorder differential equations do not reflect this characteristic of memory. Therefore, it is necessary to introduce fractionalorder differential equations since fractionalorder derivative can provide a useful instrument for description of memory and hereditary properties in many different domains. What’s more, it was showed that fractionalorder differential equations could simulate many phenomena that integerorder can not [24, 25]. In fact, fractional calculus generalizes integrals and derivatives of integerorder to arbitrary orders [26–28]. It has been widely used in many different fields, such as finance [21], physics [22], biology [23], medicine [29] and so on.
Although bilinear incidence rate is mostly used in epidemic models, recently many literatures considered Holling typeII functional response as the incidence rate [2, 30–32]. In addition, our research is mainly to control the development of HBV. In fact, interferon (IFN) and lamivudine (LAM) are mainly used for treatment of HBV. They both have the effects of antivirus and antifibrosis. IFN is the most commonly used antiviral drug, it can inhibit virus’ replication with extensive, indirect and speciesspecific. LAM is a cytosine nucleoside analogue that inhibits viral reverse transcriptase.
However, few literatures use fractional differential equations to describe the transmission of HBV with drug treatment. Therefore, based on the model (2), combining the fractionalorder derivatives and the Holling typeII functional response, we propose the following improved HBV model with drug treatment:where and it represents the drug effect on HBV by IFN. and it represents the drug effect on HBV by LAM. The biological meanings of state variables and parameters are shown in Table 1.

The organization of our current paper is as follows: In Section 2, we give some properties of fractionalorder calculus. In Section 3, firstly, the existence and uniqueness of the positive solutions are proved; secondly, the sufficient conditions for the existence and the stability of two equilibriums for system (3) are obtained. In Section 4, some numerical simulations are performed to verify our theoretical results. Finally, this paper ends up with a brief discussion.
2. Basic Properties of FractionalOrder Calculus
In fractionalorder calculus, we know there are many fractionalorder integration and fractionalorder differentiation that have been defined, for example, the GrunwaldLetnikov (GL) definition, the RiemannLiouville (RL) definition and the Caputo definition. Since the initial conditions are in the same form as for the integerorder differential equations, we will adapt the Caputo’s definition in our paper. Another advantage of this definition is that applied problems require definitions of fractional derivatives, where there are clear interpretations of initial conditions, which contain , , , etc.
Definition 1 (see [26]). The Caputo fractionalorder derivative of a continuous function : can be defined aswhere , .
Lemma 2 (see [33], (generalized mean value theorem)). Assume that , are continuous functions on the interval , then we havewhere .
Lemma 3 (see [33]). Assume that , are continuous functions on the interval , and , then we have(i)If , for , then is nondecreasing for each .(ii)If , for , then is nonincreasing for each .
Lemma 4 (see [34]). Suppose that the vector function : satisfies the following conditions:(i)Function is Lebesgue measurable with respect to ;(ii)Function is continuous with respect to on ;(iii) is continuous with respect to on ;(iv), , , where , are two positive constants.Then the initial value problems (IVP) has a unique solution.
3. Analysis of System (3)
In this section, firstly, the existence and uniqueness of positive solution is proved; secondly, the basic reproduction number and the existence conditions for both equilibriums (diseasefree equilibrium and endemic equilibrium) are obtained; lastly, the conditions for the stability of both equilibriums are obtained.
3.1. The Existence and Uniqueness of Positive Solutions
Theorem 5. System (3) with any positive initial value has a unique solution, and it remains within , where
Proof. We will accomplish the proof through three steps:
Step 1. In this step, we will show that system (3) with any positive initial value has a unique solution.
Denote the right side of system (3) as vector function , then the corresponding conditions (i)(iii) of Lemma 4 are satisfied. So, we only need to prove that system (3) satisfies the fourth condition of Lemma 4.
Let , , , , , , andthen system (3) can be reduced to the following form:Denote and thenAccording to Lemma 4, system (3) with any positive initial values has a unique solution.
Step 2. We will prove that the solution of system (3) with positive initial values is always nonnegative. From system (3), we easily getFrom Lemma 3, we have , , for any . As a result, the solution of system (3) will remain in .
Step 3. Assume that for the biological justification. Adding the first two equations of system (3), we have which implies that Since for any , then we haveprovided that .
From the last equation of system (3) we getand from which we get that as time large enough.
From the results of step 2 and step 3, we know that is a positive invariant set with respect to system (3).
Since is positive invariant with respect to system (3), we only need to consider this system within in the rest of this section.
3.2. Basic Reproduction Number and the Existence of Both Equilibriums
In the process of studying disease transmission, what we most concerned about is how to eradicate or control the disease within a range. Many epidemiological models have a diseasefree equilibrium where the disease is not present in the population. Therefore, we need a threshold to estimate whether the disease is present or not in the population. This threshold is called basic reproduction number. By the next generation matrix approach given in [35], we can get the basic reproduction number of system (3) as follows: where
Corollary 6. When , we obtain . Here is called basic reproduction number, and it is a threshold to estimate whether the disease is persistent or not.
In addition, through the eigenvalue analysis method, we can easily prove that the diseasefree equilibrium is stable when ; and it is unstable when .
Let the right side of system (3) equal to zero, we obtain an algebraic equations as follows
By simple calculation, we obtain two equilibriums of the system (3), namely:
There always exists a diseasefree equilibrium ;
If , then there exists a unique endemic equilibrium , where
3.3. Stability Analysis of the Two Equilibriums
Theorem 7. The diseasefree equilibrium is locally asymptotically stable if , and it is unstable if the inequality is reversed.
Proof. The Jacobian matrix evaluated at the diseasefree equilibrium is given by and the corresponding characteristic equation at the diseasefree equilibrium isIt is easily seen that one of the eigenvalues of (22) is , and the other two eigenvalues are determined by the following quadratic equation:whereIf , then and both roots of (23) have negative real part. Thus, the diseasefree equilibrium is locally asymptotically stable.
If , then , which means that one root of (23) is positive. Thus, the diseasefree equilibrium is unstable.
Remark 8. When , if the parameters satisfy , then the diseasefree equilibrium is locally asymptotically stable, which means no drug treatments are employed.
Remark 9. If the parameters satisfy , then the diseasefree equilibrium will be unstable if there is no drug treatment. We can employ suitable drug treatment to control the disease, specifically we have the following:
If , , which means only LAM is employed to treat HBV, then we can select suitable such that . According to Theorem 7, is locally asymptotically stable.
If , , which means only IFN is employed to treat HBV, then we can select suitable such that . According to Theorem 7, is locally asymptotically stable.
If and , which means two drugs are used to treat HBV, then we can select suitable and such that . According to Theorem 7, is locally asymptotically stable.
To discuss the local stability of the endemic equilibrium , we consider the linearized system of (3) at . The Jacobian matrix at the endemic equilibrium is given by:
The corresponding characteristic equation at the endemic equilibrium iswhere
Proposition 10. The endemic equilibrium is locally asymptotic stable if all of the eigenvalues of Eq.(26) satisfy , .
Denote as the discriminant of , where
In order to discuss the stability of the endemic equilibrium , we get the following result by use of the same method as in [36].
Proposition 11. If the discriminant , then the endemic equilibrium is locally asymptotically stable if and only if the RouthHurwitz conditions are satisfied; that is to say, , , .
If the discriminant , , , , and , then the endemic equilibrium is locally asymptotically stable.
If the discriminant , , , , , then the endemic equilibrium is locally asymptotically stable.
If the discriminant , , , , then the endemic equilibrium is unstable.
About the global stability of the diseasefree equilibrium , we have the following result.
Theorem 12. The diseasefree equilibrium is globally asymptotically stable if .
Proof. Consider the following Lyapunov functionThe derivative of along a solution of system (3) is:and the last inequality is obtained since .
If , then we have . In addition, we know that the maximum invariant set for is singleton . According to the LaSalle’s invariance principle, we know that all solutions in converge to .
This completes the proof.
4. Numerical Simulations for System (3)
In previous sections, we have predicted some dynamical behavior of system (3). In this section, we will explore some numerical simulations to verify our theoretical results. The values of parameters are given in Table 1, and most of the values are taken from [37]. Since we are only concerned with the number of infected hepatocytes and the concentration of free virions, only and are shown in our illustration.
Figure 1 shows that when , if the value of is relatively big, then , and the diseasefree equilibrium is locally asymptotically stable.
(a) u1=u2=0, α=0.9
(b) u1=u2=0, α=0.90
(c) u1=u2=0, α=0.90
(d) u1=u2=0, α=0.90
Figure 2 shows that when , if the value of is relatively small, then , and the diseasefree equilibrium is unstable.
(a) u1=u2=0, α=0.95
(b) u1=u2=0, α=0.95
(c) u1=u2=0, α=0.95
(d) u1=u2=0, α=0.95
Remark 13. From Figures 1 and 2, we find that the value of cure rate is very important for the dynamics of the system. If the value of is relatively small, the disease will persist; while the disease will eradicate if the value of is relatively big. In addition, we also find that the initial values are not sensitive to the dynamical behaviors.
Figure 3 shows that if , as long as , then the diseasefree equilibrium is stable.
(a) u1=0, u2=5, α=0.95
(b) u1=0, u2=5, α=0.95
(c) u1=0, u2=5, α=0.95
(d) u1=0, u2=5, α=0.95
Figure 4 shows that if , as long as , then the diseasefree equilibrium is stable.
(a) u1=0.7, u2=0, α=0.95
(b) u1=0.7, u2=0, α=0.95
(c) u1=0.7, u2=0, α=0.95
(d) u1=0.7, u2=0, α=0.95
Figure 5 shows that if , as long as , then the diseasefree equilibrium is stable.
(a) u1=0.2, u2=5, α=0.95
(b) u1=0.2, u2=5, α=0.95
(c) u1=0.2, u2=5, α=0.95
(d) u1=0.2, u2=5, α=0.95
Remark 14. From these three figures, we found that whatever the initial values are, if we choose suitable values of and/or such that , then the diseasefree equilibrium is stable. The numerical simulation results are consistent with the Theorem 7 and Remark 9. That is to say, is a threshold to determine whether the disease is persistent or not.
From (a)(b) of Figure 6 we find that the diseasefree equilibrium is always asymptotically stable. At this case, only LAM (without IFN) is used to treat HBV.
(a) u1=0, u2=5, =0.110
(b) u1=0, u2=5, =0.110
(c) u1=0.7, u2=0, =0.768
(d) u1=0.7, u2=0, =0.768
(e) u1=0.2, u2=5, =0.088
(f) u1=0.2, u2=5, =0.088
From (c)(d) of Figure 6 we find that the diseasefree equilibrium is always asymptotically stable. At this case, only IFN (without LAM) is used to control HBV.
From (e)(f) of Figure 6 we find that the diseasefree equilibrium is always asymptotically stable. At this case, both IFN and LAM are used to control HBV.
Remark 15. From Figure 6, we can see that no matter one or two control measures are taken, the diseasefree equilibrium is always stable, as long as the values of and satisfy .
From Figures 6(a) and 6(e), we can see that as the value of increases, the peak of disease infection decreases; from Figure 6(c), we can see that the peak of disease infection increases as the value of increases. These simulation results indicate that the value of is an important factor for the peak of disease infection.
From Figures 7(a) and 7(b) we find that when , and the endemic equilibrium is always stable for three different initial values.
(a) α=1,
(b) α=1,
(c) α=0.5,
(d) α=0.5,
From Figures 7(c) and 7(d) we find that when , although , the endemic equilibrium is unstable for three different initial values.
Remark 16. From Figure 7, we find that under the condition of , the endemic equilibrium may be stable or unstable for different values of .
5. Discussion
In this article, we construct a fractionalorder model with Holling II functional response to describe the transmission of HBV. We do the following work:
We proved the existence and uniqueness of positive solutions.
The basic reproduction number and the sufficient conditions for the existence of two equilibriums are obtained.
We analyzed the stability of equilibriums are analyzed.
Some numerical simulations are performed to verify our theoretical prediction.
Remark 17. When , and , the system (3) will degenerated to the model in [20].
Remark 18. From Figures 3–5, we can see that whatever the initial value are, the diseasefree equilibrium is always asymptotically stable, as long as we take suitable control measures (i.e., we can take one or two control measures such that ).
Remark 19. By construct Lyanpunov functions, we proved the global stability of the diseasefree equilibrium for system (3). However, according to Remark 16, we find the endemic equilibrium may stable or unstable; and the value of is crucial for the dynamics of system (3).
Remark 20. In fact, HBV can be transmitted not only from virustocell, but also from celltocell. In this article, only virustocell is included in our model. We will consider both virustocell and celltocell in our future work.
In this article, the fractional derivatives of the system (3) has the same order for different variables. In our future work, we will consider the dynamic behavior of models with different fractional order derivatives.
Data Availability
The data used to support the findings of this study are included within the article.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Authors’ Contributions
Each of the authors, Ruiqing Shi, Ting Lu, and Cuihong Wang contributed to each part of this work equally and read and approved the final version of the manuscript.
Acknowledgments
This work is partly supported by Fund Program for the Scientific Activities of Selected Returned Overseas Professionals in Shanxi Province under Grant 201825.
References
 J. Wang, Y. Xiao, and R. A. Cheke, “Modelling the effects of contaminated environments on HFMD infections in mainland China,” BioSystems, vol. 140, pp. 1–7, 2016. View at: Publisher Site  Google Scholar
 H. Qi, L. Liu, and X. Meng, “Dynamics of a nonautonomous stochastic SIS epidemic model with double epidemic hypothesis,” Complexity, vol. 2017, Article ID 4861391, 14 pages, 2017. View at: Publisher Site  Google Scholar
 J. Zhang, Y. Li, Z. Jin, and H. Zhu, “Dynamics analysis of an avian influenza A (H7N9) epidemic model with vaccination and seasonality,” Complexity, vol. 2019, Article ID 4161287, p. 1, 2019. View at: Publisher Site  Google Scholar
 A. R. Carvalho and C. M. Pinto, “Energence of drugresistance in HIV dynamics under distict HAART regimes,” Communications in Nonlinear Science and Numerical Simulation, vol. 30, no. 13, pp. 207–226, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 B. Tang, Y. Xiao, S. Tang, and J. Wu, “Modelling weekly vector control against Dengue in the Guangdong Province of China,” Journal of Theoretical Biology, vol. 410, pp. 65–76, 2016. View at: Publisher Site  Google Scholar
 J. Lin, R. Xu, and X. Tian, “Threshold dynamics of an HIV1 virus model with both virustocell and celltocell transmissions, intracellular delay, and humoral immunity,” Applied Mathematics and Computation, vol. 315, pp. 516–530, 2017. View at: Publisher Site  Google Scholar
 M. A. Nowak and C. R. M. Bangham, “Population dynamics of immune responses to persistent viruses,” Science, vol. 272, no. 5258, pp. 74–79, 1996. View at: Publisher Site  Google Scholar
 M. A. Nowak and R. M. May, Virus Dynamics: Mathematical Principles of Immunology and Virology, Oxford University Press, London, UK, 2000. View at: MathSciNet
 M. A. Nowak, S. Bonhoeffer, A. M. Hill, R. Boehme, H. C. Thomas, and H. Mcdade, “Viral dynamics in hepatitis B virus infection,” Proceedings of the National Acadamy of Sciences of the United States of America, vol. 93, no. 9, pp. 4398–4402, 1996. View at: Publisher Site  Google Scholar
 S. M. Ciupe, R. M. Ribeiro, P. W. Nelson, and A. S. Perelson, “Modeling the mechanisms of acute hepatitis B virus infection,” Journal of Theoretical Biology, vol. 247, no. 1, pp. 23–35, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 K. Wang and W. Wang, “Propagation of HBV with spatial dependence,” Mathematical Biosciences, vol. 210, no. 1, pp. 78–95, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 A. U. Neumann, N. P. Lam, H. Dahari et al., “Hepatitis C viral dynamics in vivo and the antiviral efficacy of interferonα therapy,” Science, vol. 282, no. 5386, pp. 103–107, 1998. View at: Publisher Site  Google Scholar
 A. S. Perelson, E. Herrmann, F. Micol, and S. Zeuzem, “New kinetic models for the hepatitis C virus,” Hepatology, vol. 42, no. 4, pp. 749–754, 2005. View at: Publisher Site  Google Scholar
 WHO, Fact Sheet WHO/204 Revised October 2000: Hepatitis B, 2000, http://www.who.int/mediacentred/factsheets/fs204/en/.
 K. Moriyama, H. Okamoto, F. Tsuda, and M. Mayumi, “Reduced precore transcription and enhanced corepregenome transcription of hepatitis B virus DNA after replacement of the precorecore promoter with sequences associated with e antigenseronegative persistent infections,” Virology, vol. 226, no. 2, pp. 269–280, 1996. View at: Publisher Site  Google Scholar
 S. R. Lewin, R. M. Ribeiro, T. Walters et al., “Analysis of hepatitis B viral load decline under potent therapy: Complex decay profiles observed,” Hepatology, vol. 34, no. 5, pp. 1012–1020, 2001. View at: Publisher Site  Google Scholar
 L. G. Guidotti, R. Rochford, J. Chung, M. Shapiro, R. Purcell, and F. V. Chisari, “Viral clearance without destruction of infected cells during acute HBV infection,” Science, vol. 284, no. 5415, pp. 825–829, 1999. View at: Publisher Site  Google Scholar
 A. S. Perelson, “Modelling viral and immune system dynamics,” Nature Reviews Immunology, vol. 2, no. 1, pp. 28–36, 2002. View at: Publisher Site  Google Scholar
 H. Dahari, E. Shudo, R. M. Ribeiro, and A. S. Perelson, “Modeling complex decay profiles of hepatitis B virus during antiviral therapy,” Hepatology, vol. 49, no. 1, pp. 32–38, 2009. View at: Publisher Site  Google Scholar
 C. VargasDeLeón, “Stability analysis of a model for HBV infection with cure of infected cells and intracellular delay,” Applied Mathematics and Computation, vol. 219, no. 1, pp. 389–398, 2012. View at: Publisher Site  Google Scholar
 C. Huang, L. Cai, and J. Cao, “Linear control for synchronization of a fractionalorder timedelayed chaotic financial system,” Chaos, Solitons & Fractals, vol. 113, pp. 326–332, 2018. View at: Publisher Site  Google Scholar
 R. Rakkiyappan, G. Velmurugan, and J. Cao, “Stability analysis of fractionalorder complexvalued neural networks with time delays,” Chaos, Solitons & Fractals, vol. 78, pp. 297–316, 2015. View at: Publisher Site  Google Scholar
 F. A. Rihan, D. H. Abdel Rahman, and S. Lakshmanan, “A time delay model of tumourimmune system interactions: global dynamics, parameter estimation, sensitivity analysis,” Applied Mathematics and Computation, vol. 232, pp. 606–623, 2014. View at: Publisher Site  Google Scholar  MathSciNet
 C. Huang, J. Cao, M. Xiao, A. Alsaedi, and F. E. Alsaadi, “Controlling bifurcation in a delayed fractional predatorprey system with incommensurate orders,” Applied Mathematics and Computation, vol. 293, pp. 293–310, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 C. M. Pinto and A. R. Carvalho, “A latency fractional order model for HIV dynamics,” Journal of Computational and Applied Mathematics, vol. 312, pp. 240–256, 2017. View at: Publisher Site  Google Scholar  MathSciNet
 I. Podlubny, Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of their Solution and Some of their Applications, vol. 198, Academic Press, San Diego, Calif, USA, 1999. View at: MathSciNet
 Y. Li, D. Zhao, Y. Chen, I. Podlubny, and C. Zhang, “Finite energy Lyapunov function candidate for fractional order general nonlinear systems,” Communications in Nonlinear Science and Numerical Simulation, vol. 78, 2019. View at: Publisher Site  Google Scholar
 F. A. Rihan, Q. M. AlMdallal, H. J. AlSakaji, and A. Hashish, “A fractionalorder epidemic model with timedelay and nonlinear incidence rate,” Chaos, Solitons & Fractals, vol. 126, pp. 97–105, 2019. View at: Publisher Site  Google Scholar  MathSciNet
 Y. Ferdi, “Some applications of fractional order calculus to design digital filters for biomedical signal processing,” Journal of Mechanics in Medicine and Biology, vol. 12, no. 2, Article ID 12400088, 13 pages, 2012. View at: Publisher Site  Google Scholar
 L. Huo, J. Jiang, S. Gong, and B. He, “Dynamical behavior of a rumor transmission model with Hollingtype II functional response in emergency event,” Physica A: Statistical Mechanics and its Applications, vol. 450, pp. 228–240, 2016. View at: Publisher Site  Google Scholar  MathSciNet
 V. Castellanos and R. E. ChanLópez, “Existence of limit cycles in a three level trophic chain with LotkaVolterra and Holling type II functional responses,” Chaos, Solitons & Fractals, vol. 95, pp. 157–167, 2017. View at: Publisher Site  Google Scholar
 X. Meng, F. Li, and S. Gao, “Global analysis and numerical simulations of a novel stochastic ecoepidemiological model with time delay,” Applied Mathematics and Computation, vol. 339, pp. 701–726, 2018. View at: Publisher Site  Google Scholar  MathSciNet
 Z. M. Odibat and N. T. Shawagfeh, “Generalized Taylor's formula,” Applied Mathematics and Computation, vol. 186, no. 1, pp. 286–293, 2007. View at: Publisher Site  Google Scholar  MathSciNet
 S. Mao, R. Xu, and Y. Li, “A fractional order sirs model with standard incidence rate,” Journal of Beihua University (Natural Science), vol. 12, no. 4, pp. 79–382, 2012. View at: Google Scholar
 P. van den Driessche and J. Watmough, “Reproduction numbers and subthreshold endemic equilibria for compartmental models of disease transmission,” Mathematical Biosciences, vol. 180, pp. 29–48, 2002. View at: Publisher Site  Google Scholar  MathSciNet
 E. Ahmed, A. M. A. ElSayed, and H. A. A. ElSaka, “On some RouthHurwitz conditions for fractional order differential equations and their applications in Lorenz, Rssler, Chua and Chen systems,” Physics Letters A, vol. 358, pp. 1–4, 2006. View at: Google Scholar
 S. M. Salman and A. M. Yousef, “On a fractionalorder model for HBV infection with cure of infected cells,” Journal of the Egyptian Mathematical Society, vol. 25, no. 4, pp. 445–451, 2017. View at: Publisher Site  Google Scholar  MathSciNet
Copyright
Copyright © 2019 Ruiqing Shi 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.