{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\begin{array}{l}
\mathbf{if}\;n \le -48075474.504670314:\\
\;\;\;\;\frac{\frac{1}{n}}{x} + \left(\frac{\frac{\log x}{x}}{n \cdot n} - \frac{\frac{\frac{1}{2}}{x \cdot x}}{n}\right)\\
\mathbf{elif}\;n \le 55908.36668550803:\\
\;\;\;\;\sqrt[3]{{\left(1 + x\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \log \left(e^{\sqrt[3]{\left(\sqrt{{\left(1 + x\right)}^{\left(\frac{1}{n}\right)}} - \sqrt{{x}^{\left(\frac{1}{n}\right)}}\right) \cdot \left(\sqrt{{x}^{\left(\frac{1}{n}\right)}} + \sqrt{{\left(1 + x\right)}^{\left(\frac{1}{n}\right)}}\right)} \cdot \sqrt[3]{{\left(1 + x\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}}}\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{\frac{1}{n}}{x} + \left(\frac{\frac{\log x}{x}}{n \cdot n} - \frac{\frac{\frac{1}{2}}{x \cdot x}}{n}\right)\\
\end{array}double f(double x, double n) {
double r3533962 = x;
double r3533963 = 1.0;
double r3533964 = r3533962 + r3533963;
double r3533965 = n;
double r3533966 = r3533963 / r3533965;
double r3533967 = pow(r3533964, r3533966);
double r3533968 = pow(r3533962, r3533966);
double r3533969 = r3533967 - r3533968;
return r3533969;
}
double f(double x, double n) {
double r3533970 = n;
double r3533971 = -48075474.504670314;
bool r3533972 = r3533970 <= r3533971;
double r3533973 = 1.0;
double r3533974 = r3533973 / r3533970;
double r3533975 = x;
double r3533976 = r3533974 / r3533975;
double r3533977 = log(r3533975);
double r3533978 = r3533977 / r3533975;
double r3533979 = r3533970 * r3533970;
double r3533980 = r3533978 / r3533979;
double r3533981 = 0.5;
double r3533982 = r3533975 * r3533975;
double r3533983 = r3533981 / r3533982;
double r3533984 = r3533983 / r3533970;
double r3533985 = r3533980 - r3533984;
double r3533986 = r3533976 + r3533985;
double r3533987 = 55908.36668550803;
bool r3533988 = r3533970 <= r3533987;
double r3533989 = r3533973 + r3533975;
double r3533990 = pow(r3533989, r3533974);
double r3533991 = pow(r3533975, r3533974);
double r3533992 = r3533990 - r3533991;
double r3533993 = cbrt(r3533992);
double r3533994 = sqrt(r3533990);
double r3533995 = sqrt(r3533991);
double r3533996 = r3533994 - r3533995;
double r3533997 = r3533995 + r3533994;
double r3533998 = r3533996 * r3533997;
double r3533999 = cbrt(r3533998);
double r3534000 = r3533999 * r3533993;
double r3534001 = exp(r3534000);
double r3534002 = log(r3534001);
double r3534003 = r3533993 * r3534002;
double r3534004 = r3533988 ? r3534003 : r3533986;
double r3534005 = r3533972 ? r3533986 : r3534004;
return r3534005;
}



Bits error versus x



Bits error versus n
Results
if n < -48075474.504670314 or 55908.36668550803 < n Initial program 45.6
Taylor expanded around inf 33.2
Simplified32.4
if -48075474.504670314 < n < 55908.36668550803Initial program 8.2
rmApplied add-log-exp8.3
Applied add-log-exp8.3
Applied diff-log8.3
Simplified8.3
rmApplied add-cube-cbrt8.3
Applied exp-prod8.3
Applied log-pow8.3
rmApplied add-sqr-sqrt8.3
Applied add-sqr-sqrt8.3
Applied difference-of-squares8.3
Final simplification22.0
herbie shell --seed 2019158
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
(- (pow (+ x 1) (/ 1 n)) (pow x (/ 1 n))))