In this study, we study a nonlinear age-structured population models with discontinues mortality and fertility rates, motivated by the fact that different maturation period may cause the significant difference in rates. We develop a novel numerical method with two-layer boundary conditions, the linearly implicit -methods on a special mesh. With a uniform boundedness analysis of numerical solutions, the finite time convergence is proved piecewisely according to the fundamental approach for the smooth rates. For juvenile-adult models, the existence of numerical endemic equilibrium is determined by a numerical basic reproduction function, which converges to the exact one with accuracy of order 1. Moreover, it is shown that for juvenile-adult models, the global stability of the disease-free equilibrium and the local stability of the endemic equilibrium are approximately exhibited by the numerical processes. Finally, some numerical experiments on the Logistic models and tadpoles-frogs models illustrate the verification and the efficiency of our results.