A numerical method is proposed and analyzed to approximate the solution of a math-ematical model of age-dependent population dynamics with spatial diffusion, with a form of nonlinear and nonlocal system of integro-differential equations. A finite difference method along the characteristic age-time direction is presented and mixed finite elements in the spatial variable is used for the approximation. The error analysis of this integro-differential equation for mixed nonconforming finite element is present for the fully discrete scheme in this paper, and the optimal error estimates are derived. |