dsolve(x*(x-1)*diff(diff(y(x),x),x)+((a1+b1+1)*x-d1)*diff(y(x),x)+a1*b1*d1=0,y(x))