Abstract
A method is described for iteratively solving the Boltzmann equation with the angular dependence of the distribution function expanded in Legendre polynomials. Compared with earlier integral equation methods this approach is shown to have a number of advantages which can lead to an increase in computing efficiency by as much as two orders of magnitude. Numerical results for p-Ge and n-GaAs have been calculated to show the convergence properties of the present approach.