The traditional probability density evolution equations of stochastic systems are usually in high dimensions. It is very hard to obtain the solutions. Recently the development of a family of generalized density evolution equation (GDEE) provides a new possibility of tackling nonlinear stochastic systems. In the present paper, a numerical method different from the finite difference method is developed for the solution of the GDEE. In the proposed method, the formal solution is firstly obtained through the method of characteristics. Then the solution is approximated by introducing the asymptotic sequences of the Dirac δ function combined with the smart selection of representative point sets in the random parameters space. The implementation procedure of the proposed method is elaborated. Some details of the computation including the selection of the parameters are discussed. The rationality and effectiveness of the proposed method is verified by some examples. Some features of the numerical results are observed.