A numerical method is proposed and analyzed to approximate a mathematical model of age-dependent population dynamics with spatial diffusion. The model takes a form of nonlinear and nonlocal system of integro-differential equations. A finite difference method along the characteristic age-time direction is considered and primal mixed finite elements are used in the spatial variable. A priori error estimates are derived for the relevant variables.