[[_homogeneous, `class G`], _dAlembert]
Book solution method
Homogeneous equation
Mathematica ✓
cpu = 1.15496 (sec), leaf count = 2949
Maple ✓
cpu = 0.132 (sec), leaf count = 133
DSolve[y[x]*Sqrt[x^2 + y[x]^2] + x*(x + Sqrt[x^2 + y[x]^2])*y'[x] == 0,y[x],x]
Mathematica raw output
{{y[x] -> -Sqrt[(8*E^(6*C[1])*x^10 + x^16 - x^8*(8*E^(12*C[1])*x^12 + 20*E^(6*C[
1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) - x^6)^3])^(1/3) + (8*E^(12
*C[1])*x^12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) - x
^6)^3])^(2/3))/(x^6*(8*E^(12*C[1])*x^12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(
6*C[1])*x^24*(E^(6*C[1]) - x^6)^3])^(1/3))]/2}, {y[x] -> Sqrt[(8*E^(6*C[1])*x^10
+ x^16 - x^8*(8*E^(12*C[1])*x^12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1]
)*x^24*(E^(6*C[1]) - x^6)^3])^(1/3) + (8*E^(12*C[1])*x^12 + 20*E^(6*C[1])*x^18 -
x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) - x^6)^3])^(2/3))/(x^6*(8*E^(12*C[1])
*x^12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) - x^6)^3]
)^(1/3))]/2}, {y[x] -> -Sqrt[((-I)*(-8*(I + Sqrt[3])*E^(6*C[1])*x^10 + (x^8 + (8
*E^(12*C[1])*x^12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1
]) - x^6)^3])^(1/3))*(-((I + Sqrt[3])*x^8) + (-I + Sqrt[3])*(8*E^(12*C[1])*x^12
+ 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) - x^6)^3])^(1/3
))))/(x^6*(8*E^(12*C[1])*x^12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^
24*(E^(6*C[1]) - x^6)^3])^(1/3))]/(2*Sqrt[2])}, {y[x] -> Sqrt[((-I)*(-8*(I + Sqr
t[3])*E^(6*C[1])*x^10 + (x^8 + (8*E^(12*C[1])*x^12 + 20*E^(6*C[1])*x^18 - x^24 +
8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) - x^6)^3])^(1/3))*(-((I + Sqrt[3])*x^8) + (-
I + Sqrt[3])*(8*E^(12*C[1])*x^12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])
*x^24*(E^(6*C[1]) - x^6)^3])^(1/3))))/(x^6*(8*E^(12*C[1])*x^12 + 20*E^(6*C[1])*x
^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) - x^6)^3])^(1/3))]/(2*Sqrt[2])},
{y[x] -> -Sqrt[(I*(-8*(-I + Sqrt[3])*E^(6*C[1])*x^10 + (x^8 + (8*E^(12*C[1])*x^
12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) - x^6)^3])^(
1/3))*(-((-I + Sqrt[3])*x^8) + (I + Sqrt[3])*(8*E^(12*C[1])*x^12 + 20*E^(6*C[1])
*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) - x^6)^3])^(1/3))))/(x^6*(8*E^
(12*C[1])*x^12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1])
- x^6)^3])^(1/3))]/(2*Sqrt[2])}, {y[x] -> Sqrt[(I*(-8*(-I + Sqrt[3])*E^(6*C[1])*
x^10 + (x^8 + (8*E^(12*C[1])*x^12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1]
)*x^24*(E^(6*C[1]) - x^6)^3])^(1/3))*(-((-I + Sqrt[3])*x^8) + (I + Sqrt[3])*(8*E
^(12*C[1])*x^12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1])
- x^6)^3])^(1/3))))/(x^6*(8*E^(12*C[1])*x^12 + 20*E^(6*C[1])*x^18 - x^24 + 8*Sq
rt[E^(6*C[1])*x^24*(E^(6*C[1]) - x^6)^3])^(1/3))]/(2*Sqrt[2])}, {y[x] -> -Sqrt[(
-8*E^(6*C[1])*x^10 + x^16 - x^8*(8*E^(12*C[1])*x^12 - 20*E^(6*C[1])*x^18 - x^24
+ 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) + x^6)^3])^(1/3) + (8*E^(12*C[1])*x^12 - 20
*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) + x^6)^3])^(2/3))/(
x^6*(8*E^(12*C[1])*x^12 - 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^
(6*C[1]) + x^6)^3])^(1/3))]/2}, {y[x] -> Sqrt[(-8*E^(6*C[1])*x^10 + x^16 - x^8*(
8*E^(12*C[1])*x^12 - 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[
1]) + x^6)^3])^(1/3) + (8*E^(12*C[1])*x^12 - 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[
E^(6*C[1])*x^24*(E^(6*C[1]) + x^6)^3])^(2/3))/(x^6*(8*E^(12*C[1])*x^12 - 20*E^(6
*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) + x^6)^3])^(1/3))]/2}, {
y[x] -> -Sqrt[(8*(1 - I*Sqrt[3])*E^(6*C[1])*x^10 - I*(x^8 + (8*E^(12*C[1])*x^12
- 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) + x^6)^3])^(1/3
))*(-((I + Sqrt[3])*x^8) + (-I + Sqrt[3])*(8*E^(12*C[1])*x^12 - 20*E^(6*C[1])*x^
18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) + x^6)^3])^(1/3)))/(x^6*(8*E^(12*
C[1])*x^12 - 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) + x^
6)^3])^(1/3))]/(2*Sqrt[2])}, {y[x] -> Sqrt[(8*(1 - I*Sqrt[3])*E^(6*C[1])*x^10 -
I*(x^8 + (8*E^(12*C[1])*x^12 - 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^2
4*(E^(6*C[1]) + x^6)^3])^(1/3))*(-((I + Sqrt[3])*x^8) + (-I + Sqrt[3])*(8*E^(12*
C[1])*x^12 - 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) + x^
6)^3])^(1/3)))/(x^6*(8*E^(12*C[1])*x^12 - 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(
6*C[1])*x^24*(E^(6*C[1]) + x^6)^3])^(1/3))]/(2*Sqrt[2])}, {y[x] -> -Sqrt[(8*(1 +
I*Sqrt[3])*E^(6*C[1])*x^10 + I*(x^8 + (8*E^(12*C[1])*x^12 - 20*E^(6*C[1])*x^18
- x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) + x^6)^3])^(1/3))*(-((-I + Sqrt[3])*
x^8) + (I + Sqrt[3])*(8*E^(12*C[1])*x^12 - 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^
(6*C[1])*x^24*(E^(6*C[1]) + x^6)^3])^(1/3)))/(x^6*(8*E^(12*C[1])*x^12 - 20*E^(6*
C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) + x^6)^3])^(1/3))]/(2*Sqr
t[2])}, {y[x] -> Sqrt[(8*(1 + I*Sqrt[3])*E^(6*C[1])*x^10 + I*(x^8 + (8*E^(12*C[1
])*x^12 - 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) + x^6)^
3])^(1/3))*(-((-I + Sqrt[3])*x^8) + (I + Sqrt[3])*(8*E^(12*C[1])*x^12 - 20*E^(6*
C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[1]) + x^6)^3])^(1/3)))/(x^6*(
8*E^(12*C[1])*x^12 - 20*E^(6*C[1])*x^18 - x^24 + 8*Sqrt[E^(6*C[1])*x^24*(E^(6*C[
1]) + x^6)^3])^(1/3))]/(2*Sqrt[2])}}
Maple raw input
dsolve(x*(x+(x^2+y(x)^2)^(1/2))*diff(y(x),x)+y(x)*(x^2+y(x)^2)^(1/2) = 0, y(x),'implicit')
Maple raw output
Int(-(_a^2+y(x)^2)^(1/2)/_a/(2*(_a^2+y(x)^2)^(1/2)+_a),_a = _b .. x)+Intat(-(x+(
_f^2+x^2)^(1/2))/_f/(2*(_f^2+x^2)^(1/2)+x)-Int(-(-2+1/(_a^2+_f^2)^(1/2)*(2*(_a^2
+_f^2)^(1/2)+_a))*_f/_a/(2*(_a^2+_f^2)^(1/2)+_a)^2,_a = _b .. x),_f = y(x))+_C1
= 0