4.13.10 (3x2y(x)2)y(x)=2xy(x)

ODE
(3x2y(x)2)y(x)=2xy(x) ODE Classification

[[_homogeneous, `class A`], _rational, _dAlembert]

Book solution method
Exact equation, integrating factor

Mathematica
cpu = 0.178843 (sec), leaf count = 458

{{y(x)13(27ec1x2+381e2c1x412e4c1x22e3c1323+23e2c127ec1x2+381e2c1x412e4c1x22e3c13ec1)},{y(x)i(3+i)27ec1x2+381e2c1x412e4c1x22e3c13623i(3i)e2c13 22/327ec1x2+381e2c1x412e4c1x22e3c13ec13},{y(x)i(3i)27ec1x2+381e2c1x412e4c1x22e3c13623+i(3+i)e2c13 22/327ec1x2+381e2c1x412e4c1x22e3c13ec13}}

Maple
cpu = 0.025 (sec), leaf count = 40

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

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

Mathematica raw output

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

Maple raw input

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

Maple raw output

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