Stability Analysis of an SIR Epidemic Model with Non-Linear Incidence Rate and Treatment

We consider a SIR epidemic model with saturated incidence rate and treatment. We show that if the basic reproduction number, R0 is less than unity and the disease free equilibrium is locally asymptotically stable. Moreover, we show that if R0 > 1, the endemic equilibrium is locally asymptotically stable. In the end, we give some numerical results to compare our model with existing model and to show the effect of the treatment term on the model.


Introduction
While mathematical modeling of infectious diseases could be traced back to 1760 when Bernoulli used mathematical models for small pox [1], the research in infectious diseases, using deterministic mathematical models actually began in 20 th century.
Bilinear and standard incidence rates have been frequently used in classical epidemic models [2].Several different incidence rates have been proposed by researchers.After a study of the cholera epidemic spread in Bari in 1973, Capasso and Serio [3] introduced a saturated incidence rate g(I)S into epidemic models.Ruan and Wang [4] studied an epidemic model with specific nonlinear incidence rate of the form 1 l h KI S I α + where l and h were positive constants and obtained lots of interesting dynamical behaviour of the model such as a limit cycle, two limit cycles and homoclimic loop etc.In 2000, van den Driessche and Watmough [5] studied an SIS epidemic model with the incidence rate of the form ( ) Xiao and Ruan [6]   considered a special nonlinear incident rate of the form Jasmine and Amirtharaj [7] on account of the effect of limited treatment resources on the control of epidemic disease incorporated a modified SIR epidemic model with generalized incidence rate.They discussed the stability analysis of the disease-free equilibrium and endemic equilibrium with a nonlinear incidence rate.Chauhan et al. [8]  > .We also look at the effect of the transmission rate on the SIR epidemic model and we include the treatment term I τ where , 0

The Basic Mathematical Model
We modify the model of Pathak et al., by extending the transmission rate to nonlinear orders and also include the treatment term.

Modified Model
The modified model is as follows: ( ) ( ) where the transmission rate 1 where ( ) ( ) ( ) , and S t I t R t represent the number of susceptible, infective and recovered individuals at time t, respectively, b is the recruitment rate of the population, d is the natural death rate of the population, k is the proportionality constant, µ is the natural recovery rate of the infective individuals, γ is the rate at which re- covered individuals lose immunity and return to the susceptible class, and α β are the parameters which measure the effects of sociological, psychological or other mechanisms and p and q are positive constants which are greater than unity.

Steady State and Local Stability of the Critical Points
In this section, we discuss the local stability of the disease-free equilibrium and endemic equilibrium of system (2).
The system (2) has a disease free of ( ) admits a unique endemic equilibrium ( ) , , E S I R Proposition 1: If R 0 < 1, then the disease free equilibrium E 0 is locally asymptotically stable.Proof: Obviously, (3) has three roots 1 d ,

Stability Analysis of the Endemic Equilibrium
Let the endemic equilibrium  A straight forward calculation yields 2 1 3 0. a a a − > Then, it follows from Routh-Hurwitz criteria that all characteristics roots have negative real parts.Thus, the endemic equilibrium is locally asymptotically stable.

Numerical Simulations
To see the dynamical behavior of system (2), we solve the system by Runge-Kutta Felhberg 45(RKF 45) method using the parameters;

Conclusion
In this paper, we have carried out the stability of the equilibrium states using some of the tested parameters from literature reviewed in this paper.The simulation is carried out using numerical software called "maple".The effect of the treatment term in the model has been investigated and it shows that treatment has a positive effect on the endemic nature of the disease.The more the treatment is applied, the faster the disease fades out.
I) is non-monotone.
discussed the stability analysis pf SIR epidemic model with and without vaccination.