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

char *name = "Maksimov and Kolovsky, Equation (4)";

double f_if(float J, float l, float K, float U) {
        float r9546 = J;
        float r9547 = l;
        float r9548 = exp(r9547);
        float r9549 = -r9547;
        float r9550 = exp(r9549);
        float r9551 = r9548 - r9550;
        float r9552 = r9546 * r9551;
        float r9553 = K;
        float r9554 = 2;
        float r9555 = r9553 / r9554;
        float r9556 = cos(r9555);
        float r9557 = r9552 * r9556;
        float r9558 = U;
        float r9559 = r9557 + r9558;
        return r9559;
}

double f_id(double J, double l, double K, double U) {
        double r9560 = J;
        double r9561 = l;
        double r9562 = exp(r9561);
        double r9563 = -r9561;
        double r9564 = exp(r9563);
        double r9565 = r9562 - r9564;
        double r9566 = r9560 * r9565;
        double r9567 = K;
        double r9568 = 2;
        double r9569 = r9567 / r9568;
        double r9570 = cos(r9569);
        double r9571 = r9566 * r9570;
        double r9572 = U;
        double r9573 = r9571 + r9572;
        return r9573;
}


double f_of(float J, float l, float K, float U) {
        float r9574 = J;
        float r9575 = 1/3;
        float r9576 = l;
        float r9577 = 3;
        float r9578 = pow(r9576, r9577);
        float r9579 = r9575 * r9578;
        float r9580 = r9574 * r9579;
        float r9581 = 1/60;
        float r9582 = 5;
        float r9583 = pow(r9576, r9582);
        float r9584 = r9581 * r9583;
        float r9585 = 2;
        float r9586 = r9585 * r9576;
        float r9587 = r9584 + r9586;
        float r9588 = r9574 * r9587;
        float r9589 = r9580 + r9588;
        float r9590 = K;
        float r9591 = r9590 / r9585;
        float r9592 = cos(r9591);
        float r9593 = r9589 * r9592;
        float r9594 = U;
        float r9595 = r9593 + r9594;
        return r9595;
}

double f_od(double J, double l, double K, double U) {
        double r9596 = J;
        double r9597 = 1/3;
        double r9598 = l;
        double r9599 = 3;
        double r9600 = pow(r9598, r9599);
        double r9601 = r9597 * r9600;
        double r9602 = r9596 * r9601;
        double r9603 = 1/60;
        double r9604 = 5;
        double r9605 = pow(r9598, r9604);
        double r9606 = r9603 * r9605;
        double r9607 = 2;
        double r9608 = r9607 * r9598;
        double r9609 = r9606 + r9608;
        double r9610 = r9596 * r9609;
        double r9611 = r9602 + r9610;
        double r9612 = K;
        double r9613 = r9612 / r9607;
        double r9614 = cos(r9613);
        double r9615 = r9611 * r9614;
        double r9616 = U;
        double r9617 = r9615 + r9616;
        return r9617;
}

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 r9618, r9619, r9620, r9621, r9622, r9623, r9624, r9625, r9626, r9627, r9628, r9629, r9630, r9631;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9618);
        mpfr_init(r9619);
        mpfr_init(r9620);
        mpfr_init(r9621);
        mpfr_init(r9622);
        mpfr_init(r9623);
        mpfr_init(r9624);
        mpfr_init(r9625);
        mpfr_init_set_str(r9626, "2", 10, MPFR_RNDN);
        mpfr_init(r9627);
        mpfr_init(r9628);
        mpfr_init(r9629);
        mpfr_init(r9630);
        mpfr_init(r9631);
}

double f_im(double J, double l, double K, double U) {
        mpfr_set_d(r9618, J, MPFR_RNDN);
        mpfr_set_d(r9619, l, MPFR_RNDN);
        mpfr_exp(r9620, r9619, MPFR_RNDN);
        mpfr_neg(r9621, r9619, MPFR_RNDN);
        mpfr_exp(r9622, r9621, MPFR_RNDN);
        mpfr_sub(r9623, r9620, r9622, MPFR_RNDN);
        mpfr_mul(r9624, r9618, r9623, MPFR_RNDN);
        mpfr_set_d(r9625, K, MPFR_RNDN);
        ;
        mpfr_div(r9627, r9625, r9626, MPFR_RNDN);
        mpfr_cos(r9628, r9627, MPFR_RNDN);
        mpfr_mul(r9629, r9624, r9628, MPFR_RNDN);
        mpfr_set_d(r9630, U, MPFR_RNDN);
        mpfr_add(r9631, r9629, r9630, MPFR_RNDN);
        return mpfr_get_d(r9631, MPFR_RNDN);
}

