A simple second-order nonstandard numerical method for a general class of dynamical systems and its applications

In this work, we consider a class of continuous-time autonomous dynamical systems that model various important phenomena and processes encountered in real-world situations. We construct a second-order nonstandard finite difference (NSFD) method that simultaneously preserves two essential properties of the dynamical systems for all finite step sizes, namely the positivity of the solutions, the set of equilibrium points and their asymptotic stability. This NSFD method is constructed based on an appropriate choice of nonstandard denominator functions and a weighted discretization of the right-hand side functions. Under easily-verified conditions, the denominator functions guarantee second-order convergence, whereas the weights ensure the dynamic consistency. By taking advantage of the specific structure of the right-hand side functions, a simple discretization is utilized instead of the nonlocal discretization approaches commonly used in previous works. This simplifies the construction of the proposed NSFD method and, in particular, makes its asymptotic stability analysis easier. As an illustration and an important application, we apply the constructed second-order NSFD method to a well-known two-stage structured species model with recruitment. Consequently, a simple second-order NSFD scheme for the considered two-stage structured species model is derived, improving upon a first-order NSFD scheme constructed in a previous work. Numerical experiments demonstrate the advantages of the second-order NSFD scheme over a standard second-order numerical method, namely, the explicit trapezoidal method. The proposed NSFD method is simple and can be applied to a broad class of dynamical system models arising in both theory and applications. Moreover, it can be readily combined with the Richardson extrapolation technique to improve its accuracy.