\[ \boxed { {\frac {{\rm d}^{2}}{{\rm d}{x}^{2}}}y \left ( x \right ) - \left ( {x}^{2}+1 \right ) y \left ( x \right ) =0} \]
Mathematica: cpu = 0.006001 (sec), leaf count = 33 \[ \left \{\left \{y(x)\to c_1 D_{-1}\left (\sqrt {2} x\right )+c_2 D_0\left (i \sqrt {2} x\right )\right \}\right \} \]
Maple: cpu = 0.031 (sec), leaf count = 23 \[ \left \{ y \left ( x \right ) ={\it \_C1}\,{{\rm e}^{{\frac {{x}^{2}}{2} }}}+{\it \_C2}\,{{\rm e}^{{\frac {{x}^{2}}{2}}}}{\it Erf} \left ( x \right ) \right \} \]