Lyapunov direct method for investigating stability of nonstandard finite difference schemes for metapopulation models

In this paper nonstandard finite difference (NSFD) schemes of two metapopulation models are constructed. The stability properties of the discrete models are investigated by the use of the Lyapunov stability theorem. As a result of this we have proved that the NSFD schemes preserve essential properties of the metapopulation models (positivity, boundedness and monotone convergence of the solutions, equilibria and their stability properties). Especially, the basic reproduction number of the continuous models is also preserved. Numerical examples confirm the obtained theoretical results of the properties of the constructed difference schemes. The method of Lyapunov functions proves to be much simpler than the standard method for studying stability of the discrete metapopulation model in our very recent paper.