 Solution

Multiplication by $x^2$ puts the eigenvalue problem into Sturm-Liouville form, $-(x^2{y}^{\prime }{)}^{\prime }=\lambda x^{2}y$, with weight $x^2$. If $u,v$ correspond to distinct eigenvalues $\lambda ,\mu$, subtract their multiplied equations and integrate to getThe upper boundary term vanishes by the Dirichlet boundary condition. For a bounded solution near zero, $x^2{y}^{\prime }$ has a limit because $(x^2{y}^{\prime }{)}^{\prime }=-\lambda x^{2}y$ is integrable. A nonzero limit would force a $1/x$ divergence in $y$, so that limit is zero; integration then gives $y^{\prime }=O(x)$. The lower boundary term therefore also vanishes, proving the weighted orthogonality.

Set $z=xy$. For $x>0$, $y^{\prime \prime }+2y^{\prime }/x=z^{\prime \prime }/x$, and finiteness of $y(0)$ gives $z(0)=0$. Thus $z^{\prime \prime }+\lambda z=0$ with zero values at zero and $\pi$. The radial Dirichlet spectrum in three dimensions givesThe zero and negative eigenvalues yield only the zero solution under these boundary conditions. Also ${\int }_0^{\pi }{x}^2{y}_n^{2}dx=\pi /2$.

For the forced equation, the same substitution gives $z^{\prime \prime }-\beta z=x$. The Fourier sine series of $x$ on $(0,\pi )$ has coefficientsExpanding $z=\sum c_n\sin nx$ gives $-(n^2+\beta )c_n=b_n$, henceThe value at zero is defined by a limit. In fact $|\sin nx/x|\le n$, so uniform convergence holds even there, by comparison with $2\sum (n^2+\beta {)}^{-1}$. An independent check comes from solving the constant-coefficient equation directly:It satisfies both endpoint conditions and the forcing for $x>0$, and its finite origin limit agrees with the series. For $\beta >0$ the homogeneous boundary value problem has no nonzero solution, so the solution is unique.