\[\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 r8484 = c0;
        double r8485 = 2.0;
        double r8486 = w;
        double r8487 = r8485 * r8486;
        double r8488 = r8484 / r8487;
        double r8489 = d;
        double r8490 = r8489 * r8489;
        double r8491 = r8484 * r8490;
        double r8492 = h;
        double r8493 = r8486 * r8492;
        double r8494 = D;
        double r8495 = r8494 * r8494;
        double r8496 = r8493 * r8495;
        double r8497 = r8491 / r8496;
        double r8498 = r8497 * r8497;
        double r8499 = M;
        double r8500 = r8499 * r8499;
        double r8501 = r8498 - r8500;
        double r8502 = sqrt(r8501);
        double r8503 = r8497 + r8502;
        double r8504 = r8488 * r8503;
        return r8504;
}

Reproduce

Please include this information when filing a bug report:

herbie shell --seed 2019305 
(FPCore (c0 w h D d M)
  :name "Henrywood and Agarwal, Equation (13)"
  :precision binary64
  (* (/ 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))))))

Backtrace

exact?: contract violation expected: number? given: 'h5LC
egg-parsed->expr/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt650
for-loop/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt706
egg-parsed->expr/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt650
for-loop/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt706
egg-parsed->expr/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt650
(unnamed)/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/simplify.rkt498
egraph-add-exprs/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt1220
egraph-run/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/eggmath.rkt550
(unnamed)/data/pavpan/nightlies/herbie/eggmath-testversion/src/core/simplify.rkt380
simplify!/data/pavpan/nightlies/herbie/eggmath-testversion/src/mainloop.rkt2230
run-improve47/data/pavpan/nightlies/herbie/eggmath-testversion/src/mainloop.rkt3450
(unnamed)/opt/racket-7.2/collects/racket/private/more-scheme.rkt26128
run/opt/racket-7.2/share/pkgs/profile-lib/main.rkt392
profile-thunk16/opt/racket-7.2/share/pkgs/profile-lib/main.rkt90
(unnamed)/opt/racket-7.2/collects/racket/private/more-scheme.rkt26128