w0 \cdot \sqrt{1 - {\left(\frac{M \cdot D}{2 \cdot d}\right)}^{2} \cdot \frac{h}{\ell}}w0 \cdot \sqrt{1 - \left(h \cdot {\left(\frac{D}{d} \cdot \frac{M}{2}\right)}^{\left(\frac{2}{2}\right)}\right) \cdot \frac{{\left(\frac{D}{d} \cdot \frac{M}{2}\right)}^{\left(\frac{2}{2}\right)}}{\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) (h * ((double) pow(((double) (((double) (D / d)) * ((double) (M / 2.0)))), ((double) (2.0 / 2.0)))))) * ((double) (((double) pow(((double) (((double) (D / d)) * ((double) (M / 2.0)))), ((double) (2.0 / 2.0)))) / 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.0
rmApplied associate-*r/10.8
Simplified10.9
rmApplied sqr-pow10.9
Applied associate-*r*9.2
Simplified9.2
rmApplied *-un-lft-identity9.2
Applied times-frac8.3
Simplified8.3
Final simplification8.3
herbie shell --seed 2020184
(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))))))