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