Stability-preserving finite-difference methods for general multi-dimensional autonomous dynamical systems