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

char *name = "Jmat.Real.dawson";

double f_if(float x) {
        float r27621 = 1;
        float r27622 = 0.1049934947;
        float r27623 = x;
        float r27624 = r27623 * r27623;
        float r27625 = r27622 * r27624;
        float r27626 = r27621 + r27625;
        float r27627 = 0.0424060604;
        float r27628 = r27624 * r27624;
        float r27629 = r27627 * r27628;
        float r27630 = r27626 + r27629;
        float r27631 = 0.0072644182;
        float r27632 = r27628 * r27624;
        float r27633 = r27631 * r27632;
        float r27634 = r27630 + r27633;
        float r27635 = 0.0005064034;
        float r27636 = r27632 * r27624;
        float r27637 = r27635 * r27636;
        float r27638 = r27634 + r27637;
        float r27639 = 0.0001789971;
        float r27640 = r27636 * r27624;
        float r27641 = r27639 * r27640;
        float r27642 = r27638 + r27641;
        float r27643 = 0.7715471019;
        float r27644 = r27643 * r27624;
        float r27645 = r27621 + r27644;
        float r27646 = 0.2909738639;
        float r27647 = r27646 * r27628;
        float r27648 = r27645 + r27647;
        float r27649 = 0.0694555761;
        float r27650 = r27649 * r27632;
        float r27651 = r27648 + r27650;
        float r27652 = 0.0140005442;
        float r27653 = r27652 * r27636;
        float r27654 = r27651 + r27653;
        float r27655 = 0.0008327945;
        float r27656 = r27655 * r27640;
        float r27657 = r27654 + r27656;
        float r27658 = 2;
        float r27659 = r27658 * r27639;
        float r27660 = r27640 * r27624;
        float r27661 = r27659 * r27660;
        float r27662 = r27657 + r27661;
        float r27663 = r27642 / r27662;
        float r27664 = r27663 * r27623;
        return r27664;
}

double f_id(double x) {
        double r27665 = 1;
        double r27666 = 0.1049934947;
        double r27667 = x;
        double r27668 = r27667 * r27667;
        double r27669 = r27666 * r27668;
        double r27670 = r27665 + r27669;
        double r27671 = 0.0424060604;
        double r27672 = r27668 * r27668;
        double r27673 = r27671 * r27672;
        double r27674 = r27670 + r27673;
        double r27675 = 0.0072644182;
        double r27676 = r27672 * r27668;
        double r27677 = r27675 * r27676;
        double r27678 = r27674 + r27677;
        double r27679 = 0.0005064034;
        double r27680 = r27676 * r27668;
        double r27681 = r27679 * r27680;
        double r27682 = r27678 + r27681;
        double r27683 = 0.0001789971;
        double r27684 = r27680 * r27668;
        double r27685 = r27683 * r27684;
        double r27686 = r27682 + r27685;
        double r27687 = 0.7715471019;
        double r27688 = r27687 * r27668;
        double r27689 = r27665 + r27688;
        double r27690 = 0.2909738639;
        double r27691 = r27690 * r27672;
        double r27692 = r27689 + r27691;
        double r27693 = 0.0694555761;
        double r27694 = r27693 * r27676;
        double r27695 = r27692 + r27694;
        double r27696 = 0.0140005442;
        double r27697 = r27696 * r27680;
        double r27698 = r27695 + r27697;
        double r27699 = 0.0008327945;
        double r27700 = r27699 * r27684;
        double r27701 = r27698 + r27700;
        double r27702 = 2;
        double r27703 = r27702 * r27683;
        double r27704 = r27684 * r27668;
        double r27705 = r27703 * r27704;
        double r27706 = r27701 + r27705;
        double r27707 = r27686 / r27706;
        double r27708 = r27707 * r27667;
        return r27708;
}


