4.15.35 x(x2+y(x)2+x)y(x)+x2+y(x)2y(x)=0

ODE
x(x2+y(x)2+x)y(x)+x2+y(x)2y(x)=0 ODE Classification

[[_homogeneous, `class G`], _dAlembert]

Book solution method
Homogeneous equation

Mathematica
cpu = 1.15496 (sec), leaf count = 2949

{{y(x)12x16+8e6c1x10x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33x8+(x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)3)2/3x6x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33},{y(x)12x16+8e6c1x10x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33x8+(x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)3)2/3x6x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33},{y(x)i((x8+x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33)((i+3)x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33(i+3)x8)8(i+3)e6c1x10)x6x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)3322},{y(x)i((x8+x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33)((i+3)x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33(i+3)x8)8(i+3)e6c1x10)x6x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)3322},{y(x)i((x8+x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33)((i+3)x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33(i+3)x8)8(i+3)e6c1x10)x6x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)3322},{y(x)i((x8+x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33)((i+3)x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)33(i+3)x8)8(i+3)e6c1x10)x6x24+20e6c1x18+8e12c1x12+8e6c1x24(e6c1x6)3322},{y(x)12x168e6c1x10x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33x8+(x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)3)2/3x6x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33},{y(x)12x168e6c1x10x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33x8+(x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)3)2/3x6x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33},{y(x)8(1i3)e6c1x10i(x8+x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33)((i+3)x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33(i+3)x8)x6x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)3322},{y(x)8(1i3)e6c1x10i(x8+x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33)((i+3)x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33(i+3)x8)x6x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)3322},{y(x)8(1+i3)e6c1x10+i(x8+x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33)((i+3)x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33(i+3)x8)x6x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)3322},{y(x)8(1+i3)e6c1x10+i(x8+x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33)((i+3)x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)33(i+3)x8)x6x2420e6c1x18+8e12c1x12+8e6c1x24(x6+e6c1)3322}}

Maple
cpu = 0.132 (sec), leaf count = 133

{_bx1_a_a2+(y(x))2(2_a2+(y(x))2+_a)1d_a+y(x)1_f(x+_f2+x2)(2_f2+x2+x)1_bx_f_a(2+1(2_a2+_f2+_a)1_a2+_f2)(2_a2+_f2+_a)2d_ad_f+_C1=0} Mathematica raw input

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