ODE
\[ a x y'(x)+x^2 y(x) \left (\text {a1}+\text {b1} x^2\right )+x^2 y''(x)=0 \] ODE Classification
[[_2nd_order, _with_linear_symmetries]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.296505 (sec), leaf count = 136
\[\left \{\left \{y(x)\to 2^{\frac {a+1}{4}} x^{\frac {1}{2} (-a-1)} \left (x^2\right )^{\frac {a+1}{4}} e^{\frac {1}{2} i \sqrt {\text {b1}} x^2} \left (c_1 U\left (\frac {1}{4} \left (a-\frac {i \text {a1}}{\sqrt {\text {b1}}}+1\right ),\frac {a+1}{2},-i \sqrt {\text {b1}} x^2\right )+c_2 L_{\frac {1}{4} \left (-a+\frac {i \text {a1}}{\sqrt {\text {b1}}}-1\right )}^{\frac {a-1}{2}}\left (-i \sqrt {\text {b1}} x^2\right )\right )\right \}\right \}\]
Maple ✓
cpu = 0.58 (sec), leaf count = 95
\[\left [y \left (x \right ) = \textit {\_C1} \,{\mathrm e}^{-\frac {i \sqrt {\mathit {b1}}\, x^{2}}{2}} \KummerM \left (\frac {i \mathit {a1} +a \sqrt {\mathit {b1}}+\sqrt {\mathit {b1}}}{4 \sqrt {\mathit {b1}}}, \frac {a}{2}+\frac {1}{2}, i \sqrt {\mathit {b1}}\, x^{2}\right )+\textit {\_C2} \,{\mathrm e}^{-\frac {i \sqrt {\mathit {b1}}\, x^{2}}{2}} \KummerU \left (\frac {i \mathit {a1} +a \sqrt {\mathit {b1}}+\sqrt {\mathit {b1}}}{4 \sqrt {\mathit {b1}}}, \frac {a}{2}+\frac {1}{2}, i \sqrt {\mathit {b1}}\, x^{2}\right )\right ]\] Mathematica raw input
DSolve[x^2*(a1 + b1*x^2)*y[x] + a*x*y'[x] + x^2*y''[x] == 0,y[x],x]
Mathematica raw output
{{y[x] -> 2^((1 + a)/4)*E^((I/2)*Sqrt[b1]*x^2)*x^((-1 - a)/2)*(x^2)^((1 + a)/4)*(C[1]*HypergeometricU[(1 + a - (I*a1)/Sqrt[b1])/4, (1 + a)/2, (-I)*Sqrt[b1]*x^2] + C[2]*LaguerreL[(-1 - a + (I*a1)/Sqrt[b1])/4, (-1 + a)/2, (-I)*Sqrt[b1]*x^2])}}
Maple raw input
dsolve(x^2*diff(diff(y(x),x),x)+a*x*diff(y(x),x)+x^2*(b1*x^2+a1)*y(x) = 0, y(x))
Maple raw output
[y(x) = _C1*exp(-1/2*I*b1^(1/2)*x^2)*KummerM(1/4*(I*a1+a*b1^(1/2)+b1^(1/2))/b1^(1/2),1/2*a+1/2,I*b1^(1/2)*x^2)+_C2*exp(-1/2*I*b1^(1/2)*x^2)*KummerU(1/4*(I*a1+a*b1^(1/2)+b1^(1/2))/b1^(1/2),1/2*a+1/2,I*b1^(1/2)*x^2)]