\left({\left(\frac{d}{h}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\ell}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left(1 - \left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{h}{\ell}\right)\left(\sqrt[3]{\left(\left({\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{1}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{\sqrt[3]{d}}{h}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left({\left(\frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - \left(\left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{\sqrt[3]{h} \cdot \sqrt[3]{h}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right) \cdot \frac{\sqrt[3]{h}}{\sqrt[3]{\ell}}\right)} \cdot \sqrt[3]{\left(\left({\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{1}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{\sqrt[3]{d}}{h}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left({\left(\frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - \left(\left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{\sqrt[3]{h} \cdot \sqrt[3]{h}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right) \cdot \frac{\sqrt[3]{h}}{\sqrt[3]{\ell}}\right)}\right) \cdot \sqrt[3]{\left(\left({\left(\frac{\sqrt[3]{d} \cdot \sqrt[3]{d}}{1}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{\sqrt[3]{d}}{h}\right)}^{\left(\frac{1}{2}\right)}\right) \cdot \left({\left(\frac{1}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)} \cdot {\left(\frac{d}{\sqrt[3]{\ell}}\right)}^{\left(\frac{1}{2}\right)}\right)\right) \cdot \left(1 - \left(\left(\frac{1}{2} \cdot {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2}\right) \cdot \frac{\sqrt[3]{h} \cdot \sqrt[3]{h}}{\sqrt[3]{\ell} \cdot \sqrt[3]{\ell}}\right) \cdot \frac{\sqrt[3]{h}}{\sqrt[3]{\ell}}\right)}double code(double d, double h, double l, double M, double D) {
return ((double) (((double) (((double) pow(((double) (d / h)), ((double) (1.0 / 2.0)))) * ((double) pow(((double) (d / l)), ((double) (1.0 / 2.0)))))) * ((double) (1.0 - ((double) (((double) (((double) (1.0 / 2.0)) * ((double) pow(((double) (((double) (M * D)) / ((double) (2.0 * d)))), 2.0)))) * ((double) (h / l))))))));
}
double code(double d, double h, double l, double M, double D) {
return ((double) (((double) (((double) cbrt(((double) (((double) (((double) (((double) pow(((double) (((double) (((double) cbrt(d)) * ((double) cbrt(d)))) / 1.0)), ((double) (1.0 / 2.0)))) * ((double) pow(((double) (((double) cbrt(d)) / h)), ((double) (1.0 / 2.0)))))) * ((double) (((double) pow(((double) (1.0 / ((double) (((double) cbrt(l)) * ((double) cbrt(l)))))), ((double) (1.0 / 2.0)))) * ((double) pow(((double) (d / ((double) cbrt(l)))), ((double) (1.0 / 2.0)))))))) * ((double) (1.0 - ((double) (((double) (((double) (((double) (1.0 / 2.0)) * ((double) pow(((double) (((double) (M * D)) / ((double) (2.0 * d)))), 2.0)))) * ((double) (((double) (((double) cbrt(h)) * ((double) cbrt(h)))) / ((double) (((double) cbrt(l)) * ((double) cbrt(l)))))))) * ((double) (((double) cbrt(h)) / ((double) cbrt(l)))))))))))) * ((double) cbrt(((double) (((double) (((double) (((double) pow(((double) (((double) (((double) cbrt(d)) * ((double) cbrt(d)))) / 1.0)), ((double) (1.0 / 2.0)))) * ((double) pow(((double) (((double) cbrt(d)) / h)), ((double) (1.0 / 2.0)))))) * ((double) (((double) pow(((double) (1.0 / ((double) (((double) cbrt(l)) * ((double) cbrt(l)))))), ((double) (1.0 / 2.0)))) * ((double) pow(((double) (d / ((double) cbrt(l)))), ((double) (1.0 / 2.0)))))))) * ((double) (1.0 - ((double) (((double) (((double) (((double) (1.0 / 2.0)) * ((double) pow(((double) (((double) (M * D)) / ((double) (2.0 * d)))), 2.0)))) * ((double) (((double) (((double) cbrt(h)) * ((double) cbrt(h)))) / ((double) (((double) cbrt(l)) * ((double) cbrt(l)))))))) * ((double) (((double) cbrt(h)) / ((double) cbrt(l)))))))))))))) * ((double) cbrt(((double) (((double) (((double) (((double) pow(((double) (((double) (((double) cbrt(d)) * ((double) cbrt(d)))) / 1.0)), ((double) (1.0 / 2.0)))) * ((double) pow(((double) (((double) cbrt(d)) / h)), ((double) (1.0 / 2.0)))))) * ((double) (((double) pow(((double) (1.0 / ((double) (((double) cbrt(l)) * ((double) cbrt(l)))))), ((double) (1.0 / 2.0)))) * ((double) pow(((double) (d / ((double) cbrt(l)))), ((double) (1.0 / 2.0)))))))) * ((double) (1.0 - ((double) (((double) (((double) (((double) (1.0 / 2.0)) * ((double) pow(((double) (((double) (M * D)) / ((double) (2.0 * d)))), 2.0)))) * ((double) (((double) (((double) cbrt(h)) * ((double) cbrt(h)))) / ((double) (((double) cbrt(l)) * ((double) cbrt(l)))))))) * ((double) (((double) cbrt(h)) / ((double) cbrt(l))))))))))))));
}



Bits error versus d



Bits error versus h



Bits error versus l



Bits error versus M



Bits error versus D
Results
Initial program 26.7
rmApplied *-un-lft-identity26.7
Applied add-cube-cbrt27.0
Applied times-frac27.0
Applied unpow-prod-down22.3
rmApplied add-cube-cbrt22.4
Applied *-un-lft-identity22.4
Applied times-frac22.4
Applied unpow-prod-down18.5
rmApplied add-cube-cbrt18.6
Applied add-cube-cbrt18.6
Applied times-frac18.6
Applied associate-*r*15.8
rmApplied add-cube-cbrt16.1
Final simplification16.1
herbie shell --seed 2020123 +o rules:numerics
(FPCore (d h l M D)
:name "Henrywood and Agarwal, Equation (12)"
:precision binary64
(* (* (pow (/ d h) (/ 1 2)) (pow (/ d l) (/ 1 2))) (- 1 (* (* (/ 1 2) (pow (/ (* M D) (* 2 d)) 2)) (/ h l)))))