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