#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 r30331 = atan2(1.0, 0.0);
        float r30332 = 2;
        float r30333 = r30331 * r30332;
        float r30334 = sqrt(r30333);
        float r30335 = z;
        float r30336 = 1;
        float r30337 = r30335 - r30336;
        float r30338 = 7;
        float r30339 = r30337 + r30338;
        float r30340 = 0.5;
        float r30341 = r30339 + r30340;
        float r30342 = r30337 + r30340;
        float r30343 = pow(r30341, r30342);
        float r30344 = r30334 * r30343;
        float r30345 = -r30341;
        float r30346 = exp(r30345);
        float r30347 = r30344 * r30346;
        float r30348 = 0.9999999999998099;
        float r30349 = 676.5203681218851;
        float r30350 = r30337 + r30336;
        float r30351 = r30349 / r30350;
        float r30352 = r30348 + r30351;
        float r30353 = -1259.1392167224028;
        float r30354 = r30337 + r30332;
        float r30355 = r30353 / r30354;
        float r30356 = r30352 + r30355;
        float r30357 = 771.3234287776531;
        float r30358 = 3;
        float r30359 = r30337 + r30358;
        float r30360 = r30357 / r30359;
        float r30361 = r30356 + r30360;
        float r30362 = -176.6150291621406;
        float r30363 = 4;
        float r30364 = r30337 + r30363;
        float r30365 = r30362 / r30364;
        float r30366 = r30361 + r30365;
        float r30367 = 12.507343278686905;
        float r30368 = 5;
        float r30369 = r30337 + r30368;
        float r30370 = r30367 / r30369;
        float r30371 = r30366 + r30370;
        float r30372 = -0.13857109526572012;
        float r30373 = 6;
        float r30374 = r30337 + r30373;
        float r30375 = r30372 / r30374;
        float r30376 = r30371 + r30375;
        float r30377 = 9.984369578019572e-06;
        float r30378 = r30377 / r30339;
        float r30379 = r30376 + r30378;
        float r30380 = 1.5056327351493116e-07;
        float r30381 = 8;
        float r30382 = r30337 + r30381;
        float r30383 = r30380 / r30382;
        float r30384 = r30379 + r30383;
        float r30385 = r30347 * r30384;
        return r30385;
}

double f_id(double z) {
        double r30386 = atan2(1.0, 0.0);
        double r30387 = 2;
        double r30388 = r30386 * r30387;
        double r30389 = sqrt(r30388);
        double r30390 = z;
        double r30391 = 1;
        double r30392 = r30390 - r30391;
        double r30393 = 7;
        double r30394 = r30392 + r30393;
        double r30395 = 0.5;
        double r30396 = r30394 + r30395;
        double r30397 = r30392 + r30395;
        double r30398 = pow(r30396, r30397);
        double r30399 = r30389 * r30398;
        double r30400 = -r30396;
        double r30401 = exp(r30400);
        double r30402 = r30399 * r30401;
        double r30403 = 0.9999999999998099;
        double r30404 = 676.5203681218851;
        double r30405 = r30392 + r30391;
        double r30406 = r30404 / r30405;
        double r30407 = r30403 + r30406;
        double r30408 = -1259.1392167224028;
        double r30409 = r30392 + r30387;
        double r30410 = r30408 / r30409;
        double r30411 = r30407 + r30410;
        double r30412 = 771.3234287776531;
        double r30413 = 3;
        double r30414 = r30392 + r30413;
        double r30415 = r30412 / r30414;
        double r30416 = r30411 + r30415;
        double r30417 = -176.6150291621406;
        double r30418 = 4;
        double r30419 = r30392 + r30418;
        double r30420 = r30417 / r30419;
        double r30421 = r30416 + r30420;
        double r30422 = 12.507343278686905;
        double r30423 = 5;
        double r30424 = r30392 + r30423;
        double r30425 = r30422 / r30424;
        double r30426 = r30421 + r30425;
        double r30427 = -0.13857109526572012;
        double r30428 = 6;
        double r30429 = r30392 + r30428;
        double r30430 = r30427 / r30429;
        double r30431 = r30426 + r30430;
        double r30432 = 9.984369578019572e-06;
        double r30433 = r30432 / r30394;
        double r30434 = r30431 + r30433;
        double r30435 = 1.5056327351493116e-07;
        double r30436 = 8;
        double r30437 = r30392 + r30436;
        double r30438 = r30435 / r30437;
        double r30439 = r30434 + r30438;
        double r30440 = r30402 * r30439;
        return r30440;
}


