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