That's a bit strange coordinate system. The usual physicist's spherical coordinates ##(r,\vartheta,\varphi)## have the domains ##r>0##, ##0<\vartheta<\pi##, and ##0\leq \varphi<2 \pi##. They are defined by the transformation of the position vector from spherical to Cartesian components
$$(x,y,z)=r (\cos \varphi \sin \vartheta,\sin \varphi \sin \vartheta,\cos \vartheta).$$
The conventional way to address a spherically symmetric problem is to write the Laplacian in spherical coordinates and then make an separation ansatz
$$\psi(r,\vartheta,\varphi)=\frac{1}{r} R(r) \Theta(\vartheta) \Phi(\varphi).$$
To take out the factor ##1/r## is for convenience. The boundary conditions are not so trivial to determine as it seems on the first glance. They boil down to the fact that without loss of generality you can assume that the wave function is a scalar under rotations and that the radial Schrödinger equation should be governed by an effective Hamiltonian that is self-adjoint on ##\mathrm{L}^2(0,\infty)## (where ##R(r)## lives). Particularly this yields ##\psi(r,\vartheta,\varphi+2 \pi)=\psi(r,\vartheta,\varphi)##. From the self-adjointness of the Hamiltonian you also get that ##R(0)=0## as a boundary condition. Further the solutions for ##\Theta## should be square integrable with the weight ##\sin \vartheta##. This leads to the set ##\Theta(\vartheta) \Phi(\varphi)=Y_{lm}(\vartheta,\varphi) \propto \exp(\mathrm{i} m \varphi)##, where for each ##l \in \{0,1,2,\ldots \}## the ##m \in \{-l,-l+1,\ldots,l-1,l \}##. The radial solution should also stay finite (go to 0) with ##r \rightarrow \infty## for the scattering (bound states). With these conditions you also get the radial solutions for the case of the hydrogen atom. They are given by the Laguerre polynomials times an exponential function (for the bound states). A very good treatment, simplifying a lot by making use of the angular-momentum algebra to determine the spherical harmonics, can be found in the textbook by Messiah, vol. I.