4.12.7 x(3y(x)+2x)y(x)=y(x)2

ODE
x(3y(x)+2x)y(x)=y(x)2 ODE Classification

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

Book solution method
Homogeneous equation

Mathematica
cpu = 0.0881313 (sec), leaf count = 413

{{y(x)13(x2323ec1x2(27ec14x2)+27ec1x2x33+323ec1x2(27ec14x2)+27ec1x2x33x)},{y(x)112(2(1+i3)x2323ec1x2(27ec14x2)+27ec1x2x33+i22/3(3+i)33ec1x2(27ec14x2)+27ec1x2x334x)},{y(x)112(2i(3+i)x2323ec1x2(27ec14x2)+27ec1x2x3322/3(1+i3)33ec1x2(27ec14x2)+27ec1x2x334x)}}

Maple
cpu = 0.019 (sec), leaf count = 31

{12ln(x+y(x)x)ln(y(x)x)ln(x)_C1=0} Mathematica raw input

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

Mathematica raw output

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

Maple raw input

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

Maple raw output

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