DSolve[Derivative[1][y][x] == (Coth[1 + x] + x^2*Coth[1 + x] - Log[-1 + x] + 2*x*Coth[1 + x]*y[x] + Coth[1 + x]*y[x]^2)/Log[-1 + x],y[x],x]