4.40.27 \(x y(x) y''(x)=a y(x) y'(x)+b^2 x y(x)^3+x y'(x)^2\)

ODE
\[ x y(x) y''(x)=a y(x) y'(x)+b^2 x y(x)^3+x y'(x)^2 \] ODE Classification

[[_2nd_order, _with_linear_symmetries]]

Book solution method
TO DO

Mathematica
cpu = 55.5149 (sec), leaf count = 0 , could not solve

DSolve[x*y[x]*Derivative[2][y][x] == b^2*x*y[x]^3 + a*y[x]*Derivative[1][y][x] + x*Derivative[1][y][x]^2, y[x], x]

Maple
cpu = 2.574 (sec), leaf count = 0 , result contains DESol or ODESolStruc

\[[]\]

Mathematica raw input

DSolve[x*y[x]*y''[x] == b^2*x*y[x]^3 + a*y[x]*y'[x] + x*y'[x]^2,y[x],x]

Mathematica raw output

DSolve[x*y[x]*Derivative[2][y][x] == b^2*x*y[x]^3 + a*y[x]*Derivative[1][y][x] +
 x*Derivative[1][y][x]^2, y[x], x]

Maple raw input

dsolve(x*y(x)*diff(diff(y(x),x),x) = x*diff(y(x),x)^2+a*y(x)*diff(y(x),x)+b^2*x*y(x)^3, y(x))

Maple raw output

[]