4.39.45 2y(x)y(x)=y(x)2+8y(x)3+4y(x)2

ODE
2y(x)y(x)=y(x)2+8y(x)3+4y(x)2 ODE Classification

[[_2nd_order, _missing_x]]

Book solution method
TO DO

Mathematica
cpu = 1.18291 (sec), leaf count = 359

{{y(x)InverseFunction[2i#1c1#1(221c1)+1c1#1(21c1+2)+1F(isinh1(c121c1+2#1)|1c1+111c1)c121c1+24#12+4#1+c1&][c2+x]},{y(x)InverseFunction[2i#1c1#1(221c1)+1c1#1(21c1+2)+1F(isinh1(c121c1+2#1)|1c1+111c1)c121c1+24#12+4#1+c1&][c2+x]}}

Maple
cpu = 0.075 (sec), leaf count = 61

{y(x)14_a3+_a_C1+4_a2d_ax_C2=0,y(x)1(4_a2+_C1+4_a)_ad_ax_C2=0} Mathematica raw input

DSolve[2*y[x]*y''[x] == 4*y[x]^2 + 8*y[x]^3 + y'[x]^2,y[x],x]

Mathematica raw output

{{y[x] -> InverseFunction[((-2*I)*EllipticF[I*ArcSinh[Sqrt[C[1]/(2 + 2*Sqrt[1 - 
C[1]])]/Sqrt[#1]], (1 + Sqrt[1 - C[1]])/(1 - Sqrt[1 - C[1]])]*Sqrt[1 + C[1]/((2 
- 2*Sqrt[1 - C[1]])*#1)]*Sqrt[1 + C[1]/((2 + 2*Sqrt[1 - C[1]])*#1)]*#1)/(Sqrt[C[
1]/(2 + 2*Sqrt[1 - C[1]])]*Sqrt[C[1] + 4*#1 + 4*#1^2]) & ][x + C[2]]}, {y[x] -> 
InverseFunction[((2*I)*EllipticF[I*ArcSinh[Sqrt[C[1]/(2 + 2*Sqrt[1 - C[1]])]/Sqr
t[#1]], (1 + Sqrt[1 - C[1]])/(1 - Sqrt[1 - C[1]])]*Sqrt[1 + C[1]/((2 - 2*Sqrt[1 
- C[1]])*#1)]*Sqrt[1 + C[1]/((2 + 2*Sqrt[1 - C[1]])*#1)]*#1)/(Sqrt[C[1]/(2 + 2*S
qrt[1 - C[1]])]*Sqrt[C[1] + 4*#1 + 4*#1^2]) & ][x + C[2]]}}

Maple raw input

dsolve(2*y(x)*diff(diff(y(x),x),x) = diff(y(x),x)^2+4*y(x)^2+8*y(x)^3, y(x),'implicit')

Maple raw output

Intat(1/(4*_a^3+_C1*_a+4*_a^2)^(1/2),_a = y(x))-x-_C2 = 0, Intat(-1/((4*_a^2+_C1
+4*_a)*_a)^(1/2),_a = y(x))-x-_C2 = 0