 Solution

Every orbital angular momentum operator annihilates a radial function: substituting ${\partial }_{i}f(r)=f^{\prime }(r)x_i/r$ makes the two terms cancel. It also annihilates $r^2$. Thus the radial factor can be carried through the angular operators. Since $x_3^2$ is independent of $x_1,x_2$, $L_3\psi =0$ for every value of $a$.

Direct differentiation givesConsequently $L^2{x}_3^2={\hslash }^2(6x_3^2-2r^2)$ and $L^2{r}^2=0$. For a nonzero admissible radial factor, comparison with $\lambda (x_3^2+a{r}^2)f(r)$ forces $\lambda =6{\hslash }^2$ and $6a=-2$. Therefore the simultaneous eigenfunction occurs atIt is the axially symmetric quadratic orbital eigenfunction, with quantum numbers $l=2$, $m=0$ and angular dependence proportional to the spherical harmonic $Y_{20}$. The subtraction of $r^2/3$ removes its rotationally invariant scalar component.