\[ \boxed { x{\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) + \left ( 2\,a{x}^{3}-1 \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) + \left ( {a}^{2}{x}^{3}+a \right ) {x}^{2}y \left ( x \right ) =0} \]
Mathematica: cpu = 1.204153 (sec), leaf count = 58 \[ \left \{\left \{y(x)\to \text {DifferentialRoot}\left (\{\unicode {f818},\unicode {f817}\}\unicode {f4a1}\left \{\left (a^2 \unicode {f817}^3+a\right ) \unicode {f818}(\unicode {f817})+\left (2 \unicode {f817}^3 a-1\right ) \unicode {f818}'(\unicode {f817})+\unicode {f817} \unicode {f818}''(\unicode {f817})=0,\unicode {f818}(1)=c_1,\unicode {f818}'(1)=c_2\right \}\right )(x)\right \}\right \} \]
Maple: cpu = 0.047 (sec), leaf count = 26 \[ \left \{ y \left ( x \right ) ={\it \_C1}\,{{\rm e}^{-{\frac {a{x}^{3}}{ 3}}}}+{\it \_C2}\,{{\rm e}^{-{\frac {a{x}^{3}}{3}}}}{x}^{2} \right \} \]