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

This paper introduces a saturated treatment function into an SIS model with saturated incidence rate. The treatment function is a continuous and differential function which describes the effect of delayed treatment when the medical condition is limited and the number of infected individuals is large. Sufficient conditions for the existence and global asymptotical stability of the disease-free and endemic equilibria are given in this paper, and the nonexistence of limit cycles is also demonstrated. A backward bifurcation is found when the capacity of the treatment is low. It indicates that we should improve the efficiency and enlarge the capacity of the treatment to control the spread of diseases. By mathematical analysis and numerical simulations, it is shown that the system undergoes Hopf bifurcation and Bogdanov–Takens bifurcation.