\[ \boxed { \left ( xy \left ( x \right ) -1 \right ) ^{2}x{\frac {\rm d}{{\rm d}x}}y \left ( x \right ) + \left ( {x}^{2} \left ( y \left ( x \right ) \right ) ^{2}+1 \right ) y \left ( x \right ) =0} \]
Mathematica: cpu = 0.044006 (sec), leaf count = 25 \[ \text {Solve}\left [x y(x)-\frac {1}{x y(x)}-2 \log (y(x))=c_1,y(x)\right ] \]
Maple: cpu = 0.125 (sec), leaf count = 34 \[ \left \{ y \left ( x \right ) ={\frac {{{\rm e}^{{\it RootOf} \left ( -2 \,{{\rm e}^{{\it \_Z}}}\ln \left ( x \right ) -{{\rm e}^{2\,{\it \_Z}}} +2\,{{\rm e}^{{\it \_Z}}}{\it \_C1}+2\,{\it \_Z}\,{{\rm e}^{{\it \_Z}} }+1 \right ) }}}{x}} \right \} \]