DSolve[Derivative[2][y][x] == -(((-q + alpha*beta*x)*y[x])/((-1 + x)*x*(-a + x))) - ((a*gamma1 - (1 + alpha + beta - delta + a*(delta + gamma1))*x + (1 + alpha + beta)*x^2)*Derivative[1][y][x])/((-1 + x)*x*(-a + x)),y[x],x]