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