4.12.4 x(2y(x)x+1)y(x)+y(x)(y(x)+2x+1)=0

ODE
x(2y(x)x+1)y(x)+y(x)(y(x)+2x+1)=0 ODE Classification

[_rational, [_Abel, `2nd type`, `class B`]]

Book solution method
Exact equation, integrating factor

Mathematica
cpu = 14.9599 (sec), leaf count = 451

{{y(x)23x27c12x2+108c13x3+(27c12x27c12x2)2+27c12x3+27c12x2+108c13x3+(27c12x27c12x2)2+27c12x3323c1+x1},{y(x)(1+i3)x22/327c12x2+108c13x3+(27c12x27c12x2)2+27c12x3(1i3)27c12x2+108c13x3+(27c12x27c12x2)2+27c12x3623c1+x1},{y(x)(1i3)x22/327c12x2+108c13x3+(27c12x27c12x2)2+27c12x3(1+i3)27c12x2+108c13x3+(27c12x27c12x2)2+27c12x3623c1+x1}}

Maple
cpu = 0.032 (sec), leaf count = 76

{79ln(189189x4+4x+8y(x))76ln(126xy(x)11+x+2y(x))+718ln(378y(x)5+5x+10y(x))7ln(1+x)9+7ln(x)18_C1=0} Mathematica raw input

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

Mathematica raw output

{{y[x] -> -1 + x - (2^(1/3)*x)/(27*x*C[1]^2 - 27*x^2*C[1]^2 + Sqrt[108*x^3*C[1]^
3 + (27*x*C[1]^2 - 27*x^2*C[1]^2)^2])^(1/3) + (27*x*C[1]^2 - 27*x^2*C[1]^2 + Sqr
t[108*x^3*C[1]^3 + (27*x*C[1]^2 - 27*x^2*C[1]^2)^2])^(1/3)/(3*2^(1/3)*C[1])}, {y
[x] -> -1 + x + ((1 + I*Sqrt[3])*x)/(2^(2/3)*(27*x*C[1]^2 - 27*x^2*C[1]^2 + Sqrt
[108*x^3*C[1]^3 + (27*x*C[1]^2 - 27*x^2*C[1]^2)^2])^(1/3)) - ((1 - I*Sqrt[3])*(2
7*x*C[1]^2 - 27*x^2*C[1]^2 + Sqrt[108*x^3*C[1]^3 + (27*x*C[1]^2 - 27*x^2*C[1]^2)
^2])^(1/3))/(6*2^(1/3)*C[1])}, {y[x] -> -1 + x + ((1 - I*Sqrt[3])*x)/(2^(2/3)*(2
7*x*C[1]^2 - 27*x^2*C[1]^2 + Sqrt[108*x^3*C[1]^3 + (27*x*C[1]^2 - 27*x^2*C[1]^2)
^2])^(1/3)) - ((1 + I*Sqrt[3])*(27*x*C[1]^2 - 27*x^2*C[1]^2 + Sqrt[108*x^3*C[1]^
3 + (27*x*C[1]^2 - 27*x^2*C[1]^2)^2])^(1/3))/(6*2^(1/3)*C[1])}}

Maple raw input

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

Maple raw output

7/9*ln((189-189*x)/(-4+4*x+8*y(x)))-7/6*ln(-126*(x-y(x)-1)/(-1+x+2*y(x)))+7/18*l
n(378*y(x)/(-5+5*x+10*y(x)))-7/9*ln(-1+x)+7/18*ln(x)-_C1 = 0