\[ \boxed { x{\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) -{\frac {\rm d}{{\rm d}x}}y \left ( x \right ) +ay \left ( x \right ) =0} \]
Mathematica: cpu = 0.026003 (sec), leaf count = 46 \[ \left \{\left \{y(x)\to 2 a c_1 x J_2\left (2 \sqrt {a} \sqrt {x}\right )-2 a c_2 x Y_2\left (2 \sqrt {a} \sqrt {x}\right )\right \}\right \} \]
Maple: cpu = 0.016 (sec), leaf count = 31 \[ \left \{ y \left ( x \right ) ={\it \_C1}\,x{{\sl J}_{2}\left (2\,\sqrt { a}\sqrt {x}\right )}+{\it \_C2}\,x{{\sl Y}_{2}\left (2\,\sqrt {a}\sqrt { x}\right )} \right \} \]