An integral relation for tensor polynomials

We prove two lemmas and one theorem that allow integrating the product of an arbitrary number of unit vectors and the Legendre polynomials over a sphere of arbitrary radius. Such integral tensor products appear in solving inhomogeneous Helmholtz equations whose right-hand side is proportional to the product of a nonfixed number of unit vectors.