dsolve(diff(diff(y(x),x),x) = -((alpha+beta+1)*x^2-(alpha+beta+1+a*(gamma1+delta)-delta)*x+a*gamma1)/x/(x-1)/(x-a)*diff(y(x),x)-(alpha*beta*x-q)/x/(x-1)/(x-a)*y(x),y(x))