double f_of(float z) {
        float r30441 = atan2(1.0, 0.0);
        float r30442 = r30441 + r30441;
        float r30443 = sqrt(r30442);
        float r30444 = 0.5;
        float r30445 = 7;
        float r30446 = r30444 + r30445;
        float r30447 = z;
        float r30448 = 1;
        float r30449 = r30447 - r30448;
        float r30450 = r30446 + r30449;
        float r30451 = r30444 + r30449;
        float r30452 = pow(r30450, r30451);
        float r30453 = exp(r30450);
        float r30454 = r30452 / r30453;
        float r30455 = r30443 * r30454;
        float r30456 = -0.13857109526572012;
        float r30457 = 6;
        float r30458 = r30448 - r30457;
        float r30459 = r30447 - r30458;
        float r30460 = r30456 / r30459;
        float r30461 = 9.984369578019572e-06;
        float r30462 = r30449 + r30445;
        float r30463 = r30461 / r30462;
        float r30464 = r30460 + r30463;
        float r30465 = 12.507343278686905;
        float r30466 = 5;
        float r30467 = r30447 + r30466;
        float r30468 = r30467 - r30448;
        float r30469 = r30465 / r30468;
        float r30470 = r30464 + r30469;
        float r30471 = -176.6150291621406;
        float r30472 = 4;
        float r30473 = r30447 + r30472;
        float r30474 = r30473 - r30448;
        float r30475 = r30471 / r30474;
        float r30476 = 676.5203681218851;
        float r30477 = 0;
        float r30478 = r30447 - r30477;
        float r30479 = r30476 / r30478;
        float r30480 = 0.9999999999998099;
        float r30481 = r30479 + r30480;
        float r30482 = r30475 + r30481;
        float r30483 = -1259.1392167224028;
        float r30484 = 2;
        float r30485 = r30448 - r30484;
        float r30486 = r30447 - r30485;
        float r30487 = r30483 / r30486;
        float r30488 = 771.3234287776531;
        float r30489 = 3;
        float r30490 = r30489 + r30447;
        float r30491 = r30490 - r30448;
        float r30492 = r30488 / r30491;
        float r30493 = r30487 + r30492;
        float r30494 = r30482 + r30493;
        float r30495 = r30470 + r30494;
        float r30496 = 1.5056327351493116e-07;
        float r30497 = 8;
        float r30498 = r30447 + r30497;
        float r30499 = r30498 - r30448;
        float r30500 = r30496 / r30499;
        float r30501 = r30495 + r30500;
        float r30502 = r30455 * r30501;
        return r30502;
}