double f_of(float x) {
        float r27709 = 0.2514179000665375;
        float r27710 = x;
        float r27711 = 4;
        float r27712 = pow(r27710, r27711);
        float r27713 = r27709 / r27712;
        float r27714 = 0.15298196345929327;
        float r27715 = 6;
        float r27716 = pow(r27710, r27715);
        float r27717 = r27714 / r27716;
        float r27718 = r27713 + r27717;
        float r27719 = 0.5;
        float r27720 = r27719 / r27710;
        float r27721 = fma(r27718, r27710, r27720);
        float r27722 = -2.9491516833193935e-06;
        bool r27723 = r27721 <= r27722;
        float r27724 = r27710 * r27710;
        float r27725 = r27724 * r27724;
        float r27726 = r27725 * r27725;
        float r27727 = 0.0001789971;
        float r27728 = r27724 * r27727;
        float r27729 = r27726 * r27728;
        float r27730 = 0.0005064034;
        float r27731 = 0.0072644182;
        float r27732 = r27731 * r27724;
        float r27733 = r27732 * r27725;
        float r27734 = 0.0424060604;
        float r27735 = 0.1049934947;
        float r27736 = 1;
        float r27737 = fma(r27724, r27735, r27736);
        float r27738 = fma(r27734, r27725, r27737);
        float r27739 = r27733 + r27738;
        float r27740 = fma(r27730, r27726, r27739);
        float r27741 = r27729 + r27740;
        float r27742 = 0.0008327945;
        float r27743 = r27742 * r27724;
        float r27744 = r27743 * r27726;
        float r27745 = 2;
        float r27746 = r27727 * r27745;
        float r27747 = r27725 * r27746;
        float r27748 = r27726 * r27747;
        float r27749 = r27744 + r27748;
        float r27750 = 0.0140005442;
        float r27751 = 0.2909738639;
        float r27752 = 0.7715471019;
        float r27753 = r27752 * r27710;
        float r27754 = fma(r27753, r27710, r27736);
        float r27755 = fma(r27751, r27725, r27754);
        float r27756 = 0.0694555761;
        float r27757 = r27724 * r27756;
        float r27758 = r27725 * r27757;
        float r27759 = r27755 + r27758;
        float r27760 = fma(r27726, r27750, r27759);
        float r27761 = r27749 + r27760;
        float r27762 = r27741 / r27761;
        float r27763 = r27762 * r27710;
        float r27764 = 8.508551442347999e-05;
        bool r27765 = r27721 <= r27764;
        float r27766 = 0.265709700396151;
        float r27767 = r27766 * r27712;
        float r27768 = r27736 + r27767;
        float r27769 = 0.6665536072;
        float r27770 = pow(r27710, r27745);
        float r27771 = r27769 * r27770;
        float r27772 = r27768 - r27771;
        float r27773 = r27772 * r27710;
        float r27774 = r27765 ? r27721 : r27773;
        float r27775 = r27723 ? r27763 : r27774;
        return r27775;
}

