#include <tgmath.h>
#include <gmp.h>
#include <mpfr.h>
#include <stdio.h>
#include <stdbool.h>

char *name = "Henrywood and Agarwal, Equation (9a)";

double f_if(float w0, float M, float D, float h, float l, float d) {
        float r9550 = w0;
        float r9551 = 1.0f;
        float r9552 = M;
        float r9553 = D;
        float r9554 = r9552 * r9553;
        float r9555 = 2.0f;
        float r9556 = d;
        float r9557 = r9555 * r9556;
        float r9558 = r9554 / r9557;
        float r9559 = pow(r9558, r9555);
        float r9560 = h;
        float r9561 = l;
        float r9562 = r9560 / r9561;
        float r9563 = r9559 * r9562;
        float r9564 = r9551 - r9563;
        float r9565 = sqrt(r9564);
        float r9566 = r9550 * r9565;
        return r9566;
}

double f_id(double w0, double M, double D, double h, double l, double d) {
        double r9567 = w0;
        double r9568 = 1.0;
        double r9569 = M;
        double r9570 = D;
        double r9571 = r9569 * r9570;
        double r9572 = 2.0;
        double r9573 = d;
        double r9574 = r9572 * r9573;
        double r9575 = r9571 / r9574;
        double r9576 = pow(r9575, r9572);
        double r9577 = h;
        double r9578 = l;
        double r9579 = r9577 / r9578;
        double r9580 = r9576 * r9579;
        double r9581 = r9568 - r9580;
        double r9582 = sqrt(r9581);
        double r9583 = r9567 * r9582;
        return r9583;
}


double f_of(float w0, float M, float D, float h, float l, float d) {
        float r9584 = 1.0f;
        float r9585 = l;
        float r9586 = r9584 / r9585;
        float r9587 = D;
        float r9588 = 2.0f;
        float r9589 = r9587 / r9588;
        float r9590 = d;
        float r9591 = M;
        float r9592 = r9590 / r9591;
        float r9593 = r9589 / r9592;
        float r9594 = cbrt(r9593);
        float r9595 = r9594 * r9594;
        float r9596 = cbrt(r9595);
        float r9597 = cbrt(r9594);
        float r9598 = r9596 * r9597;
        float r9599 = r9598 * r9593;
        float r9600 = r9591 * r9587;
        float r9601 = r9588 * r9590;
        float r9602 = r9600 / r9601;
        float r9603 = cbrt(r9602);
        float r9604 = pow(r9603, r9588);
        float r9605 = h;
        float r9606 = r9604 * r9605;
        float r9607 = r9599 * r9606;
        float r9608 = r9586 * r9607;
        float r9609 = r9584 - r9608;
        float r9610 = sqrt(r9609);
        float r9611 = w0;
        float r9612 = r9610 * r9611;
        return r9612;
}

double f_od(double w0, double M, double D, double h, double l, double d) {
        double r9613 = 1.0;
        double r9614 = l;
        double r9615 = r9613 / r9614;
        double r9616 = D;
        double r9617 = 2.0;
        double r9618 = r9616 / r9617;
        double r9619 = d;
        double r9620 = M;
        double r9621 = r9619 / r9620;
        double r9622 = r9618 / r9621;
        double r9623 = cbrt(r9622);
        double r9624 = r9623 * r9623;
        double r9625 = cbrt(r9624);
        double r9626 = cbrt(r9623);
        double r9627 = r9625 * r9626;
        double r9628 = r9627 * r9622;
        double r9629 = r9620 * r9616;
        double r9630 = r9617 * r9619;
        double r9631 = r9629 / r9630;
        double r9632 = cbrt(r9631);
        double r9633 = pow(r9632, r9617);
        double r9634 = h;
        double r9635 = r9633 * r9634;
        double r9636 = r9628 * r9635;
        double r9637 = r9615 * r9636;
        double r9638 = r9613 - r9637;
        double r9639 = sqrt(r9638);
        double r9640 = w0;
        double r9641 = r9639 * r9640;
        return r9641;
}

void mpfr_fmod2(mpfr_t r, mpfr_t n, mpfr_t d, mpfr_rnd_t rmd) {
        mpfr_fmod(r, n, d, rmd);
        if (mpfr_cmp_ui(r, 0) < 0) mpfr_add(r, r, d, rmd);
}


