4.14.26 \(y(x) \left (x^2 y(x)^2+1\right )+x (1-x y(x))^2 y'(x)=0\)

ODE
\[ y(x) \left (x^2 y(x)^2+1\right )+x (1-x y(x))^2 y'(x)=0 \] ODE Classification

[[_homogeneous, `class G`], _rational]

Book solution method
Homogeneous equation, isobaric equation

Mathematica
cpu = 0.41484 (sec), leaf count = 24

\[\text {Solve}\left [\frac {1}{x y(x)}+2 \log (y(x))+c_1=x y(x),y(x)\right ]\]

Maple
cpu = 0.344 (sec), leaf count = 34

\[\left [y \left (x \right ) = \frac {{\mathrm e}^{\RootOf \left (-2 \,{\mathrm e}^{\textit {\_Z}} \ln \left (x \right )-{\mathrm e}^{2 \textit {\_Z}}+2 \,{\mathrm e}^{\textit {\_Z}} \textit {\_C1} +2 \,{\mathrm e}^{\textit {\_Z}} \textit {\_Z} +1\right )}}{x}\right ]\] Mathematica raw input

DSolve[y[x]*(1 + x^2*y[x]^2) + x*(1 - x*y[x])^2*y'[x] == 0,y[x],x]

Mathematica raw output

Solve[C[1] + 2*Log[y[x]] + 1/(x*y[x]) == x*y[x], y[x]]

Maple raw input

dsolve(x*(1-x*y(x))^2*diff(y(x),x)+(1+x^2*y(x)^2)*y(x) = 0, y(x))

Maple raw output

[y(x) = exp(RootOf(-2*exp(_Z)*ln(x)-exp(2*_Z)+2*exp(_Z)*_C1+2*exp(_Z)*_Z+1))/x]