The Lotka–Volterra predator–prey system x′ = x − xy, y′ = − y+xy is a good differential equation system for testing numerical methods. This model gives rise to mutually periodic solutions surrounding the positive fixed point (1,1), provided the initial conditions are positive. Standard finite-difference methods produce solutions that spiral into or out of the positive fixed point. Previously, the author [Roeger, J. Diff. Equ. Appl. 12(9) (2006), pp. 937–948], generalized three different classes of nonstandard finite-difference methods that when applied to the predator–prey system produced periodic solutions. These methods preserve weighted area; they are symplectic with respect to a noncanonical structure and have the property that the computed points do not spiral. In this paper, we use a different approach. We apply the Jacobian matrix procedure to find a fourth class of nonstandard finite-difference methods. The Jacobian matrix method gives more general nonstandard methods that also produce periodic solutions for the predator–prey model. These methods also preserve the positivity property of the solutions.
Read full abstract