We present a numerical method that allows us to solve the transient Bo
ltzmann equation in p-Si with accuracy as high as in previous methods,
but up to 200 times faster than is possible by the direct method. It
consists in expanding the carrier distribution into a set of displaced
Maxwellian functions and using a least-squares-like approximation met
hod.