Issue 
Wuhan Univ. J. Nat. Sci.
Volume 27, Number 2, April 2022



Page(s)  104  114  
DOI  https://doi.org/10.1051/wujns/2022272104  
Published online  20 May 2022 
Mathematics
CLC number: O175
Stability Research of an SEIR Model with Distinct General Contact Rates and Infectious Force in Latent and Recovered Period
Department of Common Course, Anhui Xinhua University, Hefei
230088, Anhui , China
Received:
22
December
2021
In this work, an SEIR infectious model with distinct general contact rates and infectious force in latent and recovered period is established, and the stability of the model is studied using theoretical and numerical methods. First, we derive the basic reproduction number R_{0}, which determines whether the disease is extinct or not. Second, using the LaSalle’s invariance principle, we show that the diseasefree equilibrium is globally asymptotically stable and the disease always dies out when R_{0}<1. On the other hand, by RouthHurwitz criterion theory, we also prove that the diseasefree equilibrium is unstable and that the unique endemic equilibrium is locally asymptotically stable when R_{0}>1. Third, through the method of autonomous convergence theorem, we obtain that the unique endemic equilibrium is globally asymptotically stable and the disease persists at this endemic equilibrium when R_{0}>1. Finally, numerical simulations are carried out to confirm the theoretical analysis.
Key words: basic reproductive number / equilibrium / stability / SEIR epidemic model / general contact rate
Biography: MA Yanli, female, Associate professor, research direction: epidemic dynamics. Email: Lindamayanli@163.com
Foundation item: Supported by Anhui Province Natural Science Research Foundation (KJ2019A0875, KJ2021A1155, KJ2021A1197)
© Wuhan University 2022
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
0 Introduction
Epidemiology is the research hotspots in the spread of infectious diseases in order to trace the factors that cause the occurrence. Mathematical models describing the behavior of infectious population have long played an important role in better analyzing disease control and transmission patterns. To predict the development of infectious diseases between regions, the transmission behavior of infectious diseases is studied in host population through a number of epidemic models. In recent years, many mathematical models of ordinary differential equations had been proposed for the transmission dynamics of infectious diseases^{[113]}. These models provided theoretical and quantitative bases for the prevention and control of infectious diseases.
However, some epidemic diseases, such as viral hepatitis, SARS, and whooping cough, become conta gious only after they have been dormant in the popula tion for a period of time.
Therefore, more general systems than the types of SIRS and SIR need to be considered to research the role of hatch in the development of infectious diseases. It can be supposed that susceptible population first experience an incubation period after infection before becoming infectious. The most general form of an epidemic model is an SEIR model consisting of four population subclasses: Ssusceptible, Eexposed, I infected and Rrecovered. All other models are cases of the SEIR model under various parameter restrictions. If there is no immunity and hence no R class, the SEIS model is obtained, which means the average period of immunity tends to be zero.
Many epidemic models with the infectious force in the latent period have been performed. Li ^{[14]} analyzed the global dynamics of an SEIR model with varying total population size. Sun ^{[15]} discussed the global dynamics of an SEIR model with a varying total population size and vaccination. Sun and Hsiehc^{[16]} studied the global analysis of an SEIR model with varying population size and vaccination. Yuan et al ^{[17]} considered threshold dynamics in an SEIRS model with latency and temporary immunity. Zhang^{[18]} studied global asymptotic stability of a delayed SEIRS epidemic model with saturation incidence. Zhang et al ^{[19]} considered an SEIR epidemic system with nonlinear transmission rate. Trawicki^{[20]} proposed a deterministic SEIRS epidemic model for modeling vital dynamics, vaccinations, and temporary immunity. We can see that the disease models mentioned above do not consider the infectious force in the latent and recovered period. However, for infectious diseases such as hepatitis C and condyloma acuminatum, incubation and recovery periods may be contagious. Although it is difficult to measure infectivity exactly, it is generally accepted that the degree of infectivity declines with the acquisition of immunity.
The incidence of diseases plays an extremely important role in the analysis of infectious disease models. Incidence should be recorded as , where N is the total population size. In many models of epidemiology, the bilinear incidence rate βSI and the standard incidence rate are often considered^{[2022]}. The bilinear incidence rate is duce to the law of mass action. It has been pointed out that for standard incidence rate, it may be a good approximation if the number of available partners is large enough and everybody could not make more contacts than it is practically feasible^{[23, 24]}. It has been suggested by several authors that the disease transmission process may have a saturation incidence rate, generally written as , where aI measures the infectivity of the disease and measures the inhibitory effect of behavioral changes or the crowding effect of infected individuals when the number of susceptible individuals increases. Compared with bilinear and standard incidence rate, saturating incidence rate may be a better fit for our actual situation. Several authors have studied epidemic models with saturation incidence rate^{[2529]}. In this paper, we introduce the general contact rate β(N)SI into the infectious disease model, which is more extensive.
In this paper, motivated by the work of Refs.[1330], we study an SEIR epidemic model with different general contact rates, which also has the infectious force in the exposed and recovered period. In Section 1, we formulate an SEIR model with different general contact rates, which also has the infectious force in the exposed and recovered period. In Section 2, we determine the basic reproduction number, obtain the existence of equilibriums and study the global stability of diseasefree equilibrium by LaSalle’s invariance principle. In Section 3, we analyze the local stability of the unique endemic equilibrium by RouthHurwitz criterion theory. In Section 4, we prove the unique endemic equilibrium is globally asymptotically stable through the limiting equation theory and the periodic orbit stability theory. In Section 5, we carry out numerical simulations to confirm the theoretical analysis. In the last section, we give a conclusion and prospect for the research work.
1 Model Formulation
In this work, we are interested in an SEIR model with different general contact rates, which also has the infectious force in the exposed and recovered period.
The total population is assumed to be divided into four distinct epidemiological subgroups, namely susceptible, latent, infectious and recovered (removed), with sizes denoted by S(t), E(t), I(t), and R(t), respectively. The size of the total population is expressed in N(t) at time t, where . The transfer mechanism for the class S to the class E is guided by the function:where are general contact rate, which satisfies the following assumptions, for N>0, , , , .
We construct the ordinary differential equation of the SEIR epidemic model as follows:(1)where A is a constant recruitment rate of the population; , , are general contact rates in the latent, infectious and recovered period, respectively; α_{1}, α_{2} are the rates of diseasecaused death of the exposed and the infectious, respectively; d is the natural death rate of the population; k_{1}, k_{2} are the elimination rates of the exposed and the infectious, respectively; γ is the transfer rates between the exposed and the infectious; ε is the removed rate from the infective class to the recovered class. The above parameters are positive.
Summing up the four equations of system (1) and denotinghavingFrom the formula of N', we obtain . From a biological point of view, we study system (1) in the following feasible regionwhere represents the nonnegative cone and its lower dimensional faces. It can be shown that D is positively invariant with respect to (1).
2 Stability Analysis of the DiseaseFree Equilibrium
In this work, we aim to show the existence of a diseasefree equilibrium and an endemic equilibrium of system (1) by solving equations and study the global stability of the diseasefree equilibrium by RouthHurwitz criterion theory and LaSalle’s invariance principle.
Let the right side of equation (1) be zero, and by calculating formula (2), we get(2)where , , and we have the following equation about N:whereSo system (2) always has the diseasefree equilibrium in the interval . SinceLet R_{0} is called the basic reproduction number of system (1). It is easy to get that if R_{0}>1, F(0)<0, . And
Hence, the equation F(N)=0 has only a positive root in the interval , that is, system (1) has an unique endemic equilibrium , where are determined by (2). It is easy to obtain the following theorem.
Theorem 1 If R_{0}≤1, the diseasefree equilibrium E_{0} of system (1) is globally asymptotically stable. If R_{0}>1, the diseasefree equilibrium E_{0} is unstable and the solutions to system (1) starting sufficiently close to E_{0} in D move away from E_{0} except those starting on the invariant Saxis approach E_{0} along this axis.
Proof To prove this result, we consider the following Lyapunov functionCalculating the derivative of V(t) along the positive solution of system (1), we getFurthermore, V'(t)=0 only if or R_{0}=1. The largest invariant set in is the singleton {E_{0}}. When R_{0}≤1, the global asymptotical stability of the diseasefree equilibrium E_{0} follows from LaSalle’s invariance. When R_{0}>1, S sufficiently closes to A/d except E(t)=I(t)=R(t)=0, we have the relation V'(t)>0. Hence the solutions to system (1) starting sufficiently close to E_{0} in D move away from E_{0} except those starting on the invariant Saxis approach E_{0} along this axis.
3 Local Stability Analysis of the Endemic Equilibrium
In this section, we analyze the local stability of the endemic equilibrium E* of system (1) with the help of the RouthHurwitz criterion theory. We have the following result.
Theorem 2 If R_{0}>1, the endemic equilibrium E* of system (1) is locally asymptotically stable.
Proof The Jacobian matrix of system (1) at the E* iswhereThen, the characteristic equation of J (E*) at E* is of the formwhereBy calculation, we havewhereBy calculation, we have . By Routh Hurwitz stability, all the four eigenvalues have negative real parts. Thus, the endemic equilibrium E* of system (1) is locally asymptotically stable in D, when R_{0}>1.
4 Global Stability Analysis of the Endemic Equilibrium
In this work, we aim to study the global dynamics of the endemic equilibrium E* of system (1) by means of the limiting equation theory and the method of autonomous convergence theorem.
In order to study the global stability of the endemic equilibrium E* of system (1), the following results can be found in Ref. [30]. We consider the following two systems(3) (4)where , , , x satisfies conditions of local Lipschitz. If as , system (4) is the limit system of system (3).
Lemma 1 Let P be a locally asymptotically stable equilibrium of system (4) and ω be the ω limit set of a forward bounded solution x(t) of system (3). If and the solution as satisfying the initial condition , then , i.e. as .
Corollary 1 If solutions to system (3) are bounded and the equilibrium P of the limit system (4) is globally asymptotically stable, then any solution x(t) to system (3) satisfies as .
Let dt=τ, system (1) becomes(5)where , , , .
By calculation, we have the following differential equationWe use N as a variable in place of the variable S to give the following system(6)System (6) is equivalent to system (5). From biological considerations, we study system (6) in the closed set We denote by and the boundary and the interior of T in , respectively. For system (6), the global stability of the endemic equilibrium E* is considered when . Since as , we can obtain the following limit system(7)We make the change of variable , y=E, z=I, then the following (5) is equivalent to the above system(8)
Theorem 3 Consider the following system^{[31]}:(9)where T is an open set. If system (9) satisfies the following conditions: ① System (9) exists a compact absorbing set K⊂T and has a unique equilibrium P in ; ② P is locally asymptotically stable; ③ System (9) satisfies the PoincareBendixson criterion; ④ The periodic orbit of system (9) is asymptotically orbitally stable. Then the only equilibrium P is the globally asymptotically stable in T.
Theorem 4 A sufficient condition for a periodic orbit of system (9) to be asymptotically orbitally stable with asymptotic phase is that the linear system is asymptotically stable, where is the second additive compound matrix of the Jacobian matrix of f. System (9) is called the second compound system of the orbit P(t).
Lemma 2 Any periodic solution to system (8), if it exists, is asymptotically orbitally stable.
Proof Suppose that the solution is periodic of least period τ>0 such that . The periodic orbit is . We have the second compound system x'= J ^{[2]}(P)x of the differential system in the periodic solution is the following periodic linear system(10)whereSuppose that is a solution to system (10). LetFrom the condition (1) of Theorem 3, we can know there exists constant η>0 such thatfor all and . By direct calculations, we can obtain the following differential inequalities:(11) (12) (13)Using (12) and (13), we have(14)From (11) and (14), we get(15)where(16) (17)Rewrite the last two equations of system (8), then we obtain(18) (19)Substituting (18) into (16) and (19) into (17), we have(20) (21)ThusandFrom (15), we have as , and in turn that as . As a result, the second compound system (10) is asymptotically stable and the periodic solution is asymptotically orbit ally stable by Theorem 3.
Lemma 3 System (8) is uniformly persistent when R_{0}>1.
Proof Let G={P_{0}}, when R_{0}>1, the stable set G ^{ s } is just contained in the Saxis and thus in the boundary of Γ. It also implies that the stable set G ^{ s } is isolated in Γ. Then, when R_{0}>1, system (8) satisfies the conditions of Theorem 2 of Ref. [1], namely, the maximal compact invariant set G in the boundary of Γ is isolated and the stable set G ^{ s } of G is contained in the boundary of Γ. Therefore, system (8) is uniformly persistent in Γ when R_{0}>1.
Lemma 4 System (5) is competitive when R_{0}>1.
Proof Set , system (5) is replaced byFurthermore, the system has the following formwhere , I denotes the 4×4 unit matrix, C (t) is a vector function,
The offdiagonal entries in this matrix are non negative, the system as a whole is quasimonotone. Thus we can verify that system (5) is competitive with respect to the partial ordering defined by the orthant , see Refs. [29, 32].
Theorem 5 If R_{0}>1, the endemic equilibrium E* of system (1) is globally asymptotically stable.
5 Example and Numerical Simulation
In this section, we provide the numerical simulation of the model to illustrate the main theoretical results above to better explain the occurrence and development of the infectious disease. The number simulations for system (1) are shown in Figs. 13. Choose . Let A=0.6, d=0.05, ε=0.07, γ=0.15, β_{1}=0.08, β_{2}=0.15, β_{3}=0.03. We randomly choose six initial conditions , , (2,1.5,2.6, 5), , and in By computing, we derive R_{0}=0.475<1 and system (1) has a diseasefree equilibrium E_{0}(12,0,0,0). We give the trajectory plot and its tridimensional figure for system (1) in Fig.1 and Fig. 3(a). From Theorem 1, it can be known that E_{0} is globally asymptotically stable. The dynamic behavior of system (1) is shown in Fig.1.
Fig.1 Variational curves of S, E, I and R with time t when R
_{0}=0.475<1 Six initial conditions are (a) (5, 1.9, 3.2, 2.6), (b) (3.8,0.6,2,5.4), (c) (2.1, 5, 2.6, 5), (d) (3.5, 2.5, 0.8, 4.6), (e) (0.5, 4, 2.2, 3.6) and (f) (4.5, 1.6, 2.4, 3.1), respectively 
Choose . Let A=1.8, d=0.05, ε=0.07, γ=0.15, β_{1}=0.1, β_{2}=0.2, β_{3}=0.05. We randomly choose six initial conditions , , , , and in . By computing, we derive R_{0}=5.508>1 and system (1) has a diseasefree equilibrium E*=(17.54,4.621,5.766,8.074). We give the trajectory plot and its tridimensional figure for system (1) in Fig.2 and Fig. 3(b). According to Theorem 1, E* is globally asymptotically stable. The dynamic behavior of system (1) is shown in Fig.2.
Fig.2 Variational curves of S, E, I and R with time t when R
_{0}=5.508>1 Six initial conditions are (a) (11.5, 9.2, 2.7, 3.8), (b) (9.8, 6.6, 8, 11.4), (c) (8, 7.5, 8.6, 11), (d) (6.5, 10, 8.2, 9.1), (e) (10.5, 7.6, 8.4, 9.1) and (f) (5.8, 9.6, 7.4, 12), respectively 
Fig.3 The graph of the trajectory in (S, E, I)space (a) and (b) correspond with Fig.1 and Fig.2, respectively, and the six curves correspond to six initial conditions 
6 Conclusion
In this work, we have performed a complete mathematical analysis of the global stability problem at equilibrium with SEIR models with general contact rates and latent period, both theoretically and numerically. For the model, we give the expression for the basic reproduction number R_{0}. When R_{0}<1, as shown in Theorem 1, the diseasefree equilibrium is globally asymptotically stable by LaSalle’s invariance principle (see Fig. 1 and Fig. 3(a)), and the disease always disappears eventually. When R_{0}>1, Theorem 5 tells us that by limit equation theory and periodic orbit stability theory, the only endemic equilibrium is globally asymptotically stable (see Fig. 2 and Fig. 3(b)), and the disease persists at the endemic equilibrium if it exists initially. If the infectivity of the disease in the exposed and recovered period is smaller, the smaller the basic regeneration number, the more conducive to the elimination of the disease. Therefore, for the infectious diseases of the exposed and recovered period, not only the patients in the infected period, but also the patients in the exposed and recovery period should be controlled, so as to control and eliminate the spread of the disease more effectively. Finally, numerical simulations are carried out to confirm the correctness of the theoretical analysis.
Interestingly, the stability of the model equilibrium is affected by the general contact rate. We believe that our findings offer guidance in face of the disease. In addition, we would like to point out here that the model (1) leaves us a problem that we only take the general form of contact rate, a more abstract situation in the model into consideration, but the incidence of other specific forms such as standard or nonlinear remains undiscussed. We leave this for our future work.
References
 Hou J, Teng Z D. Continuous and impulsive vaccination of SEIR epidemic models with saturation incidence rates [J]. Mathematics and Computers in Simulation, 2009, 79(10): 30383054. [CrossRef] [MathSciNet] [Google Scholar]
 Eckalbar J C, Eckalbar W L. Dynamics of an epidemic model with quadratic treatment [J]. Nonlinear Analysis: Real World Applications, 2011, 12(1):320332. [Google Scholar]
 Hu Z X, Ma W B, Ruan S G. Analysis of SIR epidemic models with nonlinear incidence rate and treatment[J]. Mathematical Biosciences, 2012, 238 (1):1220. [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 Hao L J, Jiang G R, Liu S Y, et al. Global dynamics of an SIRS epidemic model with saturation incidence [J]. Biosystems, 2013, 114 (1): 5663. [CrossRef] [PubMed] [Google Scholar]
 Wang J J, Zhang J Z, Jin Z. Analysis of an SIR model with bilinear incidence rate [J]. Nonlinear Analysis: Real World Application, 2010, 11 (4): 23902402. [Google Scholar]
 Misra A K, Sharma A, Shukla J B. Stability analysis and optimal control of an epidemic model with awareness programs by media [J]. Biosystems, 2015, 138(1):5362. [CrossRef] [PubMed] [Google Scholar]
 Ji C Y, Jiang D Q, Shi N Z. The behavior of an SIR epidemic model with stochastic perturbation [J]. Stochastic Analysis and Applications, 2012, 30 (5) :755773. [CrossRef] [MathSciNet] [Google Scholar]
 Yang Q S, Jiang D Q, Shi N Z, et al. The ergodicity and extinction of stochastically perturbed SIR and SEIR epidemic models with saturated incidence [J]. Journal of Mathematical Analysis and Applications, 2012, 388 (1): 248271. [CrossRef] [MathSciNet] [Google Scholar]
 Muroya Y, Enatsu Y, Nakata Y. Global stability of a delayed SIRS epidemic model with a nonmonotonic incidence rate [J]. Journal of Mathematical Analysis and Applications, 2011, 377 (1): 114. [CrossRef] [MathSciNet] [Google Scholar]
 Kuniya T, Muroya Y. Global stability of a multigroup SIS epidemic model with varying total population size[J]. Applied Mathematics and Computation, 2015, 265 (15): 785798. [CrossRef] [MathSciNet] [Google Scholar]
 Muroya Y, Kuniya T, Wang J L. Stability analysis of a delayed multigroup SIS epidemic model with nonlinear incidence rates and patch structure [J]. Journal of Mathematical Analysis and Applications, 2015, 425(1): 415439. [CrossRef] [MathSciNet] [Google Scholar]
 Muroya Y, Li H X, Kuniya T. Complete global analysis of an SIRS epidemic model with graded cure and incomplete recovery rates [J]. Journal of Mathematical Analysis and Applications, 2014, 410(2): 719732. [CrossRef] [MathSciNet] [Google Scholar]
 Wang J Y, Xiao Y N, Cheke R A. Modelling the effects of contaminated environments on HFMD infections in mainland China [J]. Biosystems, 2016, 140(1): 17. [CrossRef] [PubMed] [Google Scholar]
 Li M Y, Graef J R, Wang L C, et al. Global dynamics of a SEIR model with varying total population size [J]. Mathematical Biosciences, 1999, 160(2): 191213. [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 Sun S L . Global dynamics of a SEIR model with a varying total population size and vaccination [J]. Int Journal of Math Analysis, 2012, 6(40): 19851995. [MathSciNet] [Google Scholar]
 Sun C G, Hsieh Y H. Global analysis of an SEIR model with varying population size and vaccination[J]. Applied Mathematical Modelling, 2010, 34 (10): 26852697. [CrossRef] [MathSciNet] [Google Scholar]
 Yuan Y, Bélair J. Threshold dynamics in an SEIRS model with latency and emporary immunity [J]. Journal of Mathematical Biology, 2014, 69(1): 875904. [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 Zhang T L, Teng Z D. Global asymptotic stability of a delayed SEIRS epidemic model with saturation incidence [J]. Chaos, Solitons and Fractals, 2008, 37(5): 14561468. [Google Scholar]
 Zhang Q L, Liu C, Zhang X. Analysis and control of an SEIR epidemic system with nonlinear transmission rate [J]. Complexity, Analysis and Control of Singular Biological Systems, 2012, 421(1): 203225. [Google Scholar]
 Trawicki M B . Deterministic SEIRS epidemic model for modeling vital dynamics, vaccinations, and temporary immunity [J]. Mathematics, 2017, 5(1):725. [CrossRef] [Google Scholar]
 Zhao Z, Chen L S, Song X Y. Impulsive vaccination of SEIR epidemic model with time delay and nonlinear incidence rate [J]. Mathematics and Computers in Simulation, 2008, 79(3): 500510. [CrossRef] [MathSciNet] [Google Scholar]
 Ma Y L, Liu J B, Li H X. Global dynamics of an SIQR model with vaccination and elimination hybrid strategies [J]. Mathematics, 2018, 6(12): 112. [Google Scholar]
 Fan M, Li M Y, Wang K. Global stability of an SEIS epidemic model with recruitment and a varying total population size [J]. Mathematical Biosciences, 2001, 170(2): 199208. [CrossRef] [MathSciNet] [PubMed] [Google Scholar]
 VargasDeLeón C . On the global stability of SIS, SIR and SIRS epidemic models with standard incidence [J]. Chaos, Solitons and Fractals, 2011, 44(12): 11061110. [Google Scholar]
 Li G H, Jin Z. Global stability of an SEI epidemic model [J]. Chaos, Solitons and Fractals, 2004, 21(4): 925931. [Google Scholar]
 Chen Q L, Teng Z D, Wang L, et al. The existence of codimensiontwo bifurcation in a discrete SIS epidemic model with standard incidence[J]. Nonlinear Dynamics, 2013, 71(1): 5573. [CrossRef] [Google Scholar]
 Ma Y L, Chu Z Q, Li H J. Global dynamics of an SEIQR model with saturation incidence rate and hybrid strategies [J]. Journal of University of Science and Technology of China, 2021, 51 (2): 153163(Ch). [Google Scholar]
 Amador J . The SEIQS stochastic epidemic model with external source of infection [J]. Applied Mathematical Modelling, 2016, 40(1920): 83528365. [CrossRef] [MathSciNet] [Google Scholar]
 Silva C M . A nonautonomous epidemic model with general incidence and isolation[J]. Mathematical Methods in the Applied Sciences, 2014, 37 (13):19741991. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Bai Z G . Global dynamics of a SEIR model with information dependent vaccination and periodically varying transmission rate[J]. Mathematical Methods in the Applied Sciences, 2015, 38(11): 24032410. [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Thieme H R . Convergence results and a poincarebendi xson trichotomy for asymptotically autonomous differential equations [J]. Journal of Mathematical Biology, 1992, 30(7): 755763. [CrossRef] [MathSciNet] [Google Scholar]
 Li T T, Xue Y K. Global stability analysis of a delayed SEIQR epidemic model with quarantine and latent [J]. Applied Mathematics, 2013, 4(10):109117. [CrossRef] [Google Scholar]
All Figures
Fig.1 Variational curves of S, E, I and R with time t when R
_{0}=0.475<1 Six initial conditions are (a) (5, 1.9, 3.2, 2.6), (b) (3.8,0.6,2,5.4), (c) (2.1, 5, 2.6, 5), (d) (3.5, 2.5, 0.8, 4.6), (e) (0.5, 4, 2.2, 3.6) and (f) (4.5, 1.6, 2.4, 3.1), respectively 

In the text 
Fig.2 Variational curves of S, E, I and R with time t when R
_{0}=5.508>1 Six initial conditions are (a) (11.5, 9.2, 2.7, 3.8), (b) (9.8, 6.6, 8, 11.4), (c) (8, 7.5, 8.6, 11), (d) (6.5, 10, 8.2, 9.1), (e) (10.5, 7.6, 8.4, 9.1) and (f) (5.8, 9.6, 7.4, 12), respectively 

In the text 
Fig.3 The graph of the trajectory in (S, E, I)space (a) and (b) correspond with Fig.1 and Fig.2, respectively, and the six curves correspond to six initial conditions 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.