Existence of multiple periodic solutions for an SIR model with seasonality

We study an SIR model with a seasonal contact rate and a staged treatment strategy, which is an extension of our previous work Z. Bai, Y. Zhou, Existence of two periodic solutions for a non-autonomous SIR epidemic model, Appl. Math. Model. 35 (2011) 382-391]. It is proved that the persistence and extinction of the disease are determined by the basic reproductive number (R₀) and a threshold parameter (R_c). It is obtained that the model exhibits two different bistable behaviors under certain conditions: the stable disease-free state coexists with a stable endemic periodic solution, and three endemic periodic solutions coexist with two of them being stable. Numerical simulations are presented to illustrate theoretical results.