A nonstandard numerical scheme of predictor-corrector type for epidemic models

In this paper we construct and develop a competitive nonstandard finite difference numerical scheme of predictor-corrector type for the classical SIR epidemic model. This proposed scheme is designed with the aim of obtaining dynamical consistency between the discrete solution and the solution of the continuous model. The nonstandard finite difference scheme with Conservation Law (NSFDCL) developed here satisfies some important properties associated with the considered SIR epidemic model, such as positivity, boundedness, monotonicity, stability and conservation of frequency of the oscillations. Numerical comparisons between the NSFDCL numerical scheme developed here and Runge-Kutta type schemes show its effectiveness.

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

[2]  Francisco J. Solis,et al.  Nonstandard discrete approximationspreserving stability properties of continuous mathematical models , 2004, Math. Comput. Model..

[3]  Juan Carlos Cortés,et al.  Non-standard numerical method for a mathematical model of RSV epidemiological transmission , 2008, Comput. Math. Appl..

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

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

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

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

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

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

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

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

[12]  JódarLucas,et al.  Nonstandard numerical methods for a mathematical model for influenza disease , 2008 .

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

[14]  J. Lambert.  Computational Methods in Ordinary Differential Equations , 1973 .

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

[16]  R. Mickens.  Applications of nonstandard finite difference schemes , 2000 .

[17]  S. Duncan,et al.  Whooping cough epidemics in London, 1701-1812: infecdon dynamics, seasonal forcing and the effects of malnutrition , 1996, Proceedings of the Royal Society of London. Series B: Biological Sciences.

[18]  Abba B. Gumel,et al.  A Non-Standard Finite-Difference Scheme for a Model of HIV Transmission and Control , 2003 .

[19]  Rafael J. Villanueva,et al.  Nonstandard numerical methods for a mathematical model for influenza disease , 2008, Math. Comput. Simul..

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

[21]  Hristo V. Kojouharov,et al.  Positive and elementary stable nonstandard numerical methods with applications to predator-prey models , 2006 .

[22]  Rafael J. Villanueva,et al.  A Nonstandard Dynamically Consistent Numerical Scheme Applied to Obesity Dynamics , 2008, J. Appl. Math..