#include #include #include #include #include char *name = "From Warwick Tucker's Validated Numerics"; double f_if() { float r30408 = 333.75f; float r30409 = 33096.0f; float r30410 = 6.0f; float r30411 = pow(r30409, r30410); float r30412 = r30408 * r30411; float r30413 = 77617.0f; float r30414 = r30413 * r30413; float r30415 = 11.0f; float r30416 = r30415 * r30414; float r30417 = r30409 * r30409; float r30418 = r30416 * r30417; float r30419 = -r30411; float r30420 = r30418 + r30419; float r30421 = -121.0f; float r30422 = 4.0f; float r30423 = pow(r30409, r30422); float r30424 = r30421 * r30423; float r30425 = r30420 + r30424; float r30426 = -2.0f; float r30427 = r30425 + r30426; float r30428 = r30414 * r30427; float r30429 = r30412 + r30428; float r30430 = 5.5f; float r30431 = 8.0f; float r30432 = pow(r30409, r30431); float r30433 = r30430 * r30432; float r30434 = r30429 + r30433; float r30435 = 2.0f; float r30436 = r30435 * r30409; float r30437 = r30413 / r30436; float r30438 = r30434 + r30437; return r30438; } double f_id() { double r30439 = 333.75; double r30440 = 33096.0; double r30441 = 6.0; double r30442 = pow(r30440, r30441); double r30443 = r30439 * r30442; double r30444 = 77617.0; double r30445 = r30444 * r30444; double r30446 = 11.0; double r30447 = r30446 * r30445; double r30448 = r30440 * r30440; double r30449 = r30447 * r30448; double r30450 = -r30442; double r30451 = r30449 + r30450; double r30452 = -121.0; double r30453 = 4.0; double r30454 = pow(r30440, r30453); double r30455 = r30452 * r30454; double r30456 = r30451 + r30455; double r30457 = -2.0; double r30458 = r30456 + r30457; double r30459 = r30445 * r30458; double r30460 = r30443 + r30459; double r30461 = 5.5; double r30462 = 8.0; double r30463 = pow(r30440, r30462); double r30464 = r30461 * r30463; double r30465 = r30460 + r30464; double r30466 = 2.0; double r30467 = r30466 * r30440; double r30468 = r30444 / r30467; double r30469 = r30465 + r30468; return r30469; } double f_of() { float r30470 = 333.75f; float r30471 = 33096.0f; float r30472 = 6.0f; float r30473 = pow(r30471, r30472); float r30474 = r30470 * r30473; float r30475 = 77617.0f; float r30476 = r30475 * r30475; float r30477 = 11.0f; float r30478 = r30477 * r30476; float r30479 = r30471 * r30471; float r30480 = r30478 * r30479; float r30481 = -r30473; float r30482 = r30480 + r30481; float r30483 = -121.0f; float r30484 = 4.0f; float r30485 = pow(r30471, r30484); float r30486 = r30483 * r30485; float r30487 = r30482 + r30486; float r30488 = -2.0f; float r30489 = r30487 + r30488; float r30490 = r30476 * r30489; float r30491 = r30474 + r30490; float r30492 = r30491 * r30491; float r30493 = 5.5f; float r30494 = 8.0f; float r30495 = pow(r30471, r30494); float r30496 = r30493 * r30495; float r30497 = r30496 * r30496; float r30498 = r30492 - r30497; float r30499 = r30491 - r30496; float r30500 = r30498 / r30499; float r30501 = r30500 * r30500; float r30502 = 2.0f; float r30503 = r30502 * r30471; float r30504 = r30475 / r30503; float r30505 = r30504 * r30504; float r30506 = r30501 - r30505; float r30507 = 3.0f; float r30508 = pow(r30488, r30507); float r30509 = r30485 * r30483; float r30510 = pow(r30509, r30507); float r30511 = r30508 + r30510; float r30512 = r30479 * r30477; float r30513 = r30512 * r30476; float r30514 = r30513 - r30481; float r30515 = r30511 * r30514; float r30516 = r30488 * r30488; float r30517 = r30509 * r30509; float r30518 = r30488 * r30509; float r30519 = r30517 - r30518; float r30520 = r30516 + r30519; float r30521 = r30513 * r30513; float r30522 = r30481 * r30481; float r30523 = r30521 - r30522; float r30524 = r30520 * r30523; float r30525 = r30515 + r30524; float r30526 = r30525 * r30476; float r30527 = r30474 * r30474; float r30528 = r30495 * r30493; float r30529 = r30528 * r30528; float r30530 = r30474 * r30528; float r30531 = r30529 - r30530; float r30532 = r30527 + r30531; float r30533 = r30526 * r30532; float r30534 = r30520 * r30514; float r30535 = pow(r30474, r30507); float r30536 = pow(r30528, r30507); float r30537 = r30535 + r30536; float r30538 = r30534 * r30537; float r30539 = r30533 + r30538; float r30540 = r30534 * r30532; float r30541 = r30539 / r30540; float r30542 = 1.0f; float r30543 = r30541 / r30542; float r30544 = r30471 + r30471; float r30545 = r30475 / r30544; float r30546 = r30543 - r30545; float r30547 = r30506 / r30546; return r30547; } double f_od() { double r30548 = 333.75; double r30549 = 33096.0; double r30550 = 6.0; double r30551 = pow(r30549, r30550); double r30552 = r30548 * r30551; double r30553 = 77617.0; double r30554 = r30553 * r30553; double r30555 = 11.0; double r30556 = r30555 * r30554; double r30557 = r30549 * r30549; double r30558 = r30556 * r30557; double r30559 = -r30551; double r30560 = r30558 + r30559; double r30561 = -121.0; double r30562 = 4.0; double r30563 = pow(r30549, r30562); double r30564 = r30561 * r30563; double r30565 = r30560 + r30564; double r30566 = -2.0; double r30567 = r30565 + r30566; double r30568 = r30554 * r30567; double r30569 = r30552 + r30568; double r30570 = r30569 * r30569; double r30571 = 5.5; double r30572 = 8.0; double r30573 = pow(r30549, r30572); double r30574 = r30571 * r30573; double r30575 = r30574 * r30574; double r30576 = r30570 - r30575; double r30577 = r30569 - r30574; double r30578 = r30576 / r30577; double r30579 = r30578 * r30578; double r30580 = 2.0; double r30581 = r30580 * r30549; double r30582 = r30553 / r30581; double r30583 = r30582 * r30582; double r30584 = r30579 - r30583; double r30585 = 3.0; double r30586 = pow(r30566, r30585); double r30587 = r30563 * r30561; double r30588 = pow(r30587, r30585); double r30589 = r30586 + r30588; double r30590 = r30557 * r30555; double r30591 = r30590 * r30554; double r30592 = r30591 - r30559; double r30593 = r30589 * r30592; double r30594 = r30566 * r30566; double r30595 = r30587 * r30587; double r30596 = r30566 * r30587; double r30597 = r30595 - r30596; double r30598 = r30594 + r30597; double r30599 = r30591 * r30591; double r30600 = r30559 * r30559; double r30601 = r30599 - r30600; double r30602 = r30598 * r30601; double r30603 = r30593 + r30602; double r30604 = r30603 * r30554; double r30605 = r30552 * r30552; double r30606 = r30573 * r30571; double r30607 = r30606 * r30606; double r30608 = r30552 * r30606; double r30609 = r30607 - r30608; double r30610 = r30605 + r30609; double r30611 = r30604 * r30610; double r30612 = r30598 * r30592; double r30613 = pow(r30552, r30585); double r30614 = pow(r30606, r30585); double r30615 = r30613 + r30614; double r30616 = r30612 * r30615; double r30617 = r30611 + r30616; double r30618 = r30612 * r30610; double r30619 = r30617 / r30618; double r30620 = 1.0; double r30621 = r30619 / r30620; double r30622 = r30549 + r30549; double r30623 = r30553 / r30622; double r30624 = r30621 - r30623; double r30625 = r30584 / r30624; return r30625; } 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 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; void setup_mpfr_f_im() { mpfr_set_default_prec(144); mpfr_init_set_str(r30626, "333.75", 10, MPFR_RNDN); mpfr_init_set_str(r30627, "33096", 10, MPFR_RNDN); mpfr_init_set_str(r30628, "6", 10, MPFR_RNDN); mpfr_init(r30629); mpfr_init(r30630); mpfr_init_set_str(r30631, "77617", 10, MPFR_RNDN); mpfr_init(r30632); mpfr_init_set_str(r30633, "11", 10, MPFR_RNDN); mpfr_init(r30634); mpfr_init(r30635); mpfr_init(r30636); mpfr_init(r30637); mpfr_init(r30638); mpfr_init_set_str(r30639, "-121", 10, MPFR_RNDN); mpfr_init_set_str(r30640, "4", 10, MPFR_RNDN); mpfr_init(r30641); mpfr_init(r30642); mpfr_init(r30643); mpfr_init_set_str(r30644, "-2", 10, MPFR_RNDN); mpfr_init(r30645); mpfr_init(r30646); mpfr_init(r30647); mpfr_init_set_str(r30648, "5.5", 10, MPFR_RNDN); mpfr_init_set_str(r30649, "8", 10, MPFR_RNDN); mpfr_init(r30650); mpfr_init(r30651); mpfr_init(r30652); mpfr_init_set_str(r30653, "2", 10, MPFR_RNDN); mpfr_init(r30654); mpfr_init(r30655); mpfr_init(r30656); } double f_im() { ; ; ; mpfr_pow(r30629, r30627, r30628, MPFR_RNDN); mpfr_mul(r30630, r30626, r30629, MPFR_RNDN); ; mpfr_sqr(r30632, r30631, MPFR_RNDN); ; mpfr_mul(r30634, r30633, r30632, MPFR_RNDN); mpfr_sqr(r30635, r30627, MPFR_RNDN); mpfr_mul(r30636, r30634, r30635, MPFR_RNDN); mpfr_neg(r30637, r30629, MPFR_RNDN); mpfr_add(r30638, r30636, r30637, MPFR_RNDN); ; ; mpfr_pow(r30641, r30627, r30640, MPFR_RNDN); mpfr_mul(r30642, r30639, r30641, MPFR_RNDN); mpfr_add(r30643, r30638, r30642, MPFR_RNDN); ; mpfr_add(r30645, r30643, r30644, MPFR_RNDN); mpfr_mul(r30646, r30632, r30645, MPFR_RNDN); mpfr_add(r30647, r30630, r30646, MPFR_RNDN); ; ; mpfr_pow(r30650, r30627, r30649, MPFR_RNDN); mpfr_mul(r30651, r30648, r30650, MPFR_RNDN); mpfr_add(r30652, r30647, r30651, MPFR_RNDN); ; mpfr_mul(r30654, r30653, r30627, MPFR_RNDN); mpfr_div(r30655, r30631, r30654, MPFR_RNDN); mpfr_add(r30656, r30652, r30655, MPFR_RNDN); return mpfr_get_d(r30656, MPFR_RNDN); } static mpfr_t 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, 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; void setup_mpfr_f_fm() { mpfr_set_default_prec(144); mpfr_init_set_str(r30657, "333.75", 10, MPFR_RNDN); mpfr_init_set_str(r30658, "33096", 10, MPFR_RNDN); mpfr_init_set_str(r30659, "6", 10, MPFR_RNDN); mpfr_init(r30660); mpfr_init(r30661); mpfr_init_set_str(r30662, "77617", 10, MPFR_RNDN); mpfr_init(r30663); mpfr_init_set_str(r30664, "11", 10, MPFR_RNDN); mpfr_init(r30665); mpfr_init(r30666); mpfr_init(r30667); mpfr_init(r30668); mpfr_init(r30669); mpfr_init_set_str(r30670, "-121", 10, MPFR_RNDN); mpfr_init_set_str(r30671, "4", 10, MPFR_RNDN); mpfr_init(r30672); mpfr_init(r30673); mpfr_init(r30674); mpfr_init_set_str(r30675, "-2", 10, MPFR_RNDN); mpfr_init(r30676); mpfr_init(r30677); mpfr_init(r30678); mpfr_init(r30679); mpfr_init_set_str(r30680, "5.5", 10, MPFR_RNDN); mpfr_init_set_str(r30681, "8", 10, MPFR_RNDN); mpfr_init(r30682); mpfr_init(r30683); mpfr_init(r30684); mpfr_init(r30685); mpfr_init(r30686); mpfr_init(r30687); mpfr_init(r30688); mpfr_init_set_str(r30689, "2", 10, MPFR_RNDN); mpfr_init(r30690); mpfr_init(r30691); mpfr_init(r30692); mpfr_init(r30693); mpfr_init_set_str(r30694, "3", 10, MPFR_RNDN); mpfr_init(r30695); mpfr_init(r30696); mpfr_init(r30697); mpfr_init(r30698); mpfr_init(r30699); mpfr_init(r30700); mpfr_init(r30701); mpfr_init(r30702); mpfr_init(r30703); mpfr_init(r30704); mpfr_init(r30705); mpfr_init(r30706); mpfr_init(r30707); mpfr_init(r30708); mpfr_init(r30709); mpfr_init(r30710); mpfr_init(r30711); mpfr_init(r30712); mpfr_init(r30713); mpfr_init(r30714); mpfr_init(r30715); mpfr_init(r30716); mpfr_init(r30717); mpfr_init(r30718); mpfr_init(r30719); mpfr_init(r30720); mpfr_init(r30721); mpfr_init(r30722); mpfr_init(r30723); mpfr_init(r30724); mpfr_init(r30725); mpfr_init(r30726); mpfr_init(r30727); mpfr_init(r30728); mpfr_init_set_str(r30729, "1", 10, MPFR_RNDN); mpfr_init(r30730); mpfr_init(r30731); mpfr_init(r30732); mpfr_init(r30733); mpfr_init(r30734); } double f_fm() { ; ; ; mpfr_pow(r30660, r30658, r30659, MPFR_RNDN); mpfr_mul(r30661, r30657, r30660, MPFR_RNDN); ; mpfr_sqr(r30663, r30662, MPFR_RNDN); ; mpfr_mul(r30665, r30664, r30663, MPFR_RNDN); mpfr_sqr(r30666, r30658, MPFR_RNDN); mpfr_mul(r30667, r30665, r30666, MPFR_RNDN); mpfr_neg(r30668, r30660, MPFR_RNDN); mpfr_add(r30669, r30667, r30668, MPFR_RNDN); ; ; mpfr_pow(r30672, r30658, r30671, MPFR_RNDN); mpfr_mul(r30673, r30670, r30672, MPFR_RNDN); mpfr_add(r30674, r30669, r30673, MPFR_RNDN); ; mpfr_add(r30676, r30674, r30675, MPFR_RNDN); mpfr_mul(r30677, r30663, r30676, MPFR_RNDN); mpfr_add(r30678, r30661, r30677, MPFR_RNDN); mpfr_sqr(r30679, r30678, MPFR_RNDN); ; ; mpfr_pow(r30682, r30658, r30681, MPFR_RNDN); mpfr_mul(r30683, r30680, r30682, MPFR_RNDN); mpfr_sqr(r30684, r30683, MPFR_RNDN); mpfr_sub(r30685, r30679, r30684, MPFR_RNDN); mpfr_sub(r30686, r30678, r30683, MPFR_RNDN); mpfr_div(r30687, r30685, r30686, MPFR_RNDN); mpfr_sqr(r30688, r30687, MPFR_RNDN); ; mpfr_mul(r30690, r30689, r30658, MPFR_RNDN); mpfr_div(r30691, r30662, r30690, MPFR_RNDN); mpfr_sqr(r30692, r30691, MPFR_RNDN); mpfr_sub(r30693, r30688, r30692, MPFR_RNDN); ; mpfr_pow(r30695, r30675, r30694, MPFR_RNDN); mpfr_mul(r30696, r30672, r30670, MPFR_RNDN); mpfr_pow(r30697, r30696, r30694, MPFR_RNDN); mpfr_add(r30698, r30695, r30697, MPFR_RNDN); mpfr_mul(r30699, r30666, r30664, MPFR_RNDN); mpfr_mul(r30700, r30699, r30663, MPFR_RNDN); mpfr_sub(r30701, r30700, r30668, MPFR_RNDN); mpfr_mul(r30702, r30698, r30701, MPFR_RNDN); mpfr_sqr(r30703, r30675, MPFR_RNDN); mpfr_sqr(r30704, r30696, MPFR_RNDN); mpfr_mul(r30705, r30675, r30696, MPFR_RNDN); mpfr_sub(r30706, r30704, r30705, MPFR_RNDN); mpfr_add(r30707, r30703, r30706, MPFR_RNDN); mpfr_sqr(r30708, r30700, MPFR_RNDN); mpfr_sqr(r30709, r30668, MPFR_RNDN); mpfr_sub(r30710, r30708, r30709, MPFR_RNDN); mpfr_mul(r30711, r30707, r30710, MPFR_RNDN); mpfr_add(r30712, r30702, r30711, MPFR_RNDN); mpfr_mul(r30713, r30712, r30663, MPFR_RNDN); mpfr_sqr(r30714, r30661, MPFR_RNDN); mpfr_mul(r30715, r30682, r30680, MPFR_RNDN); mpfr_sqr(r30716, r30715, MPFR_RNDN); mpfr_mul(r30717, r30661, r30715, MPFR_RNDN); mpfr_sub(r30718, r30716, r30717, MPFR_RNDN); mpfr_add(r30719, r30714, r30718, MPFR_RNDN); mpfr_mul(r30720, r30713, r30719, MPFR_RNDN); mpfr_mul(r30721, r30707, r30701, MPFR_RNDN); mpfr_pow(r30722, r30661, r30694, MPFR_RNDN); mpfr_pow(r30723, r30715, r30694, MPFR_RNDN); mpfr_add(r30724, r30722, r30723, MPFR_RNDN); mpfr_mul(r30725, r30721, r30724, MPFR_RNDN); mpfr_add(r30726, r30720, r30725, MPFR_RNDN); mpfr_mul(r30727, r30721, r30719, MPFR_RNDN); mpfr_div(r30728, r30726, r30727, MPFR_RNDN); ; mpfr_div(r30730, r30728, r30729, MPFR_RNDN); mpfr_add(r30731, r30658, r30658, MPFR_RNDN); mpfr_div(r30732, r30662, r30731, MPFR_RNDN); mpfr_sub(r30733, r30730, r30732, MPFR_RNDN); mpfr_div(r30734, r30693, r30733, MPFR_RNDN); return mpfr_get_d(r30734, MPFR_RNDN); } static mpfr_t r30735, r30736, r30737, r30738, r30739, r30740, r30741, r30742, r30743, r30744, r30745, r30746, r30747, r30748, r30749, r30750, r30751, r30752, r30753, r30754, r30755, r30756, r30757, r30758, r30759, r30760, r30761, r30762, r30763, r30764, r30765, r30766, r30767, r30768, r30769, r30770, r30771, r30772, r30773, r30774, r30775, r30776, r30777, r30778, r30779, r30780, r30781, r30782, r30783, r30784, r30785, r30786, r30787, r30788, r30789, r30790, r30791, r30792, r30793, r30794, r30795, r30796, r30797, r30798, r30799, r30800, r30801, r30802, r30803, r30804, r30805, r30806, r30807, r30808, r30809, r30810, r30811, r30812; void setup_mpfr_f_dm() { mpfr_set_default_prec(144); mpfr_init_set_str(r30735, "333.75", 10, MPFR_RNDN); mpfr_init_set_str(r30736, "33096", 10, MPFR_RNDN); mpfr_init_set_str(r30737, "6", 10, MPFR_RNDN); mpfr_init(r30738); mpfr_init(r30739); mpfr_init_set_str(r30740, "77617", 10, MPFR_RNDN); mpfr_init(r30741); mpfr_init_set_str(r30742, "11", 10, MPFR_RNDN); mpfr_init(r30743); mpfr_init(r30744); mpfr_init(r30745); mpfr_init(r30746); mpfr_init(r30747); mpfr_init_set_str(r30748, "-121", 10, MPFR_RNDN); mpfr_init_set_str(r30749, "4", 10, MPFR_RNDN); mpfr_init(r30750); mpfr_init(r30751); mpfr_init(r30752); mpfr_init_set_str(r30753, "-2", 10, MPFR_RNDN); mpfr_init(r30754); mpfr_init(r30755); mpfr_init(r30756); mpfr_init(r30757); mpfr_init_set_str(r30758, "5.5", 10, MPFR_RNDN); mpfr_init_set_str(r30759, "8", 10, MPFR_RNDN); mpfr_init(r30760); mpfr_init(r30761); mpfr_init(r30762); mpfr_init(r30763); mpfr_init(r30764); mpfr_init(r30765); mpfr_init(r30766); mpfr_init_set_str(r30767, "2", 10, MPFR_RNDN); mpfr_init(r30768); mpfr_init(r30769); mpfr_init(r30770); mpfr_init(r30771); mpfr_init_set_str(r30772, "3", 10, MPFR_RNDN); mpfr_init(r30773); mpfr_init(r30774); mpfr_init(r30775); mpfr_init(r30776); mpfr_init(r30777); mpfr_init(r30778); mpfr_init(r30779); mpfr_init(r30780); mpfr_init(r30781); mpfr_init(r30782); mpfr_init(r30783); mpfr_init(r30784); mpfr_init(r30785); mpfr_init(r30786); mpfr_init(r30787); mpfr_init(r30788); mpfr_init(r30789); mpfr_init(r30790); mpfr_init(r30791); mpfr_init(r30792); mpfr_init(r30793); mpfr_init(r30794); mpfr_init(r30795); mpfr_init(r30796); mpfr_init(r30797); mpfr_init(r30798); mpfr_init(r30799); mpfr_init(r30800); mpfr_init(r30801); mpfr_init(r30802); mpfr_init(r30803); mpfr_init(r30804); mpfr_init(r30805); mpfr_init(r30806); mpfr_init_set_str(r30807, "1", 10, MPFR_RNDN); mpfr_init(r30808); mpfr_init(r30809); mpfr_init(r30810); mpfr_init(r30811); mpfr_init(r30812); } double f_dm() { ; ; ; mpfr_pow(r30738, r30736, r30737, MPFR_RNDN); mpfr_mul(r30739, r30735, r30738, MPFR_RNDN); ; mpfr_sqr(r30741, r30740, MPFR_RNDN); ; mpfr_mul(r30743, r30742, r30741, MPFR_RNDN); mpfr_sqr(r30744, r30736, MPFR_RNDN); mpfr_mul(r30745, r30743, r30744, MPFR_RNDN); mpfr_neg(r30746, r30738, MPFR_RNDN); mpfr_add(r30747, r30745, r30746, MPFR_RNDN); ; ; mpfr_pow(r30750, r30736, r30749, MPFR_RNDN); mpfr_mul(r30751, r30748, r30750, MPFR_RNDN); mpfr_add(r30752, r30747, r30751, MPFR_RNDN); ; mpfr_add(r30754, r30752, r30753, MPFR_RNDN); mpfr_mul(r30755, r30741, r30754, MPFR_RNDN); mpfr_add(r30756, r30739, r30755, MPFR_RNDN); mpfr_sqr(r30757, r30756, MPFR_RNDN); ; ; mpfr_pow(r30760, r30736, r30759, MPFR_RNDN); mpfr_mul(r30761, r30758, r30760, MPFR_RNDN); mpfr_sqr(r30762, r30761, MPFR_RNDN); mpfr_sub(r30763, r30757, r30762, MPFR_RNDN); mpfr_sub(r30764, r30756, r30761, MPFR_RNDN); mpfr_div(r30765, r30763, r30764, MPFR_RNDN); mpfr_sqr(r30766, r30765, MPFR_RNDN); ; mpfr_mul(r30768, r30767, r30736, MPFR_RNDN); mpfr_div(r30769, r30740, r30768, MPFR_RNDN); mpfr_sqr(r30770, r30769, MPFR_RNDN); mpfr_sub(r30771, r30766, r30770, MPFR_RNDN); ; mpfr_pow(r30773, r30753, r30772, MPFR_RNDN); mpfr_mul(r30774, r30750, r30748, MPFR_RNDN); mpfr_pow(r30775, r30774, r30772, MPFR_RNDN); mpfr_add(r30776, r30773, r30775, MPFR_RNDN); mpfr_mul(r30777, r30744, r30742, MPFR_RNDN); mpfr_mul(r30778, r30777, r30741, MPFR_RNDN); mpfr_sub(r30779, r30778, r30746, MPFR_RNDN); mpfr_mul(r30780, r30776, r30779, MPFR_RNDN); mpfr_sqr(r30781, r30753, MPFR_RNDN); mpfr_sqr(r30782, r30774, MPFR_RNDN); mpfr_mul(r30783, r30753, r30774, MPFR_RNDN); mpfr_sub(r30784, r30782, r30783, MPFR_RNDN); mpfr_add(r30785, r30781, r30784, MPFR_RNDN); mpfr_sqr(r30786, r30778, MPFR_RNDN); mpfr_sqr(r30787, r30746, MPFR_RNDN); mpfr_sub(r30788, r30786, r30787, MPFR_RNDN); mpfr_mul(r30789, r30785, r30788, MPFR_RNDN); mpfr_add(r30790, r30780, r30789, MPFR_RNDN); mpfr_mul(r30791, r30790, r30741, MPFR_RNDN); mpfr_sqr(r30792, r30739, MPFR_RNDN); mpfr_mul(r30793, r30760, r30758, MPFR_RNDN); mpfr_sqr(r30794, r30793, MPFR_RNDN); mpfr_mul(r30795, r30739, r30793, MPFR_RNDN); mpfr_sub(r30796, r30794, r30795, MPFR_RNDN); mpfr_add(r30797, r30792, r30796, MPFR_RNDN); mpfr_mul(r30798, r30791, r30797, MPFR_RNDN); mpfr_mul(r30799, r30785, r30779, MPFR_RNDN); mpfr_pow(r30800, r30739, r30772, MPFR_RNDN); mpfr_pow(r30801, r30793, r30772, MPFR_RNDN); mpfr_add(r30802, r30800, r30801, MPFR_RNDN); mpfr_mul(r30803, r30799, r30802, MPFR_RNDN); mpfr_add(r30804, r30798, r30803, MPFR_RNDN); mpfr_mul(r30805, r30799, r30797, MPFR_RNDN); mpfr_div(r30806, r30804, r30805, MPFR_RNDN); ; mpfr_div(r30808, r30806, r30807, MPFR_RNDN); mpfr_add(r30809, r30736, r30736, MPFR_RNDN); mpfr_div(r30810, r30740, r30809, MPFR_RNDN); mpfr_sub(r30811, r30808, r30810, MPFR_RNDN); mpfr_div(r30812, r30771, r30811, MPFR_RNDN); return mpfr_get_d(r30812, MPFR_RNDN); }