\[ \boxed { 2\, \left ( {x}^{2}+1 \right ) \left ( {\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) \right ) ^{2}-x \left ( {\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) \right ) \left ( x+4\,{\frac {\rm d}{{\rm d}x}}y \left ( x \right ) \right ) +2\, \left ( x+{\frac {\rm d}{{\rm d}x}}y \left ( x \right ) \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) -2\,y \left ( x \right ) =0} \]
Mathematica: cpu = 0.012502 (sec), leaf count = 32 \[ \left \{\left \{y(x)\to -\frac {1}{2} \sqrt {c_2-c_1^2} x^2+c_1 x+c_2\right \}\right \} \]
Maple: cpu = 1.045 (sec), leaf count = 20 \[ \left \{ y \left ( x \right ) ={\frac {{\it \_C1}\,{x}^{2}}{2}}+{\it \_C2}\,x+{{\it \_C1}}^{2}+{{\it \_C2}}^{2} \right \} \]