\[ \left (x^2 y(x)-1\right ) y'(x)+x y(x)^2-1=0 \] ✓ Mathematica : cpu = 0.138568 (sec), leaf count = 60
DSolve[-1 + x*y[x]^2 + (-1 + x^2*y[x])*Derivative[1][y][x] == 0,y[x],x]
\[\left \{\left \{y(x)\to \frac {1}{x^2}-\frac {\sqrt {2 x^3+c_1 x^2+1}}{x^2}\right \},\left \{y(x)\to \frac {1}{x^2}+\frac {\sqrt {2 x^3+c_1 x^2+1}}{x^2}\right \}\right \}\] ✓ Maple : cpu = 0.017 (sec), leaf count = 51
dsolve((x^2*y(x)-1)*diff(y(x),x)+x*y(x)^2-1 = 0,y(x))
\[y \left (x \right ) = \frac {1+\sqrt {2 x^{3}-2 x^{2} c_{1}+1}}{x^{2}}\]