\[ y'(x)=\frac {x (x-y(x))^2 (y(x)+x)^2}{y(x)} \] ✓ Mathematica : cpu = 0.148255 (sec), leaf count = 126
\[\left \{\left \{y(x)\to -\frac {\sqrt {x^2+x^2 e^{2 x^2+4 c_1}-e^{2 x^2+4 c_1}+1}}{\sqrt {1+e^{2 x^2+4 c_1}}}\right \},\left \{y(x)\to \frac {\sqrt {x^2+x^2 e^{2 x^2+4 c_1}-e^{2 x^2+4 c_1}+1}}{\sqrt {1+e^{2 x^2+4 c_1}}}\right \}\right \}\] ✓ Maple : cpu = 0.25 (sec), leaf count = 192
\[ \left \{ y \left ( x \right ) ={\sqrt { \left ( \left ( {x}^{2}+1 \right ) {{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}+2 \right ) }{2}}}}+{\it \_C1}\, \left ( {x}^{2}-1 \right ) {{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}-2 \right ) }{2}}}} \right ) \left ( {\it \_C1}\,{{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}-2 \right ) }{2}}}}+{{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}+2 \right ) }{2}}}} \right ) } \left ( {\it \_C1}\,{{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}-2 \right ) }{2}}}}+{{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}+2 \right ) }{2}}}} \right ) ^{-1}},y \left ( x \right ) =-{\sqrt { \left ( \left ( {x}^{2}+1 \right ) {{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}+2 \right ) }{2}}}}+{\it \_C1}\, \left ( {x}^{2}-1 \right ) {{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}-2 \right ) }{2}}}} \right ) \left ( {\it \_C1}\,{{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}-2 \right ) }{2}}}}+{{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}+2 \right ) }{2}}}} \right ) } \left ( {\it \_C1}\,{{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}-2 \right ) }{2}}}}+{{\rm e}^{-{\frac {{x}^{2} \left ( {x}^{2}+2 \right ) }{2}}}} \right ) ^{-1}} \right \} \]