Lommel function

The Lommel differential equation, named after Eugen von Lommel, is an inhomogeneous form of the Bessel differential equation:

${\displaystyle z^{2}{\frac {d^{2}y}{dz^{2}}}+z{\frac {dy}{dz}}+(z^{2}-\nu ^{2})y=z^{\mu +1}.}$

Solutions are given by the Lommel functions s_μ,ν(z) and S_μ,ν(z), introduced by Eugen von Lommel (1880),

${\displaystyle s_{\mu ,\nu }(z)={\frac {\pi }{2}}\left[Y_{\nu }(z)\!\int _{0}^{z}\!\!x^{\mu }J_{\nu }(x)\,dx-J_{\nu }(z)\!\int _{0}^{z}\!\!x^{\mu }Y_{\nu }(x)\,dx\right],}$

${\displaystyle S_{\mu ,\nu }(z)=s_{\mu ,\nu }(z)+2^{\mu -1}\Gamma \left({\frac {\mu +\nu +1}{2}}\right)\Gamma \left({\frac {\mu -\nu +1}{2}}\right)\left(\sin \left[(\mu -\nu ){\frac {\pi }{2}}\right]J_{\nu }(z)-\cos \left[(\mu -\nu ){\frac {\pi }{2}}\right]Y_{\nu }(z)\right),}$

where J_ν(z) is a Bessel function of the first kind and Y_ν(z) a Bessel function of the second kind.