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

char *name = "Migdal et al, Equation (51)";

double f_if(float k, float n) {
        float r9091 = 1.0f;
        float r9092 = k;
        float r9093 = sqrt(r9092);
        float r9094 = r9091 / r9093;
        float r9095 = 2.0f;
        float r9096 = atan2(1.0, 0.0);
        float r9097 = r9095 * r9096;
        float r9098 = n;
        float r9099 = r9097 * r9098;
        float r9100 = r9091 - r9092;
        float r9101 = r9100 / r9095;
        float r9102 = pow(r9099, r9101);
        float r9103 = r9094 * r9102;
        return r9103;
}

double f_id(double k, double n) {
        double r9104 = 1.0;
        double r9105 = k;
        double r9106 = sqrt(r9105);
        double r9107 = r9104 / r9106;
        double r9108 = 2.0;
        double r9109 = atan2(1.0, 0.0);
        double r9110 = r9108 * r9109;
        double r9111 = n;
        double r9112 = r9110 * r9111;
        double r9113 = r9104 - r9105;
        double r9114 = r9113 / r9108;
        double r9115 = pow(r9112, r9114);
        double r9116 = r9107 * r9115;
        return r9116;
}


double f_of(float k, float n) {
        float r9117 = n;
        float r9118 = 2.0f;
        float r9119 = r9117 * r9118;
        float r9120 = atan2(1.0, 0.0);
        float r9121 = r9119 * r9120;
        float r9122 = 1.0f;
        float r9123 = r9122 / r9118;
        float r9124 = pow(r9121, r9123);
        float r9125 = k;
        float r9126 = r9125 / r9118;
        float r9127 = pow(r9119, r9126);
        float r9128 = pow(r9120, r9126);
        float r9129 = r9127 * r9128;
        float r9130 = r9124 / r9129;
        float r9131 = sqrt(r9125);
        float r9132 = r9130 / r9131;
        return r9132;
}

double f_od(double k, double n) {
        double r9133 = n;
        double r9134 = 2.0;
        double r9135 = r9133 * r9134;
        double r9136 = atan2(1.0, 0.0);
        double r9137 = r9135 * r9136;
        double r9138 = 1.0;
        double r9139 = r9138 / r9134;
        double r9140 = pow(r9137, r9139);
        double r9141 = k;
        double r9142 = r9141 / r9134;
        double r9143 = pow(r9135, r9142);
        double r9144 = pow(r9136, r9142);
        double r9145 = r9143 * r9144;
        double r9146 = r9140 / r9145;
        double r9147 = sqrt(r9141);
        double r9148 = r9146 / r9147;
        return r9148;
}

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 r9149, r9150, r9151, r9152, r9153, r9154, r9155, r9156, r9157, r9158, r9159, r9160, r9161;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init_set_str(r9149, "1", 10, MPFR_RNDN);
        mpfr_init(r9150);
        mpfr_init(r9151);
        mpfr_init(r9152);
        mpfr_init_set_str(r9153, "2", 10, MPFR_RNDN);
        mpfr_init(r9154);
        mpfr_init(r9155);
        mpfr_init(r9156);
        mpfr_init(r9157);
        mpfr_init(r9158);
        mpfr_init(r9159);
        mpfr_init(r9160);
        mpfr_init(r9161);
}

double f_im(double k, double n) {
        ;
        mpfr_set_d(r9150, k, MPFR_RNDN);
        mpfr_sqrt(r9151, r9150, MPFR_RNDN);
        mpfr_div(r9152, r9149, r9151, MPFR_RNDN);
        ;
        mpfr_const_pi(r9154, MPFR_RNDN);
        mpfr_mul(r9155, r9153, r9154, MPFR_RNDN);
        mpfr_set_d(r9156, n, MPFR_RNDN);
        mpfr_mul(r9157, r9155, r9156, MPFR_RNDN);
        mpfr_sub(r9158, r9149, r9150, MPFR_RNDN);
        mpfr_div(r9159, r9158, r9153, MPFR_RNDN);
        mpfr_pow(r9160, r9157, r9159, MPFR_RNDN);
        mpfr_mul(r9161, r9152, r9160, MPFR_RNDN);
        return mpfr_get_d(r9161, MPFR_RNDN);
}

