\[ \boxed { {\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) =-{\frac { \left ( 2\,{x}^{2}+1 \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) }{{x}^{3}}}+{\frac {y \left ( x \right ) }{{x}^{4}}}=0} \]
Mathematica: cpu = 0.024003 (sec), leaf count = 50 \[ \left \{\left \{y(x)\to c_1 e^{\frac {1}{2 x^2}}-\sqrt {\frac {\pi }{2}} c_2 e^{\frac {1}{2 x^2}} \text {erf}\left (\frac {1}{\sqrt {2} x}\right )\right \}\right \} \]
Maple: cpu = 0.031 (sec), leaf count = 30 \[ \left \{ y \left ( x \right ) ={\it \_C1}\,{{\rm e}^{{\frac {1}{2\,{x}^{ 2}}}}}+{\it \_C2}\,{{\rm e}^{{\frac {1}{2\,{x}^{2}}}}}{\it Erf} \left ( {\frac {\sqrt {2}}{2\,x}} \right ) \right \} \]