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