Numerical analysis of a linearly backward Euler method with truncated Wiener process for a stochastic SIS model

The paper deals with the numerical positivity, convergence and dynamical behaviors (including extinction and persistence) for stochastic SIS model. Compared with the existing numerical methods, a linearly backward Euler method with truncated Wiener process is introduced with a less computational cost and a better dynamic behavior. We discuss the numerical positivity by the truncated Wiener process, which is the basis for the investigation of convergence and dynamic behavior. The numerical dynamical behaviors (extinction and persistence) are obtained by an exponential representation for the nonlinear stochastic stability function and the large number theorem for martingale, which reproduces the existing theoretical results of exact solution. Finally, numerical examples are given to validate our numerical results for stochastic SIS model.