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

char *name = "Jmat.Real.gamma, branch z greater than 0.5";

double f_if(float z) {
        float r30567 = atan2(1.0, 0.0);
        float r30568 = 2;
        float r30569 = r30567 * r30568;
        float r30570 = sqrt(r30569);
        float r30571 = z;
        float r30572 = 1;
        float r30573 = r30571 - r30572;
        float r30574 = 7;
        float r30575 = r30573 + r30574;
        float r30576 = 0.5;
        float r30577 = r30575 + r30576;
        float r30578 = r30573 + r30576;
        float r30579 = pow(r30577, r30578);
        float r30580 = r30570 * r30579;
        float r30581 = -r30577;
        float r30582 = exp(r30581);
        float r30583 = r30580 * r30582;
        float r30584 = 0.9999999999998099;
        float r30585 = 676.5203681218851;
        float r30586 = r30573 + r30572;
        float r30587 = r30585 / r30586;
        float r30588 = r30584 + r30587;
        float r30589 = -1259.1392167224028;
        float r30590 = r30573 + r30568;
        float r30591 = r30589 / r30590;
        float r30592 = r30588 + r30591;
        float r30593 = 771.3234287776531;
        float r30594 = 3;
        float r30595 = r30573 + r30594;
        float r30596 = r30593 / r30595;
        float r30597 = r30592 + r30596;
        float r30598 = -176.6150291621406;
        float r30599 = 4;
        float r30600 = r30573 + r30599;
        float r30601 = r30598 / r30600;
        float r30602 = r30597 + r30601;
        float r30603 = 12.507343278686905;
        float r30604 = 5;
        float r30605 = r30573 + r30604;
        float r30606 = r30603 / r30605;
        float r30607 = r30602 + r30606;
        float r30608 = -0.13857109526572012;
        float r30609 = 6;
        float r30610 = r30573 + r30609;
        float r30611 = r30608 / r30610;
        float r30612 = r30607 + r30611;
        float r30613 = 9.984369578019572e-06;
        float r30614 = r30613 / r30575;
        float r30615 = r30612 + r30614;
        float r30616 = 1.5056327351493116e-07;
        float r30617 = 8;
        float r30618 = r30573 + r30617;
        float r30619 = r30616 / r30618;
        float r30620 = r30615 + r30619;
        float r30621 = r30583 * r30620;
        return r30621;
}

double f_id(double z) {
        double r30622 = atan2(1.0, 0.0);
        double r30623 = 2;
        double r30624 = r30622 * r30623;
        double r30625 = sqrt(r30624);
        double r30626 = z;
        double r30627 = 1;
        double r30628 = r30626 - r30627;
        double r30629 = 7;
        double r30630 = r30628 + r30629;
        double r30631 = 0.5;
        double r30632 = r30630 + r30631;
        double r30633 = r30628 + r30631;
        double r30634 = pow(r30632, r30633);
        double r30635 = r30625 * r30634;
        double r30636 = -r30632;
        double r30637 = exp(r30636);
        double r30638 = r30635 * r30637;
        double r30639 = 0.9999999999998099;
        double r30640 = 676.5203681218851;
        double r30641 = r30628 + r30627;
        double r30642 = r30640 / r30641;
        double r30643 = r30639 + r30642;
        double r30644 = -1259.1392167224028;
        double r30645 = r30628 + r30623;
        double r30646 = r30644 / r30645;
        double r30647 = r30643 + r30646;
        double r30648 = 771.3234287776531;
        double r30649 = 3;
        double r30650 = r30628 + r30649;
        double r30651 = r30648 / r30650;
        double r30652 = r30647 + r30651;
        double r30653 = -176.6150291621406;
        double r30654 = 4;
        double r30655 = r30628 + r30654;
        double r30656 = r30653 / r30655;
        double r30657 = r30652 + r30656;
        double r30658 = 12.507343278686905;
        double r30659 = 5;
        double r30660 = r30628 + r30659;
        double r30661 = r30658 / r30660;
        double r30662 = r30657 + r30661;
        double r30663 = -0.13857109526572012;
        double r30664 = 6;
        double r30665 = r30628 + r30664;
        double r30666 = r30663 / r30665;
        double r30667 = r30662 + r30666;
        double r30668 = 9.984369578019572e-06;
        double r30669 = r30668 / r30630;
        double r30670 = r30667 + r30669;
        double r30671 = 1.5056327351493116e-07;
        double r30672 = 8;
        double r30673 = r30628 + r30672;
        double r30674 = r30671 / r30673;
        double r30675 = r30670 + r30674;
        double r30676 = r30638 * r30675;
        return r30676;
}


