ODE
\[ y'(x)=y(x)^2 \left (a e^x+y(x)\right ) \] ODE Classification
[_Abel]
Book solution method
Abel ODE, First kind
Mathematica ✓
cpu = 1.23794 (sec), leaf count = 73
\[\text {Solve}\left [-i a e^x=\frac {2 e^{-\frac {\left (a e^x y(x)+1\right )^2}{2 y(x)^2}}}{2 c_1-i \sqrt {2 \pi } \text {erf}\left (\frac {a e^x y(x)+1}{\sqrt {2} y(x)}\right )},y(x)\right ]\]
Maple ✓
cpu = 0.098 (sec), leaf count = 50
\[\left [\textit {\_C1} +\frac {{\mathrm e}^{-\frac {\left (a \,{\mathrm e}^{x}+\frac {1}{y \left (x \right )}\right )^{2}}{2}} {\mathrm e}^{-x}}{a}+\frac {\erf \left (\frac {\left (a \,{\mathrm e}^{x}+\frac {1}{y \left (x \right )}\right ) \sqrt {2}}{2}\right ) \sqrt {2}\, \sqrt {\pi }}{2} = 0\right ]\] Mathematica raw input
DSolve[y'[x] == y[x]^2*(a*E^x + y[x]),y[x],x]
Mathematica raw output
Solve[(-I)*a*E^x == 2/(E^((1 + a*E^x*y[x])^2/(2*y[x]^2))*(2*C[1] - I*Sqrt[2*Pi]*Erf[(1 + a*E^x*y[x])/(Sqrt[2]*y[x])])), y[x]]
Maple raw input
dsolve(diff(y(x),x) = (a*exp(x)+y(x))*y(x)^2, y(x))
Maple raw output
[_C1+exp(-1/2*(a*exp(x)+1/y(x))^2)/a/exp(x)+1/2*erf(1/2*(a*exp(x)+1/y(x))*2^(1/2))*2^(1/2)*Pi^(1/2) = 0]