{\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.249355031536442258007987837764574123867 \cdot 10^{-7}:\\
\;\;\;\;\left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \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)}\right) \cdot \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)}\\
\mathbf{elif}\;\frac{1}{n} \le 2.865758198043198947866337655820956342723 \cdot 10^{-5}:\\
\;\;\;\;\frac{1 \cdot \log x}{\left(n \cdot n\right) \cdot x} + \left(\frac{1}{x \cdot n} - \frac{\frac{0.5}{n}}{x \cdot x}\right)\\
\mathbf{else}:\\
\;\;\;\;\left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \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)}\right) \cdot \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)}\\
\end{array}double f(double x, double n) {
double r2381603 = x;
double r2381604 = 1.0;
double r2381605 = r2381603 + r2381604;
double r2381606 = n;
double r2381607 = r2381604 / r2381606;
double r2381608 = pow(r2381605, r2381607);
double r2381609 = pow(r2381603, r2381607);
double r2381610 = r2381608 - r2381609;
return r2381610;
}
double f(double x, double n) {
double r2381611 = 1.0;
double r2381612 = n;
double r2381613 = r2381611 / r2381612;
double r2381614 = -8.249355031536442e-07;
bool r2381615 = r2381613 <= r2381614;
double r2381616 = x;
double r2381617 = r2381616 + r2381611;
double r2381618 = pow(r2381617, r2381613);
double r2381619 = pow(r2381616, r2381613);
double r2381620 = r2381618 - r2381619;
double r2381621 = cbrt(r2381620);
double r2381622 = r2381621 * r2381621;
double r2381623 = r2381621 * r2381622;
double r2381624 = cbrt(r2381623);
double r2381625 = r2381621 * r2381624;
double r2381626 = r2381625 * r2381624;
double r2381627 = 2.865758198043199e-05;
bool r2381628 = r2381613 <= r2381627;
double r2381629 = log(r2381616);
double r2381630 = r2381611 * r2381629;
double r2381631 = r2381612 * r2381612;
double r2381632 = r2381631 * r2381616;
double r2381633 = r2381630 / r2381632;
double r2381634 = r2381616 * r2381612;
double r2381635 = r2381611 / r2381634;
double r2381636 = 0.5;
double r2381637 = r2381636 / r2381612;
double r2381638 = r2381616 * r2381616;
double r2381639 = r2381637 / r2381638;
double r2381640 = r2381635 - r2381639;
double r2381641 = r2381633 + r2381640;
double r2381642 = r2381628 ? r2381641 : r2381626;
double r2381643 = r2381615 ? r2381626 : r2381642;
return r2381643;
}



Bits error versus x



Bits error versus n
Results
if (/ 1.0 n) < -8.249355031536442e-07 or 2.865758198043199e-05 < (/ 1.0 n) Initial program 8.4
rmApplied add-cube-cbrt8.4
rmApplied add-cube-cbrt8.4
rmApplied add-cube-cbrt8.4
if -8.249355031536442e-07 < (/ 1.0 n) < 2.865758198043199e-05Initial program 44.9
Taylor expanded around inf 32.4
Simplified32.4
Final simplification22.3
herbie shell --seed 2019171
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
(- (pow (+ x 1.0) (/ 1.0 n)) (pow x (/ 1.0 n))))