On the nonstandard numerical discretization of SIR epidemic model with a saturated incidence rate and vaccination

Abstract

Recently, Hoang and Egbelowo (Boletin de la Sociedad Matemàtica Mexicana, 2020) proposed a nonstandard finite difference scheme (NSFD) to get a discrete SIR epidemic model with saturated incidence rate and constant vaccination. The discrete model was derived by discretizing the right-hand sides of the system locally and the first order derivative is approximated by the generalized forward difference method but with a restrictive denominator function. Their analysis showed that the NSFD scheme is dynamically-consistent only for relatively small time-step sizes. In this paper, we propose and analyze an alternative NSFD scheme by applying nonlocal approximation and choosing the denominator function such that the proposed scheme preserves the boundedness of solutions. It is verified that the proposed discrete model is dynamically-consistent with the corresponding continuous model for all time-step size. The analytical results have been confirmed by some numerical simulations. We also show numerically that the proposed NSFD scheme is superior to the Euler method and the NSFD method proposed by Hoang and Egbelowo (2020).

Keywords

• sir epidemic model
• saturated incidence rate
• local and global stability analysis
• lyapunov function
• dynamically-consistent discretization