double f_od(double z) {
        double r30503 = atan2(1.0, 0.0);
        double r30504 = r30503 + r30503;
        double r30505 = sqrt(r30504);
        double r30506 = 0.5;
        double r30507 = 7;
        double r30508 = r30506 + r30507;
        double r30509 = z;
        double r30510 = 1;
        double r30511 = r30509 - r30510;
        double r30512 = r30508 + r30511;
        double r30513 = r30506 + r30511;
        double r30514 = pow(r30512, r30513);
        double r30515 = exp(r30512);
        double r30516 = r30514 / r30515;
        double r30517 = r30505 * r30516;
        double r30518 = -0.13857109526572012;
        double r30519 = 6;
        double r30520 = r30510 - r30519;
        double r30521 = r30509 - r30520;
        double r30522 = r30518 / r30521;
        double r30523 = 9.984369578019572e-06;
        double r30524 = r30511 + r30507;
        double r30525 = r30523 / r30524;
        double r30526 = r30522 + r30525;
        double r30527 = 12.507343278686905;
        double r30528 = 5;
        double r30529 = r30509 + r30528;
        double r30530 = r30529 - r30510;
        double r30531 = r30527 / r30530;
        double r30532 = r30526 + r30531;
        double r30533 = -176.6150291621406;
        double r30534 = 4;
        double r30535 = r30509 + r30534;
        double r30536 = r30535 - r30510;
        double r30537 = r30533 / r30536;
        double r30538 = 676.5203681218851;
        double r30539 = 0;
        double r30540 = r30509 - r30539;
        double r30541 = r30538 / r30540;
        double r30542 = 0.9999999999998099;
        double r30543 = r30541 + r30542;
        double r30544 = r30537 + r30543;
        double r30545 = -1259.1392167224028;
        double r30546 = 2;
        double r30547 = r30510 - r30546;
        double r30548 = r30509 - r30547;
        double r30549 = r30545 / r30548;
        double r30550 = 771.3234287776531;
        double r30551 = 3;
        double r30552 = r30551 + r30509;
        double r30553 = r30552 - r30510;
        double r30554 = r30550 / r30553;
        double r30555 = r30549 + r30554;
        double r30556 = r30544 + r30555;
        double r30557 = r30532 + r30556;
        double r30558 = 1.5056327351493116e-07;
        double r30559 = 8;
        double r30560 = r30509 + r30559;
        double r30561 = r30560 - r30510;
        double r30562 = r30558 / r30561;
        double r30563 = r30557 + r30562;
        double r30564 = r30517 * r30563;
        return r30564;
}

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 r30565, r30566, r30567, r30568, r30569, r30570, r30571, r30572, r30573, r30574, r30575, r30576, r30577, r30578, r30579, r30580, r30581, r30582, r30583, r30584, r30585, r30586, r30587, r30588, r30589, r30590, r30591, r30592, r30593, r30594, r30595, r30596, r30597, r30598, r30599, r30600, r30601, r30602, r30603, r30604, r30605, r30606, r30607, r30608, r30609, r30610, r30611, r30612, r30613, r30614, r30615, r30616, r30617, r30618, r30619;

void setup_mpfr_f_im() {
        mpfr_set_default_prec(1424);
        mpfr_init(r30565);
        mpfr_init_set_str(r30566, "2", 10, MPFR_RNDN);
        mpfr_init(r30567);
        mpfr_init(r30568);
        mpfr_init(r30569);
        mpfr_init_set_str(r30570, "1", 10, MPFR_RNDN);
        mpfr_init(r30571);
        mpfr_init_set_str(r30572, "7", 10, MPFR_RNDN);
        mpfr_init(r30573);
        mpfr_init_set_str(r30574, "0.5", 10, MPFR_RNDN);
        mpfr_init(r30575);
        mpfr_init(r30576);
        mpfr_init(r30577);
        mpfr_init(r30578);
        mpfr_init(r30579);
        mpfr_init(r30580);
        mpfr_init(r30581);
        mpfr_init_set_str(r30582, "0.9999999999998099", 10, MPFR_RNDN);
        mpfr_init_set_str(r30583, "676.5203681218851", 10, MPFR_RNDN);
        mpfr_init(r30584);
        mpfr_init(r30585);
        mpfr_init(r30586);
        mpfr_init_set_str(r30587, "-1259.1392167224028", 10, MPFR_RNDN);
        mpfr_init(r30588);
        mpfr_init(r30589);
        mpfr_init(r30590);
        mpfr_init_set_str(r30591, "771.3234287776531", 10, MPFR_RNDN);
        mpfr_init_set_str(r30592, "3", 10, MPFR_RNDN);
        mpfr_init(r30593);
        mpfr_init(r30594);
        mpfr_init(r30595);
        mpfr_init_set_str(r30596, "-176.6150291621406", 10, MPFR_RNDN);
        mpfr_init_set_str(r30597, "4", 10, MPFR_RNDN);
        mpfr_init(r30598);
        mpfr_init(r30599);
        mpfr_init(r30600);
        mpfr_init_set_str(r30601, "12.507343278686905", 10, MPFR_RNDN);
        mpfr_init_set_str(r30602, "5", 10, MPFR_RNDN);
        mpfr_init(r30603);
        mpfr_init(r30604);
        mpfr_init(r30605);
        mpfr_init_set_str(r30606, "-0.13857109526572012", 10, MPFR_RNDN);
        mpfr_init_set_str(r30607, "6", 10, MPFR_RNDN);
        mpfr_init(r30608);
        mpfr_init(r30609);
        mpfr_init(r30610);
        mpfr_init_set_str(r30611, "9.984369578019572e-06", 10, MPFR_RNDN);
        mpfr_init(r30612);
        mpfr_init(r30613);
        mpfr_init_set_str(r30614, "1.5056327351493116e-07", 10, MPFR_RNDN);
        mpfr_init_set_str(r30615, "8", 10, MPFR_RNDN);
        mpfr_init(r30616);
        mpfr_init(r30617);
        mpfr_init(r30618);
        mpfr_init(r30619);
}