double f_of(float z) {
        float r30677 = atan2(1.0, 0.0);
        float r30678 = r30677 + r30677;
        float r30679 = sqrt(r30678);
        float r30680 = 0.5;
        float r30681 = 7;
        float r30682 = r30680 + r30681;
        float r30683 = z;
        float r30684 = 1;
        float r30685 = r30683 - r30684;
        float r30686 = r30682 + r30685;
        float r30687 = r30680 + r30685;
        float r30688 = pow(r30686, r30687);
        float r30689 = exp(r30686);
        float r30690 = r30688 / r30689;
        float r30691 = r30679 * r30690;
        float r30692 = -0.13857109526572012;
        float r30693 = 6;
        float r30694 = r30684 - r30693;
        float r30695 = r30683 - r30694;
        float r30696 = r30692 / r30695;
        float r30697 = 9.984369578019572e-06;
        float r30698 = r30685 + r30681;
        float r30699 = r30697 / r30698;
        float r30700 = r30696 + r30699;
        float r30701 = 12.507343278686905;
        float r30702 = 5;
        float r30703 = r30683 + r30702;
        float r30704 = r30703 - r30684;
        float r30705 = r30701 / r30704;
        float r30706 = r30700 + r30705;
        float r30707 = -176.6150291621406;
        float r30708 = 4;
        float r30709 = r30683 + r30708;
        float r30710 = r30709 - r30684;
        float r30711 = r30707 / r30710;
        float r30712 = 676.5203681218851;
        float r30713 = 0;
        float r30714 = r30683 - r30713;
        float r30715 = r30712 / r30714;
        float r30716 = 0.9999999999998099;
        float r30717 = r30715 + r30716;
        float r30718 = r30711 + r30717;
        float r30719 = -1259.1392167224028;
        float r30720 = 2;
        float r30721 = r30684 - r30720;
        float r30722 = r30683 - r30721;
        float r30723 = r30719 / r30722;
        float r30724 = 771.3234287776531;
        float r30725 = 3;
        float r30726 = r30725 + r30683;
        float r30727 = r30726 - r30684;
        float r30728 = r30724 / r30727;
        float r30729 = r30723 + r30728;
        float r30730 = r30718 + r30729;
        float r30731 = r30706 + r30730;
        float r30732 = 1.5056327351493116e-07;
        float r30733 = 8;
        float r30734 = r30683 + r30733;
        float r30735 = r30734 - r30684;
        float r30736 = r30732 / r30735;
        float r30737 = r30731 + r30736;
        float r30738 = r30691 * r30737;
        return r30738;
}