double f_od(double x) {
        double r27776 = 0.2514179000665375;
        double r27777 = x;
        double r27778 = 4;
        double r27779 = pow(r27777, r27778);
        double r27780 = r27776 / r27779;
        double r27781 = 0.15298196345929327;
        double r27782 = 6;
        double r27783 = pow(r27777, r27782);
        double r27784 = r27781 / r27783;
        double r27785 = r27780 + r27784;
        double r27786 = 0.5;
        double r27787 = r27786 / r27777;
        double r27788 = fma(r27785, r27777, r27787);
        double r27789 = -2.9491516833193935e-06;
        bool r27790 = r27788 <= r27789;
        double r27791 = r27777 * r27777;
        double r27792 = r27791 * r27791;
        double r27793 = r27792 * r27792;
        double r27794 = 0.0001789971;
        double r27795 = r27791 * r27794;
        double r27796 = r27793 * r27795;
        double r27797 = 0.0005064034;
        double r27798 = 0.0072644182;
        double r27799 = r27798 * r27791;
        double r27800 = r27799 * r27792;
        double r27801 = 0.0424060604;
        double r27802 = 0.1049934947;
        double r27803 = 1;
        double r27804 = fma(r27791, r27802, r27803);
        double r27805 = fma(r27801, r27792, r27804);
        double r27806 = r27800 + r27805;
        double r27807 = fma(r27797, r27793, r27806);
        double r27808 = r27796 + r27807;
        double r27809 = 0.0008327945;
        double r27810 = r27809 * r27791;
        double r27811 = r27810 * r27793;
        double r27812 = 2;
        double r27813 = r27794 * r27812;
        double r27814 = r27792 * r27813;
        double r27815 = r27793 * r27814;
        double r27816 = r27811 + r27815;
        double r27817 = 0.0140005442;
        double r27818 = 0.2909738639;
        double r27819 = 0.7715471019;
        double r27820 = r27819 * r27777;
        double r27821 = fma(r27820, r27777, r27803);
        double r27822 = fma(r27818, r27792, r27821);
        double r27823 = 0.0694555761;
        double r27824 = r27791 * r27823;
        double r27825 = r27792 * r27824;
        double r27826 = r27822 + r27825;
        double r27827 = fma(r27793, r27817, r27826);
        double r27828 = r27816 + r27827;
        double r27829 = r27808 / r27828;
        double r27830 = r27829 * r27777;
        double r27831 = 8.508551442347999e-05;
        bool r27832 = r27788 <= r27831;
        double r27833 = 0.265709700396151;
        double r27834 = r27833 * r27779;
        double r27835 = r27803 + r27834;
        double r27836 = 0.6665536072;
        double r27837 = pow(r27777, r27812);
        double r27838 = r27836 * r27837;
        double r27839 = r27835 - r27838;
        double r27840 = r27839 * r27777;
        double r27841 = r27832 ? r27788 : r27840;
        double r27842 = r27790 ? r27830 : r27841;
        return r27842;
}

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 r27843, r27844, r27845, r27846, r27847, r27848, r27849, r27850, r27851, r27852, r27853, r27854, r27855, r27856, r27857, r27858, r27859, r27860, r27861, r27862, r27863, r27864, r27865, r27866, r27867, r27868, r27869, r27870, r27871, r27872, r27873, r27874, r27875, r27876, r27877, r27878, r27879, r27880, r27881, r27882, r27883, r27884, r27885, r27886;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27843, "1", 10, MPFR_RNDN);
        mpfr_init_set_str(r27844, "0.1049934947", 10, MPFR_RNDN);
        mpfr_init(r27845);
        mpfr_init(r27846);
        mpfr_init(r27847);
        mpfr_init(r27848);
        mpfr_init_set_str(r27849, "0.0424060604", 10, MPFR_RNDN);
        mpfr_init(r27850);
        mpfr_init(r27851);
        mpfr_init(r27852);
        mpfr_init_set_str(r27853, "0.0072644182", 10, MPFR_RNDN);
        mpfr_init(r27854);
        mpfr_init(r27855);
        mpfr_init(r27856);
        mpfr_init_set_str(r27857, "0.0005064034", 10, MPFR_RNDN);
        mpfr_init(r27858);
        mpfr_init(r27859);
        mpfr_init(r27860);
        mpfr_init_set_str(r27861, "0.0001789971", 10, MPFR_RNDN);
        mpfr_init(r27862);
        mpfr_init(r27863);
        mpfr_init(r27864);
        mpfr_init_set_str(r27865, "0.7715471019", 10, MPFR_RNDN);
        mpfr_init(r27866);
        mpfr_init(r27867);
        mpfr_init_set_str(r27868, "0.2909738639", 10, MPFR_RNDN);
        mpfr_init(r27869);
        mpfr_init(r27870);
        mpfr_init_set_str(r27871, "0.0694555761", 10, MPFR_RNDN);
        mpfr_init(r27872);
        mpfr_init(r27873);
        mpfr_init_set_str(r27874, "0.0140005442", 10, MPFR_RNDN);
        mpfr_init(r27875);
        mpfr_init(r27876);
        mpfr_init_set_str(r27877, "0.0008327945", 10, MPFR_RNDN);
        mpfr_init(r27878);
        mpfr_init(r27879);
        mpfr_init_set_str(r27880, "2", 10, MPFR_RNDN);
        mpfr_init(r27881);
        mpfr_init(r27882);
        mpfr_init(r27883);
        mpfr_init(r27884);
        mpfr_init(r27885);
        mpfr_init(r27886);
}

