\[ \boxed { \left ( 2\,x-1 \right ) {\frac {{\rm d}^{3}}{{\rm d}{x}^{3}}}y \left ( x \right ) + \left ( x+4 \right ) {\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) +2\,{\frac {\rm d}{{\rm d}x}}y \left ( x \right ) =0} \]
Mathematica: cpu = 114.051983 (sec), leaf count = 86 \[ \left \{\left \{y(x)\to \int _1^x \left (c_2 e^{-\frac {K[1]}{2}} L_{-\frac {1}{4}}^{\frac {5}{4}}\left (\frac {K[1]}{2}-\frac {1}{4}\right )+\frac {c_1 e^{-\frac {K[1]}{2}} \left (\frac {1}{4 \left (\frac {K[1]}{2}-\frac {1}{4}\right )}+1\right )}{\sqrt [4]{\frac {K[1]}{2}-\frac {1}{4}}}\right ) \, dK[1]+c_3\right \}\right \} \]
Maple: cpu = 0.078 (sec), leaf count = 38 \[ \left \{ y \left ( x \right ) ={1 \left ( {\it \_C3}+\int \!{(2\,x{\it \_C1}+{\it \_C2}){{\rm e}^{{\frac {x}{2}}}} \left ( 2\,x-1 \right ) ^{-{ \frac {3}{4}}}}\,{\rm d}x \right ) {{\rm e}^{-{\frac {x}{2}}}}{\frac {1 }{\sqrt [4]{2\,x-1}}}} \right \} \]