A Nonstandard Dynamically Consistent Numerical Scheme Applied to Obesity Dynamics

The obesity epidemic is considered a health concern of paramount importance in modern society. In this work, a nonstandard finite difference scheme has been developed with the aim to solve numerically a mathematical model for obesity population dynamics. This interacting population model represented as a system of coupled nonlinear ordinary differential equations is used to analyze, understand, and predict the dynamics of obesity populations. The construction of the proposed discrete scheme is developed such that it is dynamically consistent with the original differential equations model. Since the total population in this mathematical model is assumed constant, the proposed scheme has been constructed to satisfy the associated conservation law and positivity condition. Numerical comparisons between the competitive nonstandard scheme developed here and Euler's method show the effectiveness of the proposed nonstandard numerical scheme. Numerical examples show that the nonstandard difference scheme methodology is a good option to solve numerically different mathematical models where essential properties of the populations need to be satisfied in order to simulate the real world.

[1]  Edward H. Twizell,et al.  An unconditionally convergent discretization of the SEIR model , 2002, Math. Comput. Simul..

[2]  Abba B. Gumel.  A competitive numerical method for a chemotherapy model of two HIV subtypes , 2002, Appl. Math. Comput..

[3]  Edward H. Twizell,et al.  An unconditionally convergent finite-difference scheme for the SIR model , 2003, Appl. Math. Comput..

[4]  AnguelovRoumen,et al.  Nonstandard finite difference method by nonlocal approximation , 2003 .

[5]  Ronald E. Mickens,et al.  Calculation of denominator functions for nonstandard finite difference schemes for differential equations satisfying a positivity condition , 2007 .

[6]  R. Mickens.  Nonstandard Finite Difference Models of Differential Equations , 1993 .

[7]  Hristo V. Kojouharov,et al.  Combined nonstandard numerical methods for ODEs with polynomial right-hand sides , 2006, Math. Comput. Simul..

[8]  J. Lambert.  Numerical Methods for Ordinary Differential Systems: The Initial Value Problem , 1991 .

[9]  F. Brauer,et al.  Mathematical Models in Population Biology and Epidemiology , 2001 .

[10]  Ronald E. Mickens,et al.  Dynamic consistency: a fundamental principle for constructing nonstandard finite difference schemes for differential equations , 2005 .

[11]  Ronald E Mickens,et al.  Numerical integration of population models satisfying conservation laws: NSFD methods , 2007, Journal of biological dynamics.

[12]  Kailash C. Patidar,et al.  CONTRIBUTIONS TO THE THEORY OF NON-STANDARD FINITE DIFFERENCE METHODS AND APPLICATIONS TO SINGULAR PERTURBATION PROBLEMS , 2008 .

[13]  Roumen Anguelov,et al.  Nonstandard finite difference method by nonlocal approximation , 2003, Math. Comput. Simul..

[14]  Roumen Anguelov,et al.  Contributions to the mathematics of the nonstandard finite difference method and applications , 2001 .

[15]  Lucas Jódar,et al.  Modeling dynamics of infant obesity in the region of Valencia, Spain , 2008, Comput. Math. Appl..

[16]  Murray E. Alexander,et al.  A non-standard numerical scheme for a generalized Gause-type predator–prey model , 2004 .

[17]  Hristo V. Kojouharov,et al.  Nonstandard finite-difference schemes for general two-dimensional autonomous dynamical systems , 2005, Appl. Math. Lett..

[18]  Murray E. Alexander,et al.  A Positivity-preserving Mickens-type Discretization of an Epidemic Model , 2003 .

[19]  K. Patidar.  On the use of nonstandard finite difference methods† , 2005 .