((x1*exp(x1/gamma) + x2*exp(x2/gamma) + x3*exp(x3/gamma))/(exp(x1/gamma) + exp(x2/gamma) + exp(x3/gamma)) - (x1*exp(-x1/gamma) + x2*exp(-x2/gamma) + x3*exp(-x3/gamma))/(exp(-x1/gamma) + exp(-x2/gamma) + exp(-x3/gamma)))

Specification

?
\[\left(\left(\left(1 \leq x1 \land x1 \leq 1.79 \cdot 10^{+308}\right) \land \left(0 \leq \gamma \land \gamma \leq 1\right)\right) \land \left(1 \leq x2 \land x2 \leq 1.79 \cdot 10^{+308}\right)\right) \land \left(1 \leq x3 \land x3 \leq 1.79 \cdot 10^{+308}\right)\]
\[\begin{array}{l} \\ \begin{array}{l} t_0 := e^{\frac{x1}{\gamma}}\\ t_1 := e^{\frac{-x1}{\gamma}}\\ t_2 := e^{\frac{x3}{\gamma}}\\ t_3 := e^{\frac{-x3}{\gamma}}\\ t_4 := e^{\frac{x2}{\gamma}}\\ t_5 := e^{\frac{-x2}{\gamma}}\\ \frac{\left(x1 \cdot t\_0 + x2 \cdot t\_4\right) + x3 \cdot t\_2}{\left(t\_0 + t\_4\right) + t\_2} - \frac{\left(x1 \cdot t\_1 + x2 \cdot t\_5\right) + x3 \cdot t\_3}{\left(t\_1 + t\_5\right) + t\_3} \end{array} \end{array} \]
(FPCore (x1 gamma x2 x3)
 :precision binary64
 (let* ((t_0 (exp (/ x1 gamma)))
        (t_1 (exp (/ (- x1) gamma)))
        (t_2 (exp (/ x3 gamma)))
        (t_3 (exp (/ (- x3) gamma)))
        (t_4 (exp (/ x2 gamma)))
        (t_5 (exp (/ (- x2) gamma))))
   (-
    (/ (+ (+ (* x1 t_0) (* x2 t_4)) (* x3 t_2)) (+ (+ t_0 t_4) t_2))
    (/ (+ (+ (* x1 t_1) (* x2 t_5)) (* x3 t_3)) (+ (+ t_1 t_5) t_3)))))
double code(double x1, double gamma, double x2, double x3) {
	double t_0 = exp((x1 / gamma));
	double t_1 = exp((-x1 / gamma));
	double t_2 = exp((x3 / gamma));
	double t_3 = exp((-x3 / gamma));
	double t_4 = exp((x2 / gamma));
	double t_5 = exp((-x2 / gamma));
	return ((((x1 * t_0) + (x2 * t_4)) + (x3 * t_2)) / ((t_0 + t_4) + t_2)) - ((((x1 * t_1) + (x2 * t_5)) + (x3 * t_3)) / ((t_1 + t_5) + t_3));
}

real(8) function code(x1, gamma, x2, x3)
    real(8), intent (in) :: x1
    real(8), intent (in) :: gamma
    real(8), intent (in) :: x2
    real(8), intent (in) :: x3
    real(8) :: t_0
    real(8) :: t_1
    real(8) :: t_2
    real(8) :: t_3
    real(8) :: t_4
    real(8) :: t_5
    t_0 = exp((x1 / gamma))
    t_1 = exp((-x1 / gamma))
    t_2 = exp((x3 / gamma))
    t_3 = exp((-x3 / gamma))
    t_4 = exp((x2 / gamma))
    t_5 = exp((-x2 / gamma))
    code = ((((x1 * t_0) + (x2 * t_4)) + (x3 * t_2)) / ((t_0 + t_4) + t_2)) - ((((x1 * t_1) + (x2 * t_5)) + (x3 * t_3)) / ((t_1 + t_5) + t_3))
end function

public static double code(double x1, double gamma, double x2, double x3) {
	double t_0 = Math.exp((x1 / gamma));
	double t_1 = Math.exp((-x1 / gamma));
	double t_2 = Math.exp((x3 / gamma));
	double t_3 = Math.exp((-x3 / gamma));
	double t_4 = Math.exp((x2 / gamma));
	double t_5 = Math.exp((-x2 / gamma));
	return ((((x1 * t_0) + (x2 * t_4)) + (x3 * t_2)) / ((t_0 + t_4) + t_2)) - ((((x1 * t_1) + (x2 * t_5)) + (x3 * t_3)) / ((t_1 + t_5) + t_3));
}