static mpfr_t r9632, r9633, r9634, r9635, r9636, r9637, r9638, r9639, r9640, r9641, r9642, r9643, r9644, r9645, r9646, r9647, r9648, r9649, r9650, r9651, r9652, r9653;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9632);
        mpfr_init_set_str(r9633, "1/3", 10, MPFR_RNDN);
        mpfr_init(r9634);
        mpfr_init_set_str(r9635, "3", 10, MPFR_RNDN);
        mpfr_init(r9636);
        mpfr_init(r9637);
        mpfr_init(r9638);
        mpfr_init_set_str(r9639, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r9640, "5", 10, MPFR_RNDN);
        mpfr_init(r9641);
        mpfr_init(r9642);
        mpfr_init_set_str(r9643, "2", 10, MPFR_RNDN);
        mpfr_init(r9644);
        mpfr_init(r9645);
        mpfr_init(r9646);
        mpfr_init(r9647);
        mpfr_init(r9648);
        mpfr_init(r9649);
        mpfr_init(r9650);
        mpfr_init(r9651);
        mpfr_init(r9652);
        mpfr_init(r9653);
}

double f_fm(double J, double l, double K, double U) {
        mpfr_set_d(r9632, J, MPFR_RNDN);
        ;
        mpfr_set_d(r9634, l, MPFR_RNDN);
        ;
        mpfr_pow(r9636, r9634, r9635, MPFR_RNDN);
        mpfr_mul(r9637, r9633, r9636, MPFR_RNDN);
        mpfr_mul(r9638, r9632, r9637, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r9641, r9634, r9640, MPFR_RNDN);
        mpfr_mul(r9642, r9639, r9641, MPFR_RNDN);
        ;
        mpfr_mul(r9644, r9643, r9634, MPFR_RNDN);
        mpfr_add(r9645, r9642, r9644, MPFR_RNDN);
        mpfr_mul(r9646, r9632, r9645, MPFR_RNDN);
        mpfr_add(r9647, r9638, r9646, MPFR_RNDN);
        mpfr_set_d(r9648, K, MPFR_RNDN);
        mpfr_div(r9649, r9648, r9643, MPFR_RNDN);
        mpfr_cos(r9650, r9649, MPFR_RNDN);
        mpfr_mul(r9651, r9647, r9650, MPFR_RNDN);
        mpfr_set_d(r9652, U, MPFR_RNDN);
        mpfr_add(r9653, r9651, r9652, MPFR_RNDN);
        return mpfr_get_d(r9653, MPFR_RNDN);
}

static mpfr_t r9654, r9655, r9656, r9657, r9658, r9659, r9660, r9661, r9662, r9663, r9664, r9665, r9666, r9667, r9668, r9669, r9670, r9671, r9672, r9673, r9674, r9675;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9654);
        mpfr_init_set_str(r9655, "1/3", 10, MPFR_RNDN);
        mpfr_init(r9656);
        mpfr_init_set_str(r9657, "3", 10, MPFR_RNDN);
        mpfr_init(r9658);
        mpfr_init(r9659);
        mpfr_init(r9660);
        mpfr_init_set_str(r9661, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r9662, "5", 10, MPFR_RNDN);
        mpfr_init(r9663);
        mpfr_init(r9664);
        mpfr_init_set_str(r9665, "2", 10, MPFR_RNDN);
        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);
}

double f_dm(double J, double l, double K, double U) {
        mpfr_set_d(r9654, J, MPFR_RNDN);
        ;
        mpfr_set_d(r9656, l, MPFR_RNDN);
        ;
        mpfr_pow(r9658, r9656, r9657, MPFR_RNDN);
        mpfr_mul(r9659, r9655, r9658, MPFR_RNDN);
        mpfr_mul(r9660, r9654, r9659, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r9663, r9656, r9662, MPFR_RNDN);
        mpfr_mul(r9664, r9661, r9663, MPFR_RNDN);
        ;
        mpfr_mul(r9666, r9665, r9656, MPFR_RNDN);
        mpfr_add(r9667, r9664, r9666, MPFR_RNDN);
        mpfr_mul(r9668, r9654, r9667, MPFR_RNDN);
        mpfr_add(r9669, r9660, r9668, MPFR_RNDN);
        mpfr_set_d(r9670, K, MPFR_RNDN);
        mpfr_div(r9671, r9670, r9665, MPFR_RNDN);
        mpfr_cos(r9672, r9671, MPFR_RNDN);
        mpfr_mul(r9673, r9669, r9672, MPFR_RNDN);
        mpfr_set_d(r9674, U, MPFR_RNDN);
        mpfr_add(r9675, r9673, r9674, MPFR_RNDN);
        return mpfr_get_d(r9675, MPFR_RNDN);
}