static mpfr_t r9642, r9643, r9644, r9645, r9646, r9647, r9648, r9649, r9650, r9651, r9652, r9653, r9654, r9655, r9656, r9657, r9658;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r9642);
        mpfr_init_set_str(r9643, "1", 10, MPFR_RNDN);
        mpfr_init(r9644);
        mpfr_init(r9645);
        mpfr_init(r9646);
        mpfr_init_set_str(r9647, "2", 10, MPFR_RNDN);
        mpfr_init(r9648);
        mpfr_init(r9649);
        mpfr_init(r9650);
        mpfr_init(r9651);
        mpfr_init(r9652);
        mpfr_init(r9653);
        mpfr_init(r9654);
        mpfr_init(r9655);
        mpfr_init(r9656);
        mpfr_init(r9657);
        mpfr_init(r9658);
}

double f_im(double w0, double M, double D, double h, double l, double d) {
        mpfr_set_d(r9642, w0, MPFR_RNDN);
        ;
        mpfr_set_d(r9644, M, MPFR_RNDN);
        mpfr_set_d(r9645, D, MPFR_RNDN);
        mpfr_mul(r9646, r9644, r9645, MPFR_RNDN);
        ;
        mpfr_set_d(r9648, d, MPFR_RNDN);
        mpfr_mul(r9649, r9647, r9648, MPFR_RNDN);
        mpfr_div(r9650, r9646, r9649, MPFR_RNDN);
        mpfr_pow(r9651, r9650, r9647, MPFR_RNDN);
        mpfr_set_d(r9652, h, MPFR_RNDN);
        mpfr_set_d(r9653, l, MPFR_RNDN);
        mpfr_div(r9654, r9652, r9653, MPFR_RNDN);
        mpfr_mul(r9655, r9651, r9654, MPFR_RNDN);
        mpfr_sub(r9656, r9643, r9655, MPFR_RNDN);
        mpfr_sqrt(r9657, r9656, MPFR_RNDN);
        mpfr_mul(r9658, r9642, r9657, MPFR_RNDN);
        return mpfr_get_d(r9658, MPFR_RNDN);
}

static mpfr_t r9659, r9660, r9661, r9662, r9663, r9664, r9665, r9666, r9667, r9668, r9669, r9670, r9671, r9672, r9673, r9674, r9675, r9676, r9677, r9678, r9679, r9680, r9681, r9682, r9683, r9684, r9685, r9686, r9687;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r9659, "1", 10, MPFR_RNDN);
        mpfr_init(r9660);
        mpfr_init(r9661);
        mpfr_init(r9662);
        mpfr_init_set_str(r9663, "2", 10, MPFR_RNDN);
        mpfr_init(r9664);
        mpfr_init(r9665);
        mpfr_init(r9666);
        mpfr_init(r9667);
        mpfr_init(r9668);
        mpfr_init(r9669);
        mpfr_init(r9670);
        mpfr_init(r9671);
        mpfr_init(r9672);
        mpfr_init(r9673);
        mpfr_init(r9674);
        mpfr_init(r9675);
        mpfr_init(r9676);
        mpfr_init(r9677);
        mpfr_init(r9678);
        mpfr_init(r9679);
        mpfr_init(r9680);
        mpfr_init(r9681);
        mpfr_init(r9682);
        mpfr_init(r9683);
        mpfr_init(r9684);
        mpfr_init(r9685);
        mpfr_init(r9686);
        mpfr_init(r9687);
}

double f_fm(double w0, double M, double D, double h, double l, double d) {
        ;
        mpfr_set_d(r9660, l, MPFR_RNDN);
        mpfr_div(r9661, r9659, r9660, MPFR_RNDN);
        mpfr_set_d(r9662, D, MPFR_RNDN);
        ;
        mpfr_div(r9664, r9662, r9663, MPFR_RNDN);
        mpfr_set_d(r9665, d, MPFR_RNDN);
        mpfr_set_d(r9666, M, MPFR_RNDN);
        mpfr_div(r9667, r9665, r9666, MPFR_RNDN);
        mpfr_div(r9668, r9664, r9667, MPFR_RNDN);
        mpfr_cbrt(r9669, r9668, MPFR_RNDN);
        mpfr_mul(r9670, r9669, r9669, MPFR_RNDN);
        mpfr_cbrt(r9671, r9670, MPFR_RNDN);
        mpfr_cbrt(r9672, r9669, MPFR_RNDN);
        mpfr_mul(r9673, r9671, r9672, MPFR_RNDN);
        mpfr_mul(r9674, r9673, r9668, MPFR_RNDN);
        mpfr_mul(r9675, r9666, r9662, MPFR_RNDN);
        mpfr_mul(r9676, r9663, r9665, MPFR_RNDN);
        mpfr_div(r9677, r9675, r9676, MPFR_RNDN);
        mpfr_cbrt(r9678, r9677, MPFR_RNDN);
        mpfr_pow(r9679, r9678, r9663, MPFR_RNDN);
        mpfr_set_d(r9680, h, MPFR_RNDN);
        mpfr_mul(r9681, r9679, r9680, MPFR_RNDN);
        mpfr_mul(r9682, r9674, r9681, MPFR_RNDN);
        mpfr_mul(r9683, r9661, r9682, MPFR_RNDN);
        mpfr_sub(r9684, r9659, r9683, MPFR_RNDN);
        mpfr_sqrt(r9685, r9684, MPFR_RNDN);
        mpfr_set_d(r9686, w0, MPFR_RNDN);
        mpfr_mul(r9687, r9685, r9686, MPFR_RNDN);
        return mpfr_get_d(r9687, MPFR_RNDN);
}

