\[ a y(x)^2+x^2 y'(x)+x y(x)^3=0 \] ✓ Mathematica : cpu = 0.645525 (sec), leaf count = 78
\[\text {Solve}\left [-\frac {i a}{x}=\frac {2 e^{\frac {1}{2} \left (-\frac {i a}{x}-\frac {i}{y(x)}\right )^2}}{\sqrt {2 \pi } \text {erfi}\left (\frac {-\frac {i a}{x}-\frac {i}{y(x)}}{\sqrt {2}}\right )+2 c_1},y(x)\right ]\] ✓ Maple : cpu = 0.134 (sec), leaf count = 84
\[ \left \{ {\frac {1}{2} \left ( a\sqrt {\pi }\sqrt {2}{\it Erf} \left ( {\frac {\sqrt {2} \left ( ay \left ( x \right ) +x \right ) }{2\,xy \left ( x \right ) }} \right ) {{\rm e}^{{\frac { \left ( ay \left ( x \right ) +x \right ) ^{2}}{2\,{x}^{2} \left ( y \left ( x \right ) \right ) ^{2}}}}}+2\,x \right ) {{\rm e}^{-{\frac { \left ( \left ( a-x \right ) y \left ( x \right ) +x \right ) \left ( \left ( x+a \right ) y \left ( x \right ) +x \right ) }{2\,{x}^{2} \left ( y \left ( x \right ) \right ) ^{2}}}}}}+{\it \_C1}=0 \right \} \]