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 = (0.125 * alpha) / i;
double t_1 = (-0.0625 * (beta + alpha)) / i;
double t_2 = t_1 - (-0.0625 + t_1);
double t_3 = (beta + i) / alpha;
double t_4 = i * (beta + i);
double t_5 = (alpha + beta) + (2.0 * i);
double tmp;
if (alpha <= 9e+83) {
tmp = t_2;
} else if (alpha <= 3.2e+114) {
double tmp_1;
if (i != 0.0) {
tmp_1 = (beta + i) / (1.0 / i);
} else {
tmp_1 = t_4;
}
tmp = tmp_1 / ((t_5 * t_5) - 1.0);
} else if (alpha <= 4.9e+148) {
tmp = t_2;
} else if (alpha <= 1.3e+226) {
tmp = (t_3 * i) / ((alpha + fma(2.0, i, beta)) + -1.0);
} else if (alpha <= 5.2e+243) {
tmp = (0.0625 + t_0) - t_0;
} else if (i != 0.0) {
tmp = t_3 / (alpha / i);
} else {
tmp = t_4 / (alpha * alpha);
}
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(0.125 * alpha) / i)
t_1 = Float64(Float64(-0.0625 * Float64(beta + alpha)) / i)
t_2 = Float64(t_1 - Float64(-0.0625 + t_1))
t_3 = Float64(Float64(beta + i) / alpha)
t_4 = Float64(i * Float64(beta + i))
t_5 = Float64(Float64(alpha + beta) + Float64(2.0 * i))
tmp = 0.0
if (alpha <= 9e+83)
tmp = t_2;
elseif (alpha <= 3.2e+114)
tmp_1 = 0.0
if (i != 0.0)
tmp_1 = Float64(Float64(beta + i) / Float64(1.0 / i));
else
tmp_1 = t_4;
end
tmp = Float64(tmp_1 / Float64(Float64(t_5 * t_5) - 1.0));
elseif (alpha <= 4.9e+148)
tmp = t_2;
elseif (alpha <= 1.3e+226)
tmp = Float64(Float64(t_3 * i) / Float64(Float64(alpha + fma(2.0, i, beta)) + -1.0));
elseif (alpha <= 5.2e+243)
tmp = Float64(Float64(0.0625 + t_0) - t_0);
elseif (i != 0.0)
tmp = Float64(t_3 / Float64(alpha / i));
else
tmp = Float64(t_4 / Float64(alpha * alpha));
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[(0.125 * alpha), $MachinePrecision] / i), $MachinePrecision]}, Block[{t$95$1 = N[(N[(-0.0625 * N[(beta + alpha), $MachinePrecision]), $MachinePrecision] / i), $MachinePrecision]}, Block[{t$95$2 = N[(t$95$1 - N[(-0.0625 + t$95$1), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$3 = N[(N[(beta + i), $MachinePrecision] / alpha), $MachinePrecision]}, Block[{t$95$4 = N[(i * N[(beta + i), $MachinePrecision]), $MachinePrecision]}, Block[{t$95$5 = N[(N[(alpha + beta), $MachinePrecision] + N[(2.0 * i), $MachinePrecision]), $MachinePrecision]}, If[LessEqual[alpha, 9e+83], t$95$2, If[LessEqual[alpha, 3.2e+114], N[(If[Unequal[i, 0.0], N[(N[(beta + i), $MachinePrecision] / N[(1.0 / i), $MachinePrecision]), $MachinePrecision], t$95$4] / N[(N[(t$95$5 * t$95$5), $MachinePrecision] - 1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[alpha, 4.9e+148], t$95$2, If[LessEqual[alpha, 1.3e+226], N[(N[(t$95$3 * i), $MachinePrecision] / N[(N[(alpha + N[(2.0 * i + beta), $MachinePrecision]), $MachinePrecision] + -1.0), $MachinePrecision]), $MachinePrecision], If[LessEqual[alpha, 5.2e+243], N[(N[(0.0625 + t$95$0), $MachinePrecision] - t$95$0), $MachinePrecision], If[Unequal[i, 0.0], N[(t$95$3 / N[(alpha / i), $MachinePrecision]), $MachinePrecision], N[(t$95$4 / N[(alpha * alpha), $MachinePrecision]), $MachinePrecision]]]]]]]]]]]]]