DSolve[Derivative[2][y][x] == -((b*y[x])/(x^2*(a + x^2))) - ((a + 2*x^2)*Derivative[1][y][x])/(x*(a + x^2)),y[x],x]