Stochastic age-dependent population equations, one of the important classes of hybrid systems are studied. In general most equations of stochastic age-dependent population do not have explicit solutions. Thus numerical approximation schemes are invaluable tools for exploring their properties. The main purpose of this paper is to develop a numerical scheme and show the convergence of the numerical approximation solution to the analytic solution. In the last section a numerical example is given.