double f_od(double z) {
        double r30739 = atan2(1.0, 0.0);
        double r30740 = r30739 + r30739;
        double r30741 = sqrt(r30740);
        double r30742 = 0.5;
        double r30743 = 7;
        double r30744 = r30742 + r30743;
        double r30745 = z;
        double r30746 = 1;
        double r30747 = r30745 - r30746;
        double r30748 = r30744 + r30747;
        double r30749 = r30742 + r30747;
        double r30750 = pow(r30748, r30749);
        double r30751 = exp(r30748);
        double r30752 = r30750 / r30751;
        double r30753 = r30741 * r30752;
        double r30754 = -0.13857109526572012;
        double r30755 = 6;
        double r30756 = r30746 - r30755;
        double r30757 = r30745 - r30756;
        double r30758 = r30754 / r30757;
        double r30759 = 9.984369578019572e-06;
        double r30760 = r30747 + r30743;
        double r30761 = r30759 / r30760;
        double r30762 = r30758 + r30761;
        double r30763 = 12.507343278686905;
        double r30764 = 5;
        double r30765 = r30745 + r30764;
        double r30766 = r30765 - r30746;
        double r30767 = r30763 / r30766;
        double r30768 = r30762 + r30767;
        double r30769 = -176.6150291621406;
        double r30770 = 4;
        double r30771 = r30745 + r30770;
        double r30772 = r30771 - r30746;
        double r30773 = r30769 / r30772;
        double r30774 = 676.5203681218851;
        double r30775 = 0;
        double r30776 = r30745 - r30775;
        double r30777 = r30774 / r30776;
        double r30778 = 0.9999999999998099;
        double r30779 = r30777 + r30778;
        double r30780 = r30773 + r30779;
        double r30781 = -1259.1392167224028;
        double r30782 = 2;
        double r30783 = r30746 - r30782;
        double r30784 = r30745 - r30783;
        double r30785 = r30781 / r30784;
        double r30786 = 771.3234287776531;
        double r30787 = 3;
        double r30788 = r30787 + r30745;
        double r30789 = r30788 - r30746;
        double r30790 = r30786 / r30789;
        double r30791 = r30785 + r30790;
        double r30792 = r30780 + r30791;
        double r30793 = r30768 + r30792;
        double r30794 = 1.5056327351493116e-07;
        double r30795 = 8;
        double r30796 = r30745 + r30795;
        double r30797 = r30796 - r30746;
        double r30798 = r30794 / r30797;
        double r30799 = r30793 + r30798;
        double r30800 = r30753 * r30799;
        return r30800;
}

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 r30801, r30802, r30803, r30804, r30805, r30806, r30807, r30808, r30809, r30810, r30811, r30812, r30813, r30814, r30815, r30816, r30817, r30818, r30819, r30820, r30821, r30822, r30823, r30824, r30825, r30826, r30827, r30828, r30829, r30830, r30831, r30832, r30833, r30834, r30835, r30836, r30837, r30838, r30839, r30840, r30841, r30842, r30843, r30844, r30845, r30846, r30847, r30848, r30849, r30850, r30851, r30852, r30853, r30854, r30855;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r30801);
        mpfr_init_set_str(r30802, "2", 10, MPFR_RNDN);
        mpfr_init(r30803);
        mpfr_init(r30804);
        mpfr_init(r30805);
        mpfr_init_set_str(r30806, "1", 10, MPFR_RNDN);
        mpfr_init(r30807);
        mpfr_init_set_str(r30808, "7", 10, MPFR_RNDN);
        mpfr_init(r30809);
        mpfr_init_set_str(r30810, "0.5", 10, MPFR_RNDN);
        mpfr_init(r30811);
        mpfr_init(r30812);
        mpfr_init(r30813);
        mpfr_init(r30814);
        mpfr_init(r30815);
        mpfr_init(r30816);
        mpfr_init(r30817);
        mpfr_init_set_str(r30818, "0.9999999999998099", 10, MPFR_RNDN);
        mpfr_init_set_str(r30819, "676.5203681218851", 10, MPFR_RNDN);
        mpfr_init(r30820);
        mpfr_init(r30821);
        mpfr_init(r30822);
        mpfr_init_set_str(r30823, "-1259.1392167224028", 10, MPFR_RNDN);
        mpfr_init(r30824);
        mpfr_init(r30825);
        mpfr_init(r30826);
        mpfr_init_set_str(r30827, "771.3234287776531", 10, MPFR_RNDN);
        mpfr_init_set_str(r30828, "3", 10, MPFR_RNDN);
        mpfr_init(r30829);
        mpfr_init(r30830);
        mpfr_init(r30831);
        mpfr_init_set_str(r30832, "-176.6150291621406", 10, MPFR_RNDN);
        mpfr_init_set_str(r30833, "4", 10, MPFR_RNDN);
        mpfr_init(r30834);
        mpfr_init(r30835);
        mpfr_init(r30836);
        mpfr_init_set_str(r30837, "12.507343278686905", 10, MPFR_RNDN);
        mpfr_init_set_str(r30838, "5", 10, MPFR_RNDN);
        mpfr_init(r30839);
        mpfr_init(r30840);
        mpfr_init(r30841);
        mpfr_init_set_str(r30842, "-0.13857109526572012", 10, MPFR_RNDN);
        mpfr_init_set_str(r30843, "6", 10, MPFR_RNDN);
        mpfr_init(r30844);
        mpfr_init(r30845);
        mpfr_init(r30846);
        mpfr_init_set_str(r30847, "9.984369578019572e-06", 10, MPFR_RNDN);
        mpfr_init(r30848);
        mpfr_init(r30849);
        mpfr_init_set_str(r30850, "1.5056327351493116e-07", 10, MPFR_RNDN);
        mpfr_init_set_str(r30851, "8", 10, MPFR_RNDN);
        mpfr_init(r30852);
        mpfr_init(r30853);
        mpfr_init(r30854);
        mpfr_init(r30855);
}

