Stability Analysis of SIQS Epidemic Model with Saturated Incidence Rate

A SIQS epidemic model with saturated incidence rate is studied. Two equilibrium points exist for the system, disease-free and endemic equilibrium. The stability of the disease-free equilibrium and endemic equilibrium exists when the basic reproduction number R0, is less or greater than unity respectively. The global stability of the disease-free and endemic equilibrium is proved using Lyapunov functions and Poincare-Bendixson theorem plus Dulac’s criterion respectively.


Introduction
The isolation and treatment of symptomatic individuals coupled with the quarantining of individuals that have a high risk of having been infected, constitute two commonly used epidemic control measures.Mass quarantine can inflict significant social, psychological and economic costs without resulting in the detection of many infected individuals.Day et al. [1], Hethcote et al. [2] considered SIQS and SIQR epidemic models with three forms of incidence, which include the bilinear, standard and quarantined-adjusted incidences.
Feng and Thieme [3] considered SEIQR models with arbitrarily distributed periods of infection, including quarantine and a general incidence assumed that all infected individuals go through the quarantine stage and investigated the model dynamics.Settapat and Wirawah [4] discussed the SIQ epidemic model with constant immigration.Yang et al. [5] also studied an SIQ epidemic model with isolation and nonlinear incidence rate.
El-Marouf and Alihaby [6] studied the equilibrium points and their local stability for SIQ and SIQR epidemic models with three forms of incidence rates.They also studied the global stability of the equilibrium by constructing the new forms of Lyapunov functions.
Gbadamosi and Adebimpe [7] investigated an SIQ epidemic model with nonlinear incidence rate.They introduced the concept that describes the present and past states of the disease.
We extended the work of Gbadamosi and Adebimpe [7] to include the rates at which individuals recover and return to susceptible compartment from compartments I and Q respectively and we apply Lyapunov functions and Poincare-Bendixson theorem plus Dulac's criterion to prove the global stability of disease-free and endemic equilibria respectively.

The Model
The model that governs a system of differential equation is presented as follows: Subject to initial conditions ( ) ( ) ( ) The parameters with their descriptions are presented in We can now write

Equilibria
The system (1) has always the disease-free equilibrium at ( ) ( ) , , E S I Q =

Local Stability
In this section, we discussed the local stability of the disease-free equilibrium and endemic equilibrium for the system (1).
We state and prove the following results: Theorem 1: At 0 E , the disease-free equilibrium of the system (1) is locally asymptotically stable when 0 1 R < .
Proof: The Jacobian matrix at the point 0 E through linearization is given by By finding the eigenvalues, we have the following s λ : λ to be negative ( ) R < , the disease-free equilibrium is locally asymptotically stable.Theorem 3.1: The system (1) is locally asymptotically stable at * E if 0 1 R > , otherwise unstable.Proof: At the endemic equilibrium * E , the Jacobian matrix of the system (1) is given by: The characteristic equation of the Jacobian matrix * J is given by a a a > , by Routh Hurwitz criterion, all the eigenvalues of the system (1) has negative real part.There- fore, the endemic equilibrium of the system (1) at * E is locally asymptotically stable.

Global Stability
In this section, we study the global stability of the disease-free equilibrium and endemic equilibrium by Lyapunov function and Poincare-Bendixson theorem respectively.Theorem 3: (Dulac's Criterion) Consider the following general nonlinear autonomous system of de ( ) ( ), If the exists a function it ( ) is not identically zero and does not change sign in E, the system (*) has no close orbit lying entirely in E. if A is an annular region contained in E on which ( ) does not change sign, then there is at most one limit cycle of the system (*) in A.
Theorem 4: (The Poincare-Bendixson Theorem):Suppose that ( ) where E is an open subset of R n and that the system (*) has a rejecting Γ contained in a compact subset f of E. assume that the system (*) has only one unique equilibrium point x 0 in f, then one of the following possibilities holds.
(a) ( ) w Γ is a graphic Theorem 5: The disease-free equilibrium of the model ( 1) is globally asymptotically stable if 0 1 R < Proof: To prove this result, we construct the following Lyapunov function ( ) ( ) where 1 2

3
, and u u u are positive constants to be determined later.Differentiating equation (3) with respect to t, we obtain ( ) ( ) ( ) After rearrangements, we get Thus, the disease-free equilibrium of the system (1) is globally asymptotically stable if 0 1 R < In the next theorem, we present the global stability of the endemic equilibrium of the system (1) at * E Theorem 6: The endemic equilibrium * E of the system (1) is globally asymptotically stable if 0 1 R > .Proof: In order to prove the result, we use Dulac plus Poincare Bendixson theorem as follow ( ) 1 , , where 0, 0, 0.
Hence, by Dulac's criterion, there is no closed orbit in the first quadrant.Therefore, the endemic equilibrium is globally asymptotically stable.

Discussion of Results
The mathematical and stability analysis of SIQS epidemic model with saturated incidence rate and temporary immunity has been presented.We investigated the local stability of the disease-free equilibrium and endemic equilibrium using the basic reproduction number, 0 R .We observed that, when 0 1 R < , the disease-free equili- brium is stable at 0 E locally and endemic equilibrium is unstable which means there is tendency for the disease to die out in the long run.We proved the global stability of the disease free equilibrium and endemic equilibrium of the model using Lyapunov function and Dulac's criterion plus Poincare-Bendixson theorem respectively.

Table 1 .
The t S t S→∞ ≤ ≤And the second equation gives

Table 1 .
Descriptions of parameters.
δ The rate constant for individuals leaving the compartment I for the quarantine compartment Q α Disease-related death rate constant in compartments I and Q γ Rate at which individuals recover and return to Susceptible (S) from compartment I ε Rate at which individuals recover and return to Susceptible (S) from compartment Q