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