{\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}:\\
\;\;\;\;\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]{\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(\mathsf{fma}\left(1, \sqrt{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}}, -{\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{\frac{1}{n}}{2}\right)} \cdot {\left(\sqrt[3]{x}\right)}^{\left(\frac{\frac{1}{n}}{2}\right)}\right) + \mathsf{fma}\left(-{\left(\sqrt[3]{x}\right)}^{\left(\frac{\frac{1}{n}}{2}\right)}, {\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{\frac{1}{n}}{2}\right)}, {\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{\frac{1}{n}}{2}\right)} \cdot {\left(\sqrt[3]{x}\right)}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)\right) \cdot \left(\sqrt{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)}} + {x}^{\left(\frac{\frac{1}{n}}{2}\right)}\right)}\right)}\right)\\
\mathbf{elif}\;\frac{1}{n} \le 2.865758198043198947866337655820956342723 \cdot 10^{-5}:\\
\;\;\;\;\frac{-0.5}{\left(x \cdot n\right) \cdot x} + \left(\frac{1}{x \cdot n} + \frac{\log x \cdot 1}{n \cdot \left(x \cdot n\right)}\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(\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)}\right) \cdot \sqrt[3]{\log \left(e^{{\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 r2722830 = x;
double r2722831 = 1.0;
double r2722832 = r2722830 + r2722831;
double r2722833 = n;
double r2722834 = r2722831 / r2722833;
double r2722835 = pow(r2722832, r2722834);
double r2722836 = pow(r2722830, r2722834);
double r2722837 = r2722835 - r2722836;
return r2722837;
}
double f(double x, double n) {
double r2722838 = 1.0;
double r2722839 = n;
double r2722840 = r2722838 / r2722839;
double r2722841 = -8.249355031536442e-07;
bool r2722842 = r2722840 <= r2722841;
double r2722843 = x;
double r2722844 = r2722843 + r2722838;
double r2722845 = pow(r2722844, r2722840);
double r2722846 = pow(r2722843, r2722840);
double r2722847 = r2722845 - r2722846;
double r2722848 = cbrt(r2722847);
double r2722849 = 1.0;
double r2722850 = sqrt(r2722845);
double r2722851 = cbrt(r2722843);
double r2722852 = r2722851 * r2722851;
double r2722853 = 2.0;
double r2722854 = r2722840 / r2722853;
double r2722855 = pow(r2722852, r2722854);
double r2722856 = pow(r2722851, r2722854);
double r2722857 = r2722855 * r2722856;
double r2722858 = -r2722857;
double r2722859 = fma(r2722849, r2722850, r2722858);
double r2722860 = -r2722856;
double r2722861 = fma(r2722860, r2722855, r2722857);
double r2722862 = r2722859 + r2722861;
double r2722863 = pow(r2722843, r2722854);
double r2722864 = r2722850 + r2722863;
double r2722865 = r2722862 * r2722864;
double r2722866 = cbrt(r2722865);
double r2722867 = r2722848 * r2722866;
double r2722868 = r2722848 * r2722867;
double r2722869 = cbrt(r2722868);
double r2722870 = r2722848 * r2722869;
double r2722871 = r2722848 * r2722870;
double r2722872 = 2.865758198043199e-05;
bool r2722873 = r2722840 <= r2722872;
double r2722874 = 0.5;
double r2722875 = -r2722874;
double r2722876 = r2722843 * r2722839;
double r2722877 = r2722876 * r2722843;
double r2722878 = r2722875 / r2722877;
double r2722879 = r2722838 / r2722876;
double r2722880 = log(r2722843);
double r2722881 = r2722880 * r2722838;
double r2722882 = r2722839 * r2722876;
double r2722883 = r2722881 / r2722882;
double r2722884 = r2722879 + r2722883;
double r2722885 = r2722878 + r2722884;
double r2722886 = r2722850 - r2722863;
double r2722887 = r2722886 * r2722864;
double r2722888 = cbrt(r2722887);
double r2722889 = r2722848 * r2722888;
double r2722890 = r2722848 * r2722889;
double r2722891 = cbrt(r2722890);
double r2722892 = r2722848 * r2722891;
double r2722893 = exp(r2722847);
double r2722894 = log(r2722893);
double r2722895 = cbrt(r2722894);
double r2722896 = r2722892 * r2722895;
double r2722897 = r2722873 ? r2722885 : r2722896;
double r2722898 = r2722842 ? r2722871 : r2722897;
return r2722898;
}



Bits error versus x



Bits error versus n
if (/ 1.0 n) < -8.249355031536442e-07Initial program 0.8
rmApplied add-cube-cbrt0.8
rmApplied add-cube-cbrt0.8
rmApplied sqr-pow0.8
Applied add-sqr-sqrt0.8
Applied difference-of-squares0.8
rmApplied add-cube-cbrt0.8
Applied unpow-prod-down0.8
Applied *-un-lft-identity0.8
Applied prod-diff0.8
if -8.249355031536442e-07 < (/ 1.0 n) < 2.865758198043199e-05Initial program 44.9
Taylor expanded around inf 32.4
Simplified32.3
if 2.865758198043199e-05 < (/ 1.0 n) Initial program 24.9
rmApplied add-cube-cbrt25.0
rmApplied add-cube-cbrt25.0
rmApplied sqr-pow25.0
Applied add-sqr-sqrt25.0
Applied difference-of-squares25.0
rmApplied add-log-exp25.0
Final simplification22.3
herbie shell --seed 2019171 +o rules:numerics
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
(- (pow (+ x 1.0) (/ 1.0 n)) (pow x (/ 1.0 n))))