ODE
\[ 16 y(x)^2 y'(x)^3+2 x y'(x)-y(x)=0 \] ODE Classification
[[_1st_order, _with_linear_symmetries]]
Book solution method
Change of variable
Mathematica ✓
cpu = 0.178115 (sec), leaf count = 20
\[\left \{\left \{y(x)\to \sqrt {c_1 \left (x+2 c_1{}^2\right )}\right \}\right \}\]
Maple ✓
cpu = 0.789 (sec), leaf count = 107
\[\left [y \left (x \right ) = -\frac {2^{\frac {1}{4}} 3^{\frac {1}{4}} \left (-x^{3}\right )^{\frac {1}{4}}}{3}, y \left (x \right ) = \frac {2^{\frac {1}{4}} 3^{\frac {1}{4}} \left (-x^{3}\right )^{\frac {1}{4}}}{3}, y \left (x \right ) = -\frac {i 2^{\frac {1}{4}} 3^{\frac {1}{4}} \left (-x^{3}\right )^{\frac {1}{4}}}{3}, y \left (x \right ) = \frac {i 2^{\frac {1}{4}} 3^{\frac {1}{4}} \left (-x^{3}\right )^{\frac {1}{4}}}{3}, y \left (x \right ) = \sqrt {16 \textit {\_C1}^{3}+2 x \textit {\_C1}}, y \left (x \right ) = -\sqrt {16 \textit {\_C1}^{3}+2 x \textit {\_C1}}\right ]\] Mathematica raw input
DSolve[-y[x] + 2*x*y'[x] + 16*y[x]^2*y'[x]^3 == 0,y[x],x]
Mathematica raw output
{{y[x] -> Sqrt[C[1]*(x + 2*C[1]^2)]}}
Maple raw input
dsolve(16*y(x)^2*diff(y(x),x)^3+2*x*diff(y(x),x)-y(x) = 0, y(x))
Maple raw output
[y(x) = -1/3*2^(1/4)*3^(1/4)*(-x^3)^(1/4), y(x) = 1/3*2^(1/4)*3^(1/4)*(-x^3)^(1/4), y(x) = -1/3*I*2^(1/4)*3^(1/4)*(-x^3)^(1/4), y(x) = 1/3*I*2^(1/4)*3^(1/4)*(-x^3)^(1/4), y(x) = (16*_C1^3+2*_C1*x)^(1/2), y(x) = -(16*_C1^3+2*_C1*x)^(1/2)]