New approximate formulas for evaluating integrals over the Brillouin zone are reported in this paper. These formulas are formally similar to the Houston type, and can be obtained readily from the Gaussian quadrature formula without solving the matrix equation. Their accuracy is better than that of the Houston’s many-term expansion formulas. An application is made to the problems of the lattice vibration to show their feasibility.