On positive solutions in a phytoplankton–nutrient model

Abstract

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. [2]. When used in a default setting this code smoothly produces a negative steady state solution, which is not realistic in this application.