static mpfr_t r9688, r9689, r9690, r9691, r9692, r9693, r9694, r9695, r9696, r9697, r9698, r9699, r9700, r9701, r9702, r9703, r9704, r9705, r9706, r9707, r9708, r9709, r9710, r9711, r9712, r9713, r9714, r9715, r9716;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r9688, "1", 10, MPFR_RNDN);
        mpfr_init(r9689);
        mpfr_init(r9690);
        mpfr_init(r9691);
        mpfr_init_set_str(r9692, "2", 10, MPFR_RNDN);
        mpfr_init(r9693);
        mpfr_init(r9694);
        mpfr_init(r9695);
        mpfr_init(r9696);
        mpfr_init(r9697);
        mpfr_init(r9698);
        mpfr_init(r9699);
        mpfr_init(r9700);
        mpfr_init(r9701);
        mpfr_init(r9702);
        mpfr_init(r9703);
        mpfr_init(r9704);
        mpfr_init(r9705);
        mpfr_init(r9706);
        mpfr_init(r9707);
        mpfr_init(r9708);
        mpfr_init(r9709);
        mpfr_init(r9710);
        mpfr_init(r9711);
        mpfr_init(r9712);
        mpfr_init(r9713);
        mpfr_init(r9714);
        mpfr_init(r9715);
        mpfr_init(r9716);
}

double f_dm(double w0, double M, double D, double h, double l, double d) {
        ;
        mpfr_set_d(r9689, l, MPFR_RNDN);
        mpfr_div(r9690, r9688, r9689, MPFR_RNDN);
        mpfr_set_d(r9691, D, MPFR_RNDN);
        ;
        mpfr_div(r9693, r9691, r9692, MPFR_RNDN);
        mpfr_set_d(r9694, d, MPFR_RNDN);
        mpfr_set_d(r9695, M, MPFR_RNDN);
        mpfr_div(r9696, r9694, r9695, MPFR_RNDN);
        mpfr_div(r9697, r9693, r9696, MPFR_RNDN);
        mpfr_cbrt(r9698, r9697, MPFR_RNDN);
        mpfr_mul(r9699, r9698, r9698, MPFR_RNDN);
        mpfr_cbrt(r9700, r9699, MPFR_RNDN);
        mpfr_cbrt(r9701, r9698, MPFR_RNDN);
        mpfr_mul(r9702, r9700, r9701, MPFR_RNDN);
        mpfr_mul(r9703, r9702, r9697, MPFR_RNDN);
        mpfr_mul(r9704, r9695, r9691, MPFR_RNDN);
        mpfr_mul(r9705, r9692, r9694, MPFR_RNDN);
        mpfr_div(r9706, r9704, r9705, MPFR_RNDN);
        mpfr_cbrt(r9707, r9706, MPFR_RNDN);
        mpfr_pow(r9708, r9707, r9692, MPFR_RNDN);
        mpfr_set_d(r9709, h, MPFR_RNDN);
        mpfr_mul(r9710, r9708, r9709, MPFR_RNDN);
        mpfr_mul(r9711, r9703, r9710, MPFR_RNDN);
        mpfr_mul(r9712, r9690, r9711, MPFR_RNDN);
        mpfr_sub(r9713, r9688, r9712, MPFR_RNDN);
        mpfr_sqrt(r9714, r9713, MPFR_RNDN);
        mpfr_set_d(r9715, w0, MPFR_RNDN);
        mpfr_mul(r9716, r9714, r9715, MPFR_RNDN);
        return mpfr_get_d(r9716, MPFR_RNDN);
}