Abstract
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.
Disclosure statement
No potential conflict of interest was reported by the authors.