\[ \boxed { 2\,xy \left ( x \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) - \left ( y \left ( x \right ) \right ) ^{2}+ax=0} \]
Mathematica: cpu = 0.008501 (sec), leaf count = 41 \[ \left \{\left \{y(x)\to -\sqrt {c_1 x-a x \log (x)}\right \},\left \{y(x)\to \sqrt {c_1 x-a x \log (x)}\right \}\right \} \]
Maple: cpu = 0.016 (sec), leaf count = 33 \[ \left \{ y \left ( x \right ) =\sqrt {-ax\ln \left ( x \right ) +{\it \_C1}\,x},y \left ( x \right ) =-\sqrt {-ax\ln \left ( x \right ) +{\it \_C1}\,x} \right \} \]