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