4.10.31 (x3y(x))y(x)y(x)+3x+4=0

ODE
(x3y(x))y(x)y(x)+3x+4=0 ODE Classification

[[_homogeneous, `class C`], _rational, [_Abel, `2nd type`, `class A`]]

Book solution method
Equation linear in the variables, y(x)=f(X1X2)

Mathematica
cpu = 0.0688683 (sec), leaf count = 781

{{y(x)13(x1Root[#16(16e12c1+1024x6+9216x5+34560x4+69120x3+77760x2+46656x+11664)+#14(384x42304x35184x25184x1944)+#13(64x3+288x2+432x+216)+#12(36x2+108x+81)+#1(12x18)+1&,1])},{y(x)13(x1Root[#16(16e12c1+1024x6+9216x5+34560x4+69120x3+77760x2+46656x+11664)+#14(384x42304x35184x25184x1944)+#13(64x3+288x2+432x+216)+#12(36x2+108x+81)+#1(12x18)+1&,2])},{y(x)13(x1Root[#16(16e12c1+1024x6+9216x5+34560x4+69120x3+77760x2+46656x+11664)+#14(384x42304x35184x25184x1944)+#13(64x3+288x2+432x+216)+#12(36x2+108x+81)+#1(12x18)+1&,3])},{y(x)13(x1Root[#16(16e12c1+1024x6+9216x5+34560x4+69120x3+77760x2+46656x+11664)+#14(384x42304x35184x25184x1944)+#13(64x3+288x2+432x+216)+#12(36x2+108x+81)+#1(12x18)+1&,4])},{y(x)13(x1Root[#16(16e12c1+1024x6+9216x5+34560x4+69120x3+77760x2+46656x+11664)+#14(384x42304x35184x25184x1944)+#13(64x3+288x2+432x+216)+#12(36x2+108x+81)+#1(12x18)+1&,5])},{y(x)13(x1Root[#16(16e12c1+1024x6+9216x5+34560x4+69120x3+77760x2+46656x+11664)+#14(384x42304x35184x25184x1944)+#13(64x3+288x2+432x+216)+#12(36x2+108x+81)+#1(12x18)+1&,6])}}

Maple
cpu = 0.024 (sec), leaf count = 55

{23ln(2y(x)42x3+2x)13ln(2y(x)+2+2x3+2x)ln(3+2x)_C1=0} Mathematica raw input

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

Mathematica raw output

{{y[x] -> (x - Root[1 + (-18 - 12*x)*#1 + (81 + 108*x + 36*x^2)*#1^2 + (216 + 43
2*x + 288*x^2 + 64*x^3)*#1^3 + (-1944 - 5184*x - 5184*x^2 - 2304*x^3 - 384*x^4)*
#1^4 + (11664 + 16*E^(12*C[1]) + 46656*x + 77760*x^2 + 69120*x^3 + 34560*x^4 + 9
216*x^5 + 1024*x^6)*#1^6 & , 1]^(-1))/3}, {y[x] -> (x - Root[1 + (-18 - 12*x)*#1
 + (81 + 108*x + 36*x^2)*#1^2 + (216 + 432*x + 288*x^2 + 64*x^3)*#1^3 + (-1944 -
 5184*x - 5184*x^2 - 2304*x^3 - 384*x^4)*#1^4 + (11664 + 16*E^(12*C[1]) + 46656*
x + 77760*x^2 + 69120*x^3 + 34560*x^4 + 9216*x^5 + 1024*x^6)*#1^6 & , 2]^(-1))/3
}, {y[x] -> (x - Root[1 + (-18 - 12*x)*#1 + (81 + 108*x + 36*x^2)*#1^2 + (216 + 
432*x + 288*x^2 + 64*x^3)*#1^3 + (-1944 - 5184*x - 5184*x^2 - 2304*x^3 - 384*x^4
)*#1^4 + (11664 + 16*E^(12*C[1]) + 46656*x + 77760*x^2 + 69120*x^3 + 34560*x^4 +
 9216*x^5 + 1024*x^6)*#1^6 & , 3]^(-1))/3}, {y[x] -> (x - Root[1 + (-18 - 12*x)*
#1 + (81 + 108*x + 36*x^2)*#1^2 + (216 + 432*x + 288*x^2 + 64*x^3)*#1^3 + (-1944
 - 5184*x - 5184*x^2 - 2304*x^3 - 384*x^4)*#1^4 + (11664 + 16*E^(12*C[1]) + 4665
6*x + 77760*x^2 + 69120*x^3 + 34560*x^4 + 9216*x^5 + 1024*x^6)*#1^6 & , 4]^(-1))
/3}, {y[x] -> (x - Root[1 + (-18 - 12*x)*#1 + (81 + 108*x + 36*x^2)*#1^2 + (216 
+ 432*x + 288*x^2 + 64*x^3)*#1^3 + (-1944 - 5184*x - 5184*x^2 - 2304*x^3 - 384*x
^4)*#1^4 + (11664 + 16*E^(12*C[1]) + 46656*x + 77760*x^2 + 69120*x^3 + 34560*x^4
 + 9216*x^5 + 1024*x^6)*#1^6 & , 5]^(-1))/3}, {y[x] -> (x - Root[1 + (-18 - 12*x
)*#1 + (81 + 108*x + 36*x^2)*#1^2 + (216 + 432*x + 288*x^2 + 64*x^3)*#1^3 + (-19
44 - 5184*x - 5184*x^2 - 2304*x^3 - 384*x^4)*#1^4 + (11664 + 16*E^(12*C[1]) + 46
656*x + 77760*x^2 + 69120*x^3 + 34560*x^4 + 9216*x^5 + 1024*x^6)*#1^6 & , 6]^(-1
))/3}}

Maple raw input

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

Maple raw output

-2/3*ln((-2*y(x)-4-2*x)/(3+2*x))-1/3*ln((-2*y(x)+2+2*x)/(3+2*x))-ln(3+2*x)-_C1 =
 0