\[{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}
\]
↓
\[\begin{array}{l}
\mathbf{if}\;x \leq 5.7:\\
\;\;\;\;\mathsf{fma}\left(0.5, \frac{{\left(\mathsf{log1p}\left(x\right)\right)}^{2}}{n \cdot n}, \frac{0.16666666666666666 \cdot \left({\left(\mathsf{log1p}\left(x\right)\right)}^{3} - {\log x}^{3}\right)}{{n}^{3}}\right) + \left(\frac{{\log x}^{2}}{n \cdot n} \cdot -0.5 - \frac{\mathsf{fma}\left(-1, \mathsf{log1p}\left(x\right), \log x\right)}{n}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot \frac{{\left(e^{-0.5}\right)}^{\left(-\frac{\log x}{n}\right)}}{x \cdot n}\right) \cdot \left({\left(x + 1\right)}^{\left(\frac{0.5}{n}\right)} + {x}^{\left(\frac{0.5}{n}\right)}\right)\\
\end{array}
\]
double code(double x, double n) {
return pow((x + 1.0), (1.0 / n)) - pow(x, (1.0 / n));
}
↓
double code(double x, double n) {
double tmp;
if (x <= 5.7) {
tmp = fma(0.5, (pow(log1p(x), 2.0) / (n * n)), ((0.16666666666666666 * (pow(log1p(x), 3.0) - pow(log(x), 3.0))) / pow(n, 3.0))) + (((pow(log(x), 2.0) / (n * n)) * -0.5) - (fma(-1.0, log1p(x), log(x)) / n));
} else {
tmp = (0.5 * (pow(exp(-0.5), -(log(x) / n)) / (x * n))) * (pow((x + 1.0), (0.5 / n)) + pow(x, (0.5 / n)));
}
return tmp;
}
function code(x, n)
return Float64((Float64(x + 1.0) ^ Float64(1.0 / n)) - (x ^ Float64(1.0 / n)))
end
↓
function code(x, n)
tmp = 0.0
if (x <= 5.7)
tmp = Float64(fma(0.5, Float64((log1p(x) ^ 2.0) / Float64(n * n)), Float64(Float64(0.16666666666666666 * Float64((log1p(x) ^ 3.0) - (log(x) ^ 3.0))) / (n ^ 3.0))) + Float64(Float64(Float64((log(x) ^ 2.0) / Float64(n * n)) * -0.5) - Float64(fma(-1.0, log1p(x), log(x)) / n)));
else
tmp = Float64(Float64(0.5 * Float64((exp(-0.5) ^ Float64(-Float64(log(x) / n))) / Float64(x * n))) * Float64((Float64(x + 1.0) ^ Float64(0.5 / n)) + (x ^ Float64(0.5 / n))));
end
return tmp
end
code[x_, n_] := N[(N[Power[N[(x + 1.0), $MachinePrecision], N[(1.0 / n), $MachinePrecision]], $MachinePrecision] - N[Power[x, N[(1.0 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]
↓
code[x_, n_] := If[LessEqual[x, 5.7], N[(N[(0.5 * N[(N[Power[N[Log[1 + x], $MachinePrecision], 2.0], $MachinePrecision] / N[(n * n), $MachinePrecision]), $MachinePrecision] + N[(N[(0.16666666666666666 * N[(N[Power[N[Log[1 + x], $MachinePrecision], 3.0], $MachinePrecision] - N[Power[N[Log[x], $MachinePrecision], 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] / N[Power[n, 3.0], $MachinePrecision]), $MachinePrecision]), $MachinePrecision] + N[(N[(N[(N[Power[N[Log[x], $MachinePrecision], 2.0], $MachinePrecision] / N[(n * n), $MachinePrecision]), $MachinePrecision] * -0.5), $MachinePrecision] - N[(N[(-1.0 * N[Log[1 + x], $MachinePrecision] + N[Log[x], $MachinePrecision]), $MachinePrecision] / n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision], N[(N[(0.5 * N[(N[Power[N[Exp[-0.5], $MachinePrecision], (-N[(N[Log[x], $MachinePrecision] / n), $MachinePrecision])], $MachinePrecision] / N[(x * n), $MachinePrecision]), $MachinePrecision]), $MachinePrecision] * N[(N[Power[N[(x + 1.0), $MachinePrecision], N[(0.5 / n), $MachinePrecision]], $MachinePrecision] + N[Power[x, N[(0.5 / n), $MachinePrecision]], $MachinePrecision]), $MachinePrecision]), $MachinePrecision]]
{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}
↓
\begin{array}{l}
\mathbf{if}\;x \leq 5.7:\\
\;\;\;\;\mathsf{fma}\left(0.5, \frac{{\left(\mathsf{log1p}\left(x\right)\right)}^{2}}{n \cdot n}, \frac{0.16666666666666666 \cdot \left({\left(\mathsf{log1p}\left(x\right)\right)}^{3} - {\log x}^{3}\right)}{{n}^{3}}\right) + \left(\frac{{\log x}^{2}}{n \cdot n} \cdot -0.5 - \frac{\mathsf{fma}\left(-1, \mathsf{log1p}\left(x\right), \log x\right)}{n}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(0.5 \cdot \frac{{\left(e^{-0.5}\right)}^{\left(-\frac{\log x}{n}\right)}}{x \cdot n}\right) \cdot \left({\left(x + 1\right)}^{\left(\frac{0.5}{n}\right)} + {x}^{\left(\frac{0.5}{n}\right)}\right)\\
\end{array}