We consider a model of population dynamics whose mortality function is unbounded. We approximate the solution of the model using a discontinuous Galerkin finite element for the age variable and a backward Euler for the time variable. We present several numerical examples. It is experimentally shown that the scheme converges at the rate of $h^{3/2}$ in the case of piecewise linear polynomial space.