dsolve(diff(y(x),x) = -y(x)*(ln(x-1)+coth(1+x)*x-coth(1+x)*x^2*y(x))/x/ln(x-1),y(x))