Abstract
Houston's method of integration over the unit sphere of functions of cubic symmetry is extended by developing integration formulas using 9 and 15 values of the integrand. The 3-, 6-, 9-, and 15-term Houston formulas are used to calculate Debye Θ0 values for a number of cubic crystals. The convergence of the results is excellent except for the most anisotropic crystals.