{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\begin{array}{l}
\mathbf{if}\;\frac{1}{n} \le -3.2155947364724913 \cdot 10^{-22}:\\
\;\;\;\;\mathsf{fma}\left(-{\left(\sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)}, {\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)}, {\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)} \cdot {\left(\sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)}\right) + \sqrt[3]{\left(\left(\mathsf{fma}\left(-{\left(\sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)}, {\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)}, {\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)} \cdot {\left(\sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)}\right) + \mathsf{fma}\left(1, {\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(-{\left(\sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)}\right)\right)\right) \cdot \left({\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]{x}\right)}^{\left(\frac{1}{n}\right)}\right)\right) \cdot \left(\mathsf{fma}\left(-{\left(\sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)}, {\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)}, {\left(\sqrt[3]{x} \cdot \sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)} \cdot {\left(\sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)}\right) + \mathsf{fma}\left(1, {\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(-{\left(\sqrt[3]{x}\right)}^{\left(\frac{1}{n}\right)}\right)\right)\right)}\\
\mathbf{elif}\;\frac{1}{n} \le 4.2228249751738446 \cdot 10^{-17}:\\
\;\;\;\;\left(\frac{\log x}{x \cdot \left(n \cdot n\right)} + \frac{1}{x \cdot n}\right) - \frac{\frac{\frac{1}{2}}{n}}{x \cdot x}\\
\mathbf{else}:\\
\;\;\;\;e^{\frac{\mathsf{log1p}\left(x\right)}{n}} - {x}^{\left(\frac{1}{n}\right)}\\
\end{array}double f(double x, double n) {
double r3659827 = x;
double r3659828 = 1.0;
double r3659829 = r3659827 + r3659828;
double r3659830 = n;
double r3659831 = r3659828 / r3659830;
double r3659832 = pow(r3659829, r3659831);
double r3659833 = pow(r3659827, r3659831);
double r3659834 = r3659832 - r3659833;
return r3659834;
}
double f(double x, double n) {
double r3659835 = 1.0;
double r3659836 = n;
double r3659837 = r3659835 / r3659836;
double r3659838 = -3.2155947364724913e-22;
bool r3659839 = r3659837 <= r3659838;
double r3659840 = x;
double r3659841 = cbrt(r3659840);
double r3659842 = pow(r3659841, r3659837);
double r3659843 = -r3659842;
double r3659844 = r3659841 * r3659841;
double r3659845 = pow(r3659844, r3659837);
double r3659846 = r3659845 * r3659842;
double r3659847 = fma(r3659843, r3659845, r3659846);
double r3659848 = r3659840 + r3659835;
double r3659849 = pow(r3659848, r3659837);
double r3659850 = r3659845 * r3659843;
double r3659851 = fma(r3659835, r3659849, r3659850);
double r3659852 = r3659847 + r3659851;
double r3659853 = r3659849 - r3659846;
double r3659854 = r3659852 * r3659853;
double r3659855 = r3659854 * r3659852;
double r3659856 = cbrt(r3659855);
double r3659857 = r3659847 + r3659856;
double r3659858 = 4.2228249751738446e-17;
bool r3659859 = r3659837 <= r3659858;
double r3659860 = log(r3659840);
double r3659861 = r3659836 * r3659836;
double r3659862 = r3659840 * r3659861;
double r3659863 = r3659860 / r3659862;
double r3659864 = r3659840 * r3659836;
double r3659865 = r3659835 / r3659864;
double r3659866 = r3659863 + r3659865;
double r3659867 = 0.5;
double r3659868 = r3659867 / r3659836;
double r3659869 = r3659840 * r3659840;
double r3659870 = r3659868 / r3659869;
double r3659871 = r3659866 - r3659870;
double r3659872 = log1p(r3659840);
double r3659873 = r3659872 / r3659836;
double r3659874 = exp(r3659873);
double r3659875 = pow(r3659840, r3659837);
double r3659876 = r3659874 - r3659875;
double r3659877 = r3659859 ? r3659871 : r3659876;
double r3659878 = r3659839 ? r3659857 : r3659877;
return r3659878;
}



Bits error versus x



Bits error versus n
if (/ 1 n) < -3.2155947364724913e-22Initial program 2.7
rmApplied add-cube-cbrt2.7
Applied unpow-prod-down2.7
Applied sqr-pow2.7
Applied prod-diff2.8
rmApplied add-cbrt-cube2.8
Simplified2.7
rmApplied *-un-lft-identity2.7
Applied prod-diff2.7
rmApplied *-un-lft-identity2.7
Applied prod-diff2.7
if -3.2155947364724913e-22 < (/ 1 n) < 4.2228249751738446e-17Initial program 44.8
Taylor expanded around inf 32.2
Simplified32.2
if 4.2228249751738446e-17 < (/ 1 n) Initial program 26.6
rmApplied add-exp-log26.6
Applied pow-exp26.6
Simplified4.8
Final simplification19.6
herbie shell --seed 2019163 +o rules:numerics
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
(- (pow (+ x 1) (/ 1 n)) (pow x (/ 1 n))))