\[ \boxed { x{\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) - \left ( {x}^{2}-x \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) + \left ( x-1 \right ) y \left ( x \right ) =0} \]
Mathematica: cpu = 10.521836 (sec), leaf count = 40 \[ \left \{\left \{y(x)\to c_2 x \int _1^x \frac {e^{\frac {K[1]^2}{2}-K[1]}}{K[1]^2} \, dK[1]+c_1 x\right \}\right \} \]
Maple: cpu = 0.016 (sec), leaf count = 26 \[ \left \{ y \left ( x \right ) = \left ( \int \!{\frac {1}{{x}^{2}}{ {\rm e}^{{\frac {{x}^{2}}{2}}-x}}}\,{\rm d}x{\it \_C1}+{\it \_C2} \right ) x \right \} \]