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.

The s function can also be written as^[1]

${\displaystyle s_{\mu ,\nu }(z)={\frac {z^{\mu +1}}{(\mu -\nu +1)(\mu +\nu +1)}}{}_{1}F_{2}(1;{\frac {\mu }{2}}-{\frac {\nu }{2}}+{\frac {3}{2}},{\frac {\mu }{2}}+{\frac {\nu }{2}}+{\frac {3}{2}};-{\frac {z^{2}}{4}}),}$

where _pF_q is a generalized hypergeometric function.