Discretization method and numerical calculations of wigner function to introduce the influence of spatially varying effective mass as well as to reduce the error involved in the conventional discretization model are presented. Using this new discrete model, the transient responses of resonant-tunneling-diode are analyzed.