double f_im(double x) {
        ;
        ;
        mpfr_set_d(r27845, x, MPFR_RNDN);
        mpfr_mul(r27846, r27845, r27845, MPFR_RNDN);
        mpfr_mul(r27847, r27844, r27846, MPFR_RNDN);
        mpfr_add(r27848, r27843, r27847, MPFR_RNDN);
        ;
        mpfr_mul(r27850, r27846, r27846, MPFR_RNDN);
        mpfr_mul(r27851, r27849, r27850, MPFR_RNDN);
        mpfr_add(r27852, r27848, r27851, MPFR_RNDN);
        ;
        mpfr_mul(r27854, r27850, r27846, MPFR_RNDN);
        mpfr_mul(r27855, r27853, r27854, MPFR_RNDN);
        mpfr_add(r27856, r27852, r27855, MPFR_RNDN);
        ;
        mpfr_mul(r27858, r27854, r27846, MPFR_RNDN);
        mpfr_mul(r27859, r27857, r27858, MPFR_RNDN);
        mpfr_add(r27860, r27856, r27859, MPFR_RNDN);
        ;
        mpfr_mul(r27862, r27858, r27846, MPFR_RNDN);
        mpfr_mul(r27863, r27861, r27862, MPFR_RNDN);
        mpfr_add(r27864, r27860, r27863, MPFR_RNDN);
        ;
        mpfr_mul(r27866, r27865, r27846, MPFR_RNDN);
        mpfr_add(r27867, r27843, r27866, MPFR_RNDN);
        ;
        mpfr_mul(r27869, r27868, r27850, MPFR_RNDN);
        mpfr_add(r27870, r27867, r27869, MPFR_RNDN);
        ;
        mpfr_mul(r27872, r27871, r27854, MPFR_RNDN);
        mpfr_add(r27873, r27870, r27872, MPFR_RNDN);
        ;
        mpfr_mul(r27875, r27874, r27858, MPFR_RNDN);
        mpfr_add(r27876, r27873, r27875, MPFR_RNDN);
        ;
        mpfr_mul(r27878, r27877, r27862, MPFR_RNDN);
        mpfr_add(r27879, r27876, r27878, MPFR_RNDN);
        ;
        mpfr_mul(r27881, r27880, r27861, MPFR_RNDN);
        mpfr_mul(r27882, r27862, r27846, MPFR_RNDN);
        mpfr_mul(r27883, r27881, r27882, MPFR_RNDN);
        mpfr_add(r27884, r27879, r27883, MPFR_RNDN);
        mpfr_div(r27885, r27864, r27884, MPFR_RNDN);
        mpfr_mul(r27886, r27885, r27845, MPFR_RNDN);
        return mpfr_get_d(r27886, MPFR_RNDN);
}

