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