A NONSTANDARD FINITE DIFFERENCE METHOD APPLIED TO A MATHEMATICAL CHOLERA MODEL

In this paper, we aim to construct a nonstandard finite difference (NSFD) scheme to solve numerically a mathematical model for cholera epidemic dynamics. We first show that if the basic reproduction number is less than unity, the disease-free equilibrium (DFE) is locally asymptotically stable. Moreover, we mainly establish the global stability analysis of the DFE and endemic equilibrium by using suitable Lyapunov functionals regardless of the time step size. Finally, numerical simulations with different time step sizes and initial conditions are carried out and comparisons are made with other well-known methods to illustrate the main theoretical results.