In the mathematical modeling of population growth, and in particular of bacterial growth, parameters are either measured directly or determined by curve fitting. These parameters have large variability that depends on the experimental method and its inherent error, on differences in the actual population sample size used, as well as other factors that are difficult to account for. In this work the parameters that appear in the Monod kinetics growth model are considered random variables with specified distributions. A stochastic spectral representation of the parameters is used, together with the polynomial chaos method, to obtain a system of differential equations, which is integrated numerically to obtain the evolution of the mean and higher-order moments with respect to time. (C) 2009 Elsevier Ltd. All rights reserved.