\[ \boxed { {x}^{2}{\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) + \left ( 3\,x-1 \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) +y \left ( x \right ) =0} \]
Mathematica: cpu = 0.028004 (sec), leaf count = 37 \[ \left \{\left \{y(x)\to \frac {c_1 e^{-1/x}}{x}-\frac {c_2 e^{-1/x} \text {Ei}\left (\frac {1}{x}\right )}{x}\right \}\right \} \]
Maple: cpu = 0.016 (sec), leaf count = 25 \[ \left \{ y \left ( x \right ) ={\frac {1}{x{{\rm e}^{{x}^{-1}}}} \left ( {\it \_C1}\,{\it Ei} \left ( 1,-{x}^{-1} \right ) +{\it \_C2} \right ) } \right \} \]