The concern here is with Gauss-type quadrature rules that are exact for a m
ixture of polynomials and rational functions, the latter being selected so
as to simulate poles that may be present in the integrand. The underlying t
heory is presented as well as methods for constructing such rational Gauss
formulae. Relevant computer routines are provided and applied to a number o
f examples, including Fermi-Dirac and Bose-Einstein integrals of interest i
n solid state physics.