Radial modified Helmholtz equation

For a radial function in three dimensions, $(\Delta -1)y=g(x)$ becomes $y^{\prime \prime }+2y^{\prime }/x-y=g(x)$. Setting $u=xy$ gives $u^{\prime \prime }-u=xg(x)$, reducing it to a constant-coefficient linear ordinary differential equation. If $g=e^{-3x}/x^2$, a decaying particular solution obtained by variation of parameters isThe small-argument expansion of the exponential integral givesThe undetermined decaying homogeneous coefficient $\alpha$ is fixed by a boundary condition or by a matched asymptotic expansion.