{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\begin{array}{l}
\mathbf{if}\;n \le \frac{-5927677248097903}{268435456} \lor \neg \left(n \le \frac{6854405603415341}{65536}\right):\\
\;\;\;\;\left(\frac{\frac{1}{n}}{x} - \frac{\frac{\frac{1}{2}}{n}}{{x}^{2}}\right) + \frac{\log x \cdot 1}{x \cdot {n}^{2}}\\
\mathbf{else}:\\
\;\;\;\;{\left(\sqrt[3]{x + 1} \cdot \sqrt[3]{x + 1}\right)}^{\left(\frac{1}{n}\right)} \cdot {\left(\sqrt[3]{x + 1}\right)}^{\left(\frac{1}{n}\right)} + \left(-{x}^{\left(\frac{1}{n}\right)}\right)\\
\end{array}double f(double x, double n) {
double r84592 = x;
double r84593 = 1.0;
double r84594 = r84592 + r84593;
double r84595 = n;
double r84596 = r84593 / r84595;
double r84597 = pow(r84594, r84596);
double r84598 = pow(r84592, r84596);
double r84599 = r84597 - r84598;
return r84599;
}
double f(double x, double n) {
double r84600 = n;
double r84601 = -5927677248097903.0;
double r84602 = 268435456.0;
double r84603 = r84601 / r84602;
bool r84604 = r84600 <= r84603;
double r84605 = 6854405603415341.0;
double r84606 = 65536.0;
double r84607 = r84605 / r84606;
bool r84608 = r84600 <= r84607;
double r84609 = !r84608;
bool r84610 = r84604 || r84609;
double r84611 = 1.0;
double r84612 = r84611 / r84600;
double r84613 = x;
double r84614 = r84612 / r84613;
double r84615 = 2.0;
double r84616 = r84611 / r84615;
double r84617 = r84616 / r84600;
double r84618 = 2.0;
double r84619 = pow(r84613, r84618);
double r84620 = r84617 / r84619;
double r84621 = r84614 - r84620;
double r84622 = log(r84613);
double r84623 = r84622 * r84611;
double r84624 = pow(r84600, r84618);
double r84625 = r84613 * r84624;
double r84626 = r84623 / r84625;
double r84627 = r84621 + r84626;
double r84628 = r84613 + r84611;
double r84629 = cbrt(r84628);
double r84630 = r84629 * r84629;
double r84631 = pow(r84630, r84612);
double r84632 = pow(r84629, r84612);
double r84633 = r84631 * r84632;
double r84634 = pow(r84613, r84612);
double r84635 = -r84634;
double r84636 = r84633 + r84635;
double r84637 = r84610 ? r84627 : r84636;
return r84637;
}



Bits error versus x



Bits error versus n
Results
if n < -22082318.544752534Initial program 45.0
Taylor expanded around inf 33.2
Simplified32.6
rmApplied add-log-exp32.7
Simplified32.7
if -22082318.544752534 < n < 104589929251.33272Initial program 8.8
rmApplied add-cube-cbrt8.8
Applied unpow-prod-down8.9
rmApplied add-cbrt-cube8.9
Simplified8.9
if 104589929251.33272 < n Initial program 44.8
Taylor expanded around inf 32.1
Simplified31.5
rmApplied add-log-exp31.6
Final simplification22.1
herbie shell --seed 2019303
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
:precision binary64
(- (pow (+ x 1) (/ 1 n)) (pow x (/ 1 n))))