\[ \boxed { \left ( {x}^{2}+1 \right ) {\frac {\rm d}{{\rm d}x}}y \left ( x \right ) +xy \left ( x \right ) -1=0} \]
Mathematica: cpu = 0.012002 (sec), leaf count = 30 \[ \left \{\left \{y(x)\to \frac {c_1}{\sqrt {x^2+1}}+\frac {\sinh ^{-1}(x)}{\sqrt {x^2+1}}\right \}\right \} \]
Maple: cpu = 0.0 (sec), leaf count = 16 \[ \left \{ y \left ( x \right ) ={({\it Arcsinh} \left ( x \right ) +{\it \_C1}){\frac {1}{\sqrt {{x}^{2}+1}}}} \right \} \]