4.14.40 2x3+3x2y(x)+(x3+x2+3xy(x)2+2y(x)3)y(x)y(x)3+y(x)2=0

ODE
2x3+3x2y(x)+(x3+x2+3xy(x)2+2y(x)3)y(x)y(x)3+y(x)2=0 ODE Classification

[_rational]

Book solution method
Exact equation, integrating factor

Mathematica
cpu = 0.0442185 (sec), leaf count = 364

{{y(x)623(c1+x)22/3(729(c1x+x3)2+108(c1+x)3+27c1x+27x3)2/36729(c1x+x3)2+108(c1+x)3+27c1x+27x33},{y(x)22/3(1i3)(729(c1x+x3)2+108(c1+x)3+27c1x+27x3)2/36i23(3i)(c1+x)12729(c1x+x3)2+108(c1+x)3+27c1x+27x33},{y(x)22/3(1+i3)(729(c1x+x3)2+108(c1+x)3+27c1x+27x3)2/3+6i23(3+i)(c1+x)12729(c1x+x3)2+108(c1+x)3+27c1x+27x33}}

Maple
cpu = 0.162 (sec), leaf count = 28

{(y(x))3+(x+_C1)y(x)+x(x2+_C1)x+y(x)=0} Mathematica raw input

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

Mathematica raw output

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

Maple raw input

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

Maple raw output

(y(x)^3+(x+_C1)*y(x)+x*(x^2+_C1))/(x+y(x)) = 0