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)) <= -7.410239730006016e-26)) {
VAR = ((double) cbrt(((double) pow(((double) cbrt(((double) pow(((double) (((double) pow(((double) (1.0 + x)), ((double) (1.0 / n)))) - ((double) pow(((double) exp(((double) (((double) log(x)) / n)))), 1.0)))), 3.0)))), 3.0))));
} else {
double VAR_1;
if ((((double) (1.0 / n)) <= 1.7465690654858e-17)) {
VAR_1 = ((double) (((double) (1.0 / ((double) (n * x)))) + ((double) (((double) (((double) (1.0 / ((double) (n * n)))) * ((double) (((double) (((double) log(1.0)) / x)) + ((double) (((double) log(x)) / x)))))) - ((double) (0.5 / ((double) (x * ((double) (n * x))))))))));
} else {
VAR_1 = ((double) (((double) cbrt(((double) pow(((double) (((double) pow(((double) (1.0 + x)), ((double) (1.0 / n)))) - ((double) pow(((double) exp(((double) (((double) log(x)) / n)))), 1.0)))), 1.5)))) * ((double) cbrt(((double) pow(((double) (((double) pow(((double) (1.0 + x)), ((double) (1.0 / n)))) - ((double) pow(((double) exp(((double) (((double) log(x)) / n)))), 1.0)))), 1.5))))));
}
VAR = VAR_1;
}
return VAR;
}



Bits error versus x



Bits error versus n
Results
if (/ 1.0 n) < -7.41023973000601582e-26Initial program 6.7
rmApplied add-exp-log6.8
Applied pow-exp6.8
Simplified6.8
rmApplied add-cbrt-cube6.8
Simplified6.8
rmApplied add-cbrt-cube6.8
Simplified6.8
if -7.41023973000601582e-26 < (/ 1.0 n) < 1.7465690654858e-17Initial program 44.7
Taylor expanded around inf 32.3
Simplified32.2
if 1.7465690654858e-17 < (/ 1.0 n) Initial program 9.3
rmApplied add-exp-log9.3
Applied pow-exp9.3
Simplified9.3
rmApplied add-cbrt-cube9.3
Simplified9.3
rmApplied sqr-pow9.3
Applied cbrt-prod9.3
Simplified9.3
Simplified9.3
Final simplification24.3
herbie shell --seed 2020191
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
:precision binary64
(- (pow (+ x 1.0) (/ 1.0 n)) (pow x (/ 1.0 n))))