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

char *name = "Expression 3, p15";

double f_if(float x) {
        float r9328 = x;
        float r9329 = r9328 * r9328;
        float r9330 = r9328 * r9329;
        float r9331 = r9330 + r9329;
        return r9331;
}

double f_id(double x) {
        double r9332 = x;
        double r9333 = r9332 * r9332;
        double r9334 = r9332 * r9333;
        double r9335 = r9334 + r9333;
        return r9335;
}


double f_of(float x) {
        float r9336 = x;
        float r9337 = r9336 * r9336;
        float r9338 = r9337 + r9336;
        float r9339 = r9336 * r9338;
        return r9339;
}

double f_od(double x) {
        double r9340 = x;
        double r9341 = r9340 * r9340;
        double r9342 = r9341 + r9340;
        double r9343 = r9340 * r9342;
        return r9343;
}

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 r9344, r9345, r9346, r9347;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init(r9344);
        mpfr_init(r9345);
        mpfr_init(r9346);
        mpfr_init(r9347);
}

double f_im(double x) {
        mpfr_set_d(r9344, x, MPFR_RNDN);
        mpfr_mul(r9345, r9344, r9344, MPFR_RNDN);
        mpfr_mul(r9346, r9344, r9345, MPFR_RNDN);
        mpfr_add(r9347, r9346, r9345, MPFR_RNDN);
        return mpfr_get_d(r9347, MPFR_RNDN);
}

static mpfr_t r9348, r9349, r9350, r9351;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9348);
        mpfr_init(r9349);
        mpfr_init(r9350);
        mpfr_init(r9351);
}

double f_fm(double x) {
        mpfr_set_d(r9348, x, MPFR_RNDN);
        mpfr_mul(r9349, r9348, r9348, MPFR_RNDN);
        mpfr_add(r9350, r9349, r9348, MPFR_RNDN);
        mpfr_mul(r9351, r9348, r9350, MPFR_RNDN);
        return mpfr_get_d(r9351, MPFR_RNDN);
}

static mpfr_t r9352, r9353, r9354, r9355;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init(r9352);
        mpfr_init(r9353);
        mpfr_init(r9354);
        mpfr_init(r9355);
}

double f_dm(double x) {
        mpfr_set_d(r9352, x, MPFR_RNDN);
        mpfr_mul(r9353, r9352, r9352, MPFR_RNDN);
        mpfr_add(r9354, r9353, r9352, MPFR_RNDN);
        mpfr_mul(r9355, r9352, r9354, MPFR_RNDN);
        return mpfr_get_d(r9355, MPFR_RNDN);
}