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

char *name = "VandenBroeck and Keller, Equation (23)";

double f_if(float F, float B, float x) {
        float r19730 = x;
        float r19731 = 1;
        float r19732 = B;
        float r19733 = tan(r19732);
        float r19734 = r19731 / r19733;
        float r19735 = r19730 * r19734;
        float r19736 = -r19735;
        float r19737 = F;
        float r19738 = sin(r19732);
        float r19739 = r19737 / r19738;
        float r19740 = r19737 * r19737;
        float r19741 = 2;
        float r19742 = r19740 + r19741;
        float r19743 = r19741 * r19730;
        float r19744 = r19742 + r19743;
        float r19745 = r19731 / r19741;
        float r19746 = -r19745;
        float r19747 = pow(r19744, r19746);
        float r19748 = r19739 * r19747;
        float r19749 = r19736 + r19748;
        return r19749;
}

double f_id(double F, double B, double x) {
        double r19750 = x;
        double r19751 = 1;
        double r19752 = B;
        double r19753 = tan(r19752);
        double r19754 = r19751 / r19753;
        double r19755 = r19750 * r19754;
        double r19756 = -r19755;
        double r19757 = F;
        double r19758 = sin(r19752);
        double r19759 = r19757 / r19758;
        double r19760 = r19757 * r19757;
        double r19761 = 2;
        double r19762 = r19760 + r19761;
        double r19763 = r19761 * r19750;
        double r19764 = r19762 + r19763;
        double r19765 = r19751 / r19761;
        double r19766 = -r19765;
        double r19767 = pow(r19764, r19766);
        double r19768 = r19759 * r19767;
        double r19769 = r19756 + r19768;
        return r19769;
}


double f_of(float F, float B, float x) {
        float r19770 = x;
        float r19771 = -r19770;
        float r19772 = B;
        float r19773 = tan(r19772);
        float r19774 = r19771 / r19773;
        float r19775 = F;
        float r19776 = r19770 + r19770;
        float r19777 = 2;
        float r19778 = r19775 * r19775;
        float r19779 = r19777 + r19778;
        float r19780 = r19776 + r19779;
        float r19781 = 1;
        float r19782 = r19781 / r19777;
        float r19783 = pow(r19780, r19782);
        float r19784 = sin(r19772);
        float r19785 = r19783 * r19784;
        float r19786 = r19775 / r19785;
        float r19787 = r19774 + r19786;
        return r19787;
}

double f_od(double F, double B, double x) {
        double r19788 = x;
        double r19789 = -r19788;
        double r19790 = B;
        double r19791 = tan(r19790);
        double r19792 = r19789 / r19791;
        double r19793 = F;
        double r19794 = r19788 + r19788;
        double r19795 = 2;
        double r19796 = r19793 * r19793;
        double r19797 = r19795 + r19796;
        double r19798 = r19794 + r19797;
        double r19799 = 1;
        double r19800 = r19799 / r19795;
        double r19801 = pow(r19798, r19800);
        double r19802 = sin(r19790);
        double r19803 = r19801 * r19802;
        double r19804 = r19793 / r19803;
        double r19805 = r19792 + r19804;
        return r19805;
}

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 r19806, r19807, r19808, r19809, r19810, r19811, r19812, r19813, r19814, r19815, r19816, r19817, r19818, r19819, r19820, r19821, r19822, r19823, r19824, r19825;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r19806);
        mpfr_init_set_str(r19807, "1", 10, MPFR_RNDN);
        mpfr_init(r19808);
        mpfr_init(r19809);
        mpfr_init(r19810);
        mpfr_init(r19811);
        mpfr_init(r19812);
        mpfr_init(r19813);
        mpfr_init(r19814);
        mpfr_init(r19815);
        mpfr_init(r19816);
        mpfr_init_set_str(r19817, "2", 10, MPFR_RNDN);
        mpfr_init(r19818);
        mpfr_init(r19819);
        mpfr_init(r19820);
        mpfr_init(r19821);
        mpfr_init(r19822);
        mpfr_init(r19823);
        mpfr_init(r19824);
        mpfr_init(r19825);
}

double f_im(double F, double B, double x) {
        mpfr_set_d(r19806, x, MPFR_RNDN);
        ;
        mpfr_set_d(r19808, B, MPFR_RNDN);
        mpfr_tan(r19809, r19808, MPFR_RNDN);
        mpfr_div(r19810, r19807, r19809, MPFR_RNDN);
        mpfr_mul(r19811, r19806, r19810, MPFR_RNDN);
        mpfr_neg(r19812, r19811, MPFR_RNDN);
        mpfr_set_d(r19813, F, MPFR_RNDN);
        mpfr_sin(r19814, r19808, MPFR_RNDN);
        mpfr_div(r19815, r19813, r19814, MPFR_RNDN);
        mpfr_mul(r19816, r19813, r19813, MPFR_RNDN);
        ;
        mpfr_add(r19818, r19816, r19817, MPFR_RNDN);
        mpfr_mul(r19819, r19817, r19806, MPFR_RNDN);
        mpfr_add(r19820, r19818, r19819, MPFR_RNDN);
        mpfr_div(r19821, r19807, r19817, MPFR_RNDN);
        mpfr_neg(r19822, r19821, MPFR_RNDN);
        mpfr_pow(r19823, r19820, r19822, MPFR_RNDN);
        mpfr_mul(r19824, r19815, r19823, MPFR_RNDN);
        mpfr_add(r19825, r19812, r19824, MPFR_RNDN);
        return mpfr_get_d(r19825, MPFR_RNDN);
}

