{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \le -4.895039802440850170454212278486686446129 \cdot 10^{-17}:\\
\;\;\;\;\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)} \cdot \left(\sqrt[3]{{\left(\sqrt[3]{\sqrt[3]{x}} \cdot \left(\sqrt[3]{\sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}\right)\right)}^{\left(\frac{1}{n}\right)}} \cdot \left(\sqrt[3]{{\left(\sqrt[3]{\sqrt[3]{x}} \cdot \left(\sqrt[3]{\sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}\right)\right)}^{\left(\frac{1}{n}\right)}} \cdot \sqrt[3]{{\left(\sqrt[3]{\sqrt[3]{x}} \cdot \left(\sqrt[3]{\sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}\right)\right)}^{\left(\frac{1}{n}\right)}}\right)\right)} \cdot \left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)} \cdot \left(\sqrt[3]{{\left(\sqrt[3]{\sqrt[3]{x}} \cdot \left(\sqrt[3]{\sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}\right)\right)}^{\left(\frac{1}{n}\right)}} \cdot \left(\sqrt[3]{{\left(\sqrt[3]{\sqrt[3]{x}} \cdot \left(\sqrt[3]{\sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}\right)\right)}^{\left(\frac{1}{n}\right)}} \cdot \sqrt[3]{{\left(\sqrt[3]{\sqrt[3]{x}} \cdot \left(\sqrt[3]{\sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}\right)\right)}^{\left(\frac{1}{n}\right)}}\right)\right)} \cdot \sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)} \cdot \left(\sqrt[3]{{\left(\sqrt[3]{\sqrt[3]{x}} \cdot \left(\sqrt[3]{\sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}\right)\right)}^{\left(\frac{1}{n}\right)}} \cdot \left(\sqrt[3]{{\left(\sqrt[3]{\sqrt[3]{x}} \cdot \left(\sqrt[3]{\sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}\right)\right)}^{\left(\frac{1}{n}\right)}} \cdot \sqrt[3]{{\left(\sqrt[3]{\sqrt[3]{x}} \cdot \left(\sqrt[3]{\sqrt[3]{x}} \cdot \sqrt[3]{\sqrt[3]{x}}\right)\right)}^{\left(\frac{1}{n}\right)}}\right)\right)}\right)\\
\mathbf{elif}\;\frac{1}{n} \le 1.260859621565214258622089794815744532231 \cdot 10^{-44}:\\
\;\;\;\;\left(\frac{\log x \cdot 1}{x \cdot \left(n \cdot n\right)} + \frac{\frac{1}{x}}{n}\right) - \frac{\frac{0.5}{n}}{x \cdot x}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \sqrt{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}}\\
\end{array}double f(double x, double n) {
double r4273959 = x;
double r4273960 = 1.0;
double r4273961 = r4273959 + r4273960;
double r4273962 = n;
double r4273963 = r4273960 / r4273962;
double r4273964 = pow(r4273961, r4273963);
double r4273965 = pow(r4273959, r4273963);
double r4273966 = r4273964 - r4273965;
return r4273966;
}
double f(double x, double n) {
double r4273967 = 1.0;
double r4273968 = n;
double r4273969 = r4273967 / r4273968;
double r4273970 = -4.89503980244085e-17;
bool r4273971 = r4273969 <= r4273970;
double r4273972 = x;
double r4273973 = r4273972 + r4273967;
double r4273974 = pow(r4273973, r4273969);
double r4273975 = cbrt(r4273972);
double r4273976 = r4273975 * r4273975;
double r4273977 = pow(r4273976, r4273969);
double r4273978 = cbrt(r4273975);
double r4273979 = r4273978 * r4273978;
double r4273980 = r4273978 * r4273979;
double r4273981 = pow(r4273980, r4273969);
double r4273982 = cbrt(r4273981);
double r4273983 = r4273982 * r4273982;
double r4273984 = r4273982 * r4273983;
double r4273985 = r4273977 * r4273984;
double r4273986 = r4273974 - r4273985;
double r4273987 = cbrt(r4273986);
double r4273988 = r4273987 * r4273987;
double r4273989 = r4273987 * r4273988;
double r4273990 = 1.2608596215652143e-44;
bool r4273991 = r4273969 <= r4273990;
double r4273992 = log(r4273972);
double r4273993 = r4273992 * r4273967;
double r4273994 = r4273968 * r4273968;
double r4273995 = r4273972 * r4273994;
double r4273996 = r4273993 / r4273995;
double r4273997 = r4273967 / r4273972;
double r4273998 = r4273997 / r4273968;
double r4273999 = r4273996 + r4273998;
double r4274000 = 0.5;
double r4274001 = r4274000 / r4273968;
double r4274002 = r4273972 * r4273972;
double r4274003 = r4274001 / r4274002;
double r4274004 = r4273999 - r4274003;
double r4274005 = pow(r4273972, r4273969);
double r4274006 = r4273974 - r4274005;
double r4274007 = sqrt(r4274006);
double r4274008 = r4274007 * r4274007;
double r4274009 = r4273991 ? r4274004 : r4274008;
double r4274010 = r4273971 ? r4273989 : r4274009;
return r4274010;
}



Bits error versus x



Bits error versus n
Results
if (/ 1.0 n) < -4.89503980244085e-17Initial program 2.2
rmApplied add-cube-cbrt2.2
Applied unpow-prod-down2.2
rmApplied add-cube-cbrt2.3
rmApplied add-cube-cbrt2.3
rmApplied add-cube-cbrt2.3
if -4.89503980244085e-17 < (/ 1.0 n) < 1.2608596215652143e-44Initial program 44.5
Taylor expanded around inf 32.2
Simplified31.6
if 1.2608596215652143e-44 < (/ 1.0 n) Initial program 30.5
rmApplied add-sqr-sqrt30.5
Final simplification22.8
herbie shell --seed 2019172
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
(- (pow (+ x 1.0) (/ 1.0 n)) (pow x (/ 1.0 n))))