#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 r10396 = 1;
        float r10397 = x;
        float r10398 = r10396 / r10397;
        float r10399 = tan(r10397);
        float r10400 = r10396 / r10399;
        float r10401 = r10398 - r10400;
        return r10401;
}

double f_id(double x) {
        double r10402 = 1;
        double r10403 = x;
        double r10404 = r10402 / r10403;
        double r10405 = tan(r10403);
        double r10406 = r10402 / r10405;
        double r10407 = r10404 - r10406;
        return r10407;
}


double f_of(float x) {
        float r10408 = 1/3;
        float r10409 = x;
        float r10410 = r10408 * r10409;
        float r10411 = 1/45;
        float r10412 = 3;
        float r10413 = pow(r10409, r10412);
        float r10414 = r10411 * r10413;
        float r10415 = 2/945;
        float r10416 = 5;
        float r10417 = pow(r10409, r10416);
        float r10418 = r10415 * r10417;
        float r10419 = r10414 + r10418;
        float r10420 = r10410 + r10419;
        return r10420;
}

double f_od(double x) {
        double r10421 = 1/3;
        double r10422 = x;
        double r10423 = r10421 * r10422;
        double r10424 = 1/45;
        double r10425 = 3;
        double r10426 = pow(r10422, r10425);
        double r10427 = r10424 * r10426;
        double r10428 = 2/945;
        double r10429 = 5;
        double r10430 = pow(r10422, r10429);
        double r10431 = r10428 * r10430;
        double r10432 = r10427 + r10431;
        double r10433 = r10423 + r10432;
        return r10433;
}

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 r10434, r10435, r10436, r10437, r10438, r10439;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r10434, "1", 10, MPFR_RNDN);
        mpfr_init(r10435);
        mpfr_init(r10436);
        mpfr_init(r10437);
        mpfr_init(r10438);
        mpfr_init(r10439);
}

double f_im(double x) {
        ;
        mpfr_set_d(r10435, x, MPFR_RNDN);
        mpfr_div(r10436, r10434, r10435, MPFR_RNDN);
        mpfr_tan(r10437, r10435, MPFR_RNDN);
        mpfr_div(r10438, r10434, r10437, MPFR_RNDN);
        mpfr_sub(r10439, r10436, r10438, MPFR_RNDN);
        return mpfr_get_d(r10439, MPFR_RNDN);
}

static mpfr_t r10440, r10441, r10442, r10443, r10444, r10445, r10446, r10447, r10448, r10449, r10450, r10451, r10452;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r10440, "1/3", 10, MPFR_RNDN);
        mpfr_init(r10441);
        mpfr_init(r10442);
        mpfr_init_set_str(r10443, "1/45", 10, MPFR_RNDN);
        mpfr_init_set_str(r10444, "3", 10, MPFR_RNDN);
        mpfr_init(r10445);
        mpfr_init(r10446);
        mpfr_init_set_str(r10447, "2/945", 10, MPFR_RNDN);
        mpfr_init_set_str(r10448, "5", 10, MPFR_RNDN);
        mpfr_init(r10449);
        mpfr_init(r10450);
        mpfr_init(r10451);
        mpfr_init(r10452);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r10441, x, MPFR_RNDN);
        mpfr_mul(r10442, r10440, r10441, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r10445, r10441, r10444, MPFR_RNDN);
        mpfr_mul(r10446, r10443, r10445, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r10449, r10441, r10448, MPFR_RNDN);
        mpfr_mul(r10450, r10447, r10449, MPFR_RNDN);
        mpfr_add(r10451, r10446, r10450, MPFR_RNDN);
        mpfr_add(r10452, r10442, r10451, MPFR_RNDN);
        return mpfr_get_d(r10452, MPFR_RNDN);
}

static mpfr_t r10453, r10454, r10455, r10456, r10457, r10458, r10459, r10460, r10461, r10462, r10463, r10464, r10465;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2384);
        mpfr_init_set_str(r10453, "1/3", 10, MPFR_RNDN);
        mpfr_init(r10454);
        mpfr_init(r10455);
        mpfr_init_set_str(r10456, "1/45", 10, MPFR_RNDN);
        mpfr_init_set_str(r10457, "3", 10, MPFR_RNDN);
        mpfr_init(r10458);
        mpfr_init(r10459);
        mpfr_init_set_str(r10460, "2/945", 10, MPFR_RNDN);
        mpfr_init_set_str(r10461, "5", 10, MPFR_RNDN);
        mpfr_init(r10462);
        mpfr_init(r10463);
        mpfr_init(r10464);
        mpfr_init(r10465);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r10454, x, MPFR_RNDN);
        mpfr_mul(r10455, r10453, r10454, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r10458, r10454, r10457, MPFR_RNDN);
        mpfr_mul(r10459, r10456, r10458, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r10462, r10454, r10461, MPFR_RNDN);
        mpfr_mul(r10463, r10460, r10462, MPFR_RNDN);
        mpfr_add(r10464, r10459, r10463, MPFR_RNDN);
        mpfr_add(r10465, r10455, r10464, MPFR_RNDN);
        return mpfr_get_d(r10465, MPFR_RNDN);
}