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