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