#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 r9324 = J;
        float r9325 = l;
        float r9326 = exp(r9325);
        float r9327 = -r9325;
        float r9328 = exp(r9327);
        float r9329 = r9326 - r9328;
        float r9330 = r9324 * r9329;
        float r9331 = K;
        float r9332 = 2;
        float r9333 = r9331 / r9332;
        float r9334 = cos(r9333);
        float r9335 = r9330 * r9334;
        float r9336 = U;
        float r9337 = r9335 + r9336;
        return r9337;
}

double f_id(double J, double l, double K, double U) {
        double r9338 = J;
        double r9339 = l;
        double r9340 = exp(r9339);
        double r9341 = -r9339;
        double r9342 = exp(r9341);
        double r9343 = r9340 - r9342;
        double r9344 = r9338 * r9343;
        double r9345 = K;
        double r9346 = 2;
        double r9347 = r9345 / r9346;
        double r9348 = cos(r9347);
        double r9349 = r9344 * r9348;
        double r9350 = U;
        double r9351 = r9349 + r9350;
        return r9351;
}


double f_of(float J, float l, float K, float U) {
        float r9352 = J;
        float r9353 = 1/3;
        float r9354 = l;
        float r9355 = 3;
        float r9356 = pow(r9354, r9355);
        float r9357 = r9353 * r9356;
        float r9358 = r9352 * r9357;
        float r9359 = 1/60;
        float r9360 = 5;
        float r9361 = pow(r9354, r9360);
        float r9362 = r9359 * r9361;
        float r9363 = 2;
        float r9364 = r9363 * r9354;
        float r9365 = r9362 + r9364;
        float r9366 = r9352 * r9365;
        float r9367 = r9358 + r9366;
        float r9368 = K;
        float r9369 = r9368 / r9363;
        float r9370 = cos(r9369);
        float r9371 = r9367 * r9370;
        float r9372 = U;
        float r9373 = r9371 + r9372;
        return r9373;
}

double f_od(double J, double l, double K, double U) {
        double r9374 = J;
        double r9375 = 1/3;
        double r9376 = l;
        double r9377 = 3;
        double r9378 = pow(r9376, r9377);
        double r9379 = r9375 * r9378;
        double r9380 = r9374 * r9379;
        double r9381 = 1/60;
        double r9382 = 5;
        double r9383 = pow(r9376, r9382);
        double r9384 = r9381 * r9383;
        double r9385 = 2;
        double r9386 = r9385 * r9376;
        double r9387 = r9384 + r9386;
        double r9388 = r9374 * r9387;
        double r9389 = r9380 + r9388;
        double r9390 = K;
        double r9391 = r9390 / r9385;
        double r9392 = cos(r9391);
        double r9393 = r9389 * r9392;
        double r9394 = U;
        double r9395 = r9393 + r9394;
        return r9395;
}

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 r9396, r9397, r9398, r9399, r9400, r9401, r9402, r9403, r9404, r9405, r9406, r9407, r9408, r9409;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9396);
        mpfr_init(r9397);
        mpfr_init(r9398);
        mpfr_init(r9399);
        mpfr_init(r9400);
        mpfr_init(r9401);
        mpfr_init(r9402);
        mpfr_init(r9403);
        mpfr_init_set_str(r9404, "2", 10, MPFR_RNDN);
        mpfr_init(r9405);
        mpfr_init(r9406);
        mpfr_init(r9407);
        mpfr_init(r9408);
        mpfr_init(r9409);
}

double f_im(double J, double l, double K, double U) {
        mpfr_set_d(r9396, J, MPFR_RNDN);
        mpfr_set_d(r9397, l, MPFR_RNDN);
        mpfr_exp(r9398, r9397, MPFR_RNDN);
        mpfr_neg(r9399, r9397, MPFR_RNDN);
        mpfr_exp(r9400, r9399, MPFR_RNDN);
        mpfr_sub(r9401, r9398, r9400, MPFR_RNDN);
        mpfr_mul(r9402, r9396, r9401, MPFR_RNDN);
        mpfr_set_d(r9403, K, MPFR_RNDN);
        ;
        mpfr_div(r9405, r9403, r9404, MPFR_RNDN);
        mpfr_cos(r9406, r9405, MPFR_RNDN);
        mpfr_mul(r9407, r9402, r9406, MPFR_RNDN);
        mpfr_set_d(r9408, U, MPFR_RNDN);
        mpfr_add(r9409, r9407, r9408, MPFR_RNDN);
        return mpfr_get_d(r9409, MPFR_RNDN);
}

