\[ \left (-x^2+2 x y(x)+y(x)^2\right ) y'(x)+x^2+2 x y(x)-y(x)^2=0 \] ✓ Mathematica : cpu = 0.224416 (sec), leaf count = 75
\[\left \{\left \{y(x)\to \frac {1}{2} \left (e^{c_1}-\sqrt {-4 x^2+4 e^{c_1} x+e^{2 c_1}}\right )\right \},\left \{y(x)\to \frac {1}{2} \left (\sqrt {-4 x^2+4 e^{c_1} x+e^{2 c_1}}+e^{c_1}\right )\right \}\right \}\] ✓ Maple : cpu = 0.053 (sec), leaf count = 55
\[\left \{y \left (x \right ) = \frac {1-\sqrt {-4 c_{1}^{2} x^{2}+4 c_{1} x +1}}{2 c_{1}}, y \left (x \right ) = \frac {1+\sqrt {-4 c_{1}^{2} x^{2}+4 c_{1} x +1}}{2 c_{1}}\right \}\]