Finite difference numerical methods are developed for the solution system in the biomedical sciences; namely, fox-rabies model. First-order methods and second-order method are developed to solve the fox-rabies equations. The fox-rabies model is extended to one-space dimension to incorporate diffusion. The reaction terms in these systems of partial differential equations contain non-linear expressions. It is seen that the numerical solutions are obtained by solving non-linear algebraic system at each time step, as opposed to solving anon-linear algebraic system which is often required when integrating non-linear partial differential equations. The numerical methods proposed for the solution of the initial-value problem for the fox-rabies model are characterized to be implicit. In each case, however, it 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 seen that the proposed methods have an identical stability properties to those of the well-known, first-order, Euler method. The proposed methods for the numerical solution of partial differential equations are seen to be economical and reliable. Error analysis for the methods, computer implementation and numerical results are discussed. The stability of the numerical method is analyzed using maximum principle analysis.
| Identifer | oai:union.ndltd.org:bl.uk/oai:ethos.bl.uk:247493 |
| Date | January 2002 |
| Creators | Abo Elrish, Mohamed Rasmy |
| Contributors | Twizell, E. H. |
| Publisher | Brunel University |
| Source Sets | Ethos UK |
| Detected Language | English |
| Type | Electronic Thesis or Dissertation |
Page generated in 0.0021 seconds