\[\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 r10357298 = c0;
        double r10357299 = 2.0;
        double r10357300 = w;
        double r10357301 = r10357299 * r10357300;
        double r10357302 = r10357298 / r10357301;
        double r10357303 = d;
        double r10357304 = r10357303 * r10357303;
        double r10357305 = r10357298 * r10357304;
        double r10357306 = h;
        double r10357307 = r10357300 * r10357306;
        double r10357308 = D;
        double r10357309 = r10357308 * r10357308;
        double r10357310 = r10357307 * r10357309;
        double r10357311 = r10357305 / r10357310;
        double r10357312 = r10357311 * r10357311;
        double r10357313 = M;
        double r10357314 = r10357313 * r10357313;
        double r10357315 = r10357312 - r10357314;
        double r10357316 = sqrt(r10357315);
        double r10357317 = r10357311 + r10357316;
        double r10357318 = r10357302 * r10357317;
        return r10357318;
}

Reproduce

herbie shell --seed 2019158 +o rules:numerics
(FPCore (c0 w h D d M)
  :name "Henrywood and Agarwal, Equation (13)"
  (* (/ c0 (* 2 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