\[ y(x) y'(x)+y(x)^2+4 x (x+1)=0 \] ✓ Mathematica : cpu = 0.0121716 (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.023 (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 \} \]