static mpfr_t r9410, r9411, r9412, r9413, r9414, r9415, r9416, r9417, r9418, r9419, r9420, r9421, r9422, r9423, r9424, r9425, r9426, r9427, r9428, r9429, r9430, r9431;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9410);
        mpfr_init_set_str(r9411, "1/3", 10, MPFR_RNDN);
        mpfr_init(r9412);
        mpfr_init_set_str(r9413, "3", 10, MPFR_RNDN);
        mpfr_init(r9414);
        mpfr_init(r9415);
        mpfr_init(r9416);
        mpfr_init_set_str(r9417, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r9418, "5", 10, MPFR_RNDN);
        mpfr_init(r9419);
        mpfr_init(r9420);
        mpfr_init_set_str(r9421, "2", 10, MPFR_RNDN);
        mpfr_init(r9422);
        mpfr_init(r9423);
        mpfr_init(r9424);
        mpfr_init(r9425);
        mpfr_init(r9426);
        mpfr_init(r9427);
        mpfr_init(r9428);
        mpfr_init(r9429);
        mpfr_init(r9430);
        mpfr_init(r9431);
}

double f_fm(double J, double l, double K, double U) {
        mpfr_set_d(r9410, J, MPFR_RNDN);
        ;
        mpfr_set_d(r9412, l, MPFR_RNDN);
        ;
        mpfr_pow(r9414, r9412, r9413, MPFR_RNDN);
        mpfr_mul(r9415, r9411, r9414, MPFR_RNDN);
        mpfr_mul(r9416, r9410, r9415, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r9419, r9412, r9418, MPFR_RNDN);
        mpfr_mul(r9420, r9417, r9419, MPFR_RNDN);
        ;
        mpfr_mul(r9422, r9421, r9412, MPFR_RNDN);
        mpfr_add(r9423, r9420, r9422, MPFR_RNDN);
        mpfr_mul(r9424, r9410, r9423, MPFR_RNDN);
        mpfr_add(r9425, r9416, r9424, MPFR_RNDN);
        mpfr_set_d(r9426, K, MPFR_RNDN);
        mpfr_div(r9427, r9426, r9421, MPFR_RNDN);
        mpfr_cos(r9428, r9427, MPFR_RNDN);
        mpfr_mul(r9429, r9425, r9428, MPFR_RNDN);
        mpfr_set_d(r9430, U, MPFR_RNDN);
        mpfr_add(r9431, r9429, r9430, MPFR_RNDN);
        return mpfr_get_d(r9431, MPFR_RNDN);
}

static mpfr_t r9432, r9433, r9434, r9435, r9436, r9437, r9438, r9439, r9440, r9441, r9442, r9443, r9444, r9445, r9446, r9447, r9448, r9449, r9450, r9451, r9452, r9453;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1360);
        mpfr_init(r9432);
        mpfr_init_set_str(r9433, "1/3", 10, MPFR_RNDN);
        mpfr_init(r9434);
        mpfr_init_set_str(r9435, "3", 10, MPFR_RNDN);
        mpfr_init(r9436);
        mpfr_init(r9437);
        mpfr_init(r9438);
        mpfr_init_set_str(r9439, "1/60", 10, MPFR_RNDN);
        mpfr_init_set_str(r9440, "5", 10, MPFR_RNDN);
        mpfr_init(r9441);
        mpfr_init(r9442);
        mpfr_init_set_str(r9443, "2", 10, MPFR_RNDN);
        mpfr_init(r9444);
        mpfr_init(r9445);
        mpfr_init(r9446);
        mpfr_init(r9447);
        mpfr_init(r9448);
        mpfr_init(r9449);
        mpfr_init(r9450);
        mpfr_init(r9451);
        mpfr_init(r9452);
        mpfr_init(r9453);
}

double f_dm(double J, double l, double K, double U) {
        mpfr_set_d(r9432, J, MPFR_RNDN);
        ;
        mpfr_set_d(r9434, l, MPFR_RNDN);
        ;
        mpfr_pow(r9436, r9434, r9435, MPFR_RNDN);
        mpfr_mul(r9437, r9433, r9436, MPFR_RNDN);
        mpfr_mul(r9438, r9432, r9437, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r9441, r9434, r9440, MPFR_RNDN);
        mpfr_mul(r9442, r9439, r9441, MPFR_RNDN);
        ;
        mpfr_mul(r9444, r9443, r9434, MPFR_RNDN);
        mpfr_add(r9445, r9442, r9444, MPFR_RNDN);
        mpfr_mul(r9446, r9432, r9445, MPFR_RNDN);
        mpfr_add(r9447, r9438, r9446, MPFR_RNDN);
        mpfr_set_d(r9448, K, MPFR_RNDN);
        mpfr_div(r9449, r9448, r9443, MPFR_RNDN);
        mpfr_cos(r9450, r9449, MPFR_RNDN);
        mpfr_mul(r9451, r9447, r9450, MPFR_RNDN);
        mpfr_set_d(r9452, U, MPFR_RNDN);
        mpfr_add(r9453, r9451, r9452, MPFR_RNDN);
        return mpfr_get_d(r9453, MPFR_RNDN);
}