{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \le -64468894.28699002:\\
\;\;\;\;\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \left(\log \left(\sqrt[3]{e^{\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)}}}} \cdot \sqrt[3]{e^{\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) + \log \left(\sqrt[3]{e^{\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)\\
\mathbf{elif}\;\frac{1}{n} \le 3.000508809945158 \cdot 10^{-29}:\\
\;\;\;\;\left(\frac{\frac{1}{n}}{x} - \frac{-\log x}{\left(n \cdot n\right) \cdot x}\right) - \frac{\frac{1}{2}}{x \cdot \left(x \cdot n\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \left(\log \left(\sqrt[3]{e^{\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)}}}} \cdot \sqrt[3]{e^{\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) + \log \left(\sqrt[3]{e^{\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)\\
\end{array}double f(double x, double n) {
double r860564 = x;
double r860565 = 1.0;
double r860566 = r860564 + r860565;
double r860567 = n;
double r860568 = r860565 / r860567;
double r860569 = pow(r860566, r860568);
double r860570 = pow(r860564, r860568);
double r860571 = r860569 - r860570;
return r860571;
}
double f(double x, double n) {
double r860572 = 1.0;
double r860573 = n;
double r860574 = r860572 / r860573;
double r860575 = -64468894.28699002;
bool r860576 = r860574 <= r860575;
double r860577 = x;
double r860578 = r860577 + r860572;
double r860579 = pow(r860578, r860574);
double r860580 = pow(r860577, r860574);
double r860581 = r860579 - r860580;
double r860582 = cbrt(r860581);
double r860583 = r860582 * r860582;
double r860584 = exp(r860583);
double r860585 = cbrt(r860584);
double r860586 = r860585 * r860585;
double r860587 = log(r860586);
double r860588 = log(r860585);
double r860589 = r860587 + r860588;
double r860590 = r860582 * r860589;
double r860591 = 3.000508809945158e-29;
bool r860592 = r860574 <= r860591;
double r860593 = r860574 / r860577;
double r860594 = log(r860577);
double r860595 = -r860594;
double r860596 = r860573 * r860573;
double r860597 = r860596 * r860577;
double r860598 = r860595 / r860597;
double r860599 = r860593 - r860598;
double r860600 = 0.5;
double r860601 = r860577 * r860573;
double r860602 = r860577 * r860601;
double r860603 = r860600 / r860602;
double r860604 = r860599 - r860603;
double r860605 = r860592 ? r860604 : r860590;
double r860606 = r860576 ? r860590 : r860605;
return r860606;
}



Bits error versus x



Bits error versus n
Results
if (/ 1 n) < -64468894.28699002 or 3.000508809945158e-29 < (/ 1 n) Initial program 10.5
rmApplied add-log-exp10.5
rmApplied add-cube-cbrt10.5
Applied exp-prod10.5
Applied log-pow10.5
rmApplied add-cube-cbrt10.8
Applied log-prod10.8
if -64468894.28699002 < (/ 1 n) < 3.000508809945158e-29Initial program 43.5
Taylor expanded around inf 31.9
Simplified31.3
Final simplification22.5
herbie shell --seed 2019128
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
(- (pow (+ x 1) (/ 1 n)) (pow x (/ 1 n))))