static mpfr_t r27887, r27888, r27889, r27890, r27891, r27892, r27893, r27894, r27895, r27896, r27897, r27898, r27899, r27900, r27901, r27902, r27903, r27904, r27905, r27906, r27907, r27908, r27909, r27910, r27911, r27912, r27913, r27914, r27915, r27916, r27917, r27918, r27919, r27920, r27921, r27922, r27923, r27924, r27925, r27926, r27927, r27928, r27929, r27930, r27931, r27932, r27933, r27934, r27935, r27936, r27937, r27938, r27939, r27940, r27941, r27942, r27943, r27944, r27945, r27946, r27947, r27948, r27949, r27950, r27951, r27952, r27953;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27887, "0.2514179000665375", 10, MPFR_RNDN);
        mpfr_init(r27888);
        mpfr_init_set_str(r27889, "4", 10, MPFR_RNDN);
        mpfr_init(r27890);
        mpfr_init(r27891);
        mpfr_init_set_str(r27892, "0.15298196345929327", 10, MPFR_RNDN);
        mpfr_init_set_str(r27893, "6", 10, MPFR_RNDN);
        mpfr_init(r27894);
        mpfr_init(r27895);
        mpfr_init(r27896);
        mpfr_init_set_str(r27897, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27898);
        mpfr_init(r27899);
        mpfr_init_set_str(r27900, "-2.9491516833193935e-06", 10, MPFR_RNDN);
        mpfr_init(r27901);
        mpfr_init(r27902);
        mpfr_init(r27903);
        mpfr_init(r27904);
        mpfr_init_set_str(r27905, "0.0001789971", 10, MPFR_RNDN);
        mpfr_init(r27906);
        mpfr_init(r27907);
        mpfr_init_set_str(r27908, "0.0005064034", 10, MPFR_RNDN);
        mpfr_init_set_str(r27909, "0.0072644182", 10, MPFR_RNDN);
        mpfr_init(r27910);
        mpfr_init(r27911);
        mpfr_init_set_str(r27912, "0.0424060604", 10, MPFR_RNDN);
        mpfr_init_set_str(r27913, "0.1049934947", 10, MPFR_RNDN);
        mpfr_init_set_str(r27914, "1", 10, MPFR_RNDN);
        mpfr_init(r27915);
        mpfr_init(r27916);
        mpfr_init(r27917);
        mpfr_init(r27918);
        mpfr_init(r27919);
        mpfr_init_set_str(r27920, "0.0008327945", 10, MPFR_RNDN);
        mpfr_init(r27921);
        mpfr_init(r27922);
        mpfr_init_set_str(r27923, "2", 10, MPFR_RNDN);
        mpfr_init(r27924);
        mpfr_init(r27925);
        mpfr_init(r27926);
        mpfr_init(r27927);
        mpfr_init_set_str(r27928, "0.0140005442", 10, MPFR_RNDN);
        mpfr_init_set_str(r27929, "0.2909738639", 10, MPFR_RNDN);
        mpfr_init_set_str(r27930, "0.7715471019", 10, MPFR_RNDN);
        mpfr_init(r27931);
        mpfr_init(r27932);
        mpfr_init(r27933);
        mpfr_init_set_str(r27934, "0.0694555761", 10, MPFR_RNDN);
        mpfr_init(r27935);
        mpfr_init(r27936);
        mpfr_init(r27937);
        mpfr_init(r27938);
        mpfr_init(r27939);
        mpfr_init(r27940);
        mpfr_init(r27941);
        mpfr_init_set_str(r27942, "8.508551442347999e-05", 10, MPFR_RNDN);
        mpfr_init(r27943);
        mpfr_init_set_str(r27944, "0.265709700396151", 10, MPFR_RNDN);
        mpfr_init(r27945);
        mpfr_init(r27946);
        mpfr_init_set_str(r27947, "0.6665536072", 10, MPFR_RNDN);
        mpfr_init(r27948);
        mpfr_init(r27949);
        mpfr_init(r27950);
        mpfr_init(r27951);
        mpfr_init(r27952);
        mpfr_init(r27953);
}

