Abstract
The transient and stationary Boltzmann equation are solved by using a Legendre-Laguerre polynomial expansion for p-Si. Its resolution technique is based on a least-square fitting algorithm. We show this method to be as accurate as the exact ones (particular or direct iterative process); the great interest lies on its fast computation, roughly 200 times faster than a monte-Carlo code for instance, for the same accuracy.