The Pauli principle is included in a deterministic Boltzmann solver based on a spherical harmonics expansion of the distribution function. The Newton-Raphson scheme is applied to solve the nonlinear Boltzmann equation, and it is found that the inclusion of the Pauli principle introduces no numerical problems, even for multi-dimensional semiconductor devices. As a numerical example, the impact of the Pauli principle is numerically investigated for a double-gate MOSFET.