On positive solutions in a phytoplankton–nutrient model

On the basis of an application from aquatic ecology, we discuss the behaviour of the widely used time integration package VODE by Brown et al. (SIAM J. Sci. Statist. Comput. 10 (1989) 1038). When used in a default setting this code smoothly produces a negative steady state solution, which is not realistic in this application.