{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\begin{array}{l}
\mathbf{if}\;n \le -13450.5463236965716:\\
\;\;\;\;\frac{\frac{1}{n}}{x} - \left(\log \left(e^{\frac{0.5}{{x}^{2} \cdot n}}\right) - \frac{\log x \cdot 1}{x \cdot {n}^{2}}\right)\\
\mathbf{elif}\;n \le 1358332.53943660879:\\
\;\;\;\;\log \left(e^{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{n}}{x} - \left(\frac{\frac{0.5}{n}}{{x}^{2}} - \frac{\frac{\log x \cdot 1}{x}}{{n}^{2}}\right)\\
\end{array}double f(double x, double n) {
double r38842 = x;
double r38843 = 1.0;
double r38844 = r38842 + r38843;
double r38845 = n;
double r38846 = r38843 / r38845;
double r38847 = pow(r38844, r38846);
double r38848 = pow(r38842, r38846);
double r38849 = r38847 - r38848;
return r38849;
}
double f(double x, double n) {
double r38850 = n;
double r38851 = -13450.546323696572;
bool r38852 = r38850 <= r38851;
double r38853 = 1.0;
double r38854 = r38853 / r38850;
double r38855 = x;
double r38856 = r38854 / r38855;
double r38857 = 0.5;
double r38858 = 2.0;
double r38859 = pow(r38855, r38858);
double r38860 = r38859 * r38850;
double r38861 = r38857 / r38860;
double r38862 = exp(r38861);
double r38863 = log(r38862);
double r38864 = log(r38855);
double r38865 = r38864 * r38853;
double r38866 = pow(r38850, r38858);
double r38867 = r38855 * r38866;
double r38868 = r38865 / r38867;
double r38869 = r38863 - r38868;
double r38870 = r38856 - r38869;
double r38871 = 1358332.5394366088;
bool r38872 = r38850 <= r38871;
double r38873 = r38855 + r38853;
double r38874 = pow(r38873, r38854);
double r38875 = pow(r38855, r38854);
double r38876 = r38874 - r38875;
double r38877 = exp(r38876);
double r38878 = log(r38877);
double r38879 = r38857 / r38850;
double r38880 = r38879 / r38859;
double r38881 = r38865 / r38855;
double r38882 = r38881 / r38866;
double r38883 = r38880 - r38882;
double r38884 = r38856 - r38883;
double r38885 = r38872 ? r38878 : r38884;
double r38886 = r38852 ? r38870 : r38885;
return r38886;
}



Bits error versus x



Bits error versus n
Results
if n < -13450.546323696572Initial program 44.8
Taylor expanded around inf 32.2
Simplified31.4
rmApplied add-log-exp31.5
Simplified31.5
if -13450.546323696572 < n < 1358332.5394366088Initial program 8.2
rmApplied add-log-exp8.5
Applied add-log-exp8.4
Applied diff-log8.4
Simplified8.4
if 1358332.5394366088 < n Initial program 44.9
Taylor expanded around inf 33.2
Simplified32.4
rmApplied associate-/r*32.4
Final simplification22.0
herbie shell --seed 2020018
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
:precision binary64
(- (pow (+ x 1) (/ 1 n)) (pow x (/ 1 n))))