dsolve(x*(x-1)*diff(diff(y(x),x),x)+((a+1)*x+b)*diff(y(x),x)=0,y(x))