{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\begin{array}{l}
\mathbf{if}\;n \le -13657.800064202142:\\
\;\;\;\;\left(\frac{\frac{-1}{2}}{\left(x \cdot x\right) \cdot n} + \frac{\frac{1}{n}}{x}\right) + \frac{\frac{\log x}{n \cdot n}}{x}\\
\mathbf{elif}\;n \le 1591665190.7901402:\\
\;\;\;\;{\left(1 + x\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}\\
\mathbf{else}:\\
\;\;\;\;\sqrt{\log \left(e^{\frac{\frac{-1}{2}}{x \cdot \left(x \cdot n\right)}}\right) + \frac{\frac{1}{x}}{n}} \cdot \sqrt{\log \left(e^{\frac{\frac{-1}{2}}{x \cdot \left(x \cdot n\right)}}\right) + \frac{\frac{1}{x}}{n}} + \frac{\frac{\log x}{n \cdot n}}{x}\\
\end{array}double f(double x, double n) {
double r1949143 = x;
double r1949144 = 1.0;
double r1949145 = r1949143 + r1949144;
double r1949146 = n;
double r1949147 = r1949144 / r1949146;
double r1949148 = pow(r1949145, r1949147);
double r1949149 = pow(r1949143, r1949147);
double r1949150 = r1949148 - r1949149;
return r1949150;
}
double f(double x, double n) {
double r1949151 = n;
double r1949152 = -13657.800064202142;
bool r1949153 = r1949151 <= r1949152;
double r1949154 = -0.5;
double r1949155 = x;
double r1949156 = r1949155 * r1949155;
double r1949157 = r1949156 * r1949151;
double r1949158 = r1949154 / r1949157;
double r1949159 = 1.0;
double r1949160 = r1949159 / r1949151;
double r1949161 = r1949160 / r1949155;
double r1949162 = r1949158 + r1949161;
double r1949163 = log(r1949155);
double r1949164 = r1949151 * r1949151;
double r1949165 = r1949163 / r1949164;
double r1949166 = r1949165 / r1949155;
double r1949167 = r1949162 + r1949166;
double r1949168 = 1591665190.7901402;
bool r1949169 = r1949151 <= r1949168;
double r1949170 = r1949159 + r1949155;
double r1949171 = pow(r1949170, r1949160);
double r1949172 = pow(r1949155, r1949160);
double r1949173 = r1949171 - r1949172;
double r1949174 = r1949155 * r1949151;
double r1949175 = r1949155 * r1949174;
double r1949176 = r1949154 / r1949175;
double r1949177 = exp(r1949176);
double r1949178 = log(r1949177);
double r1949179 = r1949159 / r1949155;
double r1949180 = r1949179 / r1949151;
double r1949181 = r1949178 + r1949180;
double r1949182 = sqrt(r1949181);
double r1949183 = r1949182 * r1949182;
double r1949184 = r1949183 + r1949166;
double r1949185 = r1949169 ? r1949173 : r1949184;
double r1949186 = r1949153 ? r1949167 : r1949185;
return r1949186;
}



Bits error versus x



Bits error versus n
Results
if n < -13657.800064202142Initial program 45.0
rmApplied *-un-lft-identity45.0
Applied pow-unpow45.0
Simplified45.0
Taylor expanded around inf 33.2
Simplified32.4
rmApplied add-log-exp32.3
Taylor expanded around -inf 33.2
Simplified32.4
if -13657.800064202142 < n < 1591665190.7901402Initial program 8.9
rmApplied *-un-lft-identity8.9
Applied pow-unpow8.9
Simplified8.9
if 1591665190.7901402 < n Initial program 45.7
rmApplied *-un-lft-identity45.7
Applied pow-unpow45.7
Simplified45.7
Taylor expanded around inf 32.3
Simplified31.6
rmApplied add-log-exp31.7
rmApplied add-sqr-sqrt31.5
Final simplification22.1
herbie shell --seed 2019144
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
(- (pow (+ x 1) (/ 1 n)) (pow x (/ 1 n))))