ODE
\[ x^4 y''(x)=x \left (x^2+2 y(x)\right ) y'(x)-4 y(x)^2 \] ODE Classification
[[_2nd_order, _with_linear_symmetries], [_2nd_order, _reducible, _mu_x_y1], [_2nd_order, _reducible, _mu_xy]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.242225 (sec), leaf count = 83
\[\left \{\left \{y(x)\to \frac {x^2 \left (\left (1-i \sqrt {-1-c_1}\right ) x^{2 i \sqrt {-1-c_1}}+\left (1+i \sqrt {-1-c_1}\right ) c_2\right )}{x^{2 i \sqrt {-1-c_1}}+c_2}\right \}\right \}\]
Maple ✓
cpu = 0.698 (sec), leaf count = 23
\[[y \left (x \right ) = \tanh \left (-\ln \left (x \right ) \textit {\_C1} +\textit {\_C2} \textit {\_C1} \right ) x^{2} \textit {\_C1} +x^{2}]\] Mathematica raw input
DSolve[x^4*y''[x] == -4*y[x]^2 + x*(x^2 + 2*y[x])*y'[x],y[x],x]
Mathematica raw output
{{y[x] -> (x^2*(x^((2*I)*Sqrt[-1 - C[1]])*(1 - I*Sqrt[-1 - C[1]]) + (1 + I*Sqrt[-1 - C[1]])*C[2]))/(x^((2*I)*Sqrt[-1 - C[1]]) + C[2])}}
Maple raw input
dsolve(x^4*diff(diff(y(x),x),x) = x*(x^2+2*y(x))*diff(y(x),x)-4*y(x)^2, y(x))
Maple raw output
[y(x) = tanh(-ln(x)*_C1+_C2*_C1)*x^2*_C1+x^2]