Bifurcation analysis of an SIRS epidemic model with standard incidence rate and saturated treatment function

An epidemic model with standard incidence rate and saturated treatment function of infectious individuals is proposed to understand the effect of the capacity for treatment of infective individuals on the disease spread. The treatment function in this paper is a continuous and differential function which exhibits the effect of delayed treatment when the rate of treatment is lower and the number of infected individuals is getting larger. It is proved that the existence and stability of the disease-free and endemic equilibria for the model are not only related to the basic reproduction number but also to the capacity for treatment of infective individuals. And a backward bifurcation is found when the capacity is not enough. By computing the first Lyapunov coefficient, we can determine the type of Hopf bifurcation, i.e., subcritical Hopf bifurcation or supercritical Hopf bifurcation. We also show that under some conditions the model undergoes Bogdanov-Takens bifurcation. Finally, numerical simulations are given to support some of the theoretical results.