double f_im(double z) {
        mpfr_const_pi(r30565, MPFR_RNDN);
        ;
        mpfr_mul(r30567, r30565, r30566, MPFR_RNDN);
        mpfr_sqrt(r30568, r30567, MPFR_RNDN);
        mpfr_set_d(r30569, z, MPFR_RNDN);
        ;
        mpfr_sub(r30571, r30569, r30570, MPFR_RNDN);
        ;
        mpfr_add(r30573, r30571, r30572, MPFR_RNDN);
        ;
        mpfr_add(r30575, r30573, r30574, MPFR_RNDN);
        mpfr_add(r30576, r30571, r30574, MPFR_RNDN);
        mpfr_pow(r30577, r30575, r30576, MPFR_RNDN);
        mpfr_mul(r30578, r30568, r30577, MPFR_RNDN);
        mpfr_neg(r30579, r30575, MPFR_RNDN);
        mpfr_exp(r30580, r30579, MPFR_RNDN);
        mpfr_mul(r30581, r30578, r30580, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30584, r30571, r30570, MPFR_RNDN);
        mpfr_div(r30585, r30583, r30584, MPFR_RNDN);
        mpfr_add(r30586, r30582, r30585, MPFR_RNDN);
        ;
        mpfr_add(r30588, r30571, r30566, MPFR_RNDN);
        mpfr_div(r30589, r30587, r30588, MPFR_RNDN);
        mpfr_add(r30590, r30586, r30589, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30593, r30571, r30592, MPFR_RNDN);
        mpfr_div(r30594, r30591, r30593, MPFR_RNDN);
        mpfr_add(r30595, r30590, r30594, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30598, r30571, r30597, MPFR_RNDN);
        mpfr_div(r30599, r30596, r30598, MPFR_RNDN);
        mpfr_add(r30600, r30595, r30599, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30603, r30571, r30602, MPFR_RNDN);
        mpfr_div(r30604, r30601, r30603, MPFR_RNDN);
        mpfr_add(r30605, r30600, r30604, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30608, r30571, r30607, MPFR_RNDN);
        mpfr_div(r30609, r30606, r30608, MPFR_RNDN);
        mpfr_add(r30610, r30605, r30609, MPFR_RNDN);
        ;
        mpfr_div(r30612, r30611, r30573, MPFR_RNDN);
        mpfr_add(r30613, r30610, r30612, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30616, r30571, r30615, MPFR_RNDN);
        mpfr_div(r30617, r30614, r30616, MPFR_RNDN);
        mpfr_add(r30618, r30613, r30617, MPFR_RNDN);
        mpfr_mul(r30619, r30581, r30618, MPFR_RNDN);
        return mpfr_get_d(r30619, MPFR_RNDN);
}

static mpfr_t r30620, r30621, r30622, r30623, r30624, r30625, r30626, r30627, r30628, r30629, r30630, r30631, r30632, r30633, r30634, r30635, r30636, r30637, r30638, r30639, r30640, r30641, r30642, r30643, r30644, r30645, r30646, r30647, r30648, r30649, r30650, r30651, r30652, r30653, r30654, r30655, r30656, r30657, r30658, r30659, r30660, r30661, r30662, r30663, r30664, r30665, r30666, r30667, r30668, r30669, r30670, r30671, r30672, r30673, r30674, r30675, r30676, r30677, r30678, r30679, r30680, r30681;

