ODE
\[ -4 x^2 y(x) y'(x)^2+y'(x)^4+16 x y(x)^2 y'(x)-16 y(x)^3=0 \] ODE Classification
[[_homogeneous, `class G`]]
Book solution method
No Missing Variables ODE, Solve for \(x\)
Mathematica ✓
cpu = 3.16828 (sec), leaf count = 321
\[\left \{\text {Solve}\left [\frac {\sqrt {x^2-4 \sqrt {y(x)}} \sqrt {y(x)} \left (4 \log \left (\sqrt {x^2-4 \sqrt {y(x)}}+x\right )-\log (y(x))\right )}{\sqrt {\left (x^2-4 \sqrt {y(x)}\right ) y(x)}}+\log (y(x))=4 c_1,y(x)\right ],\text {Solve}\left [\frac {\sqrt {x^2-4 \sqrt {y(x)}} \sqrt {y(x)} \left (\log (y(x))-4 \log \left (\sqrt {x^2-4 \sqrt {y(x)}}+x\right )\right )}{\sqrt {\left (x^2-4 \sqrt {y(x)}\right ) y(x)}}+\log (y(x))=4 c_1,y(x)\right ],\text {Solve}\left [\frac {\sqrt {x^2+4 \sqrt {y(x)}} \sqrt {y(x)} \left (4 \log \left (\sqrt {x^2+4 \sqrt {y(x)}}+x\right )-\log (y(x))\right )}{\sqrt {\left (x^2+4 \sqrt {y(x)}\right ) y(x)}}+\log (y(x))=4 c_1,y(x)\right ],\text {Solve}\left [\frac {\sqrt {x^2+4 \sqrt {y(x)}} \sqrt {y(x)} \left (\log (y(x))-4 \log \left (\sqrt {x^2+4 \sqrt {y(x)}}+x\right )\right )}{\sqrt {\left (x^2+4 \sqrt {y(x)}\right ) y(x)}}+\log (y(x))=4 c_1,y(x)\right ]\right \}\]
Maple ✓
cpu = 0.587 (sec), leaf count = 118
\[\left [y \left (x \right ) = \frac {x^{4}}{16}, \frac {\left (\sqrt {x^{2}-4 \sqrt {y \left (x \right )}}+x \right )^{-\frac {\sqrt {x^{2} y \left (x \right )-4 y \left (x \right )^{\frac {3}{2}}}}{\sqrt {x^{2}-4 \sqrt {y \left (x \right )}}\, \sqrt {y \left (x \right )}}} \left (\sqrt {x^{2}-4 \sqrt {y \left (x \right )}}-x \right )^{\frac {\sqrt {x^{2} y \left (x \right )-4 y \left (x \right )^{\frac {3}{2}}}}{\sqrt {x^{2}-4 \sqrt {y \left (x \right )}}\, \sqrt {y \left (x \right )}}}}{\sqrt {y \left (x \right )}}-\textit {\_C1} = 0\right ]\] Mathematica raw input
DSolve[-16*y[x]^3 + 16*x*y[x]^2*y'[x] - 4*x^2*y[x]*y'[x]^2 + y'[x]^4 == 0,y[x],x]
Mathematica raw output
{Solve[Log[y[x]] + ((4*Log[x + Sqrt[x^2 - 4*Sqrt[y[x]]]] - Log[y[x]])*Sqrt[x^2 - 4*Sqrt[y[x]]]*Sqrt[y[x]])/Sqrt[(x^2 - 4*Sqrt[y[x]])*y[x]] == 4*C[1], y[x]], Solve[Log[y[x]] + ((-4*Log[x + Sqrt[x^2 - 4*Sqrt[y[x]]]] + Log[y[x]])*Sqrt[x^2 - 4*Sqrt[y[x]]]*Sqrt[y[x]])/Sqrt[(x^2 - 4*Sqrt[y[x]])*y[x]] == 4*C[1], y[x]], Solve[Log[y[x]] + ((4*Log[x + Sqrt[x^2 + 4*Sqrt[y[x]]]] - Log[y[x]])*Sqrt[x^2 + 4*Sqrt[y[x]]]*Sqrt[y[x]])/Sqrt[(x^2 + 4*Sqrt[y[x]])*y[x]] == 4*C[1], y[x]], Solve[Log[y[x]] + ((-4*Log[x + Sqrt[x^2 + 4*Sqrt[y[x]]]] + Log[y[x]])*Sqrt[x^2 + 4*Sqrt[y[x]]]*Sqrt[y[x]])/Sqrt[(x^2 + 4*Sqrt[y[x]])*y[x]] == 4*C[1], y[x]]}
Maple raw input
dsolve(diff(y(x),x)^4-4*x^2*y(x)*diff(y(x),x)^2+16*x*y(x)^2*diff(y(x),x)-16*y(x)^3 = 0, y(x))
Maple raw output
[y(x) = 1/16*x^4, 1/y(x)^(1/2)/(((x^2-4*y(x)^(1/2))^(1/2)+x)^(1/(x^2-4*y(x)^(1/2))^(1/2)/y(x)^(1/2)*(x^2*y(x)-4*y(x)^(3/2))^(1/2)))*((x^2-4*y(x)^(1/2))^(1/2)-x)^(1/(x^2-4*y(x)^(1/2))^(1/2)/y(x)^(1/2)*(x^2*y(x)-4*y(x)^(3/2))^(1/2))-_C1 = 0]