In this paper, we construct explicit nonstandard Runge-Kutta (ENRK) methods which have higher accuracy order and preserve two important properties of autonomous dynamical systems, namely, the positivity and linear stability. These methods are based on the classical explicit Runge-Kutta methods, where instead of the usual h in the formulas there stands a function ϕ(h). It is proved that the constructed methods preserve the accuracy order of the original Runge-Kutta methods. The numerical simulations confirm the validity of the obtained theoretical results. difference schemes; Positive nonstandard finite difference methods; Elementary stable. Our first objective is to construct difference schemes preserving the linear stability of the equilibrium points of System (1) for all finite step-size h > 0. These schemes are called also elementary stable [3,7,8]. It should be emphasized that standard finite difference schemes cannot preserve properties of the differential equations for any step-sizes h > 0, including the linear stability. Mickens called this phenomenon numerical instability [25].The construction of elementary stable difference schemes play especially important role in numerical solution of differential equations and numerical simulation of nonlinear dynamical systems. In 2005, Dimitrov and Kojouharov [7] proposed a method for constructing elementary stable NSFD methods for general two-dimensional autonomous dynamical systems. These NSFD methods are based on the explicit and implicit Euler and the second-order Runge-Kutta methods. Later, in 2007 these results are extented for the general n-dimensional dynamical systems, namely, NSFD schemes preserving elementary stability are constructed based on the θ-methods and the second-order Runge-Kutta methods [8]. One important action in the construction of the elementary stable NSFD schemes is the replacement of the standard denominator function ϕ(h) = h by the nonstandard denominator function, which is Email addresses: dangquanga@cic.vast.vn (Quang A Dang),