static mpfr_t r19826, r19827, r19828, r19829, r19830, r19831, r19832, r19833, r19834, r19835, r19836, r19837, r19838, r19839, r19840, r19841, r19842, r19843;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r19826);
        mpfr_init(r19827);
        mpfr_init(r19828);
        mpfr_init(r19829);
        mpfr_init(r19830);
        mpfr_init(r19831);
        mpfr_init(r19832);
        mpfr_init_set_str(r19833, "2", 10, MPFR_RNDN);
        mpfr_init(r19834);
        mpfr_init(r19835);
        mpfr_init(r19836);
        mpfr_init_set_str(r19837, "1", 10, MPFR_RNDN);
        mpfr_init(r19838);
        mpfr_init(r19839);
        mpfr_init(r19840);
        mpfr_init(r19841);
        mpfr_init(r19842);
        mpfr_init(r19843);
}

double f_fm(double F, double B, double x) {
        mpfr_set_d(r19826, x, MPFR_RNDN);
        mpfr_neg(r19827, r19826, MPFR_RNDN);
        mpfr_set_d(r19828, B, MPFR_RNDN);
        mpfr_tan(r19829, r19828, MPFR_RNDN);
        mpfr_div(r19830, r19827, r19829, MPFR_RNDN);
        mpfr_set_d(r19831, F, MPFR_RNDN);
        mpfr_add(r19832, r19826, r19826, MPFR_RNDN);
        ;
        mpfr_mul(r19834, r19831, r19831, MPFR_RNDN);
        mpfr_add(r19835, r19833, r19834, MPFR_RNDN);
        mpfr_add(r19836, r19832, r19835, MPFR_RNDN);
        ;
        mpfr_div(r19838, r19837, r19833, MPFR_RNDN);
        mpfr_pow(r19839, r19836, r19838, MPFR_RNDN);
        mpfr_sin(r19840, r19828, MPFR_RNDN);
        mpfr_mul(r19841, r19839, r19840, MPFR_RNDN);
        mpfr_div(r19842, r19831, r19841, MPFR_RNDN);
        mpfr_add(r19843, r19830, r19842, MPFR_RNDN);
        return mpfr_get_d(r19843, MPFR_RNDN);
}

static mpfr_t r19844, r19845, r19846, r19847, r19848, r19849, r19850, r19851, r19852, r19853, r19854, r19855, r19856, r19857, r19858, r19859, r19860, r19861;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r19844);
        mpfr_init(r19845);
        mpfr_init(r19846);
        mpfr_init(r19847);
        mpfr_init(r19848);
        mpfr_init(r19849);
        mpfr_init(r19850);
        mpfr_init_set_str(r19851, "2", 10, MPFR_RNDN);
        mpfr_init(r19852);
        mpfr_init(r19853);
        mpfr_init(r19854);
        mpfr_init_set_str(r19855, "1", 10, MPFR_RNDN);
        mpfr_init(r19856);
        mpfr_init(r19857);
        mpfr_init(r19858);
        mpfr_init(r19859);
        mpfr_init(r19860);
        mpfr_init(r19861);
}

double f_dm(double F, double B, double x) {
        mpfr_set_d(r19844, x, MPFR_RNDN);
        mpfr_neg(r19845, r19844, MPFR_RNDN);
        mpfr_set_d(r19846, B, MPFR_RNDN);
        mpfr_tan(r19847, r19846, MPFR_RNDN);
        mpfr_div(r19848, r19845, r19847, MPFR_RNDN);
        mpfr_set_d(r19849, F, MPFR_RNDN);
        mpfr_add(r19850, r19844, r19844, MPFR_RNDN);
        ;
        mpfr_mul(r19852, r19849, r19849, MPFR_RNDN);
        mpfr_add(r19853, r19851, r19852, MPFR_RNDN);
        mpfr_add(r19854, r19850, r19853, MPFR_RNDN);
        ;
        mpfr_div(r19856, r19855, r19851, MPFR_RNDN);
        mpfr_pow(r19857, r19854, r19856, MPFR_RNDN);
        mpfr_sin(r19858, r19846, MPFR_RNDN);
        mpfr_mul(r19859, r19857, r19858, MPFR_RNDN);
        mpfr_div(r19860, r19849, r19859, MPFR_RNDN);
        mpfr_add(r19861, r19848, r19860, MPFR_RNDN);
        return mpfr_get_d(r19861, MPFR_RNDN);
}