void setup_mpfr_f_fm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r30620);
        mpfr_init(r30621);
        mpfr_init(r30622);
        mpfr_init_set_str(r30623, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r30624, "7", 10, MPFR_RNDN);
        mpfr_init(r30625);
        mpfr_init(r30626);
        mpfr_init_set_str(r30627, "1", 10, MPFR_RNDN);
        mpfr_init(r30628);
        mpfr_init(r30629);
        mpfr_init(r30630);
        mpfr_init(r30631);
        mpfr_init(r30632);
        mpfr_init(r30633);
        mpfr_init(r30634);
        mpfr_init_set_str(r30635, "-0.13857109526572012", 10, MPFR_RNDN);
        mpfr_init_set_str(r30636, "6", 10, MPFR_RNDN);
        mpfr_init(r30637);
        mpfr_init(r30638);
        mpfr_init(r30639);
        mpfr_init_set_str(r30640, "9.984369578019572e-06", 10, MPFR_RNDN);
        mpfr_init(r30641);
        mpfr_init(r30642);
        mpfr_init(r30643);
        mpfr_init_set_str(r30644, "12.507343278686905", 10, MPFR_RNDN);
        mpfr_init_set_str(r30645, "5", 10, MPFR_RNDN);
        mpfr_init(r30646);
        mpfr_init(r30647);
        mpfr_init(r30648);
        mpfr_init(r30649);
        mpfr_init_set_str(r30650, "-176.6150291621406", 10, MPFR_RNDN);
        mpfr_init_set_str(r30651, "4", 10, MPFR_RNDN);
        mpfr_init(r30652);
        mpfr_init(r30653);
        mpfr_init(r30654);
        mpfr_init_set_str(r30655, "676.5203681218851", 10, MPFR_RNDN);
        mpfr_init_set_str(r30656, "0", 10, MPFR_RNDN);
        mpfr_init(r30657);
        mpfr_init(r30658);
        mpfr_init_set_str(r30659, "0.9999999999998099", 10, MPFR_RNDN);
        mpfr_init(r30660);
        mpfr_init(r30661);
        mpfr_init_set_str(r30662, "-1259.1392167224028", 10, MPFR_RNDN);
        mpfr_init_set_str(r30663, "2", 10, MPFR_RNDN);
        mpfr_init(r30664);
        mpfr_init(r30665);
        mpfr_init(r30666);
        mpfr_init_set_str(r30667, "771.3234287776531", 10, MPFR_RNDN);
        mpfr_init_set_str(r30668, "3", 10, MPFR_RNDN);
        mpfr_init(r30669);
        mpfr_init(r30670);
        mpfr_init(r30671);
        mpfr_init(r30672);
        mpfr_init(r30673);
        mpfr_init(r30674);
        mpfr_init_set_str(r30675, "1.5056327351493116e-07", 10, MPFR_RNDN);
        mpfr_init_set_str(r30676, "8", 10, MPFR_RNDN);
        mpfr_init(r30677);
        mpfr_init(r30678);
        mpfr_init(r30679);
        mpfr_init(r30680);
        mpfr_init(r30681);
}

double f_fm(double z) {
        mpfr_const_pi(r30620, MPFR_RNDN);
        mpfr_add(r30621, r30620, r30620, MPFR_RNDN);
        mpfr_sqrt(r30622, r30621, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30625, r30623, r30624, MPFR_RNDN);
        mpfr_set_d(r30626, z, MPFR_RNDN);
        ;
        mpfr_sub(r30628, r30626, r30627, MPFR_RNDN);
        mpfr_add(r30629, r30625, r30628, MPFR_RNDN);
        mpfr_add(r30630, r30623, r30628, MPFR_RNDN);
        mpfr_pow(r30631, r30629, r30630, MPFR_RNDN);
        mpfr_exp(r30632, r30629, MPFR_RNDN);
        mpfr_div(r30633, r30631, r30632, MPFR_RNDN);
        mpfr_mul(r30634, r30622, r30633, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30637, r30627, r30636, MPFR_RNDN);
        mpfr_sub(r30638, r30626, r30637, MPFR_RNDN);
        mpfr_div(r30639, r30635, r30638, MPFR_RNDN);
        ;
        mpfr_add(r30641, r30628, r30624, MPFR_RNDN);
        mpfr_div(r30642, r30640, r30641, MPFR_RNDN);
        mpfr_add(r30643, r30639, r30642, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30646, r30626, r30645, MPFR_RNDN);
        mpfr_sub(r30647, r30646, r30627, MPFR_RNDN);
        mpfr_div(r30648, r30644, r30647, MPFR_RNDN);
        mpfr_add(r30649, r30643, r30648, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30652, r30626, r30651, MPFR_RNDN);
        mpfr_sub(r30653, r30652, r30627, MPFR_RNDN);
        mpfr_div(r30654, r30650, r30653, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30657, r30626, r30656, MPFR_RNDN);
        mpfr_div(r30658, r30655, r30657, MPFR_RNDN);
        ;
        mpfr_add(r30660, r30658, r30659, MPFR_RNDN);
        mpfr_add(r30661, r30654, r30660, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30664, r30627, r30663, MPFR_RNDN);
        mpfr_sub(r30665, r30626, r30664, MPFR_RNDN);
        mpfr_div(r30666, r30662, r30665, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30669, r30668, r30626, MPFR_RNDN);
        mpfr_sub(r30670, r30669, r30627, MPFR_RNDN);
        mpfr_div(r30671, r30667, r30670, MPFR_RNDN);
        mpfr_add(r30672, r30666, r30671, MPFR_RNDN);
        mpfr_add(r30673, r30661, r30672, MPFR_RNDN);
        mpfr_add(r30674, r30649, r30673, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30677, r30626, r30676, MPFR_RNDN);
        mpfr_sub(r30678, r30677, r30627, MPFR_RNDN);
        mpfr_div(r30679, r30675, r30678, MPFR_RNDN);
        mpfr_add(r30680, r30674, r30679, MPFR_RNDN);
        mpfr_mul(r30681, r30634, r30680, MPFR_RNDN);
        return mpfr_get_d(r30681, MPFR_RNDN);
}

