#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 r9581 = w0;
        float r9582 = 1.0f;
        float r9583 = M;
        float r9584 = D;
        float r9585 = r9583 * r9584;
        float r9586 = 2.0f;
        float r9587 = d;
        float r9588 = r9586 * r9587;
        float r9589 = r9585 / r9588;
        float r9590 = pow(r9589, r9586);
        float r9591 = h;
        float r9592 = l;
        float r9593 = r9591 / r9592;
        float r9594 = r9590 * r9593;
        float r9595 = r9582 - r9594;
        float r9596 = sqrt(r9595);
        float r9597 = r9581 * r9596;
        return r9597;
}

double f_id(double w0, double M, double D, double h, double l, double d) {
        double r9598 = w0;
        double r9599 = 1.0;
        double r9600 = M;
        double r9601 = D;
        double r9602 = r9600 * r9601;
        double r9603 = 2.0;
        double r9604 = d;
        double r9605 = r9603 * r9604;
        double r9606 = r9602 / r9605;
        double r9607 = pow(r9606, r9603);
        double r9608 = h;
        double r9609 = l;
        double r9610 = r9608 / r9609;
        double r9611 = r9607 * r9610;
        double r9612 = r9599 - r9611;
        double r9613 = sqrt(r9612);
        double r9614 = r9598 * r9613;
        return r9614;
}


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

double f_od(double w0, double M, double D, double h, double l, double d) {
        double r9644 = 1.0;
        double r9645 = l;
        double r9646 = r9644 / r9645;
        double r9647 = D;
        double r9648 = 2.0;
        double r9649 = r9647 / r9648;
        double r9650 = d;
        double r9651 = M;
        double r9652 = r9650 / r9651;
        double r9653 = r9649 / r9652;
        double r9654 = cbrt(r9653);
        double r9655 = r9654 * r9654;
        double r9656 = cbrt(r9655);
        double r9657 = cbrt(r9654);
        double r9658 = r9656 * r9657;
        double r9659 = r9658 * r9653;
        double r9660 = r9651 * r9647;
        double r9661 = r9648 * r9650;
        double r9662 = r9660 / r9661;
        double r9663 = cbrt(r9662);
        double r9664 = pow(r9663, r9648);
        double r9665 = h;
        double r9666 = r9664 * r9665;
        double r9667 = r9659 * r9666;
        double r9668 = r9646 * r9667;
        double r9669 = r9644 - r9668;
        double r9670 = sqrt(r9669);
        double r9671 = w0;
        double r9672 = r9670 * r9671;
        return r9672;
}

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 r9673, r9674, r9675, r9676, r9677, r9678, r9679, r9680, r9681, r9682, r9683, r9684, r9685, r9686, r9687, r9688, r9689;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(336);
        mpfr_init(r9673);
        mpfr_init_set_str(r9674, "1", 10, MPFR_RNDN);
        mpfr_init(r9675);
        mpfr_init(r9676);
        mpfr_init(r9677);
        mpfr_init_set_str(r9678, "2", 10, MPFR_RNDN);
        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);
        mpfr_init(r9688);
        mpfr_init(r9689);
}

double f_im(double w0, double M, double D, double h, double l, double d) {
        mpfr_set_d(r9673, w0, MPFR_RNDN);
        ;
        mpfr_set_d(r9675, M, MPFR_RNDN);
        mpfr_set_d(r9676, D, MPFR_RNDN);
        mpfr_mul(r9677, r9675, r9676, MPFR_RNDN);
        ;
        mpfr_set_d(r9679, d, MPFR_RNDN);
        mpfr_mul(r9680, r9678, r9679, MPFR_RNDN);
        mpfr_div(r9681, r9677, r9680, MPFR_RNDN);
        mpfr_pow(r9682, r9681, r9678, MPFR_RNDN);
        mpfr_set_d(r9683, h, MPFR_RNDN);
        mpfr_set_d(r9684, l, MPFR_RNDN);
        mpfr_div(r9685, r9683, r9684, MPFR_RNDN);
        mpfr_mul(r9686, r9682, r9685, MPFR_RNDN);
        mpfr_sub(r9687, r9674, r9686, MPFR_RNDN);
        mpfr_sqrt(r9688, r9687, MPFR_RNDN);
        mpfr_mul(r9689, r9673, r9688, MPFR_RNDN);
        return mpfr_get_d(r9689, MPFR_RNDN);
}

