#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 r26065 = w0;
        float r26066 = 1;
        float r26067 = M;
        float r26068 = D;
        float r26069 = r26067 * r26068;
        float r26070 = 2;
        float r26071 = d;
        float r26072 = r26070 * r26071;
        float r26073 = r26069 / r26072;
        float r26074 = pow(r26073, r26070);
        float r26075 = h;
        float r26076 = l;
        float r26077 = r26075 / r26076;
        float r26078 = r26074 * r26077;
        float r26079 = r26066 - r26078;
        float r26080 = sqrt(r26079);
        float r26081 = r26065 * r26080;
        return r26081;
}

double f_id(double w0, double M, double D, double h, double l, double d) {
        double r26082 = w0;
        double r26083 = 1;
        double r26084 = M;
        double r26085 = D;
        double r26086 = r26084 * r26085;
        double r26087 = 2;
        double r26088 = d;
        double r26089 = r26087 * r26088;
        double r26090 = r26086 / r26089;
        double r26091 = pow(r26090, r26087);
        double r26092 = h;
        double r26093 = l;
        double r26094 = r26092 / r26093;
        double r26095 = r26091 * r26094;
        double r26096 = r26083 - r26095;
        double r26097 = sqrt(r26096);
        double r26098 = r26082 * r26097;
        return r26098;
}


double f_of(float w0, float M, float D, float h, float l, float d) {
        float r26099 = w0;
        float r26100 = 1;
        float r26101 = M;
        float r26102 = D;
        float r26103 = r26101 * r26102;
        float r26104 = 2;
        float r26105 = d;
        float r26106 = r26104 * r26105;
        float r26107 = r26103 / r26106;
        float r26108 = h;
        float r26109 = r26107 * r26108;
        float r26110 = r26107 * r26109;
        float r26111 = l;
        float r26112 = r26110 / r26111;
        float r26113 = r26100 - r26112;
        float r26114 = sqrt(r26113);
        float r26115 = r26099 * r26114;
        return r26115;
}

double f_od(double w0, double M, double D, double h, double l, double d) {
        double r26116 = w0;
        double r26117 = 1;
        double r26118 = M;
        double r26119 = D;
        double r26120 = r26118 * r26119;
        double r26121 = 2;
        double r26122 = d;
        double r26123 = r26121 * r26122;
        double r26124 = r26120 / r26123;
        double r26125 = h;
        double r26126 = r26124 * r26125;
        double r26127 = r26124 * r26126;
        double r26128 = l;
        double r26129 = r26127 / r26128;
        double r26130 = r26117 - r26129;
        double r26131 = sqrt(r26130);
        double r26132 = r26116 * r26131;
        return r26132;
}

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 r26133, r26134, r26135, r26136, r26137, r26138, r26139, r26140, r26141, r26142, r26143, r26144, r26145, r26146, r26147, r26148, r26149;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r26133);
        mpfr_init_set_str(r26134, "1", 10, MPFR_RNDN);
        mpfr_init(r26135);
        mpfr_init(r26136);
        mpfr_init(r26137);
        mpfr_init_set_str(r26138, "2", 10, MPFR_RNDN);
        mpfr_init(r26139);
        mpfr_init(r26140);
        mpfr_init(r26141);
        mpfr_init(r26142);
        mpfr_init(r26143);
        mpfr_init(r26144);
        mpfr_init(r26145);
        mpfr_init(r26146);
        mpfr_init(r26147);
        mpfr_init(r26148);
        mpfr_init(r26149);
}

double f_im(double w0, double M, double D, double h, double l, double d) {
        mpfr_set_d(r26133, w0, MPFR_RNDN);
        ;
        mpfr_set_d(r26135, M, MPFR_RNDN);
        mpfr_set_d(r26136, D, MPFR_RNDN);
        mpfr_mul(r26137, r26135, r26136, MPFR_RNDN);
        ;
        mpfr_set_d(r26139, d, MPFR_RNDN);
        mpfr_mul(r26140, r26138, r26139, MPFR_RNDN);
        mpfr_div(r26141, r26137, r26140, MPFR_RNDN);
        mpfr_pow(r26142, r26141, r26138, MPFR_RNDN);
        mpfr_set_d(r26143, h, MPFR_RNDN);
        mpfr_set_d(r26144, l, MPFR_RNDN);
        mpfr_div(r26145, r26143, r26144, MPFR_RNDN);
        mpfr_mul(r26146, r26142, r26145, MPFR_RNDN);
        mpfr_sub(r26147, r26134, r26146, MPFR_RNDN);
        mpfr_sqrt(r26148, r26147, MPFR_RNDN);
        mpfr_mul(r26149, r26133, r26148, MPFR_RNDN);
        return mpfr_get_d(r26149, MPFR_RNDN);
}

