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

char *name = "invcot (example 3.9)";

double f_if(float x) {
        float r26514 = 1;
        float r26515 = x;
        float r26516 = r26514 / r26515;
        float r26517 = tan(r26515);
        float r26518 = r26514 / r26517;
        float r26519 = r26516 - r26518;
        return r26519;
}

double f_id(double x) {
        double r26520 = 1;
        double r26521 = x;
        double r26522 = r26520 / r26521;
        double r26523 = tan(r26521);
        double r26524 = r26520 / r26523;
        double r26525 = r26522 - r26524;
        return r26525;
}


double f_of(float x) {
        float r26526 = 1/45;
        float r26527 = x;
        float r26528 = 3;
        float r26529 = pow(r26527, r26528);
        float r26530 = r26526 * r26529;
        float r26531 = 2/945;
        float r26532 = 5;
        float r26533 = pow(r26527, r26532);
        float r26534 = r26531 * r26533;
        float r26535 = 1/3;
        float r26536 = r26535 * r26527;
        float r26537 = r26534 + r26536;
        float r26538 = r26530 + r26537;
        return r26538;
}

double f_od(double x) {
        double r26539 = 1/45;
        double r26540 = x;
        double r26541 = 3;
        double r26542 = pow(r26540, r26541);
        double r26543 = r26539 * r26542;
        double r26544 = 2/945;
        double r26545 = 5;
        double r26546 = pow(r26540, r26545);
        double r26547 = r26544 * r26546;
        double r26548 = 1/3;
        double r26549 = r26548 * r26540;
        double r26550 = r26547 + r26549;
        double r26551 = r26543 + r26550;
        return r26551;
}

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 r26552, r26553, r26554, r26555, r26556, r26557;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r26552, "1", 10, MPFR_RNDN);
        mpfr_init(r26553);
        mpfr_init(r26554);
        mpfr_init(r26555);
        mpfr_init(r26556);
        mpfr_init(r26557);
}

double f_im(double x) {
        ;
        mpfr_set_d(r26553, x, MPFR_RNDN);
        mpfr_div(r26554, r26552, r26553, MPFR_RNDN);
        mpfr_tan(r26555, r26553, MPFR_RNDN);
        mpfr_div(r26556, r26552, r26555, MPFR_RNDN);
        mpfr_sub(r26557, r26554, r26556, MPFR_RNDN);
        return mpfr_get_d(r26557, MPFR_RNDN);
}

static mpfr_t r26558, r26559, r26560, r26561, r26562, r26563, r26564, r26565, r26566, r26567, r26568, r26569, r26570;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r26558, "1/45", 10, MPFR_RNDN);
        mpfr_init(r26559);
        mpfr_init_set_str(r26560, "3", 10, MPFR_RNDN);
        mpfr_init(r26561);
        mpfr_init(r26562);
        mpfr_init_set_str(r26563, "2/945", 10, MPFR_RNDN);
        mpfr_init_set_str(r26564, "5", 10, MPFR_RNDN);
        mpfr_init(r26565);
        mpfr_init(r26566);
        mpfr_init_set_str(r26567, "1/3", 10, MPFR_RNDN);
        mpfr_init(r26568);
        mpfr_init(r26569);
        mpfr_init(r26570);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r26559, x, MPFR_RNDN);
        ;
        mpfr_pow(r26561, r26559, r26560, MPFR_RNDN);
        mpfr_mul(r26562, r26558, r26561, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r26565, r26559, r26564, MPFR_RNDN);
        mpfr_mul(r26566, r26563, r26565, MPFR_RNDN);
        ;
        mpfr_mul(r26568, r26567, r26559, MPFR_RNDN);
        mpfr_add(r26569, r26566, r26568, MPFR_RNDN);
        mpfr_add(r26570, r26562, r26569, MPFR_RNDN);
        return mpfr_get_d(r26570, MPFR_RNDN);
}

static mpfr_t r26571, r26572, r26573, r26574, r26575, r26576, r26577, r26578, r26579, r26580, r26581, r26582, r26583;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r26571, "1/45", 10, MPFR_RNDN);
        mpfr_init(r26572);
        mpfr_init_set_str(r26573, "3", 10, MPFR_RNDN);
        mpfr_init(r26574);
        mpfr_init(r26575);
        mpfr_init_set_str(r26576, "2/945", 10, MPFR_RNDN);
        mpfr_init_set_str(r26577, "5", 10, MPFR_RNDN);
        mpfr_init(r26578);
        mpfr_init(r26579);
        mpfr_init_set_str(r26580, "1/3", 10, MPFR_RNDN);
        mpfr_init(r26581);
        mpfr_init(r26582);
        mpfr_init(r26583);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r26572, x, MPFR_RNDN);
        ;
        mpfr_pow(r26574, r26572, r26573, MPFR_RNDN);
        mpfr_mul(r26575, r26571, r26574, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r26578, r26572, r26577, MPFR_RNDN);
        mpfr_mul(r26579, r26576, r26578, MPFR_RNDN);
        ;
        mpfr_mul(r26581, r26580, r26572, MPFR_RNDN);
        mpfr_add(r26582, r26579, r26581, MPFR_RNDN);
        mpfr_add(r26583, r26575, r26582, MPFR_RNDN);
        return mpfr_get_d(r26583, MPFR_RNDN);
}