Background {#Sec1}
==========

Mathematical modeling is a power full tool to study the dynamic of different diseases in the real world phenomena (Zaman et al. [@CR20], [@CR21]; Mann and Roberts [@CR13]; Zhao et al. [@CR22]). Mathematician and Biologist used different widespread models to understand the spreading of transmittible maladies in the population. Hepatitis B is a globally health problem and a prominent root of expiry in the world. Billion of populates are septic with Hepatitis B virus contamination. Throughout the world 2 billion individuals have been sick with Hepatitis B virus and around 2 million individuals live chronically infected with Hepatitis B. Every year nearly 780,000 individuals expires due to chronic or acute Hepatitis B virus infection (WHO [@CR19]). Thus from the evidence it is clear that Hepatitis B virus infection is responsible for about 80 % of the primary liver cancer.

The transmission of Hepatitis B virus are very common and can be transmitted form one individual to another individual on both ways horizontally and vertically, such as transmission of blood, semen, vaginal secretions (Lavanchy [@CR8]; Lok et al. [@CR12]; McMahon [@CR15]), unprotected sexual contact, sharing of razers, blades or tooth brushes (Chang [@CR4]) and infected mother to her/his child during the time of birth, respectively. However Hepatitis B virus cannot be transmitted through water, food, hugging, kissing and causal contact such as in the work place, school etc. (McMahon [@CR15]). The mode of transmission of Hepatitis B virus and HIV are the same, but Hepatitis B virus is fifty to hundreds time more infectious (Thornley et al. [@CR16]; Ali and Zaman [@CR1]).

Infection of Hepatitis B has different phases, such as acute Hepatitis, chronic and carrier Hepatitis. Acute infection of Hepatitis B denoted the first 6 months subsequently exposed some one to Hepatitis B virus. Herein the resistant organism is frequently capable to flawless the infection from the physique then the body must be recovered in a few months complectly for some one, but for the remaining infection grow and leads to more serious stage called chronic stage or life long sickness. Chronic stage of Hepatitis B raises to the infection, that happens when Hepatitis B virus remain in an individual body and over period, the contagion can develop stern healthiness problem. Individuals with carrier Hepatitis frequently have no antiquity of acute sickness, but it can cause liver damaging which become the cause liver damaging and could also grow to liver tumor Mann and Roberts ([@CR13]).

Numerous of authors have premeditated different widespread models, which label the dynamic of communicable diseases (Zou et al. [@CR23]; Thornley et al. [@CR16]; Mann and Roberts [@CR13]). The incidence rate is one of the key concept and theatres an important role in the field of mathematical modeling. Bilinear occurrence rate $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\beta SI$$\end{document}$ frequently recycled in various epidemic problems (Fan et al. [@CR5]; Li and Ma [@CR9]; Zou et al. [@CR23]), where $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\beta$$\end{document}$ represents the contact rate, *S* represents susceptible individuals and *I* represents the infectious individuals. For the first time Capasso and Serio introduced the saturated incidence rate $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\frac{\beta SI}{1+\alpha I},$$\end{document}$ which is the generalized form of bilinear incidence rate (Capasso and Serio [@CR3]). This occurrence rate is additional sensible than the bilinear occurrence rate, especially in the case of sexually transmittible diseases, since it embraces the negotiating alteration and swarming influence of the virulent individuals and inhibits the unboundedness of the interaction rate by indicating appropriate parameters, which was recycled in numerous of epidemic problems (Kar and Jana [@CR7]; Gomes et al. [@CR6]; Liu and Yang [@CR11]).

In this article, we present an epidemic problem for the transmission dynamic of Hepatitis B virus with saturated incidence rate, which is the modified version of Zou et al. ([@CR23]). After developing the new model, we find the basic reproduction number by using the well known approach i.i. next generation matrix (NGM) approach. Furthermore, we investigate the possible equilibriums i.e. disease free and endemic equilibria and show the local asymptotic stability as well as global asymptotic stability at both equilibriums. For the local asymptotic stability, we use linearization and Routh-Hurwitz criteria, while to discuss the global stability, we then use the classic Lypanavo function theory and geometrical approach. Finally the numerical simulation are achieved by exhausting Runge-Kutta method of order fourth scheme to show the feasibility of the obtained results and the role of saturated incidence rate.

The organization of the paper is as follows. In the second section, we presents the proposed model and studied its different analysis, including basic reproduction number and equilibriums. In section "Stability analysis", we study the stability analysis and prove the local as well as global stability. Section "Numerical analysis" is devoted to numerical simulation and discussion. Finally a brief conclusion is presented in section "Conclusion".

Mathematical model and its analysis {#Sec2}
===================================

Herein this section, we presents a Hepatitis B virus transmission epidemic model. For this, we split the entire populace into seven epidemiological subclasses, susceptible *S*(*t*), latent *L*(*t*), acute infected *A*(*t*), chronic infected individuals *B*(*t*), carrier individuals *C*(*t*), recovered with permanent immunity *R*(*t*) and vaccinated *V*(*t*). Thus the compartmental mathematical model can be represented by the succeeding system of seven ordinary differential equations,$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \frac{dS(t)}{dt} & = b\xi (1-\eta C(t))+\phi V(t)-\frac{\beta S(t)A(t)}{1+\alpha C(t)} -\frac{\gamma \beta S(t)B(t)}{1+\alpha C(t)}-\frac{\zeta \beta S(t)C(t)}{1+\alpha C(t)}-(\mu _0+v)S(t),\nonumber \\ \frac{dL(t)}{dt} & = \frac{\beta S(t)A(t)}{1+\alpha C(t)}+\frac{\gamma \beta S(t)B(t)}{1+\alpha C(t)}+\frac{\zeta \beta S(t)C(t)}{1+\alpha C(t)}-(\sigma +\mu _0)L(t),\nonumber \\ \frac{dA(t)}{dt}& = \sigma L(t)-(\mu _0+\gamma _1+\psi )A(t),\nonumber \\ \frac{dB(t)}{dt}& = p\gamma _1A(t)-(\mu _0+\mu _1+\gamma _2)B(t),\nonumber \\ \frac{dC(t)}{dt}& = b \xi \eta C(t)+(1-p)\gamma _1A(t)-(\mu _0+\mu _2+\gamma _3)C(t),\nonumber \\ \frac{dR(t)}{dt}& = \psi A(t)+\gamma _2 B(t)+\gamma _3 C(t)-\mu _0R(t),\nonumber \\ \frac{dV(t)}{dt}& = b(1-\xi )+v S(t)-(\mu _0+\phi )V(t). \end{aligned}$$\end{document}$$Here, *b* represents the birth rate, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\xi$$\end{document}$ represents the birth rate without successful vaccination, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\eta$$\end{document}$ represents the proportion of prenatally infected individuals, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\phi$$\end{document}$ represents the rate of waning vaccine induced immunity, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\beta$$\end{document}$ represents the transmission rate from susceptible to infected, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\gamma$$\end{document}$, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\zeta$$\end{document}$ represent the reduced transmission rate of chronic and carrier Hepatitis B infected individuals, respectively. $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\mu _0$$\end{document}$ represents the death rate, which occur naturally, *v* represents the vaccination rate, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma$$\end{document}$ represents the moving rate from latent class to acute class, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\gamma _1$$\end{document}$ represents the moving rate from acute to chronic and carrier, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\gamma _2$$\end{document}$ represents the moving rate of chronic carrier to immune, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\gamma _3$$\end{document}$ represents the moving rate of carrier to immune, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\mu _1,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\mu _2$$\end{document}$ represents the death rate, which occur from the Hepatitis B and *p* represents the average probability of those individuals, who fails to recovers in acute stage and goes to chronic carrier.

Let "$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$N(t)$$\end{document}$" represents the entire populace, such that, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$N(t)=S(t)+L(t)+A(t)+B(t)+C(t)+R(t)+V(t).$$\end{document}$ The initial conditions with the proposed model, make sure that, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$N(t)\ge 0.$$\end{document}$ Hence the entire populace "$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$N(t)$$\end{document}$" is bounded and persist positive for $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$t>0.$$\end{document}$ Now the time derivative of *N*(*t*) yields$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \frac{dN(t)}{dt}=b-\mu _0N(t)-\mu _1B(t)-\mu _2C(t), \end{aligned}$$\end{document}$$which implies that, for large time *t*,  that is $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$t\rightarrow \infty ,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$N(t)\le \frac{b}{\mu _0}.$$\end{document}$ Thus the feasible region for our proposed model is$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \Omega =\left\{ (S(t),L(t),A(t),B(t),C(t),R(t),V(t))\in R_{+}^7,N\le \frac{b}{\mu _0}\right\} . \end{aligned}$$\end{document}$$Therefore, we study the dynamic of our proposed model (1) in the feasible region $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\Omega ,$$\end{document}$ which is an attracting set.

For obtaining the Jacobian matrix of our model, we take the reduced system, because *R* dose not appears explicitly in all others classes, so the Jacobian matrix of the reduced system (with out *R*) becomes$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} J=\left( \begin{array}{cccccc} -a_{11} &{} 0 &{} -a_{13} &{} -a_{14} &{} -a_{15} &{} \phi \\ a_{21} &{} -a_{22} &{} a_{23} &{} a_{24} &{} a_{25} &{} 0\\ 0 &{} \sigma &{} -a_{33} &{} 0 &{} 0 &{} 0\\ 0 &{} 0 &{} p\gamma _1 &{} a_{44} &{} 0 &{} 0\\ 0 &{} 0 &{} a_{53} &{} 0 &{} -a_{55} &{} 0\\ v &{} 0 &{} 0 &{} 0 &{} 0 &{} -a_{66} \end{array} \right) , \end{aligned}$$\end{document}$$where$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} a_{11} &= \mu _0+v+\frac{1}{1+\alpha C}(\beta A+\gamma \beta B+\zeta \beta C),~a_{13}=a_{23}=\frac{\beta S}{1+\alpha C},~a_{14}=a_{24}=\frac{\gamma \beta B}{1+\alpha C},\\ a_{15} &= \frac{1}{(1+\alpha C)^2}(\zeta \beta S-\alpha \beta S A-\alpha \gamma \beta S B),~a_{21}=a_{11}-(\mu _0+v),~a_{22}=\sigma +\mu _0,\\ a_{33} &= \mu _0+\gamma _1+\psi , ~a_{44}=\mu _0+\mu _1+\gamma _2,~a_{55}=b\xi \eta -(\mu _0+\mu _2+\gamma _3),~a_{66}=\mu _0+\phi . \end{aligned}$$\end{document}$$

