A Monte Carlo method for calculation of the carrier mobility in degenerate bulk semiconductors at zero electric field is presented. The method is obtained as a limiting case of an existing small-signal approach replacing the distribution function by the Fermi-Dirac distribution which is valid at zero electric field. The general form of the Boltzmann equation which takes into account the Pauli exclusion principle in the scattering term is used to derive the integral representation of a Boltzmann-like equation for a small perturbation of the distribution function. The method allows calculation of the whole mobility tensor in comparison with the one particle Monte Carlo algorithm which is traditionally used to compute low field carrier mobility.