Efficient Algorithms for Solving Systems of Ordinary Differential Equations for Ecosystem Modeling

摘要

This report presents three packages of subroutines for the numerical solution of systems of first order ordinary differential equations. The three integration methods were chosen to provide an efficient means of solving the full range of initial value problems encountered in basin ecosystem modeling. The subroutines have been designed to handle all aspects of the numerical integration as automatically as possible so that their successful use does not require a sophisticated knowledge of the techniques of numerical analysis. After a general introduction to the principles and concepts of numerical integration, the particular methods used in the codes are discussed in some detail. Full instructions for use are given, and some simple examples are presented. The complete code listings, in ANSI Fortran IV, are contained in appendices.