ODE
\[ -2 a x y'(x)+a (a+1) y(x)+x^2 y''(x)=0 \] ODE Classification
[[_2nd_order, _with_linear_symmetries]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.166191 (sec), leaf count = 105
\[\left \{\left \{y(x)\to x^{\frac {-\sqrt {\frac {1}{a^2+a}} \sqrt {a}-\sqrt {\frac {1}{a^2+a}} a^{3/2}+2 \sqrt {a+1} a+\sqrt {a+1}}{2 \sqrt {a+1}}} \left (c_2 x^{\sqrt {a} \sqrt {a+1} \sqrt {\frac {1}{a^2+a}}}+c_1\right )\right \}\right \}\]
Maple ✓
cpu = 0.014 (sec), leaf count = 17
\[[y \left (x \right ) = \textit {\_C1} \,x^{a}+\textit {\_C2} \,x^{1+a}]\] Mathematica raw input
DSolve[a*(1 + a)*y[x] - 2*a*x*y'[x] + x^2*y''[x] == 0,y[x],x]
Mathematica raw output
{{y[x] -> x^((Sqrt[1 + a] + 2*a*Sqrt[1 + a] - Sqrt[a]*Sqrt[(a + a^2)^(-1)] - a^(3/2)*Sqrt[(a + a^2)^(-1)])/(2*Sqrt[1 + a]))*(C[1] + x^(Sqrt[a]*Sqrt[1 + a]*Sqrt[(a + a^2)^(-1)])*C[2])}}
Maple raw input
dsolve(x^2*diff(diff(y(x),x),x)-2*a*x*diff(y(x),x)+a*(1+a)*y(x) = 0, y(x))
Maple raw output
[y(x) = _C1*x^a+_C2*x^(1+a)]