def code(x1, gamma, x2, x3):
	t_0 = math.exp((x1 / gamma))
	t_1 = math.exp((-x1 / gamma))
	t_2 = math.exp((x3 / gamma))
	t_3 = math.exp((-x3 / gamma))
	t_4 = math.exp((x2 / gamma))
	t_5 = math.exp((-x2 / gamma))
	return ((((x1 * t_0) + (x2 * t_4)) + (x3 * t_2)) / ((t_0 + t_4) + t_2)) - ((((x1 * t_1) + (x2 * t_5)) + (x3 * t_3)) / ((t_1 + t_5) + t_3))

function code(x1, gamma, x2, x3)
	t_0 = exp(Float64(x1 / gamma))
	t_1 = exp(Float64(Float64(-x1) / gamma))
	t_2 = exp(Float64(x3 / gamma))
	t_3 = exp(Float64(Float64(-x3) / gamma))
	t_4 = exp(Float64(x2 / gamma))
	t_5 = exp(Float64(Float64(-x2) / gamma))
	return Float64(Float64(Float64(Float64(Float64(x1 * t_0) + Float64(x2 * t_4)) + Float64(x3 * t_2)) / Float64(Float64(t_0 + t_4) + t_2)) - Float64(Float64(Float64(Float64(x1 * t_1) + Float64(x2 * t_5)) + Float64(x3 * t_3)) / Float64(Float64(t_1 + t_5) + t_3)))
end

function tmp = code(x1, gamma, x2, x3)
	t_0 = exp((x1 / gamma));
	t_1 = exp((-x1 / gamma));
	t_2 = exp((x3 / gamma));
	t_3 = exp((-x3 / gamma));
	t_4 = exp((x2 / gamma));
	t_5 = exp((-x2 / gamma));
	tmp = ((((x1 * t_0) + (x2 * t_4)) + (x3 * t_2)) / ((t_0 + t_4) + t_2)) - ((((x1 * t_1) + (x2 * t_5)) + (x3 * t_3)) / ((t_1 + t_5) + t_3));
end

code[x1_, gamma_, x2_, x3_] := Block[{t$95$0 = N[Exp[N[(x1 / gamma), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$1 = N[Exp[N[((-x1) / gamma), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$2 = N[Exp[N[(x3 / gamma), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$3 = N[Exp[N[((-x3) / gamma), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$4 = N[Exp[N[(x2 / gamma), $MachinePrecision]], $MachinePrecision]}, Block[{t$95$5 = N[Exp[N[((-x2) / gamma), $MachinePrecision]], $MachinePrecision]}, N[(N[(N[(N[(N[(x1 * t$95$0), $MachinePrecision] + N[(x2 * t$95$4), $MachinePrecision]), $MachinePrecision] + N[(x3 * t$95$2), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$0 + t$95$4), $MachinePrecision] + t$95$2), $MachinePrecision]), $MachinePrecision] - N[(N[(N[(N[(x1 * t$95$1), $MachinePrecision] + N[(x2 * t$95$5), $MachinePrecision]), $MachinePrecision] + N[(x3 * t$95$3), $MachinePrecision]), $MachinePrecision] / N[(N[(t$95$1 + t$95$5), $MachinePrecision] + t$95$3), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]]]]]]

\begin{array}{l}

\\
\begin{array}{l}
t_0 := e^{\frac{x1}{\gamma}}\\
t_1 := e^{\frac{-x1}{\gamma}}\\
t_2 := e^{\frac{x3}{\gamma}}\\
t_3 := e^{\frac{-x3}{\gamma}}\\
t_4 := e^{\frac{x2}{\gamma}}\\
t_5 := e^{\frac{-x2}{\gamma}}\\
\frac{\left(x1 \cdot t\_0 + x2 \cdot t\_4\right) + x3 \cdot t\_2}{\left(t\_0 + t\_4\right) + t\_2} - \frac{\left(x1 \cdot t\_1 + x2 \cdot t\_5\right) + x3 \cdot t\_3}{\left(t\_1 + t\_5\right) + t\_3}
\end{array}
\end{array}

Cannot sample enough valid points. (more)