{\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(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \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)}\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]{\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 \sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}}} \cdot \left(\sqrt[3]{\log \left(\frac{e^{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}}}{e^{{x}^{\left(\frac{1}{n}\right)}}}\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 r1300867 = x;
double r1300868 = 1.0;
double r1300869 = r1300867 + r1300868;
double r1300870 = n;
double r1300871 = r1300868 / r1300870;
double r1300872 = pow(r1300869, r1300871);
double r1300873 = pow(r1300867, r1300871);
double r1300874 = r1300872 - r1300873;
return r1300874;
}
double f(double x, double n) {
double r1300875 = 1.0;
double r1300876 = n;
double r1300877 = r1300875 / r1300876;
double r1300878 = -6.267281974669497e-08;
bool r1300879 = r1300877 <= r1300878;
double r1300880 = x;
double r1300881 = r1300880 + r1300875;
double r1300882 = pow(r1300881, r1300877);
double r1300883 = pow(r1300880, r1300877);
double r1300884 = r1300882 - r1300883;
double r1300885 = cbrt(r1300884);
double r1300886 = sqrt(r1300882);
double r1300887 = 2.0;
double r1300888 = r1300877 / r1300887;
double r1300889 = pow(r1300880, r1300888);
double r1300890 = r1300886 - r1300889;
double r1300891 = r1300886 + r1300889;
double r1300892 = r1300890 * r1300891;
double r1300893 = cbrt(r1300892);
double r1300894 = r1300885 * r1300893;
double r1300895 = r1300885 * r1300894;
double r1300896 = 7.385489355116676e-07;
bool r1300897 = r1300877 <= r1300896;
double r1300898 = r1300880 * r1300876;
double r1300899 = r1300875 / r1300898;
double r1300900 = log(r1300880);
double r1300901 = r1300900 / r1300876;
double r1300902 = r1300899 * r1300901;
double r1300903 = r1300902 + r1300899;
double r1300904 = 0.5;
double r1300905 = r1300880 * r1300880;
double r1300906 = r1300905 * r1300876;
double r1300907 = r1300904 / r1300906;
double r1300908 = r1300903 - r1300907;
double r1300909 = r1300885 * r1300885;
double r1300910 = r1300909 * r1300885;
double r1300911 = cbrt(r1300910);
double r1300912 = exp(r1300882);
double r1300913 = exp(r1300883);
double r1300914 = r1300912 / r1300913;
double r1300915 = log(r1300914);
double r1300916 = cbrt(r1300915);
double r1300917 = r1300916 * r1300885;
double r1300918 = r1300911 * r1300917;
double r1300919 = r1300897 ? r1300908 : r1300918;
double r1300920 = r1300879 ? r1300895 : r1300919;
return r1300920;
}



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))))