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