dv/v=xdx/(x^2 -1); lnv=xln(x-1/[x+1])/2 -1/2int ln(x-1/[x+1])dx;
lnv=1/2int (x+1)xdx/(x-1)