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

char *name = "VandenBroeck and Keller, Equation (24)";

double f_if(float B, float x) {
        float r23490 = x;
        float r23491 = 1;
        float r23492 = B;
        float r23493 = tan(r23492);
        float r23494 = r23491 / r23493;
        float r23495 = r23490 * r23494;
        float r23496 = -r23495;
        float r23497 = sin(r23492);
        float r23498 = r23491 / r23497;
        float r23499 = r23496 + r23498;
        return r23499;
}

double f_id(double B, double x) {
        double r23500 = x;
        double r23501 = 1;
        double r23502 = B;
        double r23503 = tan(r23502);
        double r23504 = r23501 / r23503;
        double r23505 = r23500 * r23504;
        double r23506 = -r23505;
        double r23507 = sin(r23502);
        double r23508 = r23501 / r23507;
        double r23509 = r23506 + r23508;
        return r23509;
}


double f_of(float B, float x) {
        float r23510 = 1;
        float r23511 = B;
        float r23512 = sin(r23511);
        float r23513 = r23510 / r23512;
        float r23514 = x;
        float r23515 = -r23514;
        float r23516 = r23515 / r23512;
        float r23517 = cos(r23511);
        float r23518 = r23516 * r23517;
        float r23519 = r23513 + r23518;
        return r23519;
}

double f_od(double B, double x) {
        double r23520 = 1;
        double r23521 = B;
        double r23522 = sin(r23521);
        double r23523 = r23520 / r23522;
        double r23524 = x;
        double r23525 = -r23524;
        double r23526 = r23525 / r23522;
        double r23527 = cos(r23521);
        double r23528 = r23526 * r23527;
        double r23529 = r23523 + r23528;
        return r23529;
}

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 r23530, r23531, r23532, r23533, r23534, r23535, r23536, r23537, r23538, r23539;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r23530);
        mpfr_init_set_str(r23531, "1", 10, MPFR_RNDN);
        mpfr_init(r23532);
        mpfr_init(r23533);
        mpfr_init(r23534);
        mpfr_init(r23535);
        mpfr_init(r23536);
        mpfr_init(r23537);
        mpfr_init(r23538);
        mpfr_init(r23539);
}

double f_im(double B, double x) {
        mpfr_set_d(r23530, x, MPFR_RNDN);
        ;
        mpfr_set_d(r23532, B, MPFR_RNDN);
        mpfr_tan(r23533, r23532, MPFR_RNDN);
        mpfr_div(r23534, r23531, r23533, MPFR_RNDN);
        mpfr_mul(r23535, r23530, r23534, MPFR_RNDN);
        mpfr_neg(r23536, r23535, MPFR_RNDN);
        mpfr_sin(r23537, r23532, MPFR_RNDN);
        mpfr_div(r23538, r23531, r23537, MPFR_RNDN);
        mpfr_add(r23539, r23536, r23538, MPFR_RNDN);
        return mpfr_get_d(r23539, MPFR_RNDN);
}

static mpfr_t r23540, r23541, r23542, r23543, r23544, r23545, r23546, r23547, r23548, r23549;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r23540, "1", 10, MPFR_RNDN);
        mpfr_init(r23541);
        mpfr_init(r23542);
        mpfr_init(r23543);
        mpfr_init(r23544);
        mpfr_init(r23545);
        mpfr_init(r23546);
        mpfr_init(r23547);
        mpfr_init(r23548);
        mpfr_init(r23549);
}

double f_fm(double B, double x) {
        ;
        mpfr_set_d(r23541, B, MPFR_RNDN);
        mpfr_sin(r23542, r23541, MPFR_RNDN);
        mpfr_div(r23543, r23540, r23542, MPFR_RNDN);
        mpfr_set_d(r23544, x, MPFR_RNDN);
        mpfr_neg(r23545, r23544, MPFR_RNDN);
        mpfr_div(r23546, r23545, r23542, MPFR_RNDN);
        mpfr_cos(r23547, r23541, MPFR_RNDN);
        mpfr_mul(r23548, r23546, r23547, MPFR_RNDN);
        mpfr_add(r23549, r23543, r23548, MPFR_RNDN);
        return mpfr_get_d(r23549, MPFR_RNDN);
}

static mpfr_t r23550, r23551, r23552, r23553, r23554, r23555, r23556, r23557, r23558, r23559;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r23550, "1", 10, MPFR_RNDN);
        mpfr_init(r23551);
        mpfr_init(r23552);
        mpfr_init(r23553);
        mpfr_init(r23554);
        mpfr_init(r23555);
        mpfr_init(r23556);
        mpfr_init(r23557);
        mpfr_init(r23558);
        mpfr_init(r23559);
}

double f_dm(double B, double x) {
        ;
        mpfr_set_d(r23551, B, MPFR_RNDN);
        mpfr_sin(r23552, r23551, MPFR_RNDN);
        mpfr_div(r23553, r23550, r23552, MPFR_RNDN);
        mpfr_set_d(r23554, x, MPFR_RNDN);
        mpfr_neg(r23555, r23554, MPFR_RNDN);
        mpfr_div(r23556, r23555, r23552, MPFR_RNDN);
        mpfr_cos(r23557, r23551, MPFR_RNDN);
        mpfr_mul(r23558, r23556, r23557, MPFR_RNDN);
        mpfr_add(r23559, r23553, r23558, MPFR_RNDN);
        return mpfr_get_d(r23559, MPFR_RNDN);
}