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



Bits error versus w0



Bits error versus M



Bits error versus D



Bits error versus h



Bits error versus l



Bits error versus d
Results
Initial program 14.1
rmApplied sqr-pow14.1
Applied associate-*l*12.6
rmApplied associate-*r/8.5
rmApplied clear-num8.5
Final simplification8.5
herbie shell --seed 2020162
(FPCore (w0 M D h l d)
:name "Henrywood and Agarwal, Equation (9a)"
:precision binary64
(* w0 (sqrt (- 1.0 (* (pow (/ (* M D) (* 2.0 d)) 2.0) (/ h l))))))