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

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

double f_if(float K, float m, float n, float M, float l) {
        float r9454 = K;
        float r9455 = m;
        float r9456 = n;
        float r9457 = r9455 + r9456;
        float r9458 = r9454 * r9457;
        float r9459 = 2;
        float r9460 = r9458 / r9459;
        float r9461 = M;
        float r9462 = r9460 - r9461;
        float r9463 = cos(r9462);
        float r9464 = r9457 / r9459;
        float r9465 = r9464 - r9461;
        float r9466 = pow(r9465, r9459);
        float r9467 = -r9466;
        float r9468 = l;
        float r9469 = r9455 - r9456;
        float r9470 = fabs(r9469);
        float r9471 = r9468 - r9470;
        float r9472 = r9467 - r9471;
        float r9473 = exp(r9472);
        float r9474 = r9463 * r9473;
        return r9474;
}

double f_id(double K, double m, double n, double M, double l) {
        double r9475 = K;
        double r9476 = m;
        double r9477 = n;
        double r9478 = r9476 + r9477;
        double r9479 = r9475 * r9478;
        double r9480 = 2;
        double r9481 = r9479 / r9480;
        double r9482 = M;
        double r9483 = r9481 - r9482;
        double r9484 = cos(r9483);
        double r9485 = r9478 / r9480;
        double r9486 = r9485 - r9482;
        double r9487 = pow(r9486, r9480);
        double r9488 = -r9487;
        double r9489 = l;
        double r9490 = r9476 - r9477;
        double r9491 = fabs(r9490);
        double r9492 = r9489 - r9491;
        double r9493 = r9488 - r9492;
        double r9494 = exp(r9493);
        double r9495 = r9484 * r9494;
        return r9495;
}


double f_of(float __attribute__((unused)) K, float m, float n, float M, float l) {
        float r9496 = n;
        float r9497 = m;
        float r9498 = r9496 + r9497;
        float r9499 = 2;
        float r9500 = r9498 / r9499;
        float r9501 = M;
        float r9502 = r9500 - r9501;
        float r9503 = pow(r9502, r9499);
        float r9504 = -r9503;
        float r9505 = l;
        float r9506 = r9497 - r9496;
        float r9507 = fabs(r9506);
        float r9508 = r9505 - r9507;
        float r9509 = r9504 - r9508;
        float r9510 = exp(r9509);
        return r9510;
}

double f_od(double __attribute__((unused)) K, double m, double n, double M, double l) {
        double r9511 = n;
        double r9512 = m;
        double r9513 = r9511 + r9512;
        double r9514 = 2;
        double r9515 = r9513 / r9514;
        double r9516 = M;
        double r9517 = r9515 - r9516;
        double r9518 = pow(r9517, r9514);
        double r9519 = -r9518;
        double r9520 = l;
        double r9521 = r9512 - r9511;
        double r9522 = fabs(r9521);
        double r9523 = r9520 - r9522;
        double r9524 = r9519 - r9523;
        double r9525 = exp(r9524);
        return r9525;
}

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 r9526, r9527, r9528, r9529, r9530, r9531, r9532, r9533, r9534, r9535, r9536, r9537, r9538, r9539, r9540, r9541, r9542, r9543, r9544, r9545, r9546;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1616);
        mpfr_init(r9526);
        mpfr_init(r9527);
        mpfr_init(r9528);
        mpfr_init(r9529);
        mpfr_init(r9530);
        mpfr_init_set_str(r9531, "2", 10, MPFR_RNDN);
        mpfr_init(r9532);
        mpfr_init(r9533);
        mpfr_init(r9534);
        mpfr_init(r9535);
        mpfr_init(r9536);
        mpfr_init(r9537);
        mpfr_init(r9538);
        mpfr_init(r9539);
        mpfr_init(r9540);
        mpfr_init(r9541);
        mpfr_init(r9542);
        mpfr_init(r9543);
        mpfr_init(r9544);
        mpfr_init(r9545);
        mpfr_init(r9546);
}

