A class of nonstandard numerical methods for autonomous dynamical systems

A nonstandard finite difference method for solving autonomous dynamical systems is constructed. The proposed numerical method is computationally efficient and easy to implement. It is designed so that it preserves positivity of solutions and the local behavior of the dynamical system near equilibria.