Numerical Modeling of Some Systems in the Biomedical Sciences

Abstract

Abstract: Finite-difference numerical methods are developed for the solution of some systems in the biomedical sciences; namely, the predator-prey model and the SEIR (Susceptible/Exposed/Infectious/Recovered) measles model. First-order methods are developed to solve the predator-prey model, and a second-order method is developed to solve the SEIR measles model. The predator-prey model is extended to a one-space dimension to incorporate diffusion. The SEIR measles model is extended to a one-space dimension to incorporate (i) diffusion, (ii) convection, and (iii) diffusion-convection. The SEIR measles model is further extended to model diffusion in two-space dimensions The reaction terms in these systems of partial differential equations contain non-linear expressions. Nevetheless, it is seen that the numerical solutions are obtained by solving a linear algebraic system at each time step, as opposed to solving a non-linear algebraic system, which is often required when integrating non-linear partial differential equations. The development of each numerical method was made in the light of experience gained in solving the system of ordinary differential equations for each system The numerical methods proposed for the solution of the initial-value problem for the predator-prey and measles models are characterized as implicit. However, in each case, it is seen that the numerical solutions are obtained explicitly. In a series of numerical experiments, in which the ordinary differential equations are solved first of all, it is seen that the proposed methods have superior stability properties to those of the well-known, first-order Euler method to which they are compared. Incorporating the proposed methods into the numerical solution of partial differential equations is seen to lead to economical and reliable methods for solving the systems