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