This paper considers the problem of computing the harmonic expansion of functions defined on the sphere. We begin by proving convolution theorems that relate the convolution of two functions on the sphere to a multiplication in the sprectral domain, as well as the multiplication of two functions on the sphere to a convolution in the spectral domain. These convolution theorems are then used to develop a sampling theorem on the sphere