\frac{a \cdot {k}^{m}}{\left(1 + 10 \cdot k\right) + k \cdot k}\begin{array}{l}
\mathbf{if}\;k \le 31303448207475440:\\
\;\;\;\;\frac{{\left(\sqrt[3]{k}\right)}^{m}}{\frac{1}{{\left(\sqrt[3]{\sqrt[3]{k}} \cdot \sqrt[3]{\sqrt[3]{k}}\right)}^{m}}} \cdot \left(\frac{{\left(\sqrt[3]{k}\right)}^{m}}{\frac{k \cdot \left(10 + k\right) + 1}{{\left(\sqrt[3]{\sqrt[3]{k}}\right)}^{m}}} \cdot a\right)\\
\mathbf{else}:\\
\;\;\;\;\frac{99 \cdot \left(\left({\left({\left(\frac{1}{k}\right)}^{\frac{-2}{3}}\right)}^{m} \cdot a\right) \cdot {\left({\left(\frac{1}{k}\right)}^{\frac{-1}{3}}\right)}^{m}\right)}{{k}^{4}} + \left(\frac{{\left({\left(\frac{1}{k}\right)}^{\frac{-2}{3}}\right)}^{m} \cdot a}{k} \cdot \frac{{\left({\left(\frac{1}{k}\right)}^{\frac{-1}{3}}\right)}^{m}}{k} - \frac{10 \cdot \left(\left({\left({\left(\frac{1}{k}\right)}^{\frac{-2}{3}}\right)}^{m} \cdot a\right) \cdot {\left({\left(\frac{1}{k}\right)}^{\frac{-1}{3}}\right)}^{m}\right)}{{k}^{3}}\right)\\
\end{array}double f(double a, double k, double m) {
double r132292 = a;
double r132293 = k;
double r132294 = m;
double r132295 = pow(r132293, r132294);
double r132296 = r132292 * r132295;
double r132297 = 1.0;
double r132298 = 10.0;
double r132299 = r132298 * r132293;
double r132300 = r132297 + r132299;
double r132301 = r132293 * r132293;
double r132302 = r132300 + r132301;
double r132303 = r132296 / r132302;
return r132303;
}
double f(double a, double k, double m) {
double r132304 = k;
double r132305 = 3.130344820747544e+16;
bool r132306 = r132304 <= r132305;
double r132307 = cbrt(r132304);
double r132308 = m;
double r132309 = pow(r132307, r132308);
double r132310 = 1.0;
double r132311 = cbrt(r132307);
double r132312 = r132311 * r132311;
double r132313 = pow(r132312, r132308);
double r132314 = r132310 / r132313;
double r132315 = r132309 / r132314;
double r132316 = 10.0;
double r132317 = r132316 + r132304;
double r132318 = r132304 * r132317;
double r132319 = 1.0;
double r132320 = r132318 + r132319;
double r132321 = pow(r132311, r132308);
double r132322 = r132320 / r132321;
double r132323 = r132309 / r132322;
double r132324 = a;
double r132325 = r132323 * r132324;
double r132326 = r132315 * r132325;
double r132327 = 99.0;
double r132328 = r132310 / r132304;
double r132329 = -0.6666666666666666;
double r132330 = pow(r132328, r132329);
double r132331 = pow(r132330, r132308);
double r132332 = r132331 * r132324;
double r132333 = -0.3333333333333333;
double r132334 = pow(r132328, r132333);
double r132335 = pow(r132334, r132308);
double r132336 = r132332 * r132335;
double r132337 = r132327 * r132336;
double r132338 = 4.0;
double r132339 = pow(r132304, r132338);
double r132340 = r132337 / r132339;
double r132341 = r132332 / r132304;
double r132342 = r132335 / r132304;
double r132343 = r132341 * r132342;
double r132344 = r132316 * r132336;
double r132345 = 3.0;
double r132346 = pow(r132304, r132345);
double r132347 = r132344 / r132346;
double r132348 = r132343 - r132347;
double r132349 = r132340 + r132348;
double r132350 = r132306 ? r132326 : r132349;
return r132350;
}



Bits error versus a



Bits error versus k



Bits error versus m
Results
if k < 3.130344820747544e+16Initial program 0.1
Simplified0.0
rmApplied add-cube-cbrt0.0
Applied unpow-prod-down0.0
Applied associate-/l*0.0
rmApplied add-cube-cbrt0.0
Applied unpow-prod-down0.0
Applied *-un-lft-identity0.0
Applied times-frac0.0
Applied unpow-prod-down0.1
Applied times-frac0.1
Applied associate-*l*0.0
if 3.130344820747544e+16 < k Initial program 5.4
Simplified5.4
rmApplied add-cube-cbrt5.4
Applied unpow-prod-down5.4
Applied associate-/l*5.4
Taylor expanded around inf 5.4
Simplified0.3
Final simplification0.1
herbie shell --seed 2019294
(FPCore (a k m)
:name "Falkner and Boettcher, Appendix A"
:precision binary64
(/ (* a (pow k m)) (+ (+ 1 (* 10 k)) (* k k))))