Non-negative solutions of ODEs

This paper discusses procedures for enforcing non-negativity in 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.