\[ \boxed { {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) ={\frac {y \left ( x \right ) \ln \left ( x \right ) x-y \left ( x \right ) +2\,{x}^{5}b+2\,{x}^{3}a \left ( y \left ( x \right ) \right ) ^{2}}{ \left ( x\ln \left ( x \right ) -1 \right ) x}}=0} \]
Mathematica: cpu = 240.336019 (sec), leaf count = 65 \[ \left \{\left \{y(x)\to \frac {\sqrt {b} x \tan \left (\sqrt {a} \sqrt {b} \int _1^x \frac {2 K[1]^3}{K[1] \log (K[1])-1} \, dK[1]+\sqrt {a} \sqrt {b} c_1\right )}{\sqrt {a}}\right \}\right \} \]
Maple: cpu = 0.047 (sec), leaf count = 45 \[ \left \{ y \left ( x \right ) ={\frac {x}{a}\tan \left ( 2\,\int \!{ \frac {{x}^{3}}{x\ln \left ( x \right ) -1}}\,{\rm d}x\sqrt {ab}+2\,{ \it \_C1}\,\sqrt {ab} \right ) \sqrt {ab}} \right \} \]