static mpfr_t r9162, r9163, r9164, r9165, r9166, r9167, r9168, r9169, r9170, r9171, r9172, r9173, r9174, r9175, r9176, r9177;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9162);
        mpfr_init_set_str(r9163, "2", 10, MPFR_RNDN);
        mpfr_init(r9164);
        mpfr_init(r9165);
        mpfr_init(r9166);
        mpfr_init_set_str(r9167, "1", 10, MPFR_RNDN);
        mpfr_init(r9168);
        mpfr_init(r9169);
        mpfr_init(r9170);
        mpfr_init(r9171);
        mpfr_init(r9172);
        mpfr_init(r9173);
        mpfr_init(r9174);
        mpfr_init(r9175);
        mpfr_init(r9176);
        mpfr_init(r9177);
}

double f_fm(double k, double n) {
        mpfr_set_d(r9162, n, MPFR_RNDN);
        ;
        mpfr_mul(r9164, r9162, r9163, MPFR_RNDN);
        mpfr_const_pi(r9165, MPFR_RNDN);
        mpfr_mul(r9166, r9164, r9165, MPFR_RNDN);
        ;
        mpfr_div(r9168, r9167, r9163, MPFR_RNDN);
        mpfr_pow(r9169, r9166, r9168, MPFR_RNDN);
        mpfr_set_d(r9170, k, MPFR_RNDN);
        mpfr_div(r9171, r9170, r9163, MPFR_RNDN);
        mpfr_pow(r9172, r9164, r9171, MPFR_RNDN);
        mpfr_pow(r9173, r9165, r9171, MPFR_RNDN);
        mpfr_mul(r9174, r9172, r9173, MPFR_RNDN);
        mpfr_div(r9175, r9169, r9174, MPFR_RNDN);
        mpfr_sqrt(r9176, r9170, MPFR_RNDN);
        mpfr_div(r9177, r9175, r9176, MPFR_RNDN);
        return mpfr_get_d(r9177, MPFR_RNDN);
}

static mpfr_t r9178, r9179, r9180, r9181, r9182, r9183, r9184, r9185, r9186, r9187, r9188, r9189, r9190, r9191, r9192, r9193;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9178);
        mpfr_init_set_str(r9179, "2", 10, MPFR_RNDN);
        mpfr_init(r9180);
        mpfr_init(r9181);
        mpfr_init(r9182);
        mpfr_init_set_str(r9183, "1", 10, MPFR_RNDN);
        mpfr_init(r9184);
        mpfr_init(r9185);
        mpfr_init(r9186);
        mpfr_init(r9187);
        mpfr_init(r9188);
        mpfr_init(r9189);
        mpfr_init(r9190);
        mpfr_init(r9191);
        mpfr_init(r9192);
        mpfr_init(r9193);
}

double f_dm(double k, double n) {
        mpfr_set_d(r9178, n, MPFR_RNDN);
        ;
        mpfr_mul(r9180, r9178, r9179, MPFR_RNDN);
        mpfr_const_pi(r9181, MPFR_RNDN);
        mpfr_mul(r9182, r9180, r9181, MPFR_RNDN);
        ;
        mpfr_div(r9184, r9183, r9179, MPFR_RNDN);
        mpfr_pow(r9185, r9182, r9184, MPFR_RNDN);
        mpfr_set_d(r9186, k, MPFR_RNDN);
        mpfr_div(r9187, r9186, r9179, MPFR_RNDN);
        mpfr_pow(r9188, r9180, r9187, MPFR_RNDN);
        mpfr_pow(r9189, r9181, r9187, MPFR_RNDN);
        mpfr_mul(r9190, r9188, r9189, MPFR_RNDN);
        mpfr_div(r9191, r9185, r9190, MPFR_RNDN);
        mpfr_sqrt(r9192, r9186, MPFR_RNDN);
        mpfr_div(r9193, r9191, r9192, MPFR_RNDN);
        return mpfr_get_d(r9193, MPFR_RNDN);
}