Combination of nonstandard schemes and Richardson's extrapolation to improve the numerical solution of population models

In this paper we combine nonstandard finite-difference (NSFD) schemes and Richardson's extrapolation method to obtain numerical solutions of two biological systems. The first biological system deals with the dynamics of phytoplankton-nutrient interaction under nutrient recycling and the second one deals with the modeling of whooping cough in the human population. Since both models requires positive solutions, the numerical solutions need to satisfy this property. In addition, it is necessary in some cases that numerical solutions reproduce correctly the dynamical behavior while in other cases it is necessary just to find the steady state. NSFD schemes can do this. In this paper Richardson's extrapolation is applied directly to the NSFD solution to increase the order of accuracy of the numerical solutions of these biological systems. Numerical results show that Richardson's extrapolation method improves accuracy.

[1]  Shigui Ruan,et al.  Persistence and coexistence in zooplankton-phytoplankton-nutrient models with instantaneous nutrient recycling , 1993 .

[2]  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.

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

[4]  T. Anglin.  Physician management of sexually abused children and adolescents. , 1984, Current problems in pediatrics.

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

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

[7]  Tamar Frankel.  [The theory and the practice...]. , 2001, Tijdschrift voor diergeneeskunde.

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

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

[10]  Lawrence F. Shampine,et al.  Global Error Estimates for Ordinary Differential Equations , 1976, TOMS.

[11]  C. Marchant,et al.  The 1993 epidemic of pertussis in Cincinnati. Resurgence of disease in a highly immunized population of children. , 1994, The New England journal of medicine.

[12]  Clarence O. E. Burg,et al.  Application of Richardson extrapolation to the numerical solution of partial differential equations , 2009 .

[13]  J. Nelson.  The changing epidemiology of pertussis in young infants. The role of adults as reservoirs of infection. , 1978, American journal of diseases of children.

[14]  J. Cherry.  The epidemiology of pertussis and pertussis immunization in the United Kingdom and the United States: a comparative study. , 1984, Current problems in pediatrics.

[15]  C. Brezinski,et al.  Extrapolation methods , 1992 .

[16]  Bob W. Kooi,et al.  A second-order, unconditionally positive, mass-conserving integration scheme for biochemical systems , 2006 .

[17]  N. Britton.  Essential Mathematical Biology , 2004 .

[18]  The 1993 epidemic of pertussis in Cincinnati. Resurgence of disease in a highly immunized population of children. , 1994 .

[19]  Lawrence F. Shampine,et al.  Algorithm 504: GERK: Global Error Estimation For Ordinary Differential Equations [D] , 1976, TOMS.

[20]  Claude Brezinski,et al.  Extrapolation methods - theory and practice , 1993, Studies in computational mathematics.

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

[22]  Dynamics of Nutrient-Phytoplankton Interaction in the Presence of Viral Infection and Periodic Nutrient Input , 2008 .

[23]  Ronald E. Mickens,et al.  Nonstandard Finite Difference Schemes for Differential Equations , 2002 .

[24]  C. Bolley,et al.  Conservation de la positivité lors de la discrétisation des problèmes d'évolution paraboliques , 1978 .