\[ y''(x)=-\frac {y(x) \left ((1-a) a x^2-b (b+x)\right )}{x^4} \] ✗ Mathematica : cpu = 0.697167 (sec), leaf count = 0 , DifferentialRoot result
\[\left \{\left \{y(x)\to \text {DifferentialRoot}\left (\{\unicode {f818},\unicode {f817}\}\unicode {f4a1}\left \{\unicode {f818}''(\unicode {f817}) \unicode {f817}^4+\left (-a^2 \unicode {f817}^2+a \unicode {f817}^2-b \unicode {f817}-b^2\right ) \unicode {f818}(\unicode {f817})=0,\unicode {f818}(1)=c_1,\unicode {f818}'(1)=c_2\right \}\right )(x)\right \}\right \}\]
✓ Maple : cpu = 0.108 (sec), leaf count = 58
\[ \left \{ y \left ( x \right ) ={{\sl I}_{a+1}\left ({\frac {b}{x}}\right )}{\it \_C1}\,b-{{\sl K}_{a+1}\left ({\frac {b}{x}}\right )}{\it \_C2}\,b+2\, \left ( {\it \_C1}\,{{\sl I}_{a}\left ({\frac {b}{x}}\right )}+{\it \_C2}\,{{\sl K}_{a}\left ({\frac {b}{x}}\right )} \right ) \left ( ax+b/2 \right ) \right \} \]