\[ y''(x)=-\frac {y(x) \left ((1-a) a x^2-b (b+x)\right )}{x^4} \] ✓ Mathematica : cpu = 0.188589 (sec), leaf count = 73
\[\left \{\left \{y(x)\to c_1 \left (2 \left (a x+\frac {b}{2}\right ) I_a\left (\frac {b}{x}\right )+b I_{a+1}\left (\frac {b}{x}\right )\right )+c_2 \left (2 \left (a x+\frac {b}{2}\right ) K_a\left (\frac {b}{x}\right )-b K_{a+1}\left (\frac {b}{x}\right )\right )\right \}\right \}\] ✓ Maple : cpu = 0.132 (sec), leaf count = 58
\[\left \{y \left (x \right ) = c_{1} b \BesselI \left (a +1, \frac {b}{x}\right )-c_{2} b \BesselK \left (a +1, \frac {b}{x}\right )+2 \left (c_{1} \BesselI \left (a , \frac {b}{x}\right )+c_{2} \BesselK \left (a , \frac {b}{x}\right )\right ) \left (a x +\frac {b}{2}\right )\right \}\]