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

char *name = "Compound Interest";

double f_if(float i, float n) {
        float r23850 = 100;
        float r23851 = 1;
        float r23852 = i;
        float r23853 = n;
        float r23854 = r23852 / r23853;
        float r23855 = r23851 + r23854;
        float r23856 = pow(r23855, r23853);
        float r23857 = r23856 - r23851;
        float r23858 = r23857 / r23854;
        float r23859 = r23850 * r23858;
        return r23859;
}

double f_id(double i, double n) {
        double r23860 = 100;
        double r23861 = 1;
        double r23862 = i;
        double r23863 = n;
        double r23864 = r23862 / r23863;
        double r23865 = r23861 + r23864;
        double r23866 = pow(r23865, r23863);
        double r23867 = r23866 - r23861;
        double r23868 = r23867 / r23864;
        double r23869 = r23860 * r23868;
        return r23869;
}


double f_of(float i, float n) {
        float r23870 = 1;
        float r23871 = 1/2;
        float r23872 = i;
        float r23873 = r23871 * r23872;
        float r23874 = r23870 + r23873;
        float r23875 = log(r23874);
        float r23876 = 93.98030941154852;
        bool r23877 = r23875 <= r23876;
        float r23878 = r23872 * r23871;
        float r23879 = r23878 + r23870;
        float r23880 = n;
        float r23881 = 100;
        float r23882 = r23880 * r23881;
        float r23883 = r23870 / r23882;
        float r23884 = r23879 / r23883;
        float r23885 = 619.7267873526765;
        bool r23886 = r23875 <= r23885;
        float r23887 = r23870 / r23880;
        float r23888 = r23887 / r23881;
        float r23889 = exp(r23888);
        float r23890 = log(r23889);
        float r23891 = r23879 / r23890;
        float r23892 = r23872 / r23880;
        float r23893 = r23870 + r23892;
        float r23894 = pow(r23893, r23880);
        float r23895 = r23894 - r23870;
        float r23896 = r23895 / r23892;
        float r23897 = r23881 * r23896;
        float r23898 = r23886 ? r23891 : r23897;
        float r23899 = r23877 ? r23884 : r23898;
        return r23899;
}

double f_od(double i, double n) {
        double r23900 = 1;
        double r23901 = 1/2;
        double r23902 = i;
        double r23903 = r23901 * r23902;
        double r23904 = r23900 + r23903;
        double r23905 = log(r23904);
        double r23906 = 93.98030941154852;
        bool r23907 = r23905 <= r23906;
        double r23908 = r23902 * r23901;
        double r23909 = r23908 + r23900;
        double r23910 = n;
        double r23911 = 100;
        double r23912 = r23910 * r23911;
        double r23913 = r23900 / r23912;
        double r23914 = r23909 / r23913;
        double r23915 = 619.7267873526765;
        bool r23916 = r23905 <= r23915;
        double r23917 = r23900 / r23910;
        double r23918 = r23917 / r23911;
        double r23919 = exp(r23918);
        double r23920 = log(r23919);
        double r23921 = r23909 / r23920;
        double r23922 = r23902 / r23910;
        double r23923 = r23900 + r23922;
        double r23924 = pow(r23923, r23910);
        double r23925 = r23924 - r23900;
        double r23926 = r23925 / r23922;
        double r23927 = r23911 * r23926;
        double r23928 = r23916 ? r23921 : r23927;
        double r23929 = r23907 ? r23914 : r23928;
        return r23929;
}

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 r23930, r23931, r23932, r23933, r23934, r23935, r23936, r23937, r23938, r23939;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(3216);
        mpfr_init_set_str(r23930, "100", 10, MPFR_RNDN);
        mpfr_init_set_str(r23931, "1", 10, MPFR_RNDN);
        mpfr_init(r23932);
        mpfr_init(r23933);
        mpfr_init(r23934);
        mpfr_init(r23935);
        mpfr_init(r23936);
        mpfr_init(r23937);
        mpfr_init(r23938);
        mpfr_init(r23939);
}