double f_im(double z) {
        mpfr_const_pi(r30801, MPFR_RNDN);
        ;
        mpfr_mul(r30803, r30801, r30802, MPFR_RNDN);
        mpfr_sqrt(r30804, r30803, MPFR_RNDN);
        mpfr_set_d(r30805, z, MPFR_RNDN);
        ;
        mpfr_sub(r30807, r30805, r30806, MPFR_RNDN);
        ;
        mpfr_add(r30809, r30807, r30808, MPFR_RNDN);
        ;
        mpfr_add(r30811, r30809, r30810, MPFR_RNDN);
        mpfr_add(r30812, r30807, r30810, MPFR_RNDN);
        mpfr_pow(r30813, r30811, r30812, MPFR_RNDN);
        mpfr_mul(r30814, r30804, r30813, MPFR_RNDN);
        mpfr_neg(r30815, r30811, MPFR_RNDN);
        mpfr_exp(r30816, r30815, MPFR_RNDN);
        mpfr_mul(r30817, r30814, r30816, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30820, r30807, r30806, MPFR_RNDN);
        mpfr_div(r30821, r30819, r30820, MPFR_RNDN);
        mpfr_add(r30822, r30818, r30821, MPFR_RNDN);
        ;
        mpfr_add(r30824, r30807, r30802, MPFR_RNDN);
        mpfr_div(r30825, r30823, r30824, MPFR_RNDN);
        mpfr_add(r30826, r30822, r30825, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30829, r30807, r30828, MPFR_RNDN);
        mpfr_div(r30830, r30827, r30829, MPFR_RNDN);
        mpfr_add(r30831, r30826, r30830, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30834, r30807, r30833, MPFR_RNDN);
        mpfr_div(r30835, r30832, r30834, MPFR_RNDN);
        mpfr_add(r30836, r30831, r30835, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30839, r30807, r30838, MPFR_RNDN);
        mpfr_div(r30840, r30837, r30839, MPFR_RNDN);
        mpfr_add(r30841, r30836, r30840, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30844, r30807, r30843, MPFR_RNDN);
        mpfr_div(r30845, r30842, r30844, MPFR_RNDN);
        mpfr_add(r30846, r30841, r30845, MPFR_RNDN);
        ;
        mpfr_div(r30848, r30847, r30809, MPFR_RNDN);
        mpfr_add(r30849, r30846, r30848, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30852, r30807, r30851, MPFR_RNDN);
        mpfr_div(r30853, r30850, r30852, MPFR_RNDN);
        mpfr_add(r30854, r30849, r30853, MPFR_RNDN);
        mpfr_mul(r30855, r30817, r30854, MPFR_RNDN);
        return mpfr_get_d(r30855, MPFR_RNDN);
}