Basic reproduction number {#Sec3}
-------------------------

In epidemiological models the inception capacity $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0$$\end{document}$ is called the basic reproduction number is a key concept, which characterizes the anticipated average amount of new contaminations created directly and circuitously by a solitary infectious, when familiarized into a entirely susceptible populace (Anderson and May [@CR2]; van den Driessche and Watmough [@CR17]). To find the this quantity for our proposed model (1), we use the method of Van Den Driessche and Watmough ([@CR18]). Let $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\chi =(L(t),A(t),B(t),C(t)),$$\end{document}$ so from the model ([1](#Equ1){ref-type=""}), we have$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \frac{d\chi }{dt}=\bar{F}-\bar{V}. \end{aligned}$$\end{document}$$In Eq. ([5](#Equ5){ref-type=""}), $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\bar{F}$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\bar{V}$$\end{document}$ are define as$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \bar{F}=\left( \begin{array}{c} \beta S(t)A(t)+\gamma \beta S(t)B(t)+\zeta \beta S(t)C(t) \\ 0 \\ 0 \\ 0 \\ \end{array} \right) ,\bar{V}=\left( \begin{array}{c} r_1\\ r_2\\ r_3\\ r_4\\ \end{array} \right) , \end{aligned}$$\end{document}$$where $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$r_1=(\sigma +\mu _0)L(t),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$r_2=(\mu _0+\gamma _1+\psi )A(t)-\sigma L(t),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$r_3=(\mu _0+\mu _1+\gamma _2)B(t)-p\gamma _1A(t)$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$r_4=(\mu _0+\mu _2+\gamma _3)C(t)-b\xi \eta C(t)-(1-p)\gamma _1 A(t).$$\end{document}$ Now we find the Jacobian of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\bar{F}$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\bar{V}$$\end{document}$ at diseases free equilibrium $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0,$$\end{document}$ thus we have$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} F =\left( \begin{array}{cccc} 0 &\quad \beta S_0 &\quad \gamma \beta S_0 &\quad \zeta \beta S_0 \\ 0 &\quad 0 &\quad 0 &\quad 0 \\ 0 &\quad 0 &\quad 0 &\quad 0 \\ 0 &\quad 0 &\quad 0 &\quad 0 \\ \end{array} \right) ,V =\left( \begin{array}{cccc} a_{11} &\quad 0 &\quad 0 &\quad 0 \\ -\sigma &\quad a_{22} &\quad 0 &\quad 0 \\ 0 &\quad -a_{32} &\quad a_{33} &\quad 0 \\ 0 &\quad -a_{41} &\quad 0 &\quad a_{44} \\ \end{array} \right) , \end{aligned}$$\end{document}$$where $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_{11}=\sigma +\mu _0,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_{22}=\mu _0+\gamma _1+\psi ,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_{32}=p\gamma _1,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_{33}=\mu _0+\mu _1+\gamma _2,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_{41}=(1-p)\gamma _1$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_{44}=\mu _0+\mu _0+\gamma _3-b\xi \eta .$$\end{document}$ Thus $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0$$\end{document}$ is the spectral radius of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\bar{K}=F V^{-1},$$\end{document}$ that is $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0=\rho (FV^{-1} )$$\end{document}$. So the basic reproduction number $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0$$\end{document}$ for our proposed model ([1](#Equ1){ref-type=""}) becomes$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} R_0=\bar{R_1}+\bar{R_2}+\bar{R_3}, \end{aligned}$$\end{document}$$where$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \bar{R_1} &= \frac{\sigma \beta S_0}{(\sigma +\mu _0)(\gamma _1+\psi +\mu _0)},\\ \bar{R_2} &= \frac{\sigma \beta \gamma \gamma _1p S_0 }{(\sigma +\mu _0)(\gamma _1+\psi +\mu _0)(\gamma _2+\mu _0+\mu _1)},\\ \bar{R_3} &= \frac{\sigma \beta \zeta \gamma _1 (1-p) S_0}{(\sigma +\mu _0)(\gamma _1+\psi +\mu _0)(\gamma _3+\mu _0+\mu _2-b \xi \eta )} . \end{aligned}$$\end{document}$$

Equilibrium analysis {#Sec4}
--------------------

To study the dynamics of our proposed model ([1](#Equ1){ref-type=""}), we need to find the equilibriums, thus the disease free equilibrium of the model ([1](#Equ1){ref-type=""}) is denoted by $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$ and define as $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0=(S_0,0,0,0,0,0,V_0),$$\end{document}$ where$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} S_0=\frac{b(\phi +\mu _0\xi )}{\mu _0(\mu _0+v+\phi )},~V_0=\frac{b(\mu _0+v-\mu _0\xi )}{\mu _0(\mu _0+v+\phi )}. \end{aligned}$$\end{document}$$Similarly endemic equilibrium is denoted by $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1$$\end{document}$ and define as $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1=(S_1,L_1,A_1,C_1,R_1,V_1),$$\end{document}$ where$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \nonumber S_1 & = S_0-\frac{1}{q_3q_5q_7}b\xi \eta \gamma _1(1-p)q_6L_1-\frac{1}{q_7}q_2q_6L_1,\\ \nonumber L_1& = \frac{q_2q^2_3q_4q^2_5q_7(R_0-1)}{\sigma (w_1(1-p) +\beta (1+p\gamma _1)q_2q_3q_4q^2_5q_6)},\\ \nonumber A_1 &= \frac{q_2q_3q_4q^2_5q_7(R_0-1)}{w_1(1-p) +\beta (1+p\gamma _1)q_2q_3q_4q^2_5q_6},\\ \nonumber B_1 &= \frac{q_2q_3q^2_5q_7p\gamma _1(R_0-1)}{w_1(1-p) +\beta (1+p\gamma _1)q_2q_3q_4q^2_5q_6},\\ C_1 &= \frac{q_2q_3q_4q_5q_7(1-p)(R_0-1)}{w_1(1-p) +\beta (1+p\gamma _1)q_2q_3q_4q^2_5q_6},\\ \nonumber R_1 & = \frac{1}{\mu _0}(\psi A_1+\gamma _2B_1+\gamma _3C_1),\\ \nonumber V_1 &= \frac{1}{q_6}(b(1-\xi )+vS_1), \end{aligned}$$\end{document}$$where $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$q_1=\mu _0+v,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$q_2=\sigma +\mu _0,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$q_3=\mu _0+\gamma _1+\psi ,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$q_4=\mu _0+\mu _1+\gamma _2,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$q_5=\mu _0+\mu _2+\gamma _3-b\xi \eta ,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$q_6=\mu _0+\phi ,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$q_7=\mu _0(\mu _0+v+\phi )$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$w_1=\sigma \beta b \xi \eta \gamma _1 q_5(q_4+\gamma _1pq_6+\gamma _1q_2q_3q_4q_5(\alpha q_7+\beta q_6)+\alpha \beta b\xi \eta \gamma ^2_1q_4q_6)$$\end{document}$.

Stability analysis {#Sec5}
==================

Herein this section, we want to study asymptotic stability of our proposed model at both equilibria. So first, we prove the local stability of our model at illness free and endemic equilibrium, then the global asymptotic stability. For the local stability, we use the linearization method, but for global dynamic stability, we want to use both the Lypanavo function theory and geometrical approach.

Local stability analysis {#Sec6}
------------------------

For the local dynamic of the proposed model at disease free and endemic equilibrium points, we state and prove the following results.

### **Theorem 1** {#FPar1}

*If*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0<1$$\end{document}$, *then the model* (1) *is locally asymptotically stable at disease free equilibrium point*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$*and if*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0>1$$\end{document}$, *then it is the unstable saddle point*.

### *Proof* {#FPar2}

The characteristic equation of the Jacobian matrix *J* ([4](#Equ4){ref-type=""}) at disease free equilibrium $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$ becomes$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} P(\lambda )=(\lambda +\mu _0)(\lambda +\mu _0+v+\phi )(\lambda ^4+a_1\lambda ^3+a_2\lambda ^2+a_3\lambda +a_4)=0, \end{aligned}$$\end{document}$$where$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} a_1 & = a_{22}+a_{33}+a_{44}+a_{55},\\ a_2 & = a_{22}a_{33}(1-R_1)+a_{44}a_{55}+a_{55}a_{33}+a_{33}a_{44}+a_{22}a_{44}+a_{22}a_{55},\\ a_3 & = a_{55}a_{33}a_{22}(1-R_1)+a_{33}a_{44}a_{55}\left( 1-\frac{a_{22}}{a_{55}}R_2\right) +a_{22}a_{44}a_{55}(1-R_2)+a_{22}a_{33}a_{44}(1-R_1),\\ a_4 & = \sigma \gamma \beta p \gamma _1S_0a_{55}+\sigma \zeta p\gamma _1 \beta S_0a_{44}+a_{22}a_{33}a_{55}-\sigma \beta S_0a_{44}b_{55}-\sigma \zeta \beta \gamma _1S_0a_{44}. \end{aligned}$$\end{document}$$If $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0<1,$$\end{document}$ we have $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$0<\bar{R_i}<1,$$\end{document}$ for $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$i=1,2,3.$$\end{document}$ So $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_j>0,$$\end{document}$ for $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$j=1,2,3,4,$$\end{document}$ and also it is easy to show that, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_1a_2a_3>a^{2}_3+a^2_2a_4.$$\end{document}$ Therefore Routh-Herwitz criteria is satisfied, that is all the roots of the characteristic polynomial $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$P(\lambda )$$\end{document}$ have negative real parts, which ensure that the disease free equilibrium point $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$ is stable, while for $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0>1,$$\end{document}$ the system have both negative and positive eigenvalues, which states that the disease free equilibrium point is a saddle point, always unstable.

### **Theorem 2** {#FPar3}

*If*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0>1$$\end{document}$, *then the model* (1) *is locally asymptotically stable at endemic equilibrium point*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1$$\end{document}$*and if*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0<1$$\end{document}$, *then it is the unstable*.

### *Proof* {#FPar4}

Using the elementary row operation, reducing the Jacobian matrix at Eq. [4](#Equ4){ref-type=""} around the endemic equilibrium $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1$$\end{document}$, we get the following echelon form$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} J_1=\left( \begin{array}{cccccc} -l_{11} &{} 0 &{} -l_{12} &{} -l_{13} &{} -l_{14} &{} \phi \\ 0 &{} -l_{22} &{} l_{23} &{} l_{24} &{} l_{25} &{} \phi l_{21} \\ 0 &{} 0 &{} l_{33} &{} l_{34} &{} 0 &{} \sigma \phi l_{21} \\ 0 &{} 0 &{} 0 &{} l_{44} &{} l_{45} &{} -\sigma \phi l_{21} \\ 0 &{} 0 &{} 0 &{} 0 &{} l_{44} &{} \sigma \phi l_{21} \\ 0 &{} 0 &{} 0 &{} 0 &{} 0 &{} l_{66} \\ \end{array} \right) , \end{aligned}$$\end{document}$$where $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{11}=a_{11},$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{12}=a_{12},$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{13}=a_{13},$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{15}=a_{15},$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{22}=a_{11}a_{22},$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{23}=a_{23}(a_{11}-a_{12}),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{24}=a_{14}(a_{11}-a_{21}),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{25}=a_{15}(a_{11}-a_{21}),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{33}=-a_{11}a_{22}a_{33}+\sigma a_{13}(a_{11}-a_{21}),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{34}=\sigma a_{14}(a_{11}-a_{21}),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{35}=\sigma a_{14}(a_{11}-a_{21}),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{44}=-\frac{1}{p\gamma _1}\sigma a_{13}a_{44}(a_{11}a_{22}a_{33}+a_{11}-a_{21}),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{45}=a_{14}(a_{11}-a_{21}),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{45}=a_{14}(a_{11}-a_{21}),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{55}=l_{44}+\sigma a_{14}(a_{11}-a_{21})$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$l_{66}=l_{44}-\frac{\sigma \phi v a_{14} a_{44}a_{53}(\mu _0+v)}{\beta p\gamma _1\zeta a_{44}+\beta a_{45}(p\gamma \gamma _1+a_{44}))}.$$\end{document}$

Thus after some simplification and little re-arrangement the eigenvalues of the Jacobian matrix $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$J_1$$\end{document}$ are given by$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \nonumber \lambda _1 & = -a_{11},~\lambda _2=-a_{22},~\lambda _3=-a_{11}a_{22}a_{33}+\sigma a_{13}(a_{11}-a_{21}),\\ \lambda _4 & = -\frac{1}{p\gamma _1}\sigma a_{13}a_{44}(a_{11}a_{22}a_{33}+a_{11}-a_{21}),\\ \lambda _5 & = -l_{44}+\sigma a_{14}(a_{11}-a_{21}),\\ \lambda _6 & = l_{44}-\frac{\sigma \phi v a_{14} a_{44}a_{53}(\mu _0+v)}{\beta p\gamma _1\zeta a_{44}+\beta a_{45}(p\gamma \gamma _1+a_{44}))}. \end{aligned}$$\end{document}$$Obviously $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\lambda _1,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\lambda _2,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\lambda _4,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\lambda _6$$\end{document}$ have negative real parts and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\lambda _3,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\lambda _5$$\end{document}$ are negative, if $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_{22}a_{33}>\sigma \beta S_1$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_{44}S_1>\gamma B_1.$$\end{document}$ Therefore for $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0>1,$$\end{document}$ the model (1) is locally asymptotically stable at endemic equilibrium point $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1,$$\end{document}$ if $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_{22}a_{33}>\sigma \beta S_1$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$a_{44}S_1>\gamma B_1.$$\end{document}$

Global stability analysis {#Sec7}
-------------------------

To show the global stability of the proposed model, using Lypanavo function theory and geometrical approach. Here, we want to prove the global stability of disease free equilibrium point $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$ by using Lypanavo function theory, while to prove the global stability at endemic equilibrium point $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1,$$\end{document}$ we use geometrical approach. Before to presents the stability results, the following important lemma due to Li and Muldowney ([@CR10]) is quoted. So, we obtain the condition at which the proposed model is globally asymptotically stable.

### **Lemma 1** {#FPar5}

*If the system*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\dot{y}=g(y),$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$g:D\rightarrow R^{n}$$\end{document}$*containing a unique equilibrium of the form*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$y^*$$\end{document}$*and there is a compact absorbing set, then this system is globally asymptotically stable around that equilibrium, if there exist a functionD*(*y*) *and a Lozinskii measure*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\ell ,$$\end{document}$*such that*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\lim _{t\rightarrow \infty }supsup\frac{1}{t}\int ^t_0\ell (B)dt<0.$$\end{document}$

Thus regarding the global stability of the proposed model at disease free and endemic equilibrium points, we have the following results.

### **Theorem 3** {#FPar6}

*If*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0<1$$\end{document}$, *then the model* ([1](#Equ1){ref-type=""}) *is globally asymptotically stable at disease free equilibrium point*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$*and unstable otherwise*.

### *Proof* {#FPar7}

To show the global stability at disease free equilibrium point $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$, we use Lypanavo function theory, so consider the following Lypanavo function, such that$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \nonumber F(t) & = \frac{1}{2}[(S-S_0)+L(t)+A(t)+B(t)+C(t)+R(t)+(V-V_0)]^2+\lambda _1S(t)+\lambda _2L(t)\\&\quad +\lambda _3A(t)+\lambda _4B(t)+\lambda _5V(t), \end{aligned}$$\end{document}$$where $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\lambda _i$$\end{document}$ for $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$i=1,2,3,4,5$$\end{document}$ are some positive constants, which will be chosen latter. After differentiating *F*(*t*) with respect to time, we get$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \nonumber \frac{dF}{dt} & = [(S-S_0)+L(t)+A(t)+B(t)+C(t)+R(t)+(V-V_0)](b-\mu _0N(t)-\mu _1B(t))\\& \quad+\lambda _1\frac{dS}{dt}+\lambda _2\frac{dL}{dt}+\lambda _3\frac{dA}{dt}+\lambda _4\frac{dB}{dt}+\lambda _5\frac{dV}{dt}. \end{aligned}$$\end{document}$$Choosing the positive parameters $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\lambda _i=q_2q_4,$$\end{document}$ for $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$i=1,2,3,5$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\lambda _4=\sigma \beta \gamma S_0$$\end{document}$ along the solution of the model ([1](#Equ1){ref-type=""}), we get$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \nonumber \frac{dF}{dt} & = [(S-S_0)+L(t)+A(t)+B(t)+C(t)+R(t)+(V-V_0)](b-\mu _0N(t)-\mu _1B(t))\\& \quad -q_2q_4(b\xi \eta C(t)-\mu _0L(t))-q_2q_3q_4(1-\bar{R_2})A(t)-\sigma \beta \gamma \gamma _1pS_0q_4B(t)\\ \nonumber& \quad +(b-\mu _0S(t)-\mu _0V(t))q_2q_4. \end{aligned}$$\end{document}$$Since $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$N(t)\le \frac{b}{\mu _0},$$\end{document}$ which implies that $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$S(t)\le \frac{b}{\mu _0}.$$\end{document}$ Thus using $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$S_0+V_0=\frac{b}{\mu _0}$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$S(t)\le \frac{b}{\mu _0}$$\end{document}$ in Eq. ([13](#Equ13){ref-type=""}), we obtain$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \nonumber \frac{dF}{dt}\le & {} -[(S-S_0)+L(t)+A(t)+B(t)+C(t)+R(t)+(V-V_0)](\mu _0(N(t)-S(t))\\&+\mu _1B(t))-q_2q_4(b\xi \eta C(t)+\mu _0L(t))-q_2q_3q_4(1-\bar{R_2})A(t)-\sigma \beta \gamma \gamma _1pS_0q_4B(t)\\ \nonumber&-\mu _0((S(t)+V(t))-(S_0+V_0)q_2q_4. \end{aligned}$$\end{document}$$If $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0<1,$$\end{document}$ we have $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$0<\bar{R_i}<1,$$\end{document}$ for $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$i=1,2,3,$$\end{document}$ therefore $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\frac{dF}{dt}$$\end{document}$ is negative, Also $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\frac{dF}{dt}=0,$$\end{document}$ if $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$S=S_0,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$L=L_0,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$A=A_0,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$B=B_0,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$C=C_0,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R=R_0$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$V=V_0,$$\end{document}$ thus the largest compact invariant set in $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\Omega$$\end{document}$ is the singleton set $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\{F_0\},$$\end{document}$ so "LaSalle's invariant principle" implies that, the illness free equilibrium point $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$ is globally asymptotically stable.

### **Theorem 4** {#FPar8}

*If*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0>1$$\end{document}$, *then the model* ([1](#Equ1){ref-type=""}) *is globally asymptotically stable at endemic equilibrium point*$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1$$\end{document}$.

### *Proof* {#FPar9}

Let $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$J_2$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$J_{3}$$\end{document}$ be the Jacobian matrix and second additive compound matrix containing only the first three equation of the model ([1](#Equ1){ref-type=""}), then we have$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} J_2=\left( \begin{array}{ccc} -a_{11} &{} 0 &{} -a_{13} \\ a_{21} &{} a_{22} &{} a_{23} \\ 0 &{} \sigma &{} -a_{33} \\ \end{array} \right) ,~~J_{3}=\left( \begin{array}{ccc} -(a_{11}+a_{22}) &{} a_{23} &{} -a_{13}\\ a_{32} &{} -(a_{11}+a_{33}) &{} a_{12}\\ -a_{31} &{} a_{21} &{} -(a_{22}+a_{33})\\ \end{array} \right) . \end{aligned}$$\end{document}$$Let us consider the function $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$P(\chi )=P(S,L,A)=diag\left\{ \frac{S}{L},\frac{S}{L},\frac{S}{L}\right\} ,$$\end{document}$ which implies that $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$P^{-1}(\chi )=diag\left\{ \frac{L}{S},\frac{L}{S},\frac{L}{S}\right\} ,$$\end{document}$ then taking the time derivative, that is $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$P_f(\chi )$$\end{document}$, we get$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} P_f(\chi )=diag\left\{ \frac{\dot{S}}{S}-\frac{S\dot{L}}{L^{2}},\frac{\dot{S}}{S}-\frac{S\dot{L}}{L^{2}},\frac{\dot{S}}{S}-\frac{S\dot{L}}{L^{2}}\right\} . \end{aligned}$$\end{document}$$Now $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$P_fP^{-1}=diag\left\{ \frac{\dot{S}}{S}-\frac{\dot{L}}{L},\frac{\dot{S}}{S}-\frac{\dot{L}}{L},\frac{\dot{S}}{S}-\frac{\dot{L}}{L}\right\}$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$PJ_{3}2P^{-1}=J_{3}.$$\end{document}$ Thus we take $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$B=P_fP^{-1}+PJ_{3}P^{-1}$$\end{document}$, which can be written as$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} B=\left( \begin{array}{cc} B_{11} &{} B_{12}\\ B_{21} &{} B_{22} \end{array} \right) , \end{aligned}$$\end{document}$$where$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} B_{11} & = \frac{\dot{S}}{S}-\frac{\dot{L}}{L}-\frac{\beta A(t)}{1+\alpha C(t)}-\frac{\gamma \beta B(t)}{1+\alpha C(t)}-\frac{\zeta B C(t)}{1+\alpha C(t)}-2\mu _0-v-\sigma ,\\ B_{12} & = \left( \begin{array}{cc} \frac{\beta S(t)}{1+\alpha C(t)}&{} \frac{\beta S(t)}{1+\alpha C(t)}\\ \end{array} \right) , ~B_{21}=\left( \begin{array}{c} \sigma \\ 0 \\ \end{array} \right) ,~ B_{22}=\left( \begin{array}{cc} x_{11} &{} 0 \\ x_{21}&{} x_{22} \\ \end{array} \right) . \end{aligned}$$\end{document}$$In the above matrix $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$x_{11}=\frac{\dot{S}}{S}-\frac{\dot{L}}{L}-\frac{\beta A(t)}{1+\alpha C(t)}-\frac{\gamma \beta B(t)}{1+\alpha C(t)}-\frac{\zeta \beta C(t)}{1+\alpha C(t)}-2\mu _0-v-\gamma _1-\psi ,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$x_{21}=\frac{\beta A(t)}{1+\alpha C(t)}+\frac{\gamma \beta B(t)}{1+\alpha C(t)}+\frac{\zeta \beta C(t)}{1+\alpha C(t)}$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$x_{22}=\frac{\dot{S}}{S}-\frac{\dot{L}}{L}-2\mu _0-\sigma -\gamma _1-\psi .$$\end{document}$ Let $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$(b_1,b_2,b_3)$$\end{document}$ be a vector in $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R^{3}$$\end{document}$ and its norm $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\Vert .\Vert$$\end{document}$ defined by$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \Vert b_1,b_2,b_3\Vert =max\{\Vert b_1\Vert ,\Vert b_2\Vert +\Vert b_3\Vert \}. \end{aligned}$$\end{document}$$Now we take the Lozinski measure $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\ell (B)$$\end{document}$ with respect to the above norm described by Martin ([@CR14]), that is $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\ell (B)\le sup\{g_1,g_2\}=sup\{\ell (B_{11})+\Vert B_{12}\Vert ,\ell (B_{22})+\Vert B_{21}\Vert \},$$\end{document}$ where $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$g_i=\ell (B_{ii})+\Vert B_{ij}\Vert$$\end{document}$ for $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$i=1,2$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$i\ne j,$$\end{document}$ which implies that$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} g_1=\ell (B_{11})+\Vert B_{12}\Vert ,~~g_2=\ell (B_{22})+\Vert B_{21}\Vert , \end{aligned}$$\end{document}$$where $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\ell (B_{11})=\frac{\dot{S}}{S}-\frac{\dot{L}}{L}-\frac{\beta A(t)}{1+\alpha C(t)}-\frac{\gamma \beta B(t)}{1+\alpha C(t)}-\frac{\zeta B C(t)}{1+\alpha C(t)}-2\mu _0-v-\sigma$$\end{document}$, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\Vert B_{12}\Vert =\frac{\beta S(t)}{1+\alpha C(t)}$$\end{document}$, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\ell (B_{22})=max\left\{ \frac{\dot{S}}{S}-\frac{\dot{L}}{L}-2\mu _0-v-\gamma _1-\psi ,\frac{\dot{S}}{S}-\frac{\dot{L}}{L}-2\mu _0-\sigma -\gamma _1-\psi \right\} = \frac{\dot{S}}{S}-\frac{\dot{L}}{L}-2\mu _0-\gamma _1-\psi -min\{v,\sigma \}$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\Vert B_{21}\Vert =max\{\sigma ,0\}=\sigma .$$\end{document}$ Therefore $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$g_1$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$g_2$$\end{document}$ becomes$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} g_1 & = \frac{\dot{S}}{S}-\frac{\dot{L}}{L}-\frac{\beta A(t)}{1+\alpha C(t)}-\frac{\gamma \beta B(t)}{1+\alpha C(t)}-\frac{\zeta \beta C(t)}{1+\alpha C(t)}-2\mu _0-min\{v,\sigma \}+\frac{\beta S(t)}{1+\alpha C(t)},\\ g_2 & = \frac{\dot{S}}{S}-\frac{\dot{L}}{L}-2\mu _0-\sigma -\gamma _1-\psi . \end{aligned}$$\end{document}$$Thus, we can write $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$g_1\le \frac{\dot{S}}{S}-2\mu _0-min\{v,\sigma \}$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$g_2\le \frac{\dot{S}}{S}-2\mu _0-\gamma _1-\psi ,$$\end{document}$ which implies that $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\ell (B)\le \left\{ \frac{\dot{S}}{S}-2\mu _0-min\{min\{v,\sigma \},\gamma _1-\psi \}\right\} .$$\end{document}$ Hence $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\ell (B)\le \frac{\dot{S}}{S}-2\mu _0.$$\end{document}$ Now integrating the Lozinski measure $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\ell (B)$$\end{document}$ with respect to *t* in the interval \[0, *t*\] and taking $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$lim_{t\rightarrow \infty },$$\end{document}$ we obtain$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} lim_{t\rightarrow \infty }sup~sup \frac{1}{t} \int ^{t}_{0}\ell (B)dS\le -2\mu _0<0. \end{aligned}$$\end{document}$$So finally, we can write$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \bar{q}=lim_{t\rightarrow \infty }sup~sup \frac{1}{t} \int ^{t}_{0}\ell (B)dS<0. \end{aligned}$$\end{document}$$Thus the subsystem, which containing first three equations of the model ([1](#Equ1){ref-type=""}) is globally asymptotically stable around its interior equilibrium $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$(S_1,L_1,A_1).$$\end{document}$ Now consider the subsystem of the model ([1](#Equ1){ref-type=""}), such that$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \nonumber \frac{dB(t)}{dt} & = p\gamma _1A(t)-(\mu _0+\mu _1+\gamma _2)B(t),\\ \nonumber \frac{dC(t)}{dt} & = b \xi \eta C(t)+(1-p)\gamma _1A(t)-(\mu _0+\mu _2+\gamma _3)C(t),\\ \frac{dR(t)}{dt} & = \psi A(t)+\gamma _2 B(t)+\gamma _3 C(t)-\mu _0R(t),\\ \nonumber \frac{dV(t)}{dt} & = b(1-\xi )+v S(t)-(\mu _0+\phi )V(t). \end{aligned}$$\end{document}$$Taking the limit system of the model ([22](#Equ22){ref-type=""}), we get$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\begin{aligned} \nonumber \frac{dB(t)}{dt} & = p\gamma _1A_1(t)-(\mu _0+\mu _1+\gamma _2)B(t),\\ \nonumber \frac{dC(t)}{dt} & = b \xi \eta C(t)+(1-p)\gamma _1A_1(t)-(\mu _0+\mu _2+\gamma _3)C(t),\\ \frac{dR(t)}{dt} & = \psi A_1(t)+\gamma _2 B_1(t)+\gamma _3 C_1(t)-\mu _0R(t),\\ \nonumber \frac{dV(t)}{dt} & = b(1-\xi )+v S_1(t)-(\mu _0+\phi )V(t). \end{aligned}$$\end{document}$$Solving system (38) and using the initial conditions *B*(0), *C*(0),  *R*(0) and *V*(0). So for large time *t* that is $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$t\rightarrow \infty$$\end{document}$, $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$B(t)\rightarrow B_1,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$C(t)\rightarrow C_1,$$\end{document}$$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R(t)\rightarrow R_1$$\end{document}$ and $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$V(t)\rightarrow V_1,$$\end{document}$ which is sufficient to prove that the endemic equilibrium point $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$E_1$$\end{document}$ is globally asymptotically stable.

Numerical analysis {#Sec8}
==================

In this section, we want to observe the dynamical behavior of our proposed model. In order to do this, we purpose numerical results by using Runge-Kutta of order 4th scheme which have used several authors for a wide range of problems consisting of ordinary differential equations. For the simulation purpose, we use different value of parameters used in the proposed model are given in the Table [1](#Tab1){ref-type="table"}. In the set of parameters some are taken from published articles, while some of the parameters are taken in such a way, that would be much more biologically feasible. Furthermore the time interval is taken 50 months with initial population for susceptible, latent, acute infected, chronically infected, carrier, recovered and vaccinated individuals as 100, 80, 60, 40, 20, 0 and 20, respectively.

In the study of biological dynamic especially in the transmission dynamic of infectious disease sensitivity analysis is very important. Because by sensitive analysis, we are able to judge the role of every parameter and thus easily, we are in the position to develop a strategy for eradicating the disease from the community. Here we carry out a sensitivity analysis of our proposed model ([1](#Equ1){ref-type=""}) with respect to the parameter $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\alpha ,$$\end{document}$ which are presented in Fig. [1](#Fig1){ref-type="fig"}.

In Fig. [1](#Fig1){ref-type="fig"}, the sensitivity analysis of susceptible, latent, acute infected with Hepatitis B, chronic infected with Hepatitis B, carrier, recovered and vaccinated are presented, which shows that, susceptible, recovered and vaccinated population are directly proportional to saturation $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\alpha ,$$\end{document}$ while latent, acute infected, chronically infected and carrier population are inversely proportional to saturation $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\alpha .$$\end{document}$ So from this brief analysis, it is observed that, saturation is also one of the factor to prevent the spreading of Hepatitis B in the world, which is one of the top three infectious diseases. Almost one third of the population are infected with Hepatitis B.

Figure [2](#Fig2){ref-type="fig"} represents the phase space diagram of susceptible, latent, recovered and susceptible, latent, vaccinated with stable endemic equilibrium and respective set of parameters given in Table  [1](#Tab1){ref-type="table"}. Figure [3](#Fig3){ref-type="fig"} represents the phase space diagram of latent, acute, chronically infected and latent, acute, recovered. Similarly Figs. [4](#Fig4){ref-type="fig"} and  [5](#Fig5){ref-type="fig"} represents the phase space diagram of carrier, recovered, vaccinated, time, recovered, vaccinated, time, latent, carrier, time, latent and recovered respectively with stable endemic equilibrium and respective set of parameters given in Table [1](#Tab1){ref-type="table"}.Table 1The value of parameters used for numerical simulationNotationParameter descriptionValueSource*b*Birth rate0.0121Zou et al. ([@CR23])$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\xi$$\end{document}$Birth rate without successful vaccination0.05Assumed$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\eta$$\end{document}$Perintally infected individuals rate0.11Zou et al. ([@CR23])$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\phi$$\end{document}$Wanning vaccine induced immunity rate0.01Assumed$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\beta$$\end{document}$Transmission rate from susceptible to infected0.95Zou et al. ([@CR23])$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\gamma$$\end{document}$Reduced transmission rate0.12Assumed$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$d_0$$\end{document}$Natural mortality rate0.0069Assumed$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\gamma _3$$\end{document}$Vaccination rate0.9Assumed$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\sigma$$\end{document}$Moving rate from latent to acute0.0012Assumed$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\gamma _1$$\end{document}$Moving rate from acute to chronic carrier0.33Assumed$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\gamma _2$$\end{document}$Moving rate from chronic carrier to immune0.09Assumed$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$d_1$$\end{document}$Hepatitis B related death rate0.00054Assumed*p*Probability of those individuals, who fails to recover in acute class0.885Zou et al. ([@CR23])$\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$\alpha$$\end{document}$Saturation0--0.9Assumed

Fig. 1The plot shows the dynamical behavior of the proposed model ([1](#Equ1){ref-type=""}) with saturation and without saturationFig. 2The plot shows the phase space diagram of susceptible, latent, recovered and susceptible, latent and vaccinated with stable endemic equilibrium $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1$$\end{document}$ and respective set of parameters given in Table [1](#Tab1){ref-type="table"}Fig. 3The plot shows the phase space diagram of latent, acute infected, chronically infected and latent, acute infected, recovered with stable endemic equilibrium $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1$$\end{document}$ and respective set of parameters given in Table [1](#Tab1){ref-type="table"}Fig. 4The plot shows the phase space diagram of carrier, recovered and vaccinated with stable endemic equilibrium $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1$$\end{document}$ and respective set of parameters given in Table [1](#Tab1){ref-type="table"}Fig. 5The plot shows the phase space diagram of time, susceptible, vaccinated, time, latent, carrier and time, latent, recovered with stable endemic equilibrium $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1$$\end{document}$ and respective set of parameters given in Table [1](#Tab1){ref-type="table"}

Conclusion {#Sec9}
==========

The contagious disease of Hepatitis B is one of the topmost three infectious diseases reported by WHO and the incidence rate is precise significant in the study of epidemic models. In this article, we developed an epidemic problem for the spreading dynamic of Hepatitis B virus by classification of different subclasses with saturated occurrence rate. This occurrence rate is more generalized from than bilinear incidence rate and seems more realistic. Therefore we divided the host population into seven epidemiological subclasses, namely, susceptible, latent, acute infected with Hepatitis B, chronically infected with Hepatitis B, carrier, recovered and vaccinated, then formulated the model with this new features. First, we find the basic reproduction number $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0.$$\end{document}$ As in epidemiological models, the mathematical models has two steady states, infected and noninfected steady states. So we obtained the two non-negative equilibria that is disease free $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$, which is always exist and locally stable for $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0<1$$\end{document}$ and endemic equilibria $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1$$\end{document}$, which exist only, if $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0>1$$\end{document}$, therefore locally stable, if $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0>1$$\end{document}$. This portents occur, since totally the flows along the axes of susceptible, vaccinated and recovered are always attractor to the point $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$, but the flow along the axes of the infected compartment that is latent, infected with acute, chronic Hepatitis B and carrier depend on the value of $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0$$\end{document}$. So if $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0<1$$\end{document}$, then the axes of infected compartments are attractor towards the illness free equilibrium $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$. But when the $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0$$\end{document}$ acroses one that is $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$R_0>1$$\end{document}$, the axes of these compartment do not attract the disease free equilibrium $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$ and repels from it. Further more to show the global stability at $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_0$$\end{document}$, we established the Lyapnovo function, while to to prove the global stability at $\documentclass[12pt]{minimal}
                \usepackage{amsmath}
                \usepackage{wasysym} 
                \usepackage{amsfonts} 
                \usepackage{amssymb} 
                \usepackage{amsbsy}
                \usepackage{mathrsfs}
                \usepackage{upgreek}
                \setlength{\oddsidemargin}{-69pt}
                \begin{document}$$F_1$$\end{document}$, we used the geometrical approach. Finally numerics of the proposed model are retrieved out to show the feasibility of the model.

All authors contributed to the formulation of the model, discussion of results, wrote and approved the paper. All authors read and approved the final manuscript.

Acknowledgements {#FPar10}
================

This work has been partially supported by Higher Education Commission (HEC) of Pakistan under the project No. 20-1983/*R and D*/*HEC*/11.

Competing interests {#FPar11}
===================

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