\[ x^2 y''(x)-a \left (y(x)^n-y(x)\right )=0 \] ✗ Mathematica : cpu = 11.9325 (sec), leaf count = 0 , could not solve
DSolve[-(a*(-y[x] + y[x]^n)) + x^2*Derivative[2][y][x] == 0, y[x], x]
✓ Maple : cpu = 0.8 (sec), leaf count = 65
\[ \left \{ y \left ( x \right ) ={\it ODESolStruc} \left ( {\it \_a},[ \left \{ {\frac {\rm d}{{\rm d}{\it \_a}}}{\it \_b} \left ( {\it \_a} \right ) = \left ( {\it \_b} \left ( {\it \_a} \right ) \right ) ^{2} \left ( {\it \_a}\,{\it \_b} \left ( {\it \_a} \right ) a-{\it \_b} \left ( {\it \_a} \right ) {{\it \_a}}^{n}a-1 \right ) \right \} , \left \{ {\it \_a}=y \left ( x \right ) ,{\it \_b} \left ( {\it \_a} \right ) ={\frac {1}{x{\frac {\rm d}{{\rm d}x}}y \left ( x \right ) }} \right \} , \left \{ x={{\rm e}^{\int \!{\it \_b} \left ( {\it \_a} \right ) \,{\rm d}{\it \_a}+{\it \_C1}}},y \left ( x \right ) ={\it \_a} \right \} ] \right ) \right \} \]