#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 r22643 = J;
        float r22644 = l;
        float r22645 = exp(r22644);
        float r22646 = -r22644;
        float r22647 = exp(r22646);
        float r22648 = r22645 - r22647;
        float r22649 = r22643 * r22648;
        float r22650 = K;
        float r22651 = 2;
        float r22652 = r22650 / r22651;
        float r22653 = cos(r22652);
        float r22654 = r22649 * r22653;
        float r22655 = U;
        float r22656 = r22654 + r22655;
        return r22656;
}

double f_id(double J, double l, double K, double U) {
        double r22657 = J;
        double r22658 = l;
        double r22659 = exp(r22658);
        double r22660 = -r22658;
        double r22661 = exp(r22660);
        double r22662 = r22659 - r22661;
        double r22663 = r22657 * r22662;
        double r22664 = K;
        double r22665 = 2;
        double r22666 = r22664 / r22665;
        double r22667 = cos(r22666);
        double r22668 = r22663 * r22667;
        double r22669 = U;
        double r22670 = r22668 + r22669;
        return r22670;
}


double f_of(float J, float l, float K, float U) {
        float r22671 = J;
        float r22672 = 1/3;
        float r22673 = l;
        float r22674 = 3;
        float r22675 = pow(r22673, r22674);
        float r22676 = r22672 * r22675;
        float r22677 = 2;
        float r22678 = r22677 * r22673;
        float r22679 = 1/60;
        float r22680 = 5;
        float r22681 = pow(r22673, r22680);
        float r22682 = r22679 * r22681;
        float r22683 = r22678 + r22682;
        float r22684 = r22676 + r22683;
        float r22685 = r22671 * r22684;
        float r22686 = K;
        float r22687 = r22686 / r22677;
        float r22688 = cos(r22687);
        float r22689 = r22685 * r22688;
        float r22690 = U;
        float r22691 = r22689 + r22690;
        return r22691;
}

double f_od(double J, double l, double K, double U) {
        double r22692 = J;
        double r22693 = 1/3;
        double r22694 = l;
        double r22695 = 3;
        double r22696 = pow(r22694, r22695);
        double r22697 = r22693 * r22696;
        double r22698 = 2;
        double r22699 = r22698 * r22694;
        double r22700 = 1/60;
        double r22701 = 5;
        double r22702 = pow(r22694, r22701);
        double r22703 = r22700 * r22702;
        double r22704 = r22699 + r22703;
        double r22705 = r22697 + r22704;
        double r22706 = r22692 * r22705;
        double r22707 = K;
        double r22708 = r22707 / r22698;
        double r22709 = cos(r22708);
        double r22710 = r22706 * r22709;
        double r22711 = U;
        double r22712 = r22710 + r22711;
        return r22712;
}

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 r22713, r22714, r22715, r22716, r22717, r22718, r22719, r22720, r22721, r22722, r22723, r22724, r22725, r22726;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r22713);
        mpfr_init(r22714);
        mpfr_init(r22715);
        mpfr_init(r22716);
        mpfr_init(r22717);
        mpfr_init(r22718);
        mpfr_init(r22719);
        mpfr_init(r22720);
        mpfr_init_set_str(r22721, "2", 10, MPFR_RNDN);
        mpfr_init(r22722);
        mpfr_init(r22723);
        mpfr_init(r22724);
        mpfr_init(r22725);
        mpfr_init(r22726);
}

double f_im(double J, double l, double K, double U) {
        mpfr_set_d(r22713, J, MPFR_RNDN);
        mpfr_set_d(r22714, l, MPFR_RNDN);
        mpfr_exp(r22715, r22714, MPFR_RNDN);
        mpfr_neg(r22716, r22714, MPFR_RNDN);
        mpfr_exp(r22717, r22716, MPFR_RNDN);
        mpfr_sub(r22718, r22715, r22717, MPFR_RNDN);
        mpfr_mul(r22719, r22713, r22718, MPFR_RNDN);
        mpfr_set_d(r22720, K, MPFR_RNDN);
        ;
        mpfr_div(r22722, r22720, r22721, MPFR_RNDN);
        mpfr_cos(r22723, r22722, MPFR_RNDN);
        mpfr_mul(r22724, r22719, r22723, MPFR_RNDN);
        mpfr_set_d(r22725, U, MPFR_RNDN);
        mpfr_add(r22726, r22724, r22725, MPFR_RNDN);
        return mpfr_get_d(r22726, MPFR_RNDN);
}

