\[ y(x) \left (a x-b^2\right )+x^2 y''(x)+2 x y'(x)=0 \] ✓ Mathematica : cpu = 0.0757834 (sec), leaf count = 103
\[\left \{\left \{y(x)\to \frac {c_1 \Gamma \left (1-\sqrt {4 b^2+1}\right ) J_{-\sqrt {4 b^2+1}}\left (2 \sqrt {a} \sqrt {x}\right )+c_2 \Gamma \left (\sqrt {4 b^2+1}+1\right ) J_{\sqrt {4 b^2+1}}\left (2 \sqrt {a} \sqrt {x}\right )}{\sqrt {a} \sqrt {x}}\right \}\right \}\]
✓ Maple : cpu = 0.029 (sec), leaf count = 49
\[ \left \{ y \left ( x \right ) ={1 \left ( {\it \_C2}\,{{\sl Y}_{\sqrt {4\,{b}^{2}+1}}\left (2\,\sqrt {a}\sqrt {x}\right )}+{\it \_C1}\,{{\sl J}_{\sqrt {4\,{b}^{2}+1}}\left (2\,\sqrt {a}\sqrt {x}\right )} \right ) {\frac {1}{\sqrt {x}}}} \right \} \]