{\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}:\\
\;\;\;\;\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(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}}\right)}\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(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}}\right)}\\
\mathbf{elif}\;\frac{1}{n} \le 2.865758198043198947866337655820956342723 \cdot 10^{-5}:\\
\;\;\;\;\left(\frac{\frac{1}{x}}{n} + \frac{1 \cdot \log x}{\left(n \cdot n\right) \cdot x}\right) - \frac{\frac{0.5}{n}}{x \cdot x}\\
\mathbf{else}:\\
\;\;\;\;\left(\sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\right)}} \cdot \sqrt[3]{\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) \cdot \sqrt[3]{{\left(x + 1\right)}^{\left(\frac{1}{n}\right)} - {x}^{\left(\frac{1}{n}\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 r2580360 = x;
double r2580361 = 1.0;
double r2580362 = r2580360 + r2580361;
double r2580363 = n;
double r2580364 = r2580361 / r2580363;
double r2580365 = pow(r2580362, r2580364);
double r2580366 = pow(r2580360, r2580364);
double r2580367 = r2580365 - r2580366;
return r2580367;
}
double f(double x, double n) {
double r2580368 = 1.0;
double r2580369 = n;
double r2580370 = r2580368 / r2580369;
double r2580371 = -8.249355031536442e-07;
bool r2580372 = r2580370 <= r2580371;
double r2580373 = x;
double r2580374 = r2580373 + r2580368;
double r2580375 = pow(r2580374, r2580370);
double r2580376 = pow(r2580373, r2580370);
double r2580377 = r2580375 - r2580376;
double r2580378 = cbrt(r2580377);
double r2580379 = r2580378 * r2580378;
double r2580380 = r2580378 * r2580379;
double r2580381 = cbrt(r2580380);
double r2580382 = r2580378 * r2580381;
double r2580383 = r2580382 * r2580381;
double r2580384 = 2.865758198043199e-05;
bool r2580385 = r2580370 <= r2580384;
double r2580386 = r2580368 / r2580373;
double r2580387 = r2580386 / r2580369;
double r2580388 = log(r2580373);
double r2580389 = r2580368 * r2580388;
double r2580390 = r2580369 * r2580369;
double r2580391 = r2580390 * r2580373;
double r2580392 = r2580389 / r2580391;
double r2580393 = r2580387 + r2580392;
double r2580394 = 0.5;
double r2580395 = r2580394 / r2580369;
double r2580396 = r2580373 * r2580373;
double r2580397 = r2580395 / r2580396;
double r2580398 = r2580393 - r2580397;
double r2580399 = sqrt(r2580375);
double r2580400 = 2.0;
double r2580401 = r2580370 / r2580400;
double r2580402 = pow(r2580373, r2580401);
double r2580403 = r2580399 + r2580402;
double r2580404 = r2580399 - r2580402;
double r2580405 = r2580403 * r2580404;
double r2580406 = cbrt(r2580405);
double r2580407 = r2580378 * r2580406;
double r2580408 = r2580407 * r2580378;
double r2580409 = cbrt(r2580408);
double r2580410 = r2580378 * r2580409;
double r2580411 = exp(r2580377);
double r2580412 = log(r2580411);
double r2580413 = cbrt(r2580412);
double r2580414 = r2580410 * r2580413;
double r2580415 = r2580385 ? r2580398 : r2580414;
double r2580416 = r2580372 ? r2580383 : r2580415;
return r2580416;
}



Bits error versus x



Bits error versus n
Results
if (/ 1.0 n) < -8.249355031536442e-07Initial program 0.8
rmApplied add-cube-cbrt0.8
rmApplied add-cube-cbrt0.8
rmApplied add-cube-cbrt0.8
if -8.249355031536442e-07 < (/ 1.0 n) < 2.865758198043199e-05Initial program 44.9
Taylor expanded around inf 32.4
Simplified31.9
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.0
herbie shell --seed 2019171
(FPCore (x n)
:name "2nthrt (problem 3.4.6)"
(- (pow (+ x 1.0) (/ 1.0 n)) (pow x (/ 1.0 n))))