\[ x \left (-a+x^2+y(x)^2\right ) y'(x)-y(x) \left (a+x^2+y(x)^2\right )=0 \] ✓ Mathematica : cpu = 0.0417931 (sec), leaf count = 65
\[\left \{\left \{y(x)\to \frac {1}{2} \left (c_1 x-\sqrt {\left (c_1^2+4\right ) x^2-4 a}\right )\right \},\left \{y(x)\to \frac {1}{2} \left (\sqrt {\left (c_1^2+4\right ) x^2-4 a}+c_1 x\right )\right \}\right \}\]
✓ Maple : cpu = 0.106 (sec), leaf count = 112
\[ \left \{ \left ( \left ( y \left ( x \right ) \right ) ^{-2}- \left ( -{x}^{2}+a \right ) ^{-1} \right ) ^{-1}=-{x\sqrt {{x}^{2}-a}{\frac {1}{\sqrt {{\it \_C1}+4\,{\frac {a}{{x}^{2}-a}}}}}}+{\frac {{x}^{2}}{2}}-{\frac {a}{2}}, \left ( \left ( y \left ( x \right ) \right ) ^{-2}- \left ( -{x}^{2}+a \right ) ^{-1} \right ) ^{-1}={x\sqrt {{x}^{2}-a}{\frac {1}{\sqrt {{\it \_C1}+4\,{\frac {a}{{x}^{2}-a}}}}}}+{\frac {{x}^{2}}{2}}-{\frac {a}{2}} \right \} \]