[[_2nd_order, _with_linear_symmetries]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.577661 (sec), leaf count = 391
Maple ✓
cpu = 0.251 (sec), leaf count = 347
DSolve[(a2 + b2*x)*y[x] + x*(a1 + b1*x)*y'[x] + x^2*(a0 + x)*y''[x] == 0,y[x],x]
Mathematica raw output
{{y[x] -> (a0^(Sqrt[a0^2 + a1^2 - 2*a0*(a1 + 2*a2)]/a0)*C[1]*Hypergeometric2F1[-
(a1 + Sqrt[a0^2 + a1^2 - 2*a0*(a1 + 2*a2)] - a0*(b1 + Sqrt[1 - 2*b1 + b1^2 - 4*b
2]))/(2*a0), -(a1 + Sqrt[a0^2 + a1^2 - 2*a0*(a1 + 2*a2)] - a0*b1 + a0*Sqrt[1 - 2
*b1 + b1^2 - 4*b2])/(2*a0), 1 - Sqrt[a0^2 + a1^2 - 2*a0*(a1 + 2*a2)]/a0, -(x/a0)
] + x^(Sqrt[a0^2 + a1^2 - 2*a0*(a1 + 2*a2)]/a0)*C[2]*Hypergeometric2F1[(-a1 + Sq
rt[a0^2 + a1^2 - 2*a0*(a1 + 2*a2)] + a0*(b1 + Sqrt[1 - 2*b1 + b1^2 - 4*b2]))/(2*
a0), (-a1 + Sqrt[a0^2 + a1^2 - 2*a0*(a1 + 2*a2)] + a0*b1 - a0*Sqrt[1 - 2*b1 + b1
^2 - 4*b2])/(2*a0), (a0 + Sqrt[a0^2 + a1^2 - 2*a0*(a1 + 2*a2)])/a0, -(x/a0)])/(a
0^((a0 - a1 + Sqrt[a0^2 + a1^2 - 2*a0*(a1 + 2*a2)])/(2*a0))*x^((-a0 + a1 + Sqrt[
a0^2 + a1^2 - 2*a0*(a1 + 2*a2)])/(2*a0)))}}
Maple raw input
dsolve(x^2*(a0+x)*diff(diff(y(x),x),x)+x*(b1*x+a1)*diff(y(x),x)+(b2*x+a2)*y(x) = 0, y(x),'implicit')
Maple raw output
y(x) = (a0+x)^((-a0*b1+a0+a1)/a0)*(x^(1/2*(-a1+a0+(a0^2+(-2*a1-4*a2)*a0+a1^2)^(1
/2))/a0)*hypergeom([1/2*((a0^2+(-2*a1-4*a2)*a0+a1^2)^(1/2)+(b1^2-2*b1-4*b2+1)^(1
/2)*a0+(-b1+2)*a0+a1)/a0, 1/2*((a0^2+(-2*a1-4*a2)*a0+a1^2)^(1/2)-(b1^2-2*b1-4*b2
+1)^(1/2)*a0+(-b1+2)*a0+a1)/a0],[1/a0*(a0+(a0^2+(-2*a1-4*a2)*a0+a1^2)^(1/2))],-1
/a0*x)*_C1+x^(1/2*(-a1+a0-(a0^2+(-2*a1-4*a2)*a0+a1^2)^(1/2))/a0)*hypergeom([-1/2
/a0*((b1^2-2*b1-4*b2+1)^(1/2)*a0+b1*a0+(a0^2+(-2*a1-4*a2)*a0+a1^2)^(1/2)-2*a0-a1
), 1/2*(-(a0^2+(-2*a1-4*a2)*a0+a1^2)^(1/2)+(b1^2-2*b1-4*b2+1)^(1/2)*a0+(-b1+2)*a
0+a1)/a0],[(a0-(a0^2+(-2*a1-4*a2)*a0+a1^2)^(1/2))/a0],-1/a0*x)*_C2)