int x/2*(1+coth(b*x/2))*dx