4.14.30 \(x^2 (a x+3 y(x))+\left (x^3+y(x)^3\right ) y'(x)=0\)

ODE
\[ x^2 (a x+3 y(x))+\left (x^3+y(x)^3\right ) y'(x)=0 \] ODE Classification

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

Book solution method
Exact equation

Mathematica
cpu = 0.159115 (sec), leaf count = 1430

\[\left \{\left \{y(x)\to \frac {\sqrt {\frac {\sqrt [3]{3} a x^4-\sqrt [3]{3} e^{4 c_1}+\left (9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}\right ){}^{2/3}}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}-\sqrt {-\frac {6 \sqrt {2} x^3}{\sqrt {\frac {\sqrt [3]{3} a x^4-\sqrt [3]{3} e^{4 c_1}+\left (9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}\right ){}^{2/3}}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}}-\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}+\frac {\sqrt [3]{3} \left (e^{4 c_1}-a x^4\right )}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}}{\sqrt {2} \sqrt [3]{3}}\right \},\left \{y(x)\to \frac {\sqrt {\frac {\sqrt [3]{3} a x^4-\sqrt [3]{3} e^{4 c_1}+\left (9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}\right ){}^{2/3}}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}+\sqrt {-\frac {6 \sqrt {2} x^3}{\sqrt {\frac {\sqrt [3]{3} a x^4-\sqrt [3]{3} e^{4 c_1}+\left (9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}\right ){}^{2/3}}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}}-\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}+\frac {\sqrt [3]{3} \left (e^{4 c_1}-a x^4\right )}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}}{\sqrt {2} \sqrt [3]{3}}\right \},\left \{y(x)\to -\frac {\sqrt {\frac {\sqrt [3]{3} a x^4-\sqrt [3]{3} e^{4 c_1}+\left (9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}\right ){}^{2/3}}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}+\sqrt {\frac {6 \sqrt {2} x^3}{\sqrt {\frac {\sqrt [3]{3} a x^4-\sqrt [3]{3} e^{4 c_1}+\left (9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}\right ){}^{2/3}}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}}-\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}+\frac {\sqrt [3]{3} \left (e^{4 c_1}-a x^4\right )}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}}{\sqrt {2} \sqrt [3]{3}}\right \},\left \{y(x)\to \frac {\sqrt {\frac {6 \sqrt {2} x^3}{\sqrt {\frac {\sqrt [3]{3} a x^4-\sqrt [3]{3} e^{4 c_1}+\left (9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}\right ){}^{2/3}}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}}-\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}+\frac {\sqrt [3]{3} \left (e^{4 c_1}-a x^4\right )}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}-\sqrt {\frac {\sqrt [3]{3} a x^4-\sqrt [3]{3} e^{4 c_1}+\left (9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}\right ){}^{2/3}}{\sqrt [3]{9 x^6+\sqrt {3} \sqrt {27 x^{12}+\left (e^{4 c_1}-a x^4\right ){}^3}}}}}{\sqrt {2} \sqrt [3]{3}}\right \}\right \}\]

Maple
cpu = 0.014 (sec), leaf count = 35

\[ \left \{ -{\frac {1}{4}\ln \left ({\frac {a{x}^{4}+4\,{x}^{3}y \relax (x ) + \left (y \relax (x ) \right ) ^{4}}{{x}^{4}}} \right ) }-\ln \relax (x ) -{\it \_C1}=0 \right \} \] Mathematica raw input

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

Mathematica raw output

{{y[x] -> (Sqrt[(-(3^(1/3)*E^(4*C[1])) + 3^(1/3)*a*x^4 + (9*x^6 + Sqrt[3]*Sqrt[2
7*x^12 + (E^(4*C[1]) - a*x^4)^3])^(2/3))/(9*x^6 + Sqrt[3]*Sqrt[27*x^12 + (E^(4*C
[1]) - a*x^4)^3])^(1/3)] - Sqrt[(3^(1/3)*(E^(4*C[1]) - a*x^4))/(9*x^6 + Sqrt[3]*
Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(1/3) - (9*x^6 + Sqrt[3]*Sqrt[27*x^12 + 
(E^(4*C[1]) - a*x^4)^3])^(1/3) - (6*Sqrt[2]*x^3)/Sqrt[(-(3^(1/3)*E^(4*C[1])) + 3
^(1/3)*a*x^4 + (9*x^6 + Sqrt[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(2/3))/(
9*x^6 + Sqrt[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(1/3)]])/(Sqrt[2]*3^(1/3
))}, {y[x] -> (Sqrt[(-(3^(1/3)*E^(4*C[1])) + 3^(1/3)*a*x^4 + (9*x^6 + Sqrt[3]*Sq
rt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(2/3))/(9*x^6 + Sqrt[3]*Sqrt[27*x^12 + (E^
(4*C[1]) - a*x^4)^3])^(1/3)] + Sqrt[(3^(1/3)*(E^(4*C[1]) - a*x^4))/(9*x^6 + Sqrt
[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(1/3) - (9*x^6 + Sqrt[3]*Sqrt[27*x^1
2 + (E^(4*C[1]) - a*x^4)^3])^(1/3) - (6*Sqrt[2]*x^3)/Sqrt[(-(3^(1/3)*E^(4*C[1]))
 + 3^(1/3)*a*x^4 + (9*x^6 + Sqrt[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(2/3
))/(9*x^6 + Sqrt[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(1/3)]])/(Sqrt[2]*3^
(1/3))}, {y[x] -> -((Sqrt[(-(3^(1/3)*E^(4*C[1])) + 3^(1/3)*a*x^4 + (9*x^6 + Sqrt
[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(2/3))/(9*x^6 + Sqrt[3]*Sqrt[27*x^12
 + (E^(4*C[1]) - a*x^4)^3])^(1/3)] + Sqrt[(3^(1/3)*(E^(4*C[1]) - a*x^4))/(9*x^6 
+ Sqrt[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(1/3) - (9*x^6 + Sqrt[3]*Sqrt[
27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(1/3) + (6*Sqrt[2]*x^3)/Sqrt[(-(3^(1/3)*E^(4*
C[1])) + 3^(1/3)*a*x^4 + (9*x^6 + Sqrt[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3]
)^(2/3))/(9*x^6 + Sqrt[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(1/3)]])/(Sqrt
[2]*3^(1/3)))}, {y[x] -> (-Sqrt[(-(3^(1/3)*E^(4*C[1])) + 3^(1/3)*a*x^4 + (9*x^6 
+ Sqrt[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(2/3))/(9*x^6 + Sqrt[3]*Sqrt[2
7*x^12 + (E^(4*C[1]) - a*x^4)^3])^(1/3)] + Sqrt[(3^(1/3)*(E^(4*C[1]) - a*x^4))/(
9*x^6 + Sqrt[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(1/3) - (9*x^6 + Sqrt[3]
*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(1/3) + (6*Sqrt[2]*x^3)/Sqrt[(-(3^(1/3)
*E^(4*C[1])) + 3^(1/3)*a*x^4 + (9*x^6 + Sqrt[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x
^4)^3])^(2/3))/(9*x^6 + Sqrt[3]*Sqrt[27*x^12 + (E^(4*C[1]) - a*x^4)^3])^(1/3)]])
/(Sqrt[2]*3^(1/3))}}

Maple raw input

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

Maple raw output

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