static mpfr_t r22727, r22728, r22729, r22730, r22731, r22732, r22733, r22734, r22735, r22736, r22737, r22738, r22739, r22740, r22741, r22742, r22743, r22744, r22745, r22746, r22747;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r22727);
        mpfr_init_set_str(r22728, "1/3", 10, MPFR_RNDN);
        mpfr_init(r22729);
        mpfr_init_set_str(r22730, "3", 10, MPFR_RNDN);
        mpfr_init(r22731);
        mpfr_init(r22732);
        mpfr_init_set_str(r22733, "2", 10, MPFR_RNDN);
        mpfr_init(r22734);
        mpfr_init_set_str(r22735, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r22736, "5", 10, MPFR_RNDN);
        mpfr_init(r22737);
        mpfr_init(r22738);
        mpfr_init(r22739);
        mpfr_init(r22740);
        mpfr_init(r22741);
        mpfr_init(r22742);
        mpfr_init(r22743);
        mpfr_init(r22744);
        mpfr_init(r22745);
        mpfr_init(r22746);
        mpfr_init(r22747);
}

double f_fm(double J, double l, double K, double U) {
        mpfr_set_d(r22727, J, MPFR_RNDN);
        ;
        mpfr_set_d(r22729, l, MPFR_RNDN);
        ;
        mpfr_pow(r22731, r22729, r22730, MPFR_RNDN);
        mpfr_mul(r22732, r22728, r22731, MPFR_RNDN);
        ;
        mpfr_mul(r22734, r22733, r22729, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22737, r22729, r22736, MPFR_RNDN);
        mpfr_mul(r22738, r22735, r22737, MPFR_RNDN);
        mpfr_add(r22739, r22734, r22738, MPFR_RNDN);
        mpfr_add(r22740, r22732, r22739, MPFR_RNDN);
        mpfr_mul(r22741, r22727, r22740, MPFR_RNDN);
        mpfr_set_d(r22742, K, MPFR_RNDN);
        mpfr_div(r22743, r22742, r22733, MPFR_RNDN);
        mpfr_cos(r22744, r22743, MPFR_RNDN);
        mpfr_mul(r22745, r22741, r22744, MPFR_RNDN);
        mpfr_set_d(r22746, U, MPFR_RNDN);
        mpfr_add(r22747, r22745, r22746, MPFR_RNDN);
        return mpfr_get_d(r22747, MPFR_RNDN);
}

static mpfr_t r22748, r22749, r22750, r22751, r22752, r22753, r22754, r22755, r22756, r22757, r22758, r22759, r22760, r22761, r22762, r22763, r22764, r22765, r22766, r22767, r22768;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r22748);
        mpfr_init_set_str(r22749, "1/3", 10, MPFR_RNDN);
        mpfr_init(r22750);
        mpfr_init_set_str(r22751, "3", 10, MPFR_RNDN);
        mpfr_init(r22752);
        mpfr_init(r22753);
        mpfr_init_set_str(r22754, "2", 10, MPFR_RNDN);
        mpfr_init(r22755);
        mpfr_init_set_str(r22756, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r22757, "5", 10, MPFR_RNDN);
        mpfr_init(r22758);
        mpfr_init(r22759);
        mpfr_init(r22760);
        mpfr_init(r22761);
        mpfr_init(r22762);
        mpfr_init(r22763);
        mpfr_init(r22764);
        mpfr_init(r22765);
        mpfr_init(r22766);
        mpfr_init(r22767);
        mpfr_init(r22768);
}

double f_dm(double J, double l, double K, double U) {
        mpfr_set_d(r22748, J, MPFR_RNDN);
        ;
        mpfr_set_d(r22750, l, MPFR_RNDN);
        ;
        mpfr_pow(r22752, r22750, r22751, MPFR_RNDN);
        mpfr_mul(r22753, r22749, r22752, MPFR_RNDN);
        ;
        mpfr_mul(r22755, r22754, r22750, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r22758, r22750, r22757, MPFR_RNDN);
        mpfr_mul(r22759, r22756, r22758, MPFR_RNDN);
        mpfr_add(r22760, r22755, r22759, MPFR_RNDN);
        mpfr_add(r22761, r22753, r22760, MPFR_RNDN);
        mpfr_mul(r22762, r22748, r22761, MPFR_RNDN);
        mpfr_set_d(r22763, K, MPFR_RNDN);
        mpfr_div(r22764, r22763, r22754, MPFR_RNDN);
        mpfr_cos(r22765, r22764, MPFR_RNDN);
        mpfr_mul(r22766, r22762, r22765, MPFR_RNDN);
        mpfr_set_d(r22767, U, MPFR_RNDN);
        mpfr_add(r22768, r22766, r22767, MPFR_RNDN);
        return mpfr_get_d(r22768, MPFR_RNDN);
}