static mpfr_t r26150, r26151, r26152, r26153, r26154, r26155, r26156, r26157, r26158, r26159, r26160, r26161, r26162, r26163, r26164, r26165, r26166;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26150);
        mpfr_init_set_str(r26151, "1", 10, MPFR_RNDN);
        mpfr_init(r26152);
        mpfr_init(r26153);
        mpfr_init(r26154);
        mpfr_init_set_str(r26155, "2", 10, MPFR_RNDN);
        mpfr_init(r26156);
        mpfr_init(r26157);
        mpfr_init(r26158);
        mpfr_init(r26159);
        mpfr_init(r26160);
        mpfr_init(r26161);
        mpfr_init(r26162);
        mpfr_init(r26163);
        mpfr_init(r26164);
        mpfr_init(r26165);
        mpfr_init(r26166);
}

double f_fm(double w0, double M, double D, double h, double l, double d) {
        mpfr_set_d(r26150, w0, MPFR_RNDN);
        ;
        mpfr_set_d(r26152, M, MPFR_RNDN);
        mpfr_set_d(r26153, D, MPFR_RNDN);
        mpfr_mul(r26154, r26152, r26153, MPFR_RNDN);
        ;
        mpfr_set_d(r26156, d, MPFR_RNDN);
        mpfr_mul(r26157, r26155, r26156, MPFR_RNDN);
        mpfr_div(r26158, r26154, r26157, MPFR_RNDN);
        mpfr_set_d(r26159, h, MPFR_RNDN);
        mpfr_mul(r26160, r26158, r26159, MPFR_RNDN);
        mpfr_mul(r26161, r26158, r26160, MPFR_RNDN);
        mpfr_set_d(r26162, l, MPFR_RNDN);
        mpfr_div(r26163, r26161, r26162, MPFR_RNDN);
        mpfr_sub(r26164, r26151, r26163, MPFR_RNDN);
        mpfr_sqrt(r26165, r26164, MPFR_RNDN);
        mpfr_mul(r26166, r26150, r26165, MPFR_RNDN);
        return mpfr_get_d(r26166, MPFR_RNDN);
}

static mpfr_t r26167, r26168, r26169, r26170, r26171, r26172, r26173, r26174, r26175, r26176, r26177, r26178, r26179, r26180, r26181, r26182, r26183;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r26167);
        mpfr_init_set_str(r26168, "1", 10, MPFR_RNDN);
        mpfr_init(r26169);
        mpfr_init(r26170);
        mpfr_init(r26171);
        mpfr_init_set_str(r26172, "2", 10, MPFR_RNDN);
        mpfr_init(r26173);
        mpfr_init(r26174);
        mpfr_init(r26175);
        mpfr_init(r26176);
        mpfr_init(r26177);
        mpfr_init(r26178);
        mpfr_init(r26179);
        mpfr_init(r26180);
        mpfr_init(r26181);
        mpfr_init(r26182);
        mpfr_init(r26183);
}

double f_dm(double w0, double M, double D, double h, double l, double d) {
        mpfr_set_d(r26167, w0, MPFR_RNDN);
        ;
        mpfr_set_d(r26169, M, MPFR_RNDN);
        mpfr_set_d(r26170, D, MPFR_RNDN);
        mpfr_mul(r26171, r26169, r26170, MPFR_RNDN);
        ;
        mpfr_set_d(r26173, d, MPFR_RNDN);
        mpfr_mul(r26174, r26172, r26173, MPFR_RNDN);
        mpfr_div(r26175, r26171, r26174, MPFR_RNDN);
        mpfr_set_d(r26176, h, MPFR_RNDN);
        mpfr_mul(r26177, r26175, r26176, MPFR_RNDN);
        mpfr_mul(r26178, r26175, r26177, MPFR_RNDN);
        mpfr_set_d(r26179, l, MPFR_RNDN);
        mpfr_div(r26180, r26178, r26179, MPFR_RNDN);
        mpfr_sub(r26181, r26168, r26180, MPFR_RNDN);
        mpfr_sqrt(r26182, r26181, MPFR_RNDN);
        mpfr_mul(r26183, r26167, r26182, MPFR_RNDN);
        return mpfr_get_d(r26183, MPFR_RNDN);
}