static mpfr_t 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, r9717, r9718;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r9690, "1", 10, MPFR_RNDN);
        mpfr_init(r9691);
        mpfr_init(r9692);
        mpfr_init(r9693);
        mpfr_init_set_str(r9694, "2", 10, MPFR_RNDN);
        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);
        mpfr_init(r9717);
        mpfr_init(r9718);
}

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

static mpfr_t r9719, r9720, r9721, r9722, r9723, r9724, r9725, r9726, r9727, r9728, r9729, r9730, r9731, r9732, r9733, r9734, r9735, r9736, r9737, r9738, r9739, r9740, r9741, r9742, r9743, r9744, r9745, r9746, r9747;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(336);
        mpfr_init_set_str(r9719, "1", 10, MPFR_RNDN);
        mpfr_init(r9720);
        mpfr_init(r9721);
        mpfr_init(r9722);
        mpfr_init_set_str(r9723, "2", 10, MPFR_RNDN);
        mpfr_init(r9724);
        mpfr_init(r9725);
        mpfr_init(r9726);
        mpfr_init(r9727);
        mpfr_init(r9728);
        mpfr_init(r9729);
        mpfr_init(r9730);
        mpfr_init(r9731);
        mpfr_init(r9732);
        mpfr_init(r9733);
        mpfr_init(r9734);
        mpfr_init(r9735);
        mpfr_init(r9736);
        mpfr_init(r9737);
        mpfr_init(r9738);
        mpfr_init(r9739);
        mpfr_init(r9740);
        mpfr_init(r9741);
        mpfr_init(r9742);
        mpfr_init(r9743);
        mpfr_init(r9744);
        mpfr_init(r9745);
        mpfr_init(r9746);
        mpfr_init(r9747);
}

double f_dm(double w0, double M, double D, double h, double l, double d) {
        ;
        mpfr_set_d(r9720, l, MPFR_RNDN);
        mpfr_div(r9721, r9719, r9720, MPFR_RNDN);
        mpfr_set_d(r9722, D, MPFR_RNDN);
        ;
        mpfr_div(r9724, r9722, r9723, MPFR_RNDN);
        mpfr_set_d(r9725, d, MPFR_RNDN);
        mpfr_set_d(r9726, M, MPFR_RNDN);
        mpfr_div(r9727, r9725, r9726, MPFR_RNDN);
        mpfr_div(r9728, r9724, r9727, MPFR_RNDN);
        mpfr_cbrt(r9729, r9728, MPFR_RNDN);
        mpfr_mul(r9730, r9729, r9729, MPFR_RNDN);
        mpfr_cbrt(r9731, r9730, MPFR_RNDN);
        mpfr_cbrt(r9732, r9729, MPFR_RNDN);
        mpfr_mul(r9733, r9731, r9732, MPFR_RNDN);
        mpfr_mul(r9734, r9733, r9728, MPFR_RNDN);
        mpfr_mul(r9735, r9726, r9722, MPFR_RNDN);
        mpfr_mul(r9736, r9723, r9725, MPFR_RNDN);
        mpfr_div(r9737, r9735, r9736, MPFR_RNDN);
        mpfr_cbrt(r9738, r9737, MPFR_RNDN);
        mpfr_pow(r9739, r9738, r9723, MPFR_RNDN);
        mpfr_set_d(r9740, h, MPFR_RNDN);
        mpfr_mul(r9741, r9739, r9740, MPFR_RNDN);
        mpfr_mul(r9742, r9734, r9741, MPFR_RNDN);
        mpfr_mul(r9743, r9721, r9742, MPFR_RNDN);
        mpfr_sub(r9744, r9719, r9743, MPFR_RNDN);
        mpfr_sqrt(r9745, r9744, MPFR_RNDN);
        mpfr_set_d(r9746, w0, MPFR_RNDN);
        mpfr_mul(r9747, r9745, r9746, MPFR_RNDN);
        return mpfr_get_d(r9747, MPFR_RNDN);
}