ODE
\[ -2 a x y'(x)+a (a+1) y(x)+x^2 y''(x)=e^x x^{a+2} \] ODE Classification
[[_2nd_order, _linear, _nonhomogeneous]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.450798 (sec), leaf count = 303
\[\left \{\left \{y(x)\to (-x)^{\frac {1}{2} \left (-\sqrt {a} \sqrt {a+1} \sqrt {\frac {1}{a^2+a}}-1\right )} x^a \left (c_1 x^{\frac {1}{2} \left (1-\sqrt {a} \sqrt {a+1} \sqrt {\frac {1}{a^2+a}}\right )} (-x)^{\frac {1}{2} \left (\sqrt {a} \sqrt {a+1} \sqrt {\frac {1}{a^2+a}}+1\right )}+c_2 \left (-x^2\right )^{\frac {1}{2} \left (\sqrt {a} \sqrt {a+1} \sqrt {\frac {1}{a^2+a}}+1\right )}+\sqrt {a} \sqrt {a+1} \sqrt {\frac {1}{a^2+a}} x (-x)^{\sqrt {a} \sqrt {a+1} \sqrt {\frac {1}{a^2+a}}} \Gamma \left (\frac {3}{2}-\frac {1}{2} \sqrt {a} \sqrt {a+1} \sqrt {\frac {1}{a^2+a}},-x\right )-\sqrt {a} \sqrt {a+1} \sqrt {\frac {1}{a^2+a}} x \Gamma \left (\frac {1}{2} \sqrt {a} \sqrt {a+1} \sqrt {\frac {1}{a^2+a}}+\frac {3}{2},-x\right )\right )\right \}\right \}\]
Maple ✓
cpu = 0.023 (sec), leaf count = 23
\[[y \left (x \right ) = x^{a} \textit {\_C2} +x^{1+a} \textit {\_C1} +{\mathrm e}^{x} x^{a}]\] Mathematica raw input
DSolve[a*(1 + a)*y[x] - 2*a*x*y'[x] + x^2*y''[x] == E^x*x^(2 + a),y[x],x]
Mathematica raw output
{{y[x] -> (-x)^((-1 - Sqrt[a]*Sqrt[1 + a]*Sqrt[(a + a^2)^(-1)])/2)*x^a*((-x)^((1 + Sqrt[a]*Sqrt[1 + a]*Sqrt[(a + a^2)^(-1)])/2)*x^((1 - Sqrt[a]*Sqrt[1 + a]*Sqrt[(a + a^2)^(-1)])/2)*C[1] + (-x^2)^((1 + Sqrt[a]*Sqrt[1 + a]*Sqrt[(a + a^2)^(-1)])/2)*C[2] + Sqrt[a]*Sqrt[1 + a]*Sqrt[(a + a^2)^(-1)]*(-x)^(Sqrt[a]*Sqrt[1 + a]*Sqrt[(a + a^2)^(-1)])*x*Gamma[3/2 - (Sqrt[a]*Sqrt[1 + a]*Sqrt[(a + a^2)^(-1)])/2, -x] - Sqrt[a]*Sqrt[1 + a]*Sqrt[(a + a^2)^(-1)]*x*Gamma[3/2 + (Sqrt[a]*Sqrt[1 + a]*Sqrt[(a + a^2)^(-1)])/2, -x])}}
Maple raw input
dsolve(x^2*diff(diff(y(x),x),x)-2*a*x*diff(y(x),x)+a*(1+a)*y(x) = exp(x)*x^(a+2), y(x))
Maple raw output
[y(x) = x^a*_C2+x^(1+a)*_C1+exp(x)*x^a]