Stability preserving non-standard finite difference scheme for a harvesting Leslie–Gower predator–prey model

A non-standard finite difference scheme for a harvesting Leslie–Gower equations is constructed. It is shown that the obtained difference system has the same dynamics as the original continuous system, such as positivity of solutions, equilibria and their local stability properties, irrespective of the size of numerical time step. To illustrate the analytical results, we present some numerical simulations.