Abstract
A method is given for the numerical calculation of energy surface integrals within the Brillouin zone like density of states, conductivity, susceptibility, dielectric function etc. The Brillouin zone is divided into tetrahedrons in which the integrand is interpolated linearly. In this way the integration can be done analytically avoiding the histogram method. Several similar methods are discussed with regard to the quotient of accuracy and effort.