\[ \boxed { \left ( {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) \right ) ^{2}-xy \left ( x \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) + \left ( y \left ( x \right ) \right ) ^{2}\ln \left ( ay \left ( x \right ) \right ) =0} \]
Mathematica: cpu = 0.255032 (sec), leaf count = 27 \[ \left \{\left \{y(x)\to \frac {e^{\frac {c_1 x}{2}-\frac {c_1^2}{4}}}{a}\right \}\right \} \]
Maple: cpu = 0.312 (sec), leaf count = 50 \[ \left \{ y \left ( x \right ) ={\frac {1}{a}{{\rm e}^{{\frac {{x}^{2}}{4 }}}}},y \left ( x \right ) ={\frac {1}{{{\rm e}^{{{\it \_C1}}^{2}}}{ {\rm e}^{{\it \_C1}\,x}}a}},y \left ( x \right ) ={\frac {{{\rm e}^{{ \it \_C1}\,x}}}{{{\rm e}^{{{\it \_C1}}^{2}}}a}} \right \} \]