\[ \boxed { {x}^{2} \left ( 1+x \right ) {\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) -x \left ( 2\,x+1 \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) + \left ( 2\,x+1 \right ) y \left ( x \right ) =0} \]
Mathematica: cpu = 0.026003 (sec), leaf count = 18 \[ \left \{\left \{y(x)\to c_1 x+c_2 x (x+\log (x))\right \}\right \} \]
Maple: cpu = 0.015 (sec), leaf count = 15 \[ \left \{ y \left ( x \right ) =x{\it \_C1}+{\it \_C2}\,x \left ( x+\ln \left ( x \right ) \right ) \right \} \]