double f_im(double K, double m, double n, double M, double l) {
        mpfr_set_d(r9526, K, MPFR_RNDN);
        mpfr_set_d(r9527, m, MPFR_RNDN);
        mpfr_set_d(r9528, n, MPFR_RNDN);
        mpfr_add(r9529, r9527, r9528, MPFR_RNDN);
        mpfr_mul(r9530, r9526, r9529, MPFR_RNDN);
        ;
        mpfr_div(r9532, r9530, r9531, MPFR_RNDN);
        mpfr_set_d(r9533, M, MPFR_RNDN);
        mpfr_sub(r9534, r9532, r9533, MPFR_RNDN);
        mpfr_cos(r9535, r9534, MPFR_RNDN);
        mpfr_div(r9536, r9529, r9531, MPFR_RNDN);
        mpfr_sub(r9537, r9536, r9533, MPFR_RNDN);
        mpfr_pow(r9538, r9537, r9531, MPFR_RNDN);
        mpfr_neg(r9539, r9538, MPFR_RNDN);
        mpfr_set_d(r9540, l, MPFR_RNDN);
        mpfr_sub(r9541, r9527, r9528, MPFR_RNDN);
        mpfr_abs(r9542, r9541, MPFR_RNDN);
        mpfr_sub(r9543, r9540, r9542, MPFR_RNDN);
        mpfr_sub(r9544, r9539, r9543, MPFR_RNDN);
        mpfr_exp(r9545, r9544, MPFR_RNDN);
        mpfr_mul(r9546, r9535, r9545, MPFR_RNDN);
        return mpfr_get_d(r9546, MPFR_RNDN);
}

static mpfr_t r9547, r9548, r9549, r9550, r9551, r9552, r9553, r9554, r9555, r9556, r9557, r9558, r9559, r9560, r9561;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1616);
        mpfr_init(r9547);
        mpfr_init(r9548);
        mpfr_init(r9549);
        mpfr_init_set_str(r9550, "2", 10, MPFR_RNDN);
        mpfr_init(r9551);
        mpfr_init(r9552);
        mpfr_init(r9553);
        mpfr_init(r9554);
        mpfr_init(r9555);
        mpfr_init(r9556);
        mpfr_init(r9557);
        mpfr_init(r9558);
        mpfr_init(r9559);
        mpfr_init(r9560);
        mpfr_init(r9561);
}

double f_fm(double __attribute__((unused)) K, double m, double n, double M, double l) {
        mpfr_set_d(r9547, n, MPFR_RNDN);
        mpfr_set_d(r9548, m, MPFR_RNDN);
        mpfr_add(r9549, r9547, r9548, MPFR_RNDN);
        ;
        mpfr_div(r9551, r9549, r9550, MPFR_RNDN);
        mpfr_set_d(r9552, M, MPFR_RNDN);
        mpfr_sub(r9553, r9551, r9552, MPFR_RNDN);
        mpfr_pow(r9554, r9553, r9550, MPFR_RNDN);
        mpfr_neg(r9555, r9554, MPFR_RNDN);
        mpfr_set_d(r9556, l, MPFR_RNDN);
        mpfr_sub(r9557, r9548, r9547, MPFR_RNDN);
        mpfr_abs(r9558, r9557, MPFR_RNDN);
        mpfr_sub(r9559, r9556, r9558, MPFR_RNDN);
        mpfr_sub(r9560, r9555, r9559, MPFR_RNDN);
        mpfr_exp(r9561, r9560, MPFR_RNDN);
        return mpfr_get_d(r9561, MPFR_RNDN);
}

static mpfr_t r9562, r9563, r9564, r9565, r9566, r9567, r9568, r9569, r9570, r9571, r9572, r9573, r9574, r9575, r9576;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1616);
        mpfr_init(r9562);
        mpfr_init(r9563);
        mpfr_init(r9564);
        mpfr_init_set_str(r9565, "2", 10, MPFR_RNDN);
        mpfr_init(r9566);
        mpfr_init(r9567);
        mpfr_init(r9568);
        mpfr_init(r9569);
        mpfr_init(r9570);
        mpfr_init(r9571);
        mpfr_init(r9572);
        mpfr_init(r9573);
        mpfr_init(r9574);
        mpfr_init(r9575);
        mpfr_init(r9576);
}

double f_dm(double __attribute__((unused)) K, double m, double n, double M, double l) {
        mpfr_set_d(r9562, n, MPFR_RNDN);
        mpfr_set_d(r9563, m, MPFR_RNDN);
        mpfr_add(r9564, r9562, r9563, MPFR_RNDN);
        ;
        mpfr_div(r9566, r9564, r9565, MPFR_RNDN);
        mpfr_set_d(r9567, M, MPFR_RNDN);
        mpfr_sub(r9568, r9566, r9567, MPFR_RNDN);
        mpfr_pow(r9569, r9568, r9565, MPFR_RNDN);
        mpfr_neg(r9570, r9569, MPFR_RNDN);
        mpfr_set_d(r9571, l, MPFR_RNDN);
        mpfr_sub(r9572, r9563, r9562, MPFR_RNDN);
        mpfr_abs(r9573, r9572, MPFR_RNDN);
        mpfr_sub(r9574, r9571, r9573, MPFR_RNDN);
        mpfr_sub(r9575, r9570, r9574, MPFR_RNDN);
        mpfr_exp(r9576, r9575, MPFR_RNDN);
        return mpfr_get_d(r9576, MPFR_RNDN);
}