double f_im(double i, double n) {
        ;
        ;
        mpfr_set_d(r23932, i, MPFR_RNDN);
        mpfr_set_d(r23933, n, MPFR_RNDN);
        mpfr_div(r23934, r23932, r23933, MPFR_RNDN);
        mpfr_add(r23935, r23931, r23934, MPFR_RNDN);
        mpfr_pow(r23936, r23935, r23933, MPFR_RNDN);
        mpfr_sub(r23937, r23936, r23931, MPFR_RNDN);
        mpfr_div(r23938, r23937, r23934, MPFR_RNDN);
        mpfr_mul(r23939, r23930, r23938, MPFR_RNDN);
        return mpfr_get_d(r23939, MPFR_RNDN);
}

static mpfr_t r23940, r23941, r23942, r23943, r23944, r23945, r23946, r23947, r23948, r23949, r23950, r23951, r23952, r23953, r23954, r23955, r23956, r23957, r23958, r23959, r23960, r23961, r23962, r23963, r23964, r23965, r23966, r23967, r23968, r23969;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(3216);
        mpfr_init_set_str(r23940, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r23941, "1/2", 10, MPFR_RNDN);
        mpfr_init(r23942);
        mpfr_init(r23943);
        mpfr_init(r23944);
        mpfr_init(r23945);
        mpfr_init_set_str(r23946, "93.98030941154852", 10, MPFR_RNDN);
        mpfr_init(r23947);
        mpfr_init(r23948);
        mpfr_init(r23949);
        mpfr_init(r23950);
        mpfr_init_set_str(r23951, "100", 10, MPFR_RNDN);
        mpfr_init(r23952);
        mpfr_init(r23953);
        mpfr_init(r23954);
        mpfr_init_set_str(r23955, "619.7267873526765", 10, MPFR_RNDN);
        mpfr_init(r23956);
        mpfr_init(r23957);
        mpfr_init(r23958);
        mpfr_init(r23959);
        mpfr_init(r23960);
        mpfr_init(r23961);
        mpfr_init(r23962);
        mpfr_init(r23963);
        mpfr_init(r23964);
        mpfr_init(r23965);
        mpfr_init(r23966);
        mpfr_init(r23967);
        mpfr_init(r23968);
        mpfr_init(r23969);
}

double f_fm(double i, double n) {
        ;
        ;
        mpfr_set_d(r23942, i, MPFR_RNDN);
        mpfr_mul(r23943, r23941, r23942, MPFR_RNDN);
        mpfr_add(r23944, r23940, r23943, MPFR_RNDN);
        mpfr_log(r23945, r23944, MPFR_RNDN);
        ;
        mpfr_set_si(r23947, mpfr_cmp(r23945, r23946) <= 0, MPFR_RNDN);
        mpfr_mul(r23948, r23942, r23941, MPFR_RNDN);
        mpfr_add(r23949, r23948, r23940, MPFR_RNDN);
        mpfr_set_d(r23950, n, MPFR_RNDN);
        ;
        mpfr_mul(r23952, r23950, r23951, MPFR_RNDN);
        mpfr_div(r23953, r23940, r23952, MPFR_RNDN);
        mpfr_div(r23954, r23949, r23953, MPFR_RNDN);
        ;
        mpfr_set_si(r23956, mpfr_cmp(r23945, r23955) <= 0, MPFR_RNDN);
        mpfr_div(r23957, r23940, r23950, MPFR_RNDN);
        mpfr_div(r23958, r23957, r23951, MPFR_RNDN);
        mpfr_exp(r23959, r23958, MPFR_RNDN);
        mpfr_log(r23960, r23959, MPFR_RNDN);
        mpfr_div(r23961, r23949, r23960, MPFR_RNDN);
        mpfr_div(r23962, r23942, r23950, MPFR_RNDN);
        mpfr_add(r23963, r23940, r23962, MPFR_RNDN);
        mpfr_pow(r23964, r23963, r23950, MPFR_RNDN);
        mpfr_sub(r23965, r23964, r23940, MPFR_RNDN);
        mpfr_div(r23966, r23965, r23962, MPFR_RNDN);
        mpfr_mul(r23967, r23951, r23966, MPFR_RNDN);
        if (mpfr_get_si(r23956, MPFR_RNDN)) { mpfr_set(r23968, r23961, MPFR_RNDN); } else { mpfr_set(r23968, r23967, MPFR_RNDN); };
        if (mpfr_get_si(r23947, MPFR_RNDN)) { mpfr_set(r23969, r23954, MPFR_RNDN); } else { mpfr_set(r23969, r23968, MPFR_RNDN); };
        return mpfr_get_d(r23969, MPFR_RNDN);
}

static mpfr_t r23970, r23971, r23972, r23973, r23974, r23975, r23976, r23977, r23978, r23979, r23980, r23981, r23982, r23983, r23984, r23985, r23986, r23987, r23988, r23989, r23990, r23991, r23992, r23993, r23994, r23995, r23996, r23997, r23998, r23999;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(3216);
        mpfr_init_set_str(r23970, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r23971, "1/2", 10, MPFR_RNDN);
        mpfr_init(r23972);
        mpfr_init(r23973);
        mpfr_init(r23974);
        mpfr_init(r23975);
        mpfr_init_set_str(r23976, "93.98030941154852", 10, MPFR_RNDN);
        mpfr_init(r23977);
        mpfr_init(r23978);
        mpfr_init(r23979);
        mpfr_init(r23980);
        mpfr_init_set_str(r23981, "100", 10, MPFR_RNDN);
        mpfr_init(r23982);
        mpfr_init(r23983);
        mpfr_init(r23984);
        mpfr_init_set_str(r23985, "619.7267873526765", 10, MPFR_RNDN);
        mpfr_init(r23986);
        mpfr_init(r23987);
        mpfr_init(r23988);
        mpfr_init(r23989);
        mpfr_init(r23990);
        mpfr_init(r23991);
        mpfr_init(r23992);
        mpfr_init(r23993);
        mpfr_init(r23994);
        mpfr_init(r23995);
        mpfr_init(r23996);
        mpfr_init(r23997);
        mpfr_init(r23998);
        mpfr_init(r23999);
}

double f_dm(double i, double n) {
        ;
        ;
        mpfr_set_d(r23972, i, MPFR_RNDN);
        mpfr_mul(r23973, r23971, r23972, MPFR_RNDN);
        mpfr_add(r23974, r23970, r23973, MPFR_RNDN);
        mpfr_log(r23975, r23974, MPFR_RNDN);
        ;
        mpfr_set_si(r23977, mpfr_cmp(r23975, r23976) <= 0, MPFR_RNDN);
        mpfr_mul(r23978, r23972, r23971, MPFR_RNDN);
        mpfr_add(r23979, r23978, r23970, MPFR_RNDN);
        mpfr_set_d(r23980, n, MPFR_RNDN);
        ;
        mpfr_mul(r23982, r23980, r23981, MPFR_RNDN);
        mpfr_div(r23983, r23970, r23982, MPFR_RNDN);
        mpfr_div(r23984, r23979, r23983, MPFR_RNDN);
        ;
        mpfr_set_si(r23986, mpfr_cmp(r23975, r23985) <= 0, MPFR_RNDN);
        mpfr_div(r23987, r23970, r23980, MPFR_RNDN);
        mpfr_div(r23988, r23987, r23981, MPFR_RNDN);
        mpfr_exp(r23989, r23988, MPFR_RNDN);
        mpfr_log(r23990, r23989, MPFR_RNDN);
        mpfr_div(r23991, r23979, r23990, MPFR_RNDN);
        mpfr_div(r23992, r23972, r23980, MPFR_RNDN);
        mpfr_add(r23993, r23970, r23992, MPFR_RNDN);
        mpfr_pow(r23994, r23993, r23980, MPFR_RNDN);
        mpfr_sub(r23995, r23994, r23970, MPFR_RNDN);
        mpfr_div(r23996, r23995, r23992, MPFR_RNDN);
        mpfr_mul(r23997, r23981, r23996, MPFR_RNDN);
        if (mpfr_get_si(r23986, MPFR_RNDN)) { mpfr_set(r23998, r23991, MPFR_RNDN); } else { mpfr_set(r23998, r23997, MPFR_RNDN); };
        if (mpfr_get_si(r23977, MPFR_RNDN)) { mpfr_set(r23999, r23984, MPFR_RNDN); } else { mpfr_set(r23999, r23998, MPFR_RNDN); };
        return mpfr_get_d(r23999, MPFR_RNDN);
}