{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \le -6.267281974669497 \cdot 10^{-08}:\\
\;\;\;\;\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \left(\sqrt[3]{\left(\sqrt{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}} + {x}^{\left(\frac{\frac{1}{n}}{2}\right)}\right) \cdot \left(\sqrt{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}} - {x}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)} \cdot \sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}}\right)\\
\mathbf{elif}\;\frac{1}{n} \le 7.385489355116676 \cdot 10^{-07}:\\
\;\;\;\;\left(\frac{1}{x \cdot n} \cdot \frac{\log x}{n} + \frac{1}{x \cdot n}\right) - \frac{\frac{1}{2}}{\left(x \cdot x\right) \cdot n}\\
\mathbf{else}:\\
\;\;\;\;\sqrt[3]{\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}}\right)} \cdot \left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \sqrt[3]{\log \left(\frac{e^{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}}}{e^{{x}^{\left(\frac{1}{n}\right)}}}\right)}\right)\\
\end{array}double f(double x, double n) {
double r1329486 = x;
double r1329487 = 1.0;
double r1329488 = r1329486 + r1329487;
double r1329489 = n;
double r1329490 = r1329487 / r1329489;
double r1329491 = pow(r1329488, r1329490);
double r1329492 = pow(r1329486, r1329490);
double r1329493 = r1329491 - r1329492;
return r1329493;
}
double f(double x, double n) {
double r1329494 = 1.0;
double r1329495 = n;
double r1329496 = r1329494 / r1329495;
double r1329497 = -6.267281974669497e-08;
bool r1329498 = r1329496 <= r1329497;
double r1329499 = x;
double r1329500 = r1329499 + r1329494;
double r1329501 = pow(r1329500, r1329496);
double r1329502 = pow(r1329499, r1329496);
double r1329503 = r1329501 - r1329502;
double r1329504 = cbrt(r1329503);
double r1329505 = sqrt(r1329501);
double r1329506 = 2.0;
double r1329507 = r1329496 / r1329506;
double r1329508 = pow(r1329499, r1329507);
double r1329509 = r1329505 + r1329508;
double r1329510 = r1329505 - r1329508;
double r1329511 = r1329509 * r1329510;
double r1329512 = cbrt(r1329511);
double r1329513 = r1329512 * r1329504;
double r1329514 = r1329504 * r1329513;
double r1329515 = 7.385489355116676e-07;
bool r1329516 = r1329496 <= r1329515;
double r1329517 = r1329499 * r1329495;
double r1329518 = r1329494 / r1329517;
double r1329519 = log(r1329499);
double r1329520 = r1329519 / r1329495;
double r1329521 = r1329518 * r1329520;
double r1329522 = r1329521 + r1329518;
double r1329523 = 0.5;
double r1329524 = r1329499 * r1329499;
double r1329525 = r1329524 * r1329495;
double r1329526 = r1329523 / r1329525;
double r1329527 = r1329522 - r1329526;
double r1329528 = r1329504 * r1329504;
double r1329529 = r1329504 * r1329528;
double r1329530 = cbrt(r1329529);
double r1329531 = exp(r1329501);
double r1329532 = exp(r1329502);
double r1329533 = r1329531 / r1329532;
double r1329534 = log(r1329533);
double r1329535 = cbrt(r1329534);
double r1329536 = r1329504 * r1329535;
double r1329537 = r1329530 * r1329536;
double r1329538 = r1329516 ? r1329527 : r1329537;
double r1329539 = r1329498 ? r1329514 : r1329538;
return r1329539;
}



Bits error versus x



Bits error versus n
Results
if (/ 1 n) < -6.267281974669497e-08Initial program 0.6
rmApplied add-cube-cbrt0.6
rmApplied sqr-pow0.6
Applied add-sqr-sqrt0.6
Applied difference-of-squares0.6
if -6.267281974669497e-08 < (/ 1 n) < 7.385489355116676e-07Initial program 44.9
Taylor expanded around inf 33.1
Simplified33.1
if 7.385489355116676e-07 < (/ 1 n) Initial program 24.9
rmApplied add-cube-cbrt24.9
rmApplied add-cube-cbrt24.9
rmApplied add-log-exp25.0
Applied add-log-exp25.0
Applied diff-log25.0
Final simplification22.4
herbie shell --seed 2019153
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
(- (pow (+ x 1) (/ 1 n)) (pow x (/ 1 n))))