Dynamically Consistent NSFD Methods for Predator-prey System

Document Type : Research Paper

Authors

1 Department of Mathematics, Faculty of Science University of Maragheh, Maragheh, Iran‎

2 Department of Mathematics, Faculty of Science University of Maragheh, Maragheh, Iran

Abstract

In this paper, we introduce two nonstandard finite difference (NSFD) methods for solving the mathematical model of the Rosenzweig-MacArthur predator-prey system. These new proposed numerical methods have important features such as positivity and elementary stability. Numerical comparisons between the proposed methods and the other methods such as second-order and forth order Runge-Kutta methods (we refer them RK2 and RK4, respectively), Euler method, and NSFD method presented in [6] indicate that the new methods have better accuracy and convergence.

Keywords

Main Subjects

[1] Dimitrov, D.T., Kojouharov, H.V., Analysis and numerical simulation of phytoplankton-nutrient systems with nutrient loss, Mathematics and Computers in Simulation, 70(1), 2005, 33-43.
[2] Dimitrov, D.T., Kojouharov, H.V., Positive and elementary stable nonstandard numerical methods with applications to predator-prey models, Journal of Computational and Applied Mathematics, 189(1-2), 2006, 98-108.
[3] Dimitrov, D.T., Kojouharov, H.V., Nonstandard numerical methods for a class of predator-prey models with predator interference, Electronic Journal of Differential Equations (EJDE), 2007, 2007, 67-75.
[4] Dumont, Y., Chiroleu, F., Domerg, C., On a temporal model for the Chikungunya disease:modeling theory and numerics, Mathematical Biosciences, 213(1), 2008, 80-91.
[5] Mehdizadeh Khalsaraei, M., Positivity of an explicit Runge-Kutta method, Ain Shams Engineering Journal, 6(4), 2015, 1217-1223.
[6] Mehdizadeh Khalsaraei, M., Khodadoosti, F., Qualitatively stability of nonstandard 2-stage explicit Runge-Kutta methods of order two, Computational Mathematics and Mathematical Physics, 56(2), 2016, 235-242.
[7] Mickens, R.E., A nonstandard finite-difference scheme for the Lotka-Volterra system, Applied Numerical Mathematics, 45(2-3), 2003, 309-314.
[8] Mickens, R.E., Jordan, P.M., A positivity-preserving nonstandard finite difference scheme for the damped wave equation, Numerical Methods for Partial Differential Equations, 20(5), 2004, 639-649.
[9] Anguelov, R., Dumontb, Y., Lubumaa, J.M.S., Shillor, M., Dynamically consistent nonstandard finite difference schemes for epidemiological models, Journal of Computational and Applied Mathematics, 255, 2014, 161-182.
[10] Ashyralyev, A., Agirseven, D., Agarwal, R.P., Stability estimates for delay parabolic differential and difference equations, Applied and Computational Mathematics, 19(2), 2020, 175-204.
[11] Gonzalez-Parra, G., Arenas, A.J., Cogollo, M.R., Numerical-analytical solutions of predator-prey models, WSEAS Transactions on Biology and Biomedicine, 10(2), 2013, 79-87.
[12] Gurski, K.F., A simple construction of nonstandard finite-difference schemes for small nonlinear systems applied to SIR models, Computers & Mathematics with Applications, 66(11), 2013, 2165-2177.
[13] Mehdizadeh Khalsaraei, M., An improvement on the positivity results for 2-stage explicit Runge-Kutta methods, Journal of Computational and Applied Mathematics, 235(1), 2010, 137-143.
[14] Mehdizadeh Khalsaraei, M., Khodadoosti, F., A new total variation diminishing implicit nonstandard finite difference scheme for conservation laws, Computational Methods for Differential Equations, 2(2), 2014, 85-92.
[15] Mehdizadeh Khalsaraei, M., Shokri Jahandizi, R., Positivity-preserving nonstandard finite difference schemes for simulation of advection-diffusion reaction equations, Computational Methods for Differential Equations, 2(4), 2014, 256-267.
[16] Mehdizadeh Khalsaraei, M., Shokri, A., Ramos, H., Heydari, S., A positive and elementary stable nonstandard explicit scheme for a mathematical model of the influenza disease, Mathematics and Computers in Simulation, 182, 2021, 397-410.
[17] Bacaer, N., Ouifki, R., Pretorius, C., Wood, R., Williams, B., Modeling the joint epidemics of TB and HIV in a south African township, Journal of Mathematical Biology, 57(4), 2008, 557-593.
[18] Brauer, F., Castillo-Chavez, C., Mathematical Models in Population Biology and Epidemiology, Springer, New York,2001.
[19] Golbabai, A., Nikan O., Molavi-Arabshahi, M., Numerical approximation of time fractional advection-dispersion model arising from solute transport in rivers, TWMS Journal of Pure and Applied Mathematics, 10(1), 2019, 117-131.
[20] Moshtaghi, N., Saadatmandi, A., Numerical solution for diffusion equations with distributed-order in time based on sinc-Legendre collocation method, Applied and Computational Mathematics, 19(3), 2020, 317-355.
[21] Murry, J.D., Mathematical Biology, Springer-Verlag, 1989.
[22] Qalandarov, A.A., Khaldjigitov, A.A., Mathematical and numerical modeling of the coupled dynamic thermoelastic problems for isotropic bodies, TWMS Journal of Pure and Applied Mathematics, 11(1), 2020, 119-126.
[23] Anguelov, R., Lubuma, J.M.S., Contributions to the mathematics of the nonstandard finite difference method and applications, Numerical Methods for Partial Differential Equations, 17(5), 2001, 518-543.
[24] Dimitrov, D.T., Kojouharov, H.V., Stability-Preserving finite-difference methods for general multi-dimensional autonomous dynamical system, International Journal of Numerical Analysis and Modeling, 4(2), 2007, 282-292.
[25] Mickens, R.E., Nonstandard Finite Difference Models of Differential Equations, World Scientific, Singapore, 1994.
[26] Portalier, S.M., Fussmann, G.F., Loreau, M., Cherif, M., The mechanics of predator–prey interactions: First principles of physics predict predator-prey size ratios, Functional Ecology, 33(2), 2019, 323-334.