\[\frac{c0}{2 \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 \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 r4734834 = c0;
        double r4734835 = 2.0;
        double r4734836 = w;
        double r4734837 = r4734835 * r4734836;
        double r4734838 = r4734834 / r4734837;
        double r4734839 = d;
        double r4734840 = r4734839 * r4734839;
        double r4734841 = r4734834 * r4734840;
        double r4734842 = h;
        double r4734843 = r4734836 * r4734842;
        double r4734844 = D;
        double r4734845 = r4734844 * r4734844;
        double r4734846 = r4734843 * r4734845;
        double r4734847 = r4734841 / r4734846;
        double r4734848 = r4734847 * r4734847;
        double r4734849 = M;
        double r4734850 = r4734849 * r4734849;
        double r4734851 = r4734848 - r4734850;
        double r4734852 = sqrt(r4734851);
        double r4734853 = r4734847 + r4734852;
        double r4734854 = r4734838 * r4734853;
        return r4734854;
}

Reproduce

herbie shell --seed 2019168 +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