static mpfr_t r30682, r30683, r30684, r30685, r30686, r30687, r30688, r30689, r30690, r30691, r30692, r30693, r30694, r30695, r30696, r30697, r30698, r30699, r30700, r30701, r30702, r30703, r30704, r30705, r30706, r30707, r30708, r30709, r30710, r30711, r30712, r30713, r30714, r30715, r30716, r30717, r30718, r30719, r30720, r30721, r30722, r30723, r30724, r30725, r30726, r30727, r30728, r30729, r30730, r30731, r30732, r30733, r30734, r30735, r30736, r30737, r30738, r30739, r30740, r30741, r30742, r30743;

void setup_mpfr_f_dm() {
        mpfr_set_default_prec(1424);
        mpfr_init(r30682);
        mpfr_init(r30683);
        mpfr_init(r30684);
        mpfr_init_set_str(r30685, "0.5", 10, MPFR_RNDN);
        mpfr_init_set_str(r30686, "7", 10, MPFR_RNDN);
        mpfr_init(r30687);
        mpfr_init(r30688);
        mpfr_init_set_str(r30689, "1", 10, MPFR_RNDN);
        mpfr_init(r30690);
        mpfr_init(r30691);
        mpfr_init(r30692);
        mpfr_init(r30693);
        mpfr_init(r30694);
        mpfr_init(r30695);
        mpfr_init(r30696);
        mpfr_init_set_str(r30697, "-0.13857109526572012", 10, MPFR_RNDN);
        mpfr_init_set_str(r30698, "6", 10, MPFR_RNDN);
        mpfr_init(r30699);
        mpfr_init(r30700);
        mpfr_init(r30701);
        mpfr_init_set_str(r30702, "9.984369578019572e-06", 10, MPFR_RNDN);
        mpfr_init(r30703);
        mpfr_init(r30704);
        mpfr_init(r30705);
        mpfr_init_set_str(r30706, "12.507343278686905", 10, MPFR_RNDN);
        mpfr_init_set_str(r30707, "5", 10, MPFR_RNDN);
        mpfr_init(r30708);
        mpfr_init(r30709);
        mpfr_init(r30710);
        mpfr_init(r30711);
        mpfr_init_set_str(r30712, "-176.6150291621406", 10, MPFR_RNDN);
        mpfr_init_set_str(r30713, "4", 10, MPFR_RNDN);
        mpfr_init(r30714);
        mpfr_init(r30715);
        mpfr_init(r30716);
        mpfr_init_set_str(r30717, "676.5203681218851", 10, MPFR_RNDN);
        mpfr_init_set_str(r30718, "0", 10, MPFR_RNDN);
        mpfr_init(r30719);
        mpfr_init(r30720);
        mpfr_init_set_str(r30721, "0.9999999999998099", 10, MPFR_RNDN);
        mpfr_init(r30722);
        mpfr_init(r30723);
        mpfr_init_set_str(r30724, "-1259.1392167224028", 10, MPFR_RNDN);
        mpfr_init_set_str(r30725, "2", 10, MPFR_RNDN);
        mpfr_init(r30726);
        mpfr_init(r30727);
        mpfr_init(r30728);
        mpfr_init_set_str(r30729, "771.3234287776531", 10, MPFR_RNDN);
        mpfr_init_set_str(r30730, "3", 10, MPFR_RNDN);
        mpfr_init(r30731);
        mpfr_init(r30732);
        mpfr_init(r30733);
        mpfr_init(r30734);
        mpfr_init(r30735);
        mpfr_init(r30736);
        mpfr_init_set_str(r30737, "1.5056327351493116e-07", 10, MPFR_RNDN);
        mpfr_init_set_str(r30738, "8", 10, MPFR_RNDN);
        mpfr_init(r30739);
        mpfr_init(r30740);
        mpfr_init(r30741);
        mpfr_init(r30742);
        mpfr_init(r30743);
}

