Non-negative solutions of ODEs

This paper discusses procedures for enforcing non-negativity ill a range of codes for solving ordinary differential equations (ODEs). This codes implement both one-step and multistep methods, all of which use continuous extensions and have event finding capabilities. Examples are given. (c) 2005 Published by Elsevier Inc.