\[\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1}
\]
↓
\[\begin{array}{l}
t_0 := \frac{\beta + i \cdot 2}{i + \beta}\\
t_1 := \mathsf{fma}\left(i, 2, \alpha + \beta\right)\\
t_2 := t_1 + -1\\
t_3 := 1 + t_1\\
\mathbf{if}\;\alpha \leq 3.1 \cdot 10^{+103}:\\
\;\;\;\;\frac{\frac{i}{t_3}}{t_0} \cdot \frac{\frac{i}{t_0}}{t_2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{i}{\frac{t_1}{i + \left(\alpha + \beta\right)}}}{t_3} \cdot \frac{\alpha + i}{t_2}\\
\end{array}
\]
double code(double alpha, double beta, double i) {
return (((i * ((alpha + beta) + i)) * ((beta * alpha) + (i * ((alpha + beta) + i)))) / (((alpha + beta) + (2.0 * i)) * ((alpha + beta) + (2.0 * i)))) / ((((alpha + beta) + (2.0 * i)) * ((alpha + beta) + (2.0 * i))) - 1.0);
}
↓
double code(double alpha, double beta, double i) {
double t_0 = (beta + (i * 2.0)) / (i + beta);
double t_1 = fma(i, 2.0, (alpha + beta));
double t_2 = t_1 + -1.0;
double t_3 = 1.0 + t_1;
double tmp;
if (alpha <= 3.1e+103) {
tmp = ((i / t_3) / t_0) * ((i / t_0) / t_2);
} else {
tmp = ((i / (t_1 / (i + (alpha + beta)))) / t_3) * ((alpha + i) / t_2);
}
return tmp;
}
function code(alpha, beta, i)
return Float64(Float64(Float64(Float64(i * Float64(Float64(alpha + beta) + i)) * Float64(Float64(beta * alpha) + Float64(i * Float64(Float64(alpha + beta) + i)))) / Float64(Float64(Float64(alpha + beta) + Float64(2.0 * i)) * Float64(Float64(alpha + beta) + Float64(2.0 * i)))) / Float64(Float64(Float64(Float64(alpha + beta) + Float64(2.0 * i)) * Float64(Float64(alpha + beta) + Float64(2.0 * i))) - 1.0))
end
↓
function code(alpha, beta, i)
t_0 = Float64(Float64(beta + Float64(i * 2.0)) / Float64(i + beta))
t_1 = fma(i, 2.0, Float64(alpha + beta))
t_2 = Float64(t_1 + -1.0)
t_3 = Float64(1.0 + t_1)
tmp = 0.0
if (alpha <= 3.1e+103)
tmp = Float64(Float64(Float64(i / t_3) / t_0) * Float64(Float64(i / t_0) / t_2));
else
tmp = Float64(Float64(Float64(i / Float64(t_1 / Float64(i + Float64(alpha + beta)))) / t_3) * Float64(Float64(alpha + i) / t_2));
end
return tmp
end
code[alpha_, beta_, i_] := N[(N[(N[(N[(i * N[(N[(alpha + beta), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision] * N[(N[(beta * alpha), $MachinePrecision] + N[(i * N[(N[(alpha + beta), $MachinePrecision] + i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision] * N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(N[(N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision] * N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision]
↓
code[alpha_, beta_, i_] := Block[{t$95$0 = N[(N[(beta + N[(i * 2.0), $MachinePrecision]), $MachinePrecision] / N[(i + beta), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$1 = N[(i * 2.0 + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 + -1.0), $MachinePrecision]}, Block[{t$95$3 = N[(1.0 + t$95$1), $MachinePrecision]}, If[LessEqual[alpha, 3.1e+103], N[(N[(N[(i / t$95$3), $MachinePrecision] / t$95$0), $MachinePrecision] * N[(N[(i / t$95$0), $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision], N[(N[(N[(i / N[(t$95$1 / N[(i + N[(alpha + beta), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / t$95$3), $MachinePrecision] * N[(N[(alpha + i), $MachinePrecision] / t$95$2), $MachinePrecision]), $MachinePrecision]]]]]]
\frac{\frac{\left(i \cdot \left(\left(\alpha + \beta\right) + i\right)\right) \cdot \left(\beta \cdot \alpha + i \cdot \left(\left(\alpha + \beta\right) + i\right)\right)}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right)}}{\left(\left(\alpha + \beta\right) + 2 \cdot i\right) \cdot \left(\left(\alpha + \beta\right) + 2 \cdot i\right) - 1}
↓
\begin{array}{l}
t_0 := \frac{\beta + i \cdot 2}{i + \beta}\\
t_1 := \mathsf{fma}\left(i, 2, \alpha + \beta\right)\\
t_2 := t_1 + -1\\
t_3 := 1 + t_1\\
\mathbf{if}\;\alpha \leq 3.1 \cdot 10^{+103}:\\
\;\;\;\;\frac{\frac{i}{t_3}}{t_0} \cdot \frac{\frac{i}{t_0}}{t_2}\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{i}{\frac{t_1}{i + \left(\alpha + \beta\right)}}}{t_3} \cdot \frac{\alpha + i}{t_2}\\
\end{array}