\[ \left (x^2+1\right ) y(x)+y'(x)-y(x)^2-2 x=0 \] ✓ Mathematica : cpu = 0.773316 (sec), leaf count = 48
\[\left \{\left \{y(x)\to \frac {e^{\frac {x^3}{3}+x}}{c_1-\int _1^x e^{\frac {K[1]^3}{3}+K[1]} \, dK[1]}+x^2+1\right \}\right \}\]
✓ Maple : cpu = 0.081 (sec), leaf count = 34
\[ \left \{ y \left ( x \right ) ={x}^{2}+1+{1{{\rm e}^{{\frac {{x}^{3}}{3}}+x}} \left ( {\it \_C1}-\int \!{{\rm e}^{{\frac {{x}^{3}}{3}}+x}}\,{\rm d}x \right ) ^{-1}} \right \} \]