In this paper, we introduce a rigorous computational approach to prove existence of rotation invariant patterns for a nonlinear Laplace-Beltrami equation posed on the 2-sphere. After changing to spherical coordinates, the problem becomes a singular second order boundary value problem (BVP) on the interval (0,pi/2] with a removable singularity at zero. The singularity is removed by solving the equation with Taylor series on (0,delta] (with delta small) while a Chebyshev series expansion is used to solve the problem on [delta,pi/2]. The two setups are incorporated in a larger zero-finding problem of the form F(a)=0 with a containing the coefficients of the Taylor and Chebyshev series. The problem F=0 is solved rigorously using a Newton-Kantorovich argument.