Abstract: | A numerical method incorporating a combination of a difference scheme and several uniform and nonuniform quadrature rules is presented. The method is designed to solve size‐structured population equations with linear growth rate and nonlinear fertility and mortality rates. A detailed analysis of the global discretization error is carried out. Examples with known exact solutions have been solved numerically using the proposed method. The computations show that the global error is of third order as predicted by the theory. © 2002 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 1–21, 2003 |