This paper presents a numerical method for solving a nonlinear age-structured population model based on a set of piecewise constant orthogonal functions. The block-pulse functions (BPFs) method is applied to determine the numerical solution of a non-classic type of partial differential equation with an integral boundary condition. BPFs duo to the simple structure can efficiently approximate the solution of systems with local or non-local boundary conditions. Numerical results reveal the accuracy of the proposed method even for the long term simulations.