\[ x^2 \left (-y'(x)\right )+y''(x)-(x+1)^2 y(x)=0 \] ✓ Mathematica : cpu = 0.0484582 (sec), leaf count = 56
\[\left \{\left \{y(x)\to c_2 e^{\frac {x^3}{3}+x} \int _1^xe^{-\frac {1}{3} K[1]^3-2 K[1]}dK[1]+c_1 e^{\frac {x^3}{3}+x}\right \}\right \}\] ✓ Maple : cpu = 0.148 (sec), leaf count = 50
\[y \relax (x ) = c_{1} \mathit {HT}\left (0, -3, 2 \,3^{\frac {1}{3}}, \frac {3^{\frac {2}{3}} x}{3}\right ) {\mathrm e}^{-x}+c_{2} \mathit {HT}\left (0, 3, 2 \,3^{\frac {1}{3}}, -\frac {3^{\frac {2}{3}} x}{3}\right ) {\mathrm e}^{\frac {x \left (x^{2}+3\right )}{3}}\]