\[ \left (x^2+1\right ) y'(x)-x \left (x^2+1\right )+x y(x)=0 \] ✓ Mathematica : cpu = 0.0121073 (sec), leaf count = 27
\[\left \{\left \{y(x)\to \frac {1}{3} \left (x^2+1\right )+\frac {c_1}{\sqrt {x^2+1}}\right \}\right \}\] ✓ Maple : cpu = 0.009 (sec), leaf count = 20
\[ \left \{ y \left ( x \right ) ={\frac {{x}^{2}}{3}}+{\frac {1}{3}}+{{\it \_C1}{\frac {1}{\sqrt {{x}^{2}+1}}}} \right \} \]