[[_Emden, _Fowler]]
Book solution method
TO DO
Mathematica ✓
cpu = 0.0112846 (sec), leaf count = 26
Maple ✓
cpu = 0.005 (sec), leaf count = 23
DSolve[13*y[x] + 5*x*y'[x] + x^2*y''[x] == 0,y[x],x]
Mathematica raw output
{{y[x] -> (C[2]*Cos[3*Log[x]] + C[1]*Sin[3*Log[x]])/x^2}}
Maple raw input
dsolve(x^2*diff(diff(y(x),x),x)+5*x*diff(y(x),x)+13*y(x) = 0, y(x),'implicit')
Maple raw output
y(x) = (_C1*sin(3*ln(x))+_C2*cos(3*ln(x)))/x^2