\[ \boxed { y \left ( x \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) + \left ( y \left ( x \right ) \right ) ^{2}+4\,x \left ( 1+x \right ) =0} \]
Mathematica: cpu = 0.010501 (sec), leaf count = 47 \[ \left \{\left \{y(x)\to -\sqrt {c_1 e^{-2 x}-4 x^2}\right \},\left \{y(x)\to \sqrt {c_1 e^{-2 x}-4 x^2}\right \}\right \} \]
Maple: cpu = 0.016 (sec), leaf count = 37 \[ \left \{ y \left ( x \right ) =\sqrt {{{\rm e}^{-2\,x}}{\it \_C1}-4\,{x} ^{2}},y \left ( x \right ) =-\sqrt {{{\rm e}^{-2\,x}}{\it \_C1}-4\,{x}^{ 2}} \right \} \]