4.40.23 \(x y(x) y''(x)+x y'(x)^2+2 y(x) y'(x)=0\)

ODE
\[ x y(x) y''(x)+x y'(x)^2+2 y(x) y'(x)=0 \] ODE Classification

[[_2nd_order, _exact, _nonlinear], _Liouville, [_2nd_order, _with_linear_symmetries], [_2nd_order, _reducible, _mu_x_y1], [_2nd_order, _reducible, _mu_xy]]

Book solution method
TO DO

Mathematica
cpu = 0.319676 (sec), leaf count = 24

\[\left \{\left \{y(x)\to \frac {c_2 \sqrt {2-c_1 x}}{\sqrt {x}}\right \}\right \}\]

Maple
cpu = 0.12 (sec), leaf count = 38

\[\left [y \left (x \right ) = \frac {\sqrt {-2 x \left (-\textit {\_C2} x +\textit {\_C1} \right )}}{x}, y \left (x \right ) = -\frac {\sqrt {-2 x \left (-\textit {\_C2} x +\textit {\_C1} \right )}}{x}\right ]\] Mathematica raw input

DSolve[2*y[x]*y'[x] + x*y'[x]^2 + x*y[x]*y''[x] == 0,y[x],x]

Mathematica raw output

{{y[x] -> (Sqrt[2 - x*C[1]]*C[2])/Sqrt[x]}}

Maple raw input

dsolve(x*y(x)*diff(diff(y(x),x),x)+x*diff(y(x),x)^2+2*y(x)*diff(y(x),x) = 0, y(x))

Maple raw output

[y(x) = 1/x*(-2*x*(-_C2*x+_C1))^(1/2), y(x) = -1/x*(-2*x*(-_C2*x+_C1))^(1/2)]