ODE
\[ y(x) y'''(x)+y(x)^3 y'(x)-y'(x) y''(x)=0 \] ODE Classification
[[_3rd_order, _missing_x], [_3rd_order, _exact, _nonlinear], [_3rd_order, _with_linear_symmetries]]
Book solution method
TO DO
Mathematica ✓
cpu = 2.7443 (sec), leaf count = 409
\[\left \{\left \{y(x)\to \text {InverseFunction}\left [-\frac {2 i \sqrt {\frac {\text {$\#$1}^2}{2 \left (\sqrt {c_2^2-c_1}-c_2\right )}+1} \sqrt {1-\frac {\text {$\#$1}^2}{2 \left (c_2+\sqrt {c_2^2-c_1}\right )}} F\left (i \sinh ^{-1}\left (\frac {\sqrt {\frac {1}{\sqrt {c_2^2-c_1}-c_2}} \text {$\#$1}}{\sqrt {2}}\right )|\frac {c_2-\sqrt {c_2^2-c_1}}{c_2+\sqrt {c_2^2-c_1}}\right )}{\sqrt {\frac {1}{\sqrt {c_2^2-c_1}-c_2}} \sqrt {-\frac {\text {$\#$1}^4}{2}+2 \text {$\#$1}^2 c_2-2 c_1}}\& \right ]\left [c_3+x\right ]\right \},\left \{y(x)\to \text {InverseFunction}\left [\frac {2 i \sqrt {\frac {\text {$\#$1}^2}{2 \left (\sqrt {c_2^2-c_1}-c_2\right )}+1} \sqrt {1-\frac {\text {$\#$1}^2}{2 \left (c_2+\sqrt {c_2^2-c_1}\right )}} F\left (i \sinh ^{-1}\left (\frac {\sqrt {\frac {1}{\sqrt {c_2^2-c_1}-c_2}} \text {$\#$1}}{\sqrt {2}}\right )|\frac {c_2-\sqrt {c_2^2-c_1}}{c_2+\sqrt {c_2^2-c_1}}\right )}{\sqrt {\frac {1}{\sqrt {c_2^2-c_1}-c_2}} \sqrt {-\frac {\text {$\#$1}^4}{2}+2 \text {$\#$1}^2 c_2-2 c_1}}\& \right ]\left [c_3+x\right ]\right \}\right \}\]
Maple ✓
cpu = 0.309 (sec), leaf count = 77
\[ \left \{ \int ^{y \left ( x \right ) }\!-2\,{\frac {1}{\sqrt {-{{\it \_a}}^{4}+4\,{\it \_C2}\,{{\it \_a}}^{2}-4\,{{\it \_C2}}^{2}+4\,{\it \_C1}}}}{d{\it \_a}}-x-{\it \_C3}=0,\int ^{y \left ( x \right ) }\!2\,{\frac {1}{\sqrt {-{{\it \_a}}^{4}+4\,{\it \_C2}\,{{\it \_a}}^{2}-4\,{{\it \_C2}}^{2}+4\,{\it \_C1}}}}{d{\it \_a}}-x-{\it \_C3}=0 \right \} \] Mathematica raw input
DSolve[y[x]^3*y'[x] - y'[x]*y''[x] + y[x]*y'''[x] == 0,y[x],x]
Mathematica raw output
{{y[x] -> InverseFunction[((-2*I)*EllipticF[I*ArcSinh[(Sqrt[(-C[2] + Sqrt[-C[1] + C[2]^2])^(-1)]*#1)/Sqrt[2]], (C[2] - Sqrt[-C[1] + C[2]^2])/(C[2] + Sqrt[-C[1] + C[2]^2])]*Sqrt[1 + #1^2/(2*(-C[2] + Sqrt[-C[1] + C[2]^2]))]*Sqrt[1 - #1^2/(2*(C[2] + Sqrt[-C[1] + C[2]^2]))])/(Sqrt[(-C[2] + Sqrt[-C[1] + C[2]^2])^(-1)]*Sqrt[-2*C[1] + 2*C[2]*#1^2 - #1^4/2]) & ][x + C[3]]}, {y[x] -> InverseFunction[((2*I)*EllipticF[I*ArcSinh[(Sqrt[(-C[2] + Sqrt[-C[1] + C[2]^2])^(-1)]*#1)/Sqrt[2]], (C[2] - Sqrt[-C[1] + C[2]^2])/(C[2] + Sqrt[-C[1] + C[2]^2])]*Sqrt[1 + #1^2/(2*(-C[2] + Sqrt[-C[1] + C[2]^2]))]*Sqrt[1 - #1^2/(2*(C[2] + Sqrt[-C[1] + C[2]^2]))])/(Sqrt[(-C[2] + Sqrt[-C[1] + C[2]^2])^(-1)]*Sqrt[-2*C[1] + 2*C[2]*#1^2 - #1^4/2]) & ][x + C[3]]}}
Maple raw input
dsolve(y(x)*diff(diff(diff(y(x),x),x),x)-diff(y(x),x)*diff(diff(y(x),x),x)+y(x)^3*diff(y(x),x) = 0, y(x),'implicit')
Maple raw output
Intat(-2/(-_a^4+4*_C2*_a^2-4*_C2^2+4*_C1)^(1/2),_a = y(x))-x-_C3 = 0, Intat(2/(-_a^4+4*_C2*_a^2-4*_C2^2+4*_C1)^(1/2),_a = y(x))-x-_C3 = 0