static mpfr_t r30856, r30857, r30858, r30859, r30860, r30861, r30862, r30863, r30864, r30865, r30866, r30867, r30868, r30869, r30870, r30871, r30872, r30873, r30874, r30875, r30876, r30877, r30878, r30879, r30880, r30881, r30882, r30883, r30884, r30885, r30886, r30887, r30888, r30889, r30890, r30891, r30892, r30893, r30894, r30895, r30896, r30897, r30898, r30899, r30900, r30901, r30902, r30903, r30904, r30905, r30906, r30907, r30908, r30909, r30910, r30911, r30912, r30913, r30914, r30915, r30916, r30917;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r30856);
        mpfr_init(r30857);
        mpfr_init(r30858);
        mpfr_init_set_str(r30859, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r30860, "7", 10, MPFR_RNDN);
        mpfr_init(r30861);
        mpfr_init(r30862);
        mpfr_init_set_str(r30863, "1", 10, MPFR_RNDN);
        mpfr_init(r30864);
        mpfr_init(r30865);
        mpfr_init(r30866);
        mpfr_init(r30867);
        mpfr_init(r30868);
        mpfr_init(r30869);
        mpfr_init(r30870);
        mpfr_init_set_str(r30871, "-0.13857109526572012", 10, MPFR_RNDN);
        mpfr_init_set_str(r30872, "6", 10, MPFR_RNDN);
        mpfr_init(r30873);
        mpfr_init(r30874);
        mpfr_init(r30875);
        mpfr_init_set_str(r30876, "9.984369578019572e-06", 10, MPFR_RNDN);
        mpfr_init(r30877);
        mpfr_init(r30878);
        mpfr_init(r30879);
        mpfr_init_set_str(r30880, "12.507343278686905", 10, MPFR_RNDN);
        mpfr_init_set_str(r30881, "5", 10, MPFR_RNDN);
        mpfr_init(r30882);
        mpfr_init(r30883);
        mpfr_init(r30884);
        mpfr_init(r30885);
        mpfr_init_set_str(r30886, "-176.6150291621406", 10, MPFR_RNDN);
        mpfr_init_set_str(r30887, "4", 10, MPFR_RNDN);
        mpfr_init(r30888);
        mpfr_init(r30889);
        mpfr_init(r30890);
        mpfr_init_set_str(r30891, "676.5203681218851", 10, MPFR_RNDN);
        mpfr_init_set_str(r30892, "0", 10, MPFR_RNDN);
        mpfr_init(r30893);
        mpfr_init(r30894);
        mpfr_init_set_str(r30895, "0.9999999999998099", 10, MPFR_RNDN);
        mpfr_init(r30896);
        mpfr_init(r30897);
        mpfr_init_set_str(r30898, "-1259.1392167224028", 10, MPFR_RNDN);
        mpfr_init_set_str(r30899, "2", 10, MPFR_RNDN);
        mpfr_init(r30900);
        mpfr_init(r30901);
        mpfr_init(r30902);
        mpfr_init_set_str(r30903, "771.3234287776531", 10, MPFR_RNDN);
        mpfr_init_set_str(r30904, "3", 10, MPFR_RNDN);
        mpfr_init(r30905);
        mpfr_init(r30906);
        mpfr_init(r30907);
        mpfr_init(r30908);
        mpfr_init(r30909);
        mpfr_init(r30910);
        mpfr_init_set_str(r30911, "1.5056327351493116e-07", 10, MPFR_RNDN);
        mpfr_init_set_str(r30912, "8", 10, MPFR_RNDN);
        mpfr_init(r30913);
        mpfr_init(r30914);
        mpfr_init(r30915);
        mpfr_init(r30916);
        mpfr_init(r30917);
}

