\[\frac{c0}{2.0 \cdot w} \cdot \left(\frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} + \sqrt{\frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} \cdot \frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} - M \cdot M}\right)\]

\frac{c0}{2.0 \cdot w} \cdot \left(\frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} + \sqrt{\frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} \cdot \frac{c0 \cdot \left(d \cdot d\right)}{\left(w \cdot h\right) \cdot \left(D \cdot D\right)} - M \cdot M}\right)

double f(double c0, double w, double h, double D, double d, double M) {
        double r4377860 = c0;
        double r4377861 = 2.0;
        double r4377862 = w;
        double r4377863 = r4377861 * r4377862;
        double r4377864 = r4377860 / r4377863;
        double r4377865 = d;
        double r4377866 = r4377865 * r4377865;
        double r4377867 = r4377860 * r4377866;
        double r4377868 = h;
        double r4377869 = r4377862 * r4377868;
        double r4377870 = D;
        double r4377871 = r4377870 * r4377870;
        double r4377872 = r4377869 * r4377871;
        double r4377873 = r4377867 / r4377872;
        double r4377874 = r4377873 * r4377873;
        double r4377875 = M;
        double r4377876 = r4377875 * r4377875;
        double r4377877 = r4377874 - r4377876;
        double r4377878 = sqrt(r4377877);
        double r4377879 = r4377873 + r4377878;
        double r4377880 = r4377864 * r4377879;
        return r4377880;
}

Reproduce

herbie shell --seed 2019165 +o rules:numerics
(FPCore (c0 w h D d M)
  :name "Henrywood and Agarwal, Equation (13)"
  (* (/ c0 (* 2.0 w)) (+ (/ (* c0 (* d d)) (* (* w h) (* D D))) (sqrt (- (* (/ (* c0 (* d d)) (* (* w h) (* D D))) (/ (* c0 (* d d)) (* (* w h) (* D D)))) (* M M))))))

Timeout in 10.0m

Reproduce