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