\[ \left (x^2+1\right ) y'(x)-x \left (x^2+1\right )+x y(x)=0 \] ✓ Mathematica : cpu = 0.0122314 (sec), leaf count = 27
\[\left \{\left \{y(x)\to \frac {c_1}{\sqrt {x^2+1}}+\frac {1}{3} \left (x^2+1\right )\right \}\right \}\]
✓ Maple : cpu = 0.012 (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 \} \]