\[ \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.0561198 (sec), leaf count = 75
\[\left \{\left \{y(x)\to \frac {1}{2} \left (e^{c_1}-\sqrt {4 e^{c_1} x+e^{2 c_1}-4 x^2}\right )\right \},\left \{y(x)\to \frac {1}{2} \left (\sqrt {4 e^{c_1} x+e^{2 c_1}-4 x^2}+e^{c_1}\right )\right \}\right \}\] ✓ Maple : cpu = 0.076 (sec), leaf count = 55
\[ \left \{ y \left ( x \right ) ={\frac {1}{2\,{\it \_C1}} \left ( 1+\sqrt {-4\,{{\it \_C1}}^{2}{x}^{2}+4\,{\it \_C1}\,x+1} \right ) },y \left ( x \right ) ={\frac {1}{2\,{\it \_C1}} \left ( 1-\sqrt {-4\,{{\it \_C1}}^{2}{x}^{2}+4\,{\it \_C1}\,x+1} \right ) } \right \} \]