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

char *name = "Rectangular parallelepiped of dimension a×b×c";

double f_if() {
        float r41015 = 2;
        float r41016 = 1;
        float r41017 = 9;
        float r41018 = r41016 / r41017;
        float r41019 = r41016 * r41018;
        float r41020 = r41018 * r41018;
        float r41021 = r41019 + r41020;
        float r41022 = r41018 * r41016;
        float r41023 = r41021 + r41022;
        float r41024 = r41015 * r41023;
        return r41024;
}

double f_id() {
        double r41025 = 2;
        double r41026 = 1;
        double r41027 = 9;
        double r41028 = r41026 / r41027;
        double r41029 = r41026 * r41028;
        double r41030 = r41028 * r41028;
        double r41031 = r41029 + r41030;
        double r41032 = r41028 * r41026;
        double r41033 = r41031 + r41032;
        double r41034 = r41025 * r41033;
        return r41034;
}


double f_of() {
        float r41035 = 2;
        float r41036 = 1;
        float r41037 = 9;
        float r41038 = r41036 / r41037;
        float r41039 = r41036 * r41038;
        float r41040 = r41038 * r41038;
        float r41041 = r41039 + r41040;
        float r41042 = r41038 * r41036;
        float r41043 = r41041 + r41042;
        float r41044 = r41035 * r41043;
        return r41044;
}

double f_od() {
        double r41045 = 2;
        double r41046 = 1;
        double r41047 = 9;
        double r41048 = r41046 / r41047;
        double r41049 = r41046 * r41048;
        double r41050 = r41048 * r41048;
        double r41051 = r41049 + r41050;
        double r41052 = r41048 * r41046;
        double r41053 = r41051 + r41052;
        double r41054 = r41045 * r41053;
        return r41054;
}

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 r41055, r41056, r41057, r41058, r41059, r41060, r41061, r41062, r41063, r41064;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r41055, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r41056, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r41057, "9", 10, MPFR_RNDN);
        mpfr_init(r41058);
        mpfr_init(r41059);
        mpfr_init(r41060);
        mpfr_init(r41061);
        mpfr_init(r41062);
        mpfr_init(r41063);
        mpfr_init(r41064);
}

double f_im() {
        ;
        ;
        ;
        mpfr_div(r41058, r41056, r41057, MPFR_RNDN);
        mpfr_mul(r41059, r41056, r41058, MPFR_RNDN);
        mpfr_mul(r41060, r41058, r41058, MPFR_RNDN);
        mpfr_add(r41061, r41059, r41060, MPFR_RNDN);
        mpfr_mul(r41062, r41058, r41056, MPFR_RNDN);
        mpfr_add(r41063, r41061, r41062, MPFR_RNDN);
        mpfr_mul(r41064, r41055, r41063, MPFR_RNDN);
        return mpfr_get_d(r41064, MPFR_RNDN);
}

static mpfr_t r41065, r41066, r41067, r41068, r41069, r41070, r41071, r41072, r41073, r41074;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r41065, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r41066, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r41067, "9", 10, MPFR_RNDN);
        mpfr_init(r41068);
        mpfr_init(r41069);
        mpfr_init(r41070);
        mpfr_init(r41071);
        mpfr_init(r41072);
        mpfr_init(r41073);
        mpfr_init(r41074);
}

double f_fm() {
        ;
        ;
        ;
        mpfr_div(r41068, r41066, r41067, MPFR_RNDN);
        mpfr_mul(r41069, r41066, r41068, MPFR_RNDN);
        mpfr_mul(r41070, r41068, r41068, MPFR_RNDN);
        mpfr_add(r41071, r41069, r41070, MPFR_RNDN);
        mpfr_mul(r41072, r41068, r41066, MPFR_RNDN);
        mpfr_add(r41073, r41071, r41072, MPFR_RNDN);
        mpfr_mul(r41074, r41065, r41073, MPFR_RNDN);
        return mpfr_get_d(r41074, MPFR_RNDN);
}

static mpfr_t r41075, r41076, r41077, r41078, r41079, r41080, r41081, r41082, r41083, r41084;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init_set_str(r41075, "2", 10, MPFR_RNDN);
        mpfr_init_set_str(r41076, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r41077, "9", 10, MPFR_RNDN);
        mpfr_init(r41078);
        mpfr_init(r41079);
        mpfr_init(r41080);
        mpfr_init(r41081);
        mpfr_init(r41082);
        mpfr_init(r41083);
        mpfr_init(r41084);
}

double f_dm() {
        ;
        ;
        ;
        mpfr_div(r41078, r41076, r41077, MPFR_RNDN);
        mpfr_mul(r41079, r41076, r41078, MPFR_RNDN);
        mpfr_mul(r41080, r41078, r41078, MPFR_RNDN);
        mpfr_add(r41081, r41079, r41080, MPFR_RNDN);
        mpfr_mul(r41082, r41078, r41076, MPFR_RNDN);
        mpfr_add(r41083, r41081, r41082, MPFR_RNDN);
        mpfr_mul(r41084, r41075, r41083, MPFR_RNDN);
        return mpfr_get_d(r41084, MPFR_RNDN);
}