double f_fm(double z) {
        mpfr_const_pi(r30856, MPFR_RNDN);
        mpfr_add(r30857, r30856, r30856, MPFR_RNDN);
        mpfr_sqrt(r30858, r30857, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30861, r30859, r30860, MPFR_RNDN);
        mpfr_set_d(r30862, z, MPFR_RNDN);
        ;
        mpfr_sub(r30864, r30862, r30863, MPFR_RNDN);
        mpfr_add(r30865, r30861, r30864, MPFR_RNDN);
        mpfr_add(r30866, r30859, r30864, MPFR_RNDN);
        mpfr_pow(r30867, r30865, r30866, MPFR_RNDN);
        mpfr_exp(r30868, r30865, MPFR_RNDN);
        mpfr_div(r30869, r30867, r30868, MPFR_RNDN);
        mpfr_mul(r30870, r30858, r30869, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30873, r30863, r30872, MPFR_RNDN);
        mpfr_sub(r30874, r30862, r30873, MPFR_RNDN);
        mpfr_div(r30875, r30871, r30874, MPFR_RNDN);
        ;
        mpfr_add(r30877, r30864, r30860, MPFR_RNDN);
        mpfr_div(r30878, r30876, r30877, MPFR_RNDN);
        mpfr_add(r30879, r30875, r30878, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30882, r30862, r30881, MPFR_RNDN);
        mpfr_sub(r30883, r30882, r30863, MPFR_RNDN);
        mpfr_div(r30884, r30880, r30883, MPFR_RNDN);
        mpfr_add(r30885, r30879, r30884, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30888, r30862, r30887, MPFR_RNDN);
        mpfr_sub(r30889, r30888, r30863, MPFR_RNDN);
        mpfr_div(r30890, r30886, r30889, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30893, r30862, r30892, MPFR_RNDN);
        mpfr_div(r30894, r30891, r30893, MPFR_RNDN);
        ;
        mpfr_add(r30896, r30894, r30895, MPFR_RNDN);
        mpfr_add(r30897, r30890, r30896, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30900, r30863, r30899, MPFR_RNDN);
        mpfr_sub(r30901, r30862, r30900, MPFR_RNDN);
        mpfr_div(r30902, r30898, r30901, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30905, r30904, r30862, MPFR_RNDN);
        mpfr_sub(r30906, r30905, r30863, MPFR_RNDN);
        mpfr_div(r30907, r30903, r30906, MPFR_RNDN);
        mpfr_add(r30908, r30902, r30907, MPFR_RNDN);
        mpfr_add(r30909, r30897, r30908, MPFR_RNDN);
        mpfr_add(r30910, r30885, r30909, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30913, r30862, r30912, MPFR_RNDN);
        mpfr_sub(r30914, r30913, r30863, MPFR_RNDN);
        mpfr_div(r30915, r30911, r30914, MPFR_RNDN);
        mpfr_add(r30916, r30910, r30915, MPFR_RNDN);
        mpfr_mul(r30917, r30870, r30916, MPFR_RNDN);
        return mpfr_get_d(r30917, MPFR_RNDN);
}

static mpfr_t r30918, r30919, r30920, r30921, r30922, r30923, r30924, r30925, r30926, r30927, r30928, r30929, r30930, r30931, r30932, r30933, r30934, r30935, r30936, r30937, r30938, r30939, r30940, r30941, r30942, r30943, r30944, r30945, r30946, r30947, r30948, r30949, r30950, r30951, r30952, r30953, r30954, r30955, r30956, r30957, r30958, r30959, r30960, r30961, r30962, r30963, r30964, r30965, r30966, r30967, r30968, r30969, r30970, r30971, r30972, r30973, r30974, r30975, r30976, r30977, r30978, r30979;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r30918);
        mpfr_init(r30919);
        mpfr_init(r30920);
        mpfr_init_set_str(r30921, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r30922, "7", 10, MPFR_RNDN);
        mpfr_init(r30923);
        mpfr_init(r30924);
        mpfr_init_set_str(r30925, "1", 10, MPFR_RNDN);
        mpfr_init(r30926);
        mpfr_init(r30927);
        mpfr_init(r30928);
        mpfr_init(r30929);
        mpfr_init(r30930);
        mpfr_init(r30931);
        mpfr_init(r30932);
        mpfr_init_set_str(r30933, "-0.13857109526572012", 10, MPFR_RNDN);
        mpfr_init_set_str(r30934, "6", 10, MPFR_RNDN);
        mpfr_init(r30935);
        mpfr_init(r30936);
        mpfr_init(r30937);
        mpfr_init_set_str(r30938, "9.984369578019572e-06", 10, MPFR_RNDN);
        mpfr_init(r30939);
        mpfr_init(r30940);
        mpfr_init(r30941);
        mpfr_init_set_str(r30942, "12.507343278686905", 10, MPFR_RNDN);
        mpfr_init_set_str(r30943, "5", 10, MPFR_RNDN);
        mpfr_init(r30944);
        mpfr_init(r30945);
        mpfr_init(r30946);
        mpfr_init(r30947);
        mpfr_init_set_str(r30948, "-176.6150291621406", 10, MPFR_RNDN);
        mpfr_init_set_str(r30949, "4", 10, MPFR_RNDN);
        mpfr_init(r30950);
        mpfr_init(r30951);
        mpfr_init(r30952);
        mpfr_init_set_str(r30953, "676.5203681218851", 10, MPFR_RNDN);
        mpfr_init_set_str(r30954, "0", 10, MPFR_RNDN);
        mpfr_init(r30955);
        mpfr_init(r30956);
        mpfr_init_set_str(r30957, "0.9999999999998099", 10, MPFR_RNDN);
        mpfr_init(r30958);
        mpfr_init(r30959);
        mpfr_init_set_str(r30960, "-1259.1392167224028", 10, MPFR_RNDN);
        mpfr_init_set_str(r30961, "2", 10, MPFR_RNDN);
        mpfr_init(r30962);
        mpfr_init(r30963);
        mpfr_init(r30964);
        mpfr_init_set_str(r30965, "771.3234287776531", 10, MPFR_RNDN);
        mpfr_init_set_str(r30966, "3", 10, MPFR_RNDN);
        mpfr_init(r30967);
        mpfr_init(r30968);
        mpfr_init(r30969);
        mpfr_init(r30970);
        mpfr_init(r30971);
        mpfr_init(r30972);
        mpfr_init_set_str(r30973, "1.5056327351493116e-07", 10, MPFR_RNDN);
        mpfr_init_set_str(r30974, "8", 10, MPFR_RNDN);
        mpfr_init(r30975);
        mpfr_init(r30976);
        mpfr_init(r30977);
        mpfr_init(r30978);
        mpfr_init(r30979);
}

