#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 r20199 = 1.0f;
        float r20200 = x;
        float r20201 = r20199 / r20200;
        float r20202 = tan(r20200);
        float r20203 = r20199 / r20202;
        float r20204 = r20201 - r20203;
        return r20204;
}

double f_id(double x) {
        double r20205 = 1.0;
        double r20206 = x;
        double r20207 = r20205 / r20206;
        double r20208 = tan(r20206);
        double r20209 = r20205 / r20208;
        double r20210 = r20207 - r20209;
        return r20210;
}


double f_of(float x) {
        float r20211 = 0.022222222222222223f;
        float r20212 = x;
        float r20213 = 3.0f;
        float r20214 = pow(r20212, r20213);
        float r20215 = r20211 * r20214;
        float r20216 = 0.0021164021164021165f;
        float r20217 = 5.0f;
        float r20218 = pow(r20212, r20217);
        float r20219 = r20216 * r20218;
        float r20220 = 0.3333333333333333f;
        float r20221 = r20220 * r20212;
        float r20222 = r20219 + r20221;
        float r20223 = r20215 + r20222;
        return r20223;
}

double f_od(double x) {
        double r20224 = 0.022222222222222223;
        double r20225 = x;
        double r20226 = 3.0;
        double r20227 = pow(r20225, r20226);
        double r20228 = r20224 * r20227;
        double r20229 = 0.0021164021164021165;
        double r20230 = 5.0;
        double r20231 = pow(r20225, r20230);
        double r20232 = r20229 * r20231;
        double r20233 = 0.3333333333333333;
        double r20234 = r20233 * r20225;
        double r20235 = r20232 + r20234;
        double r20236 = r20228 + r20235;
        return r20236;
}

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 r20237, r20238, r20239, r20240, r20241, r20242;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r20237, "1", 10, MPFR_RNDN);
        mpfr_init(r20238);
        mpfr_init(r20239);
        mpfr_init(r20240);
        mpfr_init(r20241);
        mpfr_init(r20242);
}

double f_im(double x) {
        ;
        mpfr_set_d(r20238, x, MPFR_RNDN);
        mpfr_div(r20239, r20237, r20238, MPFR_RNDN);
        mpfr_tan(r20240, r20238, MPFR_RNDN);
        mpfr_div(r20241, r20237, r20240, MPFR_RNDN);
        mpfr_sub(r20242, r20239, r20241, MPFR_RNDN);
        return mpfr_get_d(r20242, MPFR_RNDN);
}

static mpfr_t r20243, r20244, r20245, r20246, r20247, r20248, r20249, r20250, r20251, r20252, r20253, r20254, r20255;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r20243, "1/45", 10, MPFR_RNDN);
        mpfr_init(r20244);
        mpfr_init_set_str(r20245, "3", 10, MPFR_RNDN);
        mpfr_init(r20246);
        mpfr_init(r20247);
        mpfr_init_set_str(r20248, "2/945", 10, MPFR_RNDN);
        mpfr_init_set_str(r20249, "5", 10, MPFR_RNDN);
        mpfr_init(r20250);
        mpfr_init(r20251);
        mpfr_init_set_str(r20252, "1/3", 10, MPFR_RNDN);
        mpfr_init(r20253);
        mpfr_init(r20254);
        mpfr_init(r20255);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r20244, x, MPFR_RNDN);
        ;
        mpfr_pow(r20246, r20244, r20245, MPFR_RNDN);
        mpfr_mul(r20247, r20243, r20246, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r20250, r20244, r20249, MPFR_RNDN);
        mpfr_mul(r20251, r20248, r20250, MPFR_RNDN);
        ;
        mpfr_mul(r20253, r20252, r20244, MPFR_RNDN);
        mpfr_add(r20254, r20251, r20253, MPFR_RNDN);
        mpfr_add(r20255, r20247, r20254, MPFR_RNDN);
        return mpfr_get_d(r20255, MPFR_RNDN);
}

static mpfr_t r20256, r20257, r20258, r20259, r20260, r20261, r20262, r20263, r20264, r20265, r20266, r20267, r20268;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(2448);
        mpfr_init_set_str(r20256, "1/45", 10, MPFR_RNDN);
        mpfr_init(r20257);
        mpfr_init_set_str(r20258, "3", 10, MPFR_RNDN);
        mpfr_init(r20259);
        mpfr_init(r20260);
        mpfr_init_set_str(r20261, "2/945", 10, MPFR_RNDN);
        mpfr_init_set_str(r20262, "5", 10, MPFR_RNDN);
        mpfr_init(r20263);
        mpfr_init(r20264);
        mpfr_init_set_str(r20265, "1/3", 10, MPFR_RNDN);
        mpfr_init(r20266);
        mpfr_init(r20267);
        mpfr_init(r20268);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r20257, x, MPFR_RNDN);
        ;
        mpfr_pow(r20259, r20257, r20258, MPFR_RNDN);
        mpfr_mul(r20260, r20256, r20259, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r20263, r20257, r20262, MPFR_RNDN);
        mpfr_mul(r20264, r20261, r20263, MPFR_RNDN);
        ;
        mpfr_mul(r20266, r20265, r20257, MPFR_RNDN);
        mpfr_add(r20267, r20264, r20266, MPFR_RNDN);
        mpfr_add(r20268, r20260, r20267, MPFR_RNDN);
        return mpfr_get_d(r20268, MPFR_RNDN);
}