-(1/th)*log{1+((e^(-th*x)-1)*(e^(-th*y))-1)/((e^-th)-1)}