The convergence and stability of full discretization scheme for stochastic age-structured population models

Abstract

In this paper, a fully discretization scheme based on the implicit Euler method (IM) is considered for stochastic age-structured population models. The preservation of the total population with a suitable numerical boundary condition according to the biological meanings are shown. An explicit formula of the numerical basic reproductive number Rh is proposed by the technique that the numerical process is embedded into an l1(R)-valued integrable stochastic process with infinite stochastic Leslie operators. Furthermore, the convergence and connection between Rh and the stability of numerical solution is analyzed. The preservation and detection of the analytic stability through the numerical solutions are discussed for small stepsize. Finally, some numerical experiments including an infection-age model for modified SARS epidemic illustrate the verification and efficiency of our analysis.