 Solution

Expand the radial modified Helmholtz equation as $Ly=-y^{\prime \prime }-2y^{\prime }/x+y$. Setting $y=z/x$ gives $Ly=(-z^{\prime \prime }+z)/x$, soBoundedness at zero requires $A+B=0$, whereas boundedness at infinity requires $A=0$. The respective adapted homogeneous solutions are thereforeThe first has limit one at zero. Their Wronskian is $u{v}^{\prime }-u^{\prime }v=-x^{-2}$. Since the highest-derivative coefficient of $L$ is $-1$, the Green function must have first-derivative jump $-1$. The radial modified Helmholtz Green function on the half-line isIt is continuous at the source, and its derivative jump is $-1$, giving $L_{x}G=\delta (x-\xi )$. For fixed source it tends to $\xi e^{-\xi }$ as $x\downarrow 0$ and decays at infinity. This kernel is for the original operator and ordinary measure $d\xi$; multiplying the equation by $x^2$ would change the source normalization.

For the compact source, split the integral at $\xi =x$ when $x<R$. The needed antiderivatives areThuswhile for $x>R$ the result is $e^{-x}{x}^{-1}{\int }_0^R\xi \sinh \xi d\xi$. Simplifying givesThe continuous value at zero is $1-(R+1)e^{-R}$, and $y^{\prime }(0)=0$. Both $y$ and $y^{\prime }$ match at $R$, so no extra Dirac delta appears there. The only homogeneous solution bounded at both ends is zero, establishing uniqueness.