\[ \boxed { {x}^{2}{\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) -x \left ( x+4 \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) +4\,y \left ( x \right ) =0} \]
Mathematica: cpu = 0.012502 (sec), leaf count = 41 \[ \left \{\left \{y(x)\to c_2 e^x x^4-\frac {1}{6} c_1 x \left (e^x x^3 \text {Ei}(-x)+x^2-x+2\right )\right \}\right \} \]
Maple: cpu = 0.016 (sec), leaf count = 32 \[ \left \{ y \left ( x \right ) ={\it \_C1}\,{x}^{4}{{\rm e}^{x}}+{\it \_C2}\,x \left ( {x}^{3}{{\rm e}^{x}}{\it Ei} \left ( 1,x \right ) -{x}^{ 2}+x-2 \right ) \right \} \]