A method for the computation of radiation trapping with partial frequency distribution due to the Doppler effect is presented. The method is extremely efficient numerically. It computes the time evolution of the excited-state distribution, using only O(N ln N) multiplications for each time step, where N=N_rM_v.