{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \le -0.096355036664333621:\\
\;\;\;\;\sqrt[3]{{\left({\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - e^{\log x \cdot \frac{1}{n}}\right)}^{3}}\\
\mathbf{elif}\;\frac{1}{n} \le 6.8280679940979408 \cdot 10^{-9}:\\
\;\;\;\;\left(\frac{\frac{1}{n}}{x} + 1 \cdot \left(0 - \frac{\log \left(\frac{1}{x}\right)}{x \cdot {n}^{2}}\right)\right) - \frac{\frac{0.5}{n}}{{x}^{2}}\\
\mathbf{else}:\\
\;\;\;\;\sqrt[3]{e^{\left(\sqrt[3]{\log \left({\left({\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - e^{\log x \cdot \frac{1}{n}}\right)}^{3}\right)} \cdot \sqrt[3]{\log \left({\left({\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - e^{\log x \cdot \frac{1}{n}}\right)}^{3}\right)}\right) \cdot \sqrt[3]{\log \left({\left({\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - e^{\log x \cdot \frac{1}{n}}\right)}^{3}\right)}}}\\
\end{array}double code(double x, double n) {
return ((double) (((double) pow(((double) (x + 1.0)), ((double) (1.0 / n)))) - ((double) pow(x, ((double) (1.0 / n))))));
}
double code(double x, double n) {
double VAR;
if ((((double) (1.0 / n)) <= -0.09635503666433362)) {
VAR = ((double) cbrt(((double) pow(((double) (((double) pow(((double) (x + 1.0)), ((double) (1.0 / n)))) - ((double) exp(((double) (((double) log(x)) * ((double) (1.0 / n)))))))), 3.0))));
} else {
double VAR_1;
if ((((double) (1.0 / n)) <= 6.828067994097941e-09)) {
VAR_1 = ((double) (((double) (((double) (((double) (1.0 / n)) / x)) + ((double) (1.0 * ((double) (0.0 - ((double) (((double) log(((double) (1.0 / x)))) / ((double) (x * ((double) pow(n, 2.0)))))))))))) - ((double) (((double) (0.5 / n)) / ((double) pow(x, 2.0))))));
} else {
VAR_1 = ((double) cbrt(((double) exp(((double) (((double) (((double) cbrt(((double) log(((double) pow(((double) (((double) pow(((double) (x + 1.0)), ((double) (1.0 / n)))) - ((double) exp(((double) (((double) log(x)) * ((double) (1.0 / n)))))))), 3.0)))))) * ((double) cbrt(((double) log(((double) pow(((double) (((double) pow(((double) (x + 1.0)), ((double) (1.0 / n)))) - ((double) exp(((double) (((double) log(x)) * ((double) (1.0 / n)))))))), 3.0)))))))) * ((double) cbrt(((double) log(((double) pow(((double) (((double) pow(((double) (x + 1.0)), ((double) (1.0 / n)))) - ((double) exp(((double) (((double) log(x)) * ((double) (1.0 / n)))))))), 3.0))))))))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus n
Results
if (/ 1.0 n) < -0.096355036664333621Initial program 0.2
rmApplied add-exp-log0.2
Applied pow-exp0.2
rmApplied add-cbrt-cube0.3
Simplified0.3
if -0.096355036664333621 < (/ 1.0 n) < 6.8280679940979408e-9Initial program 45.5
Taylor expanded around inf 33.2
Simplified32.7
if 6.8280679940979408e-9 < (/ 1.0 n) Initial program 6.7
rmApplied add-exp-log6.7
Applied pow-exp6.7
rmApplied add-cbrt-cube6.7
Simplified6.7
rmApplied add-exp-log6.7
Applied pow-exp6.7
Simplified6.7
rmApplied add-cube-cbrt6.7
Final simplification24.1
herbie shell --seed 2020175
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
:precision binary64
(- (pow (+ x 1.0) (/ 1.0 n)) (pow x (/ 1.0 n))))