double f_dm(double z) {
        mpfr_const_pi(r30682, MPFR_RNDN);
        mpfr_add(r30683, r30682, r30682, MPFR_RNDN);
        mpfr_sqrt(r30684, r30683, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30687, r30685, r30686, MPFR_RNDN);
        mpfr_set_d(r30688, z, MPFR_RNDN);
        ;
        mpfr_sub(r30690, r30688, r30689, MPFR_RNDN);
        mpfr_add(r30691, r30687, r30690, MPFR_RNDN);
        mpfr_add(r30692, r30685, r30690, MPFR_RNDN);
        mpfr_pow(r30693, r30691, r30692, MPFR_RNDN);
        mpfr_exp(r30694, r30691, MPFR_RNDN);
        mpfr_div(r30695, r30693, r30694, MPFR_RNDN);
        mpfr_mul(r30696, r30684, r30695, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30699, r30689, r30698, MPFR_RNDN);
        mpfr_sub(r30700, r30688, r30699, MPFR_RNDN);
        mpfr_div(r30701, r30697, r30700, MPFR_RNDN);
        ;
        mpfr_add(r30703, r30690, r30686, MPFR_RNDN);
        mpfr_div(r30704, r30702, r30703, MPFR_RNDN);
        mpfr_add(r30705, r30701, r30704, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30708, r30688, r30707, MPFR_RNDN);
        mpfr_sub(r30709, r30708, r30689, MPFR_RNDN);
        mpfr_div(r30710, r30706, r30709, MPFR_RNDN);
        mpfr_add(r30711, r30705, r30710, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30714, r30688, r30713, MPFR_RNDN);
        mpfr_sub(r30715, r30714, r30689, MPFR_RNDN);
        mpfr_div(r30716, r30712, r30715, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30719, r30688, r30718, MPFR_RNDN);
        mpfr_div(r30720, r30717, r30719, MPFR_RNDN);
        ;
        mpfr_add(r30722, r30720, r30721, MPFR_RNDN);
        mpfr_add(r30723, r30716, r30722, MPFR_RNDN);
        ;
        ;
        mpfr_sub(r30726, r30689, r30725, MPFR_RNDN);
        mpfr_sub(r30727, r30688, r30726, MPFR_RNDN);
        mpfr_div(r30728, r30724, r30727, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30731, r30730, r30688, MPFR_RNDN);
        mpfr_sub(r30732, r30731, r30689, MPFR_RNDN);
        mpfr_div(r30733, r30729, r30732, MPFR_RNDN);
        mpfr_add(r30734, r30728, r30733, MPFR_RNDN);
        mpfr_add(r30735, r30723, r30734, MPFR_RNDN);
        mpfr_add(r30736, r30711, r30735, MPFR_RNDN);
        ;
        ;
        mpfr_add(r30739, r30688, r30738, MPFR_RNDN);
        mpfr_sub(r30740, r30739, r30689, MPFR_RNDN);
        mpfr_div(r30741, r30737, r30740, MPFR_RNDN);
        mpfr_add(r30742, r30736, r30741, MPFR_RNDN);
        mpfr_mul(r30743, r30696, r30742, MPFR_RNDN);
        return mpfr_get_d(r30743, MPFR_RNDN);
}

