Fully symmetric interpolatory rules for multiple integrals over hyper-spherical surfaces

Fully symmetric interpolatory integration rules are constructed for multidimensional integrals over Un, the surface of an n-dimensional hyper-sphere. Explicit formulas for the weights are given for odd degrees 3–13. The new rules are efficient and only moderately unstable. Two randomization methods are described.