[[_Emden, _Fowler]]
Book solution method
TO DO
Mathematica ✓
cpu = 1.46083 (sec), leaf count = 199
Maple ✓
cpu = 0.148 (sec), leaf count = 178
DSolve[A*y[x] + (a + b*x + c*x^2)^2*y''[x] == 0,y[x],x]
Mathematica raw output
{{y[x] -> (Sqrt[a + x*(b + c*x)]*(E^((2*Sqrt[-b^2 + 4*a*c]*Sqrt[1 - (4*A)/(b^2 -
4*a*c)]*ArcTan[(b + 2*c*x)/Sqrt[-b^2 + 4*a*c]])/Sqrt[b^2 - 4*a*c])*C[1] + C[2]/
(Sqrt[b^2 - 4*a*c]*Sqrt[1 - (4*A)/(b^2 - 4*a*c)])))/E^((Sqrt[-b^2 + 4*a*c]*Sqrt[
1 - (4*A)/(b^2 - 4*a*c)]*ArcTan[(b + 2*c*x)/Sqrt[-b^2 + 4*a*c]])/Sqrt[b^2 - 4*a*
c])}}
Maple raw input
dsolve((c*x^2+b*x+a)^2*diff(diff(y(x),x),x)+A*y(x) = 0, y(x),'implicit')
Maple raw output
y(x) = (c*x^2+b*x+a)^(1/2)*(((I*(4*a*c-b^2)^(1/2)-2*c*x-b)/(2*c*x+b+I*(4*a*c-b^2
)^(1/2)))^(-1/2*c/(-4*a*c+b^2)^(1/2)*((-4*a*c+b^2-4*A)/c^2)^(1/2))*_C2+((I*(4*a*
c-b^2)^(1/2)-2*c*x-b)/(2*c*x+b+I*(4*a*c-b^2)^(1/2)))^(1/2*c/(-4*a*c+b^2)^(1/2)*(
(-4*a*c+b^2-4*A)/c^2)^(1/2))*_C1)