4.32.38 \(y(x) \left (\sum _{m=0}^n a(m) x^m\right )+a (1-2 x) y'(x)+x (x+1) y''(x)=0\)

ODE
\[ y(x) \left (\sum _{m=0}^n a(m) x^m\right )+a (1-2 x) y'(x)+x (x+1) y''(x)=0 \] ODE Classification

[[_2nd_order, _with_linear_symmetries]]

Book solution method
TO DO

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

DSolve[Sum[x^m*a[m], {m, 0, n}]*y[x] + a*(1 - 2*x)*Derivative[1][y][x] + x*(1 + x)*Derivative[2][y][x] == 0, y[x], x]

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

\[[]\]

Mathematica raw input

DSolve[Sum[x^m*a[m], {m, 0, n}]*y[x] + a*(1 - 2*x)*y'[x] + x*(1 + x)*y''[x] == 0,y[x],x]

Mathematica raw output

DSolve[Sum[x^m*a[m], {m, 0, n}]*y[x] + a*(1 - 2*x)*Derivative[1][y][x] + x*(1 + 
x)*Derivative[2][y][x] == 0, y[x], x]

Maple raw input

dsolve(x*(x+1)*diff(diff(y(x),x),x)+a*(1-2*x)*diff(y(x),x)+sum(a(m)*x^m,m = 0 .. n)*y(x) = 0, y(x))

Maple raw output

[]