double f_dm(double z) {
        mpfr_const_pi(r30918, MPFR_RNDN);
        mpfr_add(r30919, r30918, r30918, MPFR_RNDN);
        mpfr_sqrt(r30920, r30919, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30923, r30921, r30922, MPFR_RNDN);
        mpfr_set_d(r30924, z, MPFR_RNDN);
        ;
        mpfr_sub(r30926, r30924, r30925, MPFR_RNDN);
        mpfr_add(r30927, r30923, r30926, MPFR_RNDN);
        mpfr_add(r30928, r30921, r30926, MPFR_RNDN);
        mpfr_pow(r30929, r30927, r30928, MPFR_RNDN);
        mpfr_exp(r30930, r30927, MPFR_RNDN);
        mpfr_div(r30931, r30929, r30930, MPFR_RNDN);
        mpfr_mul(r30932, r30920, r30931, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30935, r30925, r30934, MPFR_RNDN);
        mpfr_sub(r30936, r30924, r30935, MPFR_RNDN);
        mpfr_div(r30937, r30933, r30936, MPFR_RNDN);
        ;
        mpfr_add(r30939, r30926, r30922, MPFR_RNDN);
        mpfr_div(r30940, r30938, r30939, MPFR_RNDN);
        mpfr_add(r30941, r30937, r30940, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30944, r30924, r30943, MPFR_RNDN);
        mpfr_sub(r30945, r30944, r30925, MPFR_RNDN);
        mpfr_div(r30946, r30942, r30945, MPFR_RNDN);
        mpfr_add(r30947, r30941, r30946, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30950, r30924, r30949, MPFR_RNDN);
        mpfr_sub(r30951, r30950, r30925, MPFR_RNDN);
        mpfr_div(r30952, r30948, r30951, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30955, r30924, r30954, MPFR_RNDN);
        mpfr_div(r30956, r30953, r30955, MPFR_RNDN);
        ;
        mpfr_add(r30958, r30956, r30957, MPFR_RNDN);
        mpfr_add(r30959, r30952, r30958, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30962, r30925, r30961, MPFR_RNDN);
        mpfr_sub(r30963, r30924, r30962, MPFR_RNDN);
        mpfr_div(r30964, r30960, r30963, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30967, r30966, r30924, MPFR_RNDN);
        mpfr_sub(r30968, r30967, r30925, MPFR_RNDN);
        mpfr_div(r30969, r30965, r30968, MPFR_RNDN);
        mpfr_add(r30970, r30964, r30969, MPFR_RNDN);
        mpfr_add(r30971, r30959, r30970, MPFR_RNDN);
        mpfr_add(r30972, r30947, r30971, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30975, r30924, r30974, MPFR_RNDN);
        mpfr_sub(r30976, r30975, r30925, MPFR_RNDN);
        mpfr_div(r30977, r30973, r30976, MPFR_RNDN);
        mpfr_add(r30978, r30972, r30977, MPFR_RNDN);
        mpfr_mul(r30979, r30932, r30978, MPFR_RNDN);
        return mpfr_get_d(r30979, MPFR_RNDN);
}