double f_fm(double x) {
        ;
        mpfr_set_d(r27888, x, MPFR_RNDN);
        ;
        mpfr_pow(r27890, r27888, r27889, MPFR_RNDN);
        mpfr_div(r27891, r27887, r27890, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27894, r27888, r27893, MPFR_RNDN);
        mpfr_div(r27895, r27892, r27894, MPFR_RNDN);
        mpfr_add(r27896, r27891, r27895, MPFR_RNDN);
        ;
        mpfr_div(r27898, r27897, r27888, MPFR_RNDN);
        mpfr_fma(r27899, r27896, r27888, r27898, MPFR_RNDN);
        ;
        mpfr_set_si(r27901, mpfr_cmp(r27899, r27900) <= 0, MPFR_RNDN);
        mpfr_mul(r27902, r27888, r27888, MPFR_RNDN);
        mpfr_mul(r27903, r27902, r27902, MPFR_RNDN);
        mpfr_mul(r27904, r27903, r27903, MPFR_RNDN);
        ;
        mpfr_mul(r27906, r27902, r27905, MPFR_RNDN);
        mpfr_mul(r27907, r27904, r27906, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r27910, r27909, r27902, MPFR_RNDN);
        mpfr_mul(r27911, r27910, r27903, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_fma(r27915, r27902, r27913, r27914, MPFR_RNDN);
        mpfr_fma(r27916, r27912, r27903, r27915, MPFR_RNDN);
        mpfr_add(r27917, r27911, r27916, MPFR_RNDN);
        mpfr_fma(r27918, r27908, r27904, r27917, MPFR_RNDN);
        mpfr_add(r27919, r27907, r27918, MPFR_RNDN);
        ;
        mpfr_mul(r27921, r27920, r27902, MPFR_RNDN);
        mpfr_mul(r27922, r27921, r27904, MPFR_RNDN);
        ;
        mpfr_mul(r27924, r27905, r27923, MPFR_RNDN);
        mpfr_mul(r27925, r27903, r27924, MPFR_RNDN);
        mpfr_mul(r27926, r27904, r27925, MPFR_RNDN);
        mpfr_add(r27927, r27922, r27926, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_mul(r27931, r27930, r27888, MPFR_RNDN);
        mpfr_fma(r27932, r27931, r27888, r27914, MPFR_RNDN);
        mpfr_fma(r27933, r27929, r27903, r27932, MPFR_RNDN);
        ;
        mpfr_mul(r27935, r27902, r27934, MPFR_RNDN);
        mpfr_mul(r27936, r27903, r27935, MPFR_RNDN);
        mpfr_add(r27937, r27933, r27936, MPFR_RNDN);
        mpfr_fma(r27938, r27904, r27928, r27937, MPFR_RNDN);
        mpfr_add(r27939, r27927, r27938, MPFR_RNDN);
        mpfr_div(r27940, r27919, r27939, MPFR_RNDN);
        mpfr_mul(r27941, r27940, r27888, MPFR_RNDN);
        ;
        mpfr_set_si(r27943, mpfr_cmp(r27899, r27942) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r27945, r27944, r27890, MPFR_RNDN);
        mpfr_add(r27946, r27914, r27945, MPFR_RNDN);
        ;
        mpfr_pow(r27948, r27888, r27923, MPFR_RNDN);
        mpfr_mul(r27949, r27947, r27948, MPFR_RNDN);
        mpfr_sub(r27950, r27946, r27949, MPFR_RNDN);
        mpfr_mul(r27951, r27950, r27888, MPFR_RNDN);
        if (mpfr_get_si(r27943, MPFR_RNDN)) { mpfr_set(r27952, r27899, MPFR_RNDN); } else { mpfr_set(r27952, r27951, MPFR_RNDN); };
        if (mpfr_get_si(r27901, MPFR_RNDN)) { mpfr_set(r27953, r27941, MPFR_RNDN); } else { mpfr_set(r27953, r27952, MPFR_RNDN); };
        return mpfr_get_d(r27953, MPFR_RNDN);
}

static mpfr_t r27954, r27955, r27956, r27957, r27958, r27959, r27960, r27961, r27962, r27963, r27964, r27965, r27966, r27967, r27968, r27969, r27970, r27971, r27972, r27973, r27974, r27975, r27976, r27977, r27978, r27979, r27980, r27981, r27982, r27983, r27984, r27985, r27986, r27987, r27988, r27989, r27990, r27991, r27992, r27993, r27994, r27995, r27996, r27997, r27998, r27999, r28000, r28001, r28002, r28003, r28004, r28005, r28006, r28007, r28008, r28009, r28010, r28011, r28012, r28013, r28014, r28015, r28016, r28017, r28018, r28019, r28020;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(592);
        mpfr_init_set_str(r27954, "0.2514179000665375", 10, MPFR_RNDN);
        mpfr_init(r27955);
        mpfr_init_set_str(r27956, "4", 10, MPFR_RNDN);
        mpfr_init(r27957);
        mpfr_init(r27958);
        mpfr_init_set_str(r27959, "0.15298196345929327", 10, MPFR_RNDN);
        mpfr_init_set_str(r27960, "6", 10, MPFR_RNDN);
        mpfr_init(r27961);
        mpfr_init(r27962);
        mpfr_init(r27963);
        mpfr_init_set_str(r27964, "0.5", 10, MPFR_RNDN);
        mpfr_init(r27965);
        mpfr_init(r27966);
        mpfr_init_set_str(r27967, "-2.9491516833193935e-06", 10, MPFR_RNDN);
        mpfr_init(r27968);
        mpfr_init(r27969);
        mpfr_init(r27970);
        mpfr_init(r27971);
        mpfr_init_set_str(r27972, "0.0001789971", 10, MPFR_RNDN);
        mpfr_init(r27973);
        mpfr_init(r27974);
        mpfr_init_set_str(r27975, "0.0005064034", 10, MPFR_RNDN);
        mpfr_init_set_str(r27976, "0.0072644182", 10, MPFR_RNDN);
        mpfr_init(r27977);
        mpfr_init(r27978);
        mpfr_init_set_str(r27979, "0.0424060604", 10, MPFR_RNDN);
        mpfr_init_set_str(r27980, "0.1049934947", 10, MPFR_RNDN);
        mpfr_init_set_str(r27981, "1", 10, MPFR_RNDN);
        mpfr_init(r27982);
        mpfr_init(r27983);
        mpfr_init(r27984);
        mpfr_init(r27985);
        mpfr_init(r27986);
        mpfr_init_set_str(r27987, "0.0008327945", 10, MPFR_RNDN);
        mpfr_init(r27988);
        mpfr_init(r27989);
        mpfr_init_set_str(r27990, "2", 10, MPFR_RNDN);
        mpfr_init(r27991);
        mpfr_init(r27992);
        mpfr_init(r27993);
        mpfr_init(r27994);
        mpfr_init_set_str(r27995, "0.0140005442", 10, MPFR_RNDN);
        mpfr_init_set_str(r27996, "0.2909738639", 10, MPFR_RNDN);
        mpfr_init_set_str(r27997, "0.7715471019", 10, MPFR_RNDN);
        mpfr_init(r27998);
        mpfr_init(r27999);
        mpfr_init(r28000);
        mpfr_init_set_str(r28001, "0.0694555761", 10, MPFR_RNDN);
        mpfr_init(r28002);
        mpfr_init(r28003);
        mpfr_init(r28004);
        mpfr_init(r28005);
        mpfr_init(r28006);
        mpfr_init(r28007);
        mpfr_init(r28008);
        mpfr_init_set_str(r28009, "8.508551442347999e-05", 10, MPFR_RNDN);
        mpfr_init(r28010);
        mpfr_init_set_str(r28011, "0.265709700396151", 10, MPFR_RNDN);
        mpfr_init(r28012);
        mpfr_init(r28013);
        mpfr_init_set_str(r28014, "0.6665536072", 10, MPFR_RNDN);
        mpfr_init(r28015);
        mpfr_init(r28016);
        mpfr_init(r28017);
        mpfr_init(r28018);
        mpfr_init(r28019);
        mpfr_init(r28020);
}

double f_dm(double x) {
        ;
        mpfr_set_d(r27955, x, MPFR_RNDN);
        ;
        mpfr_pow(r27957, r27955, r27956, MPFR_RNDN);
        mpfr_div(r27958, r27954, r27957, MPFR_RNDN);
        ;
        ;
        mpfr_pow(r27961, r27955, r27960, MPFR_RNDN);
        mpfr_div(r27962, r27959, r27961, MPFR_RNDN);
        mpfr_add(r27963, r27958, r27962, MPFR_RNDN);
        ;
        mpfr_div(r27965, r27964, r27955, MPFR_RNDN);
        mpfr_fma(r27966, r27963, r27955, r27965, MPFR_RNDN);
        ;
        mpfr_set_si(r27968, mpfr_cmp(r27966, r27967) <= 0, MPFR_RNDN);
        mpfr_mul(r27969, r27955, r27955, MPFR_RNDN);
        mpfr_mul(r27970, r27969, r27969, MPFR_RNDN);
        mpfr_mul(r27971, r27970, r27970, MPFR_RNDN);
        ;
        mpfr_mul(r27973, r27969, r27972, MPFR_RNDN);
        mpfr_mul(r27974, r27971, r27973, MPFR_RNDN);
        ;
        ;
        mpfr_mul(r27977, r27976, r27969, MPFR_RNDN);
        mpfr_mul(r27978, r27977, r27970, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_fma(r27982, r27969, r27980, r27981, MPFR_RNDN);
        mpfr_fma(r27983, r27979, r27970, r27982, MPFR_RNDN);
        mpfr_add(r27984, r27978, r27983, MPFR_RNDN);
        mpfr_fma(r27985, r27975, r27971, r27984, MPFR_RNDN);
        mpfr_add(r27986, r27974, r27985, MPFR_RNDN);
        ;
        mpfr_mul(r27988, r27987, r27969, MPFR_RNDN);
        mpfr_mul(r27989, r27988, r27971, MPFR_RNDN);
        ;
        mpfr_mul(r27991, r27972, r27990, MPFR_RNDN);
        mpfr_mul(r27992, r27970, r27991, MPFR_RNDN);
        mpfr_mul(r27993, r27971, r27992, MPFR_RNDN);
        mpfr_add(r27994, r27989, r27993, MPFR_RNDN);
        ;
        ;
        ;
        mpfr_mul(r27998, r27997, r27955, MPFR_RNDN);
        mpfr_fma(r27999, r27998, r27955, r27981, MPFR_RNDN);
        mpfr_fma(r28000, r27996, r27970, r27999, MPFR_RNDN);
        ;
        mpfr_mul(r28002, r27969, r28001, MPFR_RNDN);
        mpfr_mul(r28003, r27970, r28002, MPFR_RNDN);
        mpfr_add(r28004, r28000, r28003, MPFR_RNDN);
        mpfr_fma(r28005, r27971, r27995, r28004, MPFR_RNDN);
        mpfr_add(r28006, r27994, r28005, MPFR_RNDN);
        mpfr_div(r28007, r27986, r28006, MPFR_RNDN);
        mpfr_mul(r28008, r28007, r27955, MPFR_RNDN);
        ;
        mpfr_set_si(r28010, mpfr_cmp(r27966, r28009) <= 0, MPFR_RNDN);
        ;
        mpfr_mul(r28012, r28011, r27957, MPFR_RNDN);
        mpfr_add(r28013, r27981, r28012, MPFR_RNDN);
        ;
        mpfr_pow(r28015, r27955, r27990, MPFR_RNDN);
        mpfr_mul(r28016, r28014, r28015, MPFR_RNDN);
        mpfr_sub(r28017, r28013, r28016, MPFR_RNDN);
        mpfr_mul(r28018, r28017, r27955, MPFR_RNDN);
        if (mpfr_get_si(r28010, MPFR_RNDN)) { mpfr_set(r28019, r27966, MPFR_RNDN); } else { mpfr_set(r28019, r28018, MPFR_RNDN); };
        if (mpfr_get_si(r27968, MPFR_RNDN)) { mpfr_set(r28020, r28008, MPFR_RNDN); } else { mpfr_set(r28020, r28019, MPFR_RNDN); };
        return mpfr_get_d(r28020, MPFR_RNDN);
}

