\[\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 := \mathsf{fma}\left(i, 2, \beta + \alpha\right)\\
\mathbf{if}\;\beta \leq 4.8 \cdot 10^{+98}:\\
\;\;\;\;\left(\frac{i}{t_0} \cdot \frac{\beta + i}{\beta + i \cdot 2}\right) \cdot 0.25\\
\mathbf{elif}\;\beta \leq 1.55 \cdot 10^{+192} \lor \neg \left(\beta \leq 9 \cdot 10^{+201}\right):\\
\;\;\;\;\frac{\frac{i}{\frac{t_0 + -1}{i + \alpha}}}{t_0 + 1}\\
\mathbf{else}:\\
\;\;\;\;\left(0.0625 + 0.125 \cdot \frac{\beta}{i}\right) + \frac{\beta}{i} \cdot -0.125\\
\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 = fma(i, 2.0, (beta + alpha));
double tmp;
if (beta <= 4.8e+98) {
tmp = ((i / t_0) * ((beta + i) / (beta + (i * 2.0)))) * 0.25;
} else if ((beta <= 1.55e+192) || !(beta <= 9e+201)) {
tmp = (i / ((t_0 + -1.0) / (i + alpha))) / (t_0 + 1.0);
} else {
tmp = (0.0625 + (0.125 * (beta / i))) + ((beta / i) * -0.125);
}
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 = fma(i, 2.0, Float64(beta + alpha))
tmp = 0.0
if (beta <= 4.8e+98)
tmp = Float64(Float64(Float64(i / t_0) * Float64(Float64(beta + i) / Float64(beta + Float64(i * 2.0)))) * 0.25);
elseif ((beta <= 1.55e+192) || !(beta <= 9e+201))
tmp = Float64(Float64(i / Float64(Float64(t_0 + -1.0) / Float64(i + alpha))) / Float64(t_0 + 1.0));
else
tmp = Float64(Float64(0.0625 + Float64(0.125 * Float64(beta / i))) + Float64(Float64(beta / i) * -0.125));
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[(i * 2.0 + N[(beta + alpha), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[beta, 4.8e+98], N[(N[(N[(i / t$95$0), $MachinePrecision] * N[(N[(beta + i), $MachinePrecision] / N[(beta + N[(i * 2.0), $MachinePrecision]), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * 0.25), $MachinePrecision], If[Or[LessEqual[beta, 1.55e+192], N[Not[LessEqual[beta, 9e+201]], $MachinePrecision]], N[(N[(i / N[(N[(t$95$0 + -1.0), $MachinePrecision] / N[(i + alpha), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[(t$95$0 + 1.0), $MachinePrecision]), $MachinePrecision], N[(N[(0.0625 + N[(0.125 * N[(beta / i), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(beta / i), $MachinePrecision] * -0.125), $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 := \mathsf{fma}\left(i, 2, \beta + \alpha\right)\\
\mathbf{if}\;\beta \leq 4.8 \cdot 10^{+98}:\\
\;\;\;\;\left(\frac{i}{t_0} \cdot \frac{\beta + i}{\beta + i \cdot 2}\right) \cdot 0.25\\
\mathbf{elif}\;\beta \leq 1.55 \cdot 10^{+192} \lor \neg \left(\beta \leq 9 \cdot 10^{+201}\right):\\
\;\;\;\;\frac{\frac{i}{\frac{t_0 + -1}{i + \alpha}}}{t_0 + 1}\\
\mathbf{else}:\\
\;\;\;\;\left(0.0625 + 0.125 \cdot \frac{\beta}{i}\right) + \frac{\beta}{i} \cdot -0.125\\
\end{array}