{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \le -301.0001581399110932579787913709878921509 \lor \neg \left(\frac{1}{n} \le 2.323875494656154788789411439162227722323 \cdot 10^{-16}\right):\\
\;\;\;\;\left(\sqrt[3]{\sqrt[3]{{\left({\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\right)}^{3}}} \cdot \sqrt[3]{\sqrt[3]{{\left({\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\right)}^{3}}}\right) \cdot \sqrt[3]{\log \left(e^{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}}\right)}\\
\mathbf{else}:\\
\;\;\;\;\frac{1}{x} \cdot \left(\frac{1}{n} - \frac{-\log x}{{n}^{2}}\right) - \frac{0.5}{{x}^{2} \cdot n}\\
\end{array}double f(double x, double n) {
double r63465 = x;
double r63466 = 1.0;
double r63467 = r63465 + r63466;
double r63468 = n;
double r63469 = r63466 / r63468;
double r63470 = pow(r63467, r63469);
double r63471 = pow(r63465, r63469);
double r63472 = r63470 - r63471;
return r63472;
}
double f(double x, double n) {
double r63473 = 1.0;
double r63474 = n;
double r63475 = r63473 / r63474;
double r63476 = -301.0001581399111;
bool r63477 = r63475 <= r63476;
double r63478 = 2.323875494656155e-16;
bool r63479 = r63475 <= r63478;
double r63480 = !r63479;
bool r63481 = r63477 || r63480;
double r63482 = x;
double r63483 = r63482 + r63473;
double r63484 = pow(r63483, r63475);
double r63485 = pow(r63482, r63475);
double r63486 = r63484 - r63485;
double r63487 = 3.0;
double r63488 = pow(r63486, r63487);
double r63489 = cbrt(r63488);
double r63490 = cbrt(r63489);
double r63491 = r63490 * r63490;
double r63492 = exp(r63486);
double r63493 = log(r63492);
double r63494 = cbrt(r63493);
double r63495 = r63491 * r63494;
double r63496 = r63473 / r63482;
double r63497 = 1.0;
double r63498 = r63497 / r63474;
double r63499 = log(r63482);
double r63500 = -r63499;
double r63501 = 2.0;
double r63502 = pow(r63474, r63501);
double r63503 = r63500 / r63502;
double r63504 = r63498 - r63503;
double r63505 = r63496 * r63504;
double r63506 = 0.5;
double r63507 = pow(r63482, r63501);
double r63508 = r63507 * r63474;
double r63509 = r63506 / r63508;
double r63510 = r63505 - r63509;
double r63511 = r63481 ? r63495 : r63510;
return r63511;
}



Bits error versus x



Bits error versus n
Results
if (/ 1.0 n) < -301.0001581399111 or 2.323875494656155e-16 < (/ 1.0 n) Initial program 8.6
rmApplied add-cube-cbrt8.6
rmApplied add-log-exp8.6
Applied add-log-exp8.6
Applied diff-log8.6
Simplified8.6
rmApplied add-cbrt-cube8.6
Simplified8.6
rmApplied add-cbrt-cube8.6
Simplified8.6
if -301.0001581399111 < (/ 1.0 n) < 2.323875494656155e-16Initial program 45.1
Taylor expanded around inf 32.8
Simplified32.2
Final simplification22.4
herbie shell --seed 2019306
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
:precision binary64
(- (pow (+ x 1) (/ 1 n)) (pow x (/ 1 n))))