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