{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \le -8.28957440978626294 \cdot 10^{-4}:\\
\;\;\;\;\left({\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)} + {x}^{\left(\frac{\frac{1}{n}}{2}\right)}\right) \cdot \left({\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)} - {x}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)\\
\mathbf{elif}\;\frac{1}{n} \le 2.7376649343181104 \cdot 10^{-11}:\\
\;\;\;\;\frac{\mathsf{fma}\left(\mathsf{fma}\left({x}^{\left(\frac{\frac{1}{n}}{2}\right)}, {x}^{\left(\frac{\frac{1}{n}}{2}\right)} - {\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)}, {\left(x + 1\right)}^{\left(2 \cdot \frac{\frac{1}{n}}{2}\right)}\right), \left(-\left({\left({\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)}^{3} + {\left({x}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)}^{3}\right)\right) \cdot \mathsf{fma}\left(0.25, \frac{\log \left(\frac{1}{x}\right)}{x \cdot {n}^{2}}, \frac{\frac{0.25}{n}}{{x}^{2}}\right), \left(\frac{\frac{0.5}{n}}{x} \cdot \left({\left({\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)}^{3} + {\left({x}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)}^{3}\right)\right) \cdot \left({x}^{\left(\frac{\frac{1}{n}}{2}\right)} \cdot \left({x}^{\left(\frac{\frac{1}{n}}{2}\right)} - {\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)}\right) + {\left(x + 1\right)}^{\left(2 \cdot \frac{\frac{1}{n}}{2}\right)}\right)\right)}{\left({\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)} \cdot {\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)} + \left({x}^{\left(\frac{\frac{1}{n}}{2}\right)} \cdot {x}^{\left(\frac{\frac{1}{n}}{2}\right)} - {\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)} \cdot {x}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)\right) \cdot \left({\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)} \cdot {\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)} + \left({x}^{\left(\frac{\frac{1}{n}}{2}\right)} \cdot {x}^{\left(\frac{\frac{1}{n}}{2}\right)} - {\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)} \cdot {x}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{{\left({\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)}^{3} + {\left({x}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)}^{3}}{\mathsf{fma}\left({x}^{\left(\frac{\frac{1}{n}}{2}\right)}, {x}^{\left(\frac{\frac{1}{n}}{2}\right)} - {\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)}, {\left(x + 1\right)}^{\left(2 \cdot \frac{\frac{1}{n}}{2}\right)}\right)} \cdot \left({\left(x + 1\right)}^{\left(\frac{\frac{1}{n}}{2}\right)} - {x}^{\left(\frac{\frac{1}{n}}{2}\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 VAR;
if (((1.0 / n) <= -0.0008289574409786263)) {
VAR = ((pow((x + 1.0), ((1.0 / n) / 2.0)) + pow(x, ((1.0 / n) / 2.0))) * (pow((x + 1.0), ((1.0 / n) / 2.0)) - pow(x, ((1.0 / n) / 2.0))));
} else {
double VAR_1;
if (((1.0 / n) <= 2.7376649343181104e-11)) {
VAR_1 = (fma(fma(pow(x, ((1.0 / n) / 2.0)), (pow(x, ((1.0 / n) / 2.0)) - pow((x + 1.0), ((1.0 / n) / 2.0))), pow((x + 1.0), (2.0 * ((1.0 / n) / 2.0)))), (-(pow(pow((x + 1.0), ((1.0 / n) / 2.0)), 3.0) + pow(pow(x, ((1.0 / n) / 2.0)), 3.0)) * fma(0.25, (log((1.0 / x)) / (x * pow(n, 2.0))), ((0.25 / n) / pow(x, 2.0)))), ((((0.5 / n) / x) * (pow(pow((x + 1.0), ((1.0 / n) / 2.0)), 3.0) + pow(pow(x, ((1.0 / n) / 2.0)), 3.0))) * ((pow(x, ((1.0 / n) / 2.0)) * (pow(x, ((1.0 / n) / 2.0)) - pow((x + 1.0), ((1.0 / n) / 2.0)))) + pow((x + 1.0), (2.0 * ((1.0 / n) / 2.0)))))) / (((pow((x + 1.0), ((1.0 / n) / 2.0)) * pow((x + 1.0), ((1.0 / n) / 2.0))) + ((pow(x, ((1.0 / n) / 2.0)) * pow(x, ((1.0 / n) / 2.0))) - (pow((x + 1.0), ((1.0 / n) / 2.0)) * pow(x, ((1.0 / n) / 2.0))))) * ((pow((x + 1.0), ((1.0 / n) / 2.0)) * pow((x + 1.0), ((1.0 / n) / 2.0))) + ((pow(x, ((1.0 / n) / 2.0)) * pow(x, ((1.0 / n) / 2.0))) - (pow((x + 1.0), ((1.0 / n) / 2.0)) * pow(x, ((1.0 / n) / 2.0)))))));
} else {
VAR_1 = (((pow(pow((x + 1.0), ((1.0 / n) / 2.0)), 3.0) + pow(pow(x, ((1.0 / n) / 2.0)), 3.0)) / fma(pow(x, ((1.0 / n) / 2.0)), (pow(x, ((1.0 / n) / 2.0)) - pow((x + 1.0), ((1.0 / n) / 2.0))), pow((x + 1.0), (2.0 * ((1.0 / n) / 2.0))))) * (pow((x + 1.0), ((1.0 / n) / 2.0)) - pow(x, ((1.0 / n) / 2.0))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus n
Results
if (/ 1.0 n) < -0.0008289574409786263Initial program 0.5
rmApplied sqr-pow0.5
Applied sqr-pow0.5
Applied difference-of-squares0.5
if -0.0008289574409786263 < (/ 1.0 n) < 2.7376649343181104e-11Initial program 44.6
rmApplied sqr-pow44.7
Applied sqr-pow44.6
Applied difference-of-squares44.6
Taylor expanded around inf 32.6
Simplified32.1
rmApplied fma-udef32.1
Applied distribute-lft-in32.1
Simplified32.0
Simplified32.0
rmApplied flip3-+32.0
Applied associate-*r/32.0
Applied flip3-+32.0
Applied distribute-neg-frac32.0
Applied associate-*l/32.0
Applied frac-add32.0
Simplified32.0
if 2.7376649343181104e-11 < (/ 1.0 n) Initial program 6.9
rmApplied sqr-pow7.0
Applied sqr-pow6.9
Applied difference-of-squares6.9
rmApplied flip3-+7.0
Simplified7.0
Final simplification23.3
herbie shell --seed 2020078 +o rules:numerics
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
:precision binary64
(- (pow (+ x 1) (/ 1 n)) (pow x (/ 1 n))))