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

char *name = "Statistics.Sample:$skurtosis from math-functions-0.1.5.2";

double f_if(float x, float y) {
        float r40456 = x;
        float r40457 = y;
        float r40458 = r40457 * r40457;
        float r40459 = r40456 / r40458;
        float r40460 = 3.0;
        float r40461 = r40459 - r40460;
        return r40461;
}

double f_id(double x, double y) {
        double r40462 = x;
        double r40463 = y;
        double r40464 = r40463 * r40463;
        double r40465 = r40462 / r40464;
        double r40466 = 3.0;
        double r40467 = r40465 - r40466;
        return r40467;
}


double f_of(float x, float y) {
        float r40468 = x;
        float r40469 = y;
        float r40470 = r40468 / r40469;
        float r40471 = r40470 / r40469;
        float r40472 = 3.0;
        float r40473 = r40471 - r40472;
        return r40473;
}

double f_od(double x, double y) {
        double r40474 = x;
        double r40475 = y;
        double r40476 = r40474 / r40475;
        double r40477 = r40476 / r40475;
        double r40478 = 3.0;
        double r40479 = r40477 - r40478;
        return r40479;
}

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 r40480, r40481, r40482, r40483, r40484, r40485;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(400);
        mpfr_init(r40480);
        mpfr_init(r40481);
        mpfr_init(r40482);
        mpfr_init(r40483);
        mpfr_init_set_str(r40484, "3.0", 10, MPFR_RNDN);
        mpfr_init(r40485);
}

double f_im(double x, double y) {
        mpfr_set_d(r40480, x, MPFR_RNDN);
        mpfr_set_d(r40481, y, MPFR_RNDN);
        mpfr_mul(r40482, r40481, r40481, MPFR_RNDN);
        mpfr_div(r40483, r40480, r40482, MPFR_RNDN);
        ;
        mpfr_sub(r40485, r40483, r40484, MPFR_RNDN);
        return mpfr_get_d(r40485, MPFR_RNDN);
}

static mpfr_t r40486, r40487, r40488, r40489, r40490, r40491;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(400);
        mpfr_init(r40486);
        mpfr_init(r40487);
        mpfr_init(r40488);
        mpfr_init(r40489);
        mpfr_init_set_str(r40490, "3.0", 10, MPFR_RNDN);
        mpfr_init(r40491);
}

double f_fm(double x, double y) {
        mpfr_set_d(r40486, x, MPFR_RNDN);
        mpfr_set_d(r40487, y, MPFR_RNDN);
        mpfr_div(r40488, r40486, r40487, MPFR_RNDN);
        mpfr_div(r40489, r40488, r40487, MPFR_RNDN);
        ;
        mpfr_sub(r40491, r40489, r40490, MPFR_RNDN);
        return mpfr_get_d(r40491, MPFR_RNDN);
}

static mpfr_t r40492, r40493, r40494, r40495, r40496, r40497;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(400);
        mpfr_init(r40492);
        mpfr_init(r40493);
        mpfr_init(r40494);
        mpfr_init(r40495);
        mpfr_init_set_str(r40496, "3.0", 10, MPFR_RNDN);
        mpfr_init(r40497);
}

double f_dm(double x, double y) {
        mpfr_set_d(r40492, x, MPFR_RNDN);
        mpfr_set_d(r40493, y, MPFR_RNDN);
        mpfr_div(r40494, r40492, r40493, MPFR_RNDN);
        mpfr_div(r40495, r40494, r40493, MPFR_RNDN);
        ;
        mpfr_sub(r40497, r40495, r40496, MPFR_RNDN);
        return mpfr_get_d(r40497, MPFR_RNDN);
}