2438 double lmu = log(muh);
2439 double Lsh = log(sh);
2440 double sh2 = pow(sh, 2);
2441 double sh3 = pow(sh, 3);
2442 double lz2 = pow(lz, 2);
2443 double lz3 = pow(lz, 3);
2444 double sqrtz = pow(z, 0.5);
2445 double z2 = pow(z, 2);
2446 double z3 = pow(z, 3);
2452 res = -24.2913 + 0.296296 * Lsh + 0.3950617283950617 * lmu * Lsh +
2453 lmu * (1.0534979423868314 - 7.111111111111111 * lz) -
2454 11.55555 * lz - 26.8464 * sh - 4.444445 * lz * sh - 86.7723 * z +
2455 3.55556 * Lsh * z - 28.77965 * lz * z - 428.313 * sh * z +
2456 16. * Lsh * sh * z - 25.4303 * lz * sh * z + 3.555555 * Lsh * lz * sh * z +
2457 6.320987654320987 * pow(lmu, 2) + 1.7777775 * z * lz2 -
2458 18.808425 * sh * z * lz2 + 0.5925925 * z * lz3 +
2459 0.5925925 * sh * z * lz3 - 0.29629625 * sh * z * pow(lz, 4) -
2460 222.769 * sh2 - 0.0987654 * Lsh * sh2 -
2461 66.186 * lz * sh2 - 1276.44 * z * sh2 +
2462 32. * Lsh * z * sh2 + 65.7645 * lz * z * sh2 +
2463 7.1111 * Lsh * lz * z * sh2 - 8. * lz2 * sh2 -
2464 85.7585 * z * lz2 * sh2 -
2465 0.5925925 * lz3 * sh2 +
2466 1.777775 * z * lz3 * sh2 -
2467 0.8888875 * z * pow(lz, 4) * sh2 - 478.485 * sh3 -
2468 0.131687 * Lsh * sh3 - 171.951 * lz * sh3 -
2469 2553.47 * z * sh3 + 49.7778 * Lsh * z * sh3 +
2470 206.4045 * lz * z * sh3 + 10.66665 * Lsh * lz * z * sh3 -
2471 16.8889 * lz2 * sh3 -
2472 194.03575 * z * lz2 * sh3 -
2473 1.777775 * lz3 * sh3 +
2474 4.14815 * z * lz3 * sh3 -
2475 1.777775 * z * pow(lz, 4) * sh3 +
2476 lmu * sh3*(-0.00125416421712718 +
2477 0.022574955908289243 * pow(z, -3)) -
2478 0.0142243 * sh3 * pow(z, -3) +
2479 lmu * sh2*(-0.008465608465608466 +
2480 0.1523809523809524 * pow(z, -2)) -
2481 0.0132191 * sh2 * pow(z, -2) +
2482 0.0288536 * sh3 * pow(z, -2) + 0.8462 * sh * pow(z, -1) +
2483 0.367901 * sh2 * pow(z, -1) +
2484 10.8601 * sh3 * pow(z, -1) +
2485 2.758375 * lz * sh3 * pow(z, -1) +
2486 0.2962975 * lz2 * sh3 * pow(z, -1) +
2487 lmu * sh * (-0.07901234567901234 +
2488 1.4222222222222223 * pow(z, -1)) +
2489 140.368 * sh * sqrtz + 842.206 * sh2 * sqrtz +
2490 sh3*(-46.7892 * pow(z, -0.5) + 1918.36 * sqrtz) -
2491 23.3946 * pow(z, 1.5) + 96.5187 * z2 -
2492 29.7586 * Lsh * z2 - 162.7315 * lz * z2 +
2493 3.555555 * Lsh * lz * z2 + 108.781 * sh * z2 -
2494 111.923 * Lsh * sh * z2 - 448.7875 * lz * sh * z2 +
2495 14.2222 * Lsh * lz * sh * z2 +
2496 23.1111 * lz2 * z2 +
2497 3.55555 * Lsh * lz2 * z2 +
2498 73.77775 * sh * lz2 * z2 +
2499 14.222225 * Lsh * sh * lz2 * z2 +
2500 2.9629625 * lz3 * z2 +
2501 8.2963 * sh * lz3 * z2 +
2502 268.098 * sh2 * z2 -
2503 244.716 * Lsh * sh2 * z2 -
2504 816.045 * lz * sh2 * z2 +
2505 42.66665 * Lsh * lz * sh2 * z2 +
2506 178.66675 * lz2 * sh2 * z2 +
2507 32. * Lsh * lz2 * sh2 * z2 +
2508 16. * lz3 * sh2 * z2 +
2509 527.368 * sh3 * z2 -
2510 421.619 * Lsh * sh3 * z2 -
2511 1252.835 * lz * sh3 * z2 +
2512 94.815 * Lsh * lz * sh3 * z2 +
2513 340.74 * lz2 * sh3 * z2 +
2514 56.889 * Lsh * lz2 * sh3 * z2 +
2515 26.074125 * lz3 * sh3 * z2 +
2516 88.3801 * z3 + 55.2172 * Lsh * z3 +
2517 85.7285 * lz * z3 + 3.160495 * Lsh * lz * z3 +
2518 437.34 * sh * z3 + 249.663 * Lsh * sh * z3 +
2519 204.405 * lz * sh * z3 + 7.1111 * Lsh * lz * sh * z3 -
2520 30.22225 * lz2 * z3 -
2521 4.74075 * Lsh * lz2 * z3 -
2522 95.4075 * sh * lz2 * z3 -
2523 21.333325 * Lsh * sh * lz2 * z3 +
2524 823.218 * sh2 * z3 +
2525 668.137 * Lsh * sh2 * z3 +
2526 298.311 * lz * sh2 * z3 -
2527 253.3325 * lz2 * sh2 * z3 -
2528 56.889 * Lsh * lz2 * sh2 * z3 +
2529 1675.61 * sh3 * z3 +
2530 1391.36 * Lsh * sh3 * z3 +
2531 440.5585 * lz * sh3 * z3 -
2532 35.55555 * Lsh * lz * sh3 * z3 -
2533 496.89 * lz2 * sh3 * z3 -
2534 118.5185 * Lsh * lz2 * sh3 * z3;
2536 double Lshb = log(1. - sh);
2537 double sqrtsh = sqrt(sh);
2538 double sqrt4sh = sqrt(4. - sh);
2539 double Li2sh = real(
Li2(sh));
2540 double Li3sh = real(
Li3(sh));
2541 double Li3shb = real(
Li3(1. - sh));
2542 double Li4sh = real(
Li4(sh));
2543 double ash = asin(sqrtsh / 2.);
2544 double cl2 =
Cl2(2. * ash);
2545 double cl3 =
Cl3(2. * ash);
2548 res += (8.547960200246003 + 6.320987654320987 * pow(lmu, 2) +
2549 0.19753086419753085 * Li2sh * (-8. + 17. * sh) *
2550 pow(-1 + sh, -1) + 0.19753086419753085 * cl2 * (2. + sh) *
2551 sqrt4sh * sqrtsh * pow(sh, -2) -
2552 lmu * (6.716049382716049 * Lsh -
2553 0.7901234567901234 * ash * (2. + sh) * sqrt4sh * sqrtsh *
2554 pow(sh, -2) + 0.3950617283950617 * (4. - 31. * sh) *
2555 pow(sh, -1)) - 0.8888888888888888 * pow(Lsh, 2) *
2556 pow(-1 + sh, -3) * pow(sh, 2) -
2557 Lsh * (-0.19753086419753085 * ash * (2. + sh) * sqrt4sh *
2558 sqrtsh * pow(sh, -2) -
2559 0.09876543209876543 * pow(-1 + sh, -2)*
2560 (-114. + 199. * sh - 67. * pow(sh, 2) +
2561 2. * Lshb * (8. - 25. * sh + 17. * pow(sh, 2)))) +
2562 0.19753086419753085 * ash * sqrt4sh * sqrtsh * pow(-1 + sh, -3) *
2563 pow(sh, -2)*(-4. + 3. * sh + 18. * pow(sh, 2) -
2564 16. * pow(sh, 3) + 5. * pow(sh, 4)) +
2565 0.5925925925925926 * pow(ash, 2) * pow(-1 + sh, -4) *
2566 pow(sh, -1)*(12. - 33. * sh + 18. * pow(sh, 2) -
2567 4. * pow(sh, 4) + pow(sh, 5)) -
2568 0.0054869684499314125 * pow(-1 + sh, -4) * pow(sh, -1)*
2569 (288. + 5. * sh * (79. + 60. * pow(3.141592653589793, 2)) -
2570 2. * (2257. + 618. * pow(3.141592653589793, 2)) *
2571 pow(sh, 2) + 6. * (1391. +
2572 306. * pow(3.141592653589793, 2)) * pow(sh, 3) -
2573 10. * (617. + 123. * pow(3.141592653589793, 2)) *
2574 pow(sh, 4) + (1655. +
2575 312. * pow(3.141592653589793, 2)) * pow(sh, 5))) *
2578 res += pow(z, 0)*(-4.256813794434612 +
2579 6.247387287113509 * (1. - sh) +
2580 4.358260258488125 * pow(1. - sh, 2) +
2581 3.2040567770853476 * pow(1. - sh, 3) +
2582 2.512816857053804 * pow(1. - sh, 4) +
2583 log(muh)*(12.81635480205537 +
2584 6.816868131135554 * (1. - sh) +
2585 3.34395932171576 * pow(1. - sh, 2) +
2586 2.241420882666335 * pow(1. - sh, 3) +
2587 1.6783963168171898 * pow(1. - sh, 4)) +
2588 6.320987654320987 * pow(log(muh), 2));
2600 res += (0.5925925925925926 * sh * pow(Lsh, 4) * pow(-1 + sh, -3) -
2601 3.5555555555555554 * Li2sh * (1. + 2. * sh) *
2602 pow(-1 + sh, -2) + 1.1851851851851851 * (-3. + sh) *
2603 pow(Lsh, 3) * pow(-1 + sh, -2) +
2604 42.666666666666664 * lmu * pow(sh, -1) -
2605 lz * (64. * Li2sh * pow(9 - 9 * sh, -1) -
2606 3.5555555555555554 * Lsh *
2607 (2. + 2. * Lshb * (-1. + sh) + sh) * pow(-1 + sh, -2) +
2608 3.5555555555555554 * pow(Lsh, 2) * pow(-1 + sh, -1) +
2609 1.1851851851851851 *
2610 (-18. + sh * (27. + pow(3.141592653589793, 2))) *
2611 pow(-1 + sh, -1) * pow(sh, -1)) -
2612 4.273980100123001 * pow(-1 + sh, -3)*
2613 (-3. + 3. * sh + 2. * pow(sh, 2)) +
2614 3.5555555555555554 * ash * sqrt4sh * sqrtsh * pow(-1 + sh, -3) *
2615 pow(sh, -1)*(2. - 11. * sh + 3. * pow(sh, 2)) -
2616 3.5555555555555554 * Li3sh * pow(-1 + sh, -3)*
2617 (-3. - 3. * sh + 4. * pow(sh, 2)) -
2618 pow(Lsh, 2)*(-7.111111111111111 * Li2sh * sh *
2620 0.5925925925925926 * pow(-1 + sh, -3)*
2621 (12. + 3. * Lshb * (-7. + 5. * sh) +
2622 2. * sh * (-6. + pow(3.141592653589793, 2)) +
2624 21.333333333333332 * pow(ash, 2) * pow(-1 + sh, -4)*
2625 (-1. + 8. * sh - 5. * pow(sh, 2) + pow(sh, 3)) -
2626 Lsh * (-34.19184080098401 * sh * pow(-1 + sh, -3) +
2627 28.444444444444443 * Li3sh * sh * pow(-1 + sh, -3) -
2628 3.5555555555555554 * Li2sh * pow(-1 + sh, -3)*
2629 (-5. + sh + 2. * pow(sh, 2)) +
2630 0.5925925925925926 * pow(-1 + sh, -3) * pow(sh, -1)*
2631 (36. - sh * (66. + 6. * Lshb +
2632 7. * pow(3.141592653589793, 2)) +
2633 (69. - 6. * Lshb + 5. * pow(3.141592653589793, 2)) *
2634 pow(sh, 2) + 3. * (-13. + 4. * Lshb) * pow(sh, 3))) +
2635 0.11851851851851852 * pow(-1 + sh, -4) * pow(sh, -1)*
2636 (-4. * (-1. + sh) * pow(3.141592653589793, 4) * pow(sh, 2) +
2637 15. * (-1. + sh)*(-56. + 156. * sh +
2638 3. * (-49. + 8. * Li4sh) * pow(sh, 2) + 47. * pow(sh, 3))\
2639 + 5. * pow(3.141592653589793, 2)*
2640 (-3. + 8. * sh - 10. * pow(sh, 2) + 2. * pow(sh, 4)))) *
2643 res += (67.82325420628487 + 57.04716984919418 * (1. - sh) +
2644 51.233916954097545 * pow(1. - sh, 2) +
2645 45.87999913302814 * pow(1. - sh, 3) +
2646 log(z)*(-16. - 17.77777777777778 * (1. - sh) -
2647 18.469135802469136 * pow(1. - sh, 2) -
2648 18.874074074074073 * pow(1. - sh, 3) -
2649 19.152592592592594 * pow(1. - sh, 4)) +
2650 41.468913753494085 * pow(1. - sh, 4) +
2651 log(muh)*(42.666666666666664 +
2652 42.666666666666664 * (1. - sh) +
2653 42.666666666666664 * pow(1. - sh, 2) +
2654 42.666666666666664 * pow(1. - sh, 3) +
2655 42.666666666666664 * pow(1. - sh, 4))) * pow(z, 1);
2660 res += -2.3703703703703702 * (-2. + sh) * pow(3.141592653589793, 2) *
2661 pow(sh, -1) * pow(z, 1.5);
2663 res += (23.39461783961922 + 46.78923567923844 * (1. - sh) +
2664 46.78923567923844 * pow(1. - sh, 2) +
2665 46.78923567923844 * pow(1. - sh, 3) +
2666 46.78923567923844 * pow(1. - sh, 4)) * pow(z, 1.5);
2671 res += (3.5555555555555554 * Li3shb * (6. + sh) * pow(-1 + sh, -3) +
2672 3.5555555555555554 * Li3sh * (-11. + 3. * sh) *
2673 pow(-1 + sh, -3) + 0.5925925925925926 * (2. + sh) *
2674 pow(Lsh, 3) * pow(-1 + sh, -3) -
2675 lmu * (-21.333333333333332 * pow(sh, -2) -
2676 42.666666666666664 * Lsh * pow(sh, -2)) -
2677 3.5555555555555554 * Li2sh * pow(-1 + sh, -3) * pow(sh, -2)*
2678 (-3. + 10. * sh - 9. * pow(sh, 2) + pow(sh, 3)) -
2679 3.5555555555555554 * cl2 * sqrt4sh * sqrtsh * pow(-1 + sh, -3) *
2680 pow(sh, -1)*(2. + 8. * pow(sh, 2) - 5. * pow(sh, 3) +
2681 pow(sh, 4)) + 3.5555555555555554 * ash * sqrt4sh * sqrtsh *
2682 pow(-1 + sh, -3) * pow(sh, -1)*
2683 (-2. - 4. * sh + 24. * pow(sh, 2) - 15. * pow(sh, 3) +
2684 3. * pow(sh, 4)) - 4.273980100123001 * pow(-1 + sh, -4) *
2685 pow(sh, -2)*(8. - 32. * sh + 57. * pow(sh, 2) -
2686 52. * pow(sh, 3) + 13. * pow(sh, 4)) -
2687 0.8888888888888888 * pow(Lsh, 2) * pow(-1 + sh, -4) *
2688 pow(sh, -2)*(-31. + 113. * sh +
2689 3. * (-53. + 2. * Lshb) * pow(sh, 2) +
2690 (87. - 8. * Lshb) * pow(sh, 3) +
2691 2. * (-8. + Lshb) * pow(sh, 4)) -
2692 Lsh * (3.5555555555555554 * Li2sh * (-7. + 2. * sh) *
2694 21.333333333333332 * pow(ash, 2) * pow(-1 + sh, -4)*
2695 (-1. - 3. * sh + pow(sh, 2)) +
2696 0.2962962962962963 * pow(-1 + sh, -3) * pow(sh, -2)*
2697 (333. - 1032. * sh +
2698 (1089. - 46. * pow(3.141592653589793, 2)) *
2700 (-153. + pow(3.141592653589793, 2)) * pow(sh, 3) +
2701 12. * Lshb * (-3. + 10. * sh - 9. * pow(sh, 2) +
2703 3.5555555555555554 * ash * sqrt4sh * sqrtsh *
2704 pow(-1 + sh, -3) * pow(sh, -1)*
2705 (2. + 8. * pow(sh, 2) - 5. * pow(sh, 3) + pow(sh, 4))) +
2706 10.666666666666666 * pow(ash, 2) * pow(-1 + sh, -4)*
2707 (4. - 2. * sh - 24. * pow(sh, 2) + 23. * pow(sh, 3) -
2708 8. * pow(sh, 4) + pow(sh, 5)) -
2709 pow(lz, 2)*(1.7777777777777777 * Lsh * (6. + sh) *
2711 0.8888888888888888 * pow(-1 + sh, -2) * pow(sh, -2)*
2712 (-28. + 54. * sh - 47. * pow(sh, 2) + 12. * pow(sh, 3) -
2713 6. * pow(sh, 4) + pow(sh, 5))) -
2714 lz * (-3.5555555555555554 * Li2sh * (6. + sh) *
2716 1.7777777777777777 * pow(Lsh, 2) * pow(-1 + sh, -2) +
2717 42.666666666666664 * lmu * pow(sh, -2) +
2718 42.666666666666664 * pow(ash, 2) * pow(-1 + sh, -4)*
2719 (-1. - 3. * sh + pow(sh, 2)) -
2720 1.7777777777777777 * Lsh * pow(-1 + sh, -3) * pow(sh, -2)*
2721 (29. - 82. * sh + 3. * (25. + 4. * Lshb) * pow(sh, 2) +
2722 2. * (-12. + Lshb) * pow(sh, 3)) +
2723 7.111111111111111 * ash * sqrt4sh * sqrtsh *
2724 pow(-1 + sh, -3) * pow(sh, -1)*
2725 (2. + 8. * pow(sh, 2) - 5. * pow(sh, 3) + pow(sh, 4)) -
2726 0.2962962962962963 * pow(-1 + sh, -4) * pow(sh, -2)*
2727 (2. * pow(3.141592653589793, 2) * pow(sh, 2)*
2728 (4. - 11. * sh + pow(sh, 2)) +
2729 3. * pow(-1 + sh, 2)*
2730 (-108. + 216. * sh - 135. * pow(sh, 2) +
2731 32. * pow(sh, 3) - 20. * pow(sh, 4) + 3. * pow(sh, 5)))
2732 ) - 0.14814814814814814 * pow(-1 + sh, -4)*
2733 (-144. * cl3 * (-1. - 3. * sh + pow(sh, 2)) +
2734 pow(sh, -2)*(pow(-1 + sh, 2)*
2735 (-232. + 560. * sh - 757. * pow(sh, 2) +
2736 336. * pow(sh, 3) - 168. * pow(sh, 4) +
2738 2. * pow(3.141592653589793, 2)*
2739 (45. - 151. * sh + 211. * pow(sh, 2) -
2740 135. * pow(sh, 3) - 14. * pow(sh, 4) +
2741 52. * pow(sh, 5) - 16. * pow(sh, 6) + 2. * pow(sh, 7)))
2744 res += pow(z, 2)*(-165.12748700917908 -
2745 411.3508120821892 * (1. - sh) -
2746 670.1268907046218 * pow(1. - sh, 2) -
2747 925.5677817546755 * pow(1. - sh, 3) +
2748 log(muh)*(21.333333333333332 -
2749 42.666666666666664 * pow(1. - sh, 2) -
2750 99.55555555555556 * pow(1. - sh, 3) -
2751 167.11111111111111 * pow(1. - sh, 4)) +
2752 log(z)*(-92.67670236481818 -
2753 123.95171074367174 * (1. - sh) -
2754 143.75841383606902 * pow(1. - sh, 2) -
2755 143.67298964675717 * pow(1. - sh, 3) +
2756 log(muh)*(-42.666666666666664 -
2757 85.33333333333333 * (1. - sh) -
2758 128. * pow(1. - sh, 2) -
2759 170.66666666666666 * pow(1. - sh, 3) -
2760 213.33333333333334 * pow(1. - sh, 4)) -
2761 130.3275572508033 * pow(1. - sh, 4)) -
2762 1177.4249452560534 * pow(1. - sh, 4) +
2763 (26.074074074074073 +
2764 49.925925925925924 * (1. - sh) +
2765 74.4 * pow(1. - sh, 2) +
2766 99.61481481481482 * pow(1. - sh, 3) +
2767 124.74074074074075 * pow(1. - sh, 4)) * pow(log(z), 2));
2772 res += 0.3160493827160494 * pow(3.141592653589793, 2) * pow(sh, -3)*
2773 (3. + 14. * sh + 3. * pow(sh, 2)) * pow(z, 2.5);
2775 res += (62.385647572317914 + 124.77129514463583 * (1. - sh) +
2776 196.51478985280144 * pow(1. - sh, 2) +
2777 277.6161316968147 * pow(1. - sh, 3) +
2778 368.0753206766757 * pow(1. - sh, 4)) * pow(z, 2.5);
2783 res += (4.7407407407407405 * Li3shb * (2. + sh) * pow(-1 + sh, -4) -
2784 lmu * (37.925925925925924 * pow(sh, -3) -
2785 56.888888888888886 * Lsh * pow(sh, -3)) -
2786 7.111111111111111 * Li3sh * pow(-1 + sh, -3) * pow(sh, -1) +
2787 2.3703703703703702 * pow(Lsh, 3) * pow(-1 + sh, -3) *
2788 pow(sh, -1) - 2.8493200667486676 * pow(-1 + sh, -4) *
2789 pow(sh, -1)*(3. + sh + 2. * pow(sh, 2)) +
2790 2.3703703703703702 * cl2 * sqrt4sh * sqrtsh * pow(-1 + sh, -3)*
2791 (4. + 15. * sh - 29. * pow(sh, 2) + 23. * pow(sh, 3) -
2792 8. * pow(sh, 4) + pow(sh, 5)) +
2793 0.7901234567901234 * ash * sqrt4sh * sqrtsh * pow(-1 + sh, -3)*
2794 (-8. - 21. * sh - 17. * pow(sh, 2) + 23. * pow(sh, 3) -
2795 8. * pow(sh, 4) + pow(sh, 5)) +
2796 0.3950617283950617 * Li2sh * pow(-1 + sh, -4) * pow(sh, -3)*
2797 (-36. + 144. * sh - 234. * pow(sh, 2) + 208. * pow(sh, 3) -
2798 73. * pow(sh, 4) + 9. * pow(sh, 5)) -
2799 Lsh * (14.222222222222221 * (2. + sh) * pow(ash, 2) *
2801 7.111111111111111 * Li2sh * pow(-1 + sh, -3) *
2802 pow(sh, -1) - 2.3703703703703702 * ash * sqrt4sh * sqrtsh *
2804 (4. + 15. * sh - 29. * pow(sh, 2) + 23. * pow(sh, 3) -
2805 8. * pow(sh, 4) + pow(sh, 5)) -
2806 0.03292181069958848 * pow(-1 + sh, -4) * pow(sh, -4)*
2807 (30. + (1508. - 432. * Lshb) * sh +
2808 2. * (-3427. + 864. * Lshb) * pow(sh, 2) +
2809 36. * (296. - 78. * Lshb +
2810 5. * pow(3.141592653589793, 2)) * pow(sh, 3) +
2811 (-6541. + 2496. * Lshb -
2812 36. * pow(3.141592653589793, 2)) * pow(sh, 4) +
2813 (1552. - 876. * Lshb +
2814 72. * pow(3.141592653589793, 2)) * pow(sh, 5) +
2815 27. * (3. + 4. * Lshb) * pow(sh, 6))) -
2816 2.3703703703703702 * pow(ash, 2) * pow(-1 + sh, -4)*
2817 (16. - 22. * sh - 171. * pow(sh, 2) + 336. * pow(sh, 3) -
2818 300. * pow(sh, 4) + 141. * pow(sh, 5) - 33. * pow(sh, 6) +
2819 3. * pow(sh, 7)) - 0.09876543209876543 * pow(Lsh, 2) *
2820 pow(-1 + sh, -5) * pow(sh, -4)*
2821 (12. + 494. * sh - 2533. * pow(sh, 2) +
2822 (4871. - 36. * Lshb) * pow(sh, 3) +
2823 (-4634. + 72. * Lshb) * pow(sh, 4) -
2824 4. * (-541. + 9. * Lshb) * pow(sh, 5) - 428. * pow(sh, 6) +
2826 pow(lz, 2)*(2.3703703703703702 * Lsh * (2. + sh) *
2828 0.09876543209876543 * pow(-1 + sh, -3) * pow(sh, -4)*
2829 (-3. - 491. * sh + 1482. * pow(sh, 2) -
2830 1467. * pow(sh, 3) + 572. * pow(sh, 4) -
2831 132. * pow(sh, 5) + 321. * pow(sh, 6) -
2832 378. * pow(sh, 7) + 222. * pow(sh, 8) -
2833 60. * pow(sh, 9) + 6. * pow(sh, 10))) -
2834 lz * (-4.7407407407407405 * Li2sh * (2. + sh) *
2836 28.444444444444443 * (2. + sh) * pow(ash, 2) *
2838 56.888888888888886 * lmu * pow(sh, -3) +
2839 3.5555555555555554 * pow(Lsh, 2) * pow(-1 + sh, -3) *
2840 pow(sh, -1) - 4.7407407407407405 * ash * sqrt4sh * sqrtsh *
2842 (4. + 15. * sh - 29. * pow(sh, 2) + 23. * pow(sh, 3) -
2843 8. * pow(sh, 4) + pow(sh, 5)) +
2844 0.19753086419753085 * Lsh * pow(-1 + sh, -4) * pow(sh, -4)*
2845 (6. + 494. * sh - 1991. * pow(sh, 2) +
2846 2886. * pow(sh, 3) -
2847 6. * (293. + 8. * Lshb) * pow(sh, 4) -
2848 6. * (-65. + 4. * Lshb) * pow(sh, 5) + 9. * pow(sh, 6)) -
2849 0.03292181069958848 * pow(-1 + sh, -3) * pow(sh, -4)*
2850 (15. + 1538. * sh - 4641. * pow(sh, 2) +
2851 4446. * pow(sh, 3) - 1361. * pow(sh, 4) +
2852 498. * pow(sh, 5) - 1215. * pow(sh, 6) +
2853 606. * pow(sh, 7) - 84. * pow(sh, 8) - 24. * pow(sh, 9) +
2854 6. * pow(sh, 10))) +
2855 0.01646090534979424 * pow(-1 + sh, -4) * pow(sh, -4)*
2856 (19. - 7836. * sh + 30701. * pow(sh, 2) -
2857 45359. * pow(sh, 3) + (31147. + 1728. * cl3) * pow(sh, 4) +
2858 12. * (-701. + 72. * cl3) * pow(sh, 5) - 343. * pow(sh, 6) -
2859 2371. * pow(sh, 7) + 4810. * pow(sh, 8) -
2860 3164. * pow(sh, 9) + 902. * pow(sh, 10) -
2862 6. * pow(3.141592653589793, 2)*
2863 (-6. - 338. * sh + 1349. * pow(sh, 2) -
2864 1798. * pow(sh, 3) + 914. * pow(sh, 4) -
2865 42. * pow(sh, 5) - 299. * pow(sh, 6) +
2866 440. * pow(sh, 7) - 400. * pow(sh, 8) +
2867 188. * pow(sh, 9) - 44. * pow(sh, 10) + 4. * pow(sh, 11)))
2870 res += pow(z, 3)*(-563.8946426367565 -
2871 1680.2345038491285 * (1. - sh) -
2872 3329.6412593857735 * pow(1. - sh, 2) -
2873 5423.45927216734 * pow(1. - sh, 3) +
2874 log(muh)*(-37.925925925925924 -
2875 170.66666666666666 * (1. - sh) -
2876 426.6666666666667 * pow(1. - sh, 2) -
2877 824.8888888888889 * pow(1. - sh, 3) -
2878 1379.5555555555557 * pow(1. - sh, 4)) -
2879 7956.241038705334 * pow(1. - sh, 4) +
2880 log(z)*(-34.97995438478488 -
2881 30.044118068027284 * (1. - sh) +
2882 62.84401734242212 * pow(1. - sh, 2) +
2883 320.43564180164736 * pow(1. - sh, 3) +
2884 log(muh)*(-56.888888888888886 -
2885 170.66666666666666 * (1. - sh) -
2886 341.3333333333333 * pow(1. - sh, 2) -
2887 568.8888888888889 * pow(1. - sh, 3) -
2888 853.3333333333334 * pow(1. - sh, 4)) +
2889 749.5923384459219 * pow(1. - sh, 4)) +
2890 (54.22222222222222 + 150.45925925925926 * (1. - sh) +
2891 299.97037037037035 * pow(1. - sh, 2) +
2892 503.0405643738977 * pow(1. - sh, 3) +
2893 755.8095238095239 * pow(1. - sh, 4)) * pow(log(z), 2));
2898 res += 0.02708994708994709 * pow(3.141592653589793, 2) * pow(sh, -5)*
2899 (15. + 108. * sh + 314. * pow(sh, 2) + 108. * pow(sh, 3) +
2900 15. * pow(sh, 4)) * pow(z, 3.5);
2902 res += (149.725554173563 + 449.17666252068904 * (1. - sh) +
2903 943.270991293447 * pow(1. - sh, 2) +
2904 1676.9262067439056 * pow(1. - sh, 3) +
2905 2699.0704810394973 * pow(1. - sh, 4)) * pow(z, 3.5);
2910 res += (-1. * lmu * (117.33333333333333 * pow(sh, -4) -
2911 128. * Lsh * pow(sh, -4)) +
2912 14.246600333743338 * pow(-1 + sh, -3) * pow(sh, -2) -
2913 11.851851851851851 * Li3sh * pow(-1 + sh, -3) * pow(sh, -2) +
2914 3.950617283950617 * pow(Lsh, 3) * pow(-1 + sh, -3) *
2915 pow(sh, -2) - 3.5555555555555554 * cl2 * sqrt4sh * sqrtsh *
2916 (-6. - 13. * sh + 20. * pow(sh, 2) - 8. * pow(sh, 3) +
2917 pow(sh, 4)) + 0.5925925925925926 * Li2sh *
2918 pow(-1 + sh, -5) * pow(sh, -4)*
2919 (54. - 268. * sh + 544. * pow(sh, 2) - 570. * pow(sh, 3) +
2920 292. * pow(sh, 4) - 58. * pow(sh, 5) - 6. * pow(sh, 6) +
2921 3. * pow(sh, 7)) - 0.5925925925925926 * ash * sqrt4sh *
2922 sqrtsh * pow(-1 + sh, -3)*
2923 (-6. - 41. * sh - 239. * pow(sh, 2) + 671. * pow(sh, 3) -
2924 680. * pow(sh, 4) + 329. * pow(sh, 5) - 77. * pow(sh, 6) +
2925 7. * pow(sh, 7)) + 10.666666666666666 * sh * pow(ash, 2) *
2926 pow(-1 + sh, -4)*(-24. - 3. * sh + 174. * pow(sh, 2) -
2927 390. * pow(sh, 3) + 411. * pow(sh, 4) - 241. * pow(sh, 5) +
2928 80. * pow(sh, 6) - 14. * pow(sh, 7) + pow(sh, 8)) -
2929 Lsh * (-11.851851851851851 * Li2sh * pow(-1 + sh, -3) *
2930 pow(sh, -2) + 3.5555555555555554 * ash * sqrt4sh * sqrtsh *
2931 (-6. - 13. * sh + 20. * pow(sh, 2) - 8. * pow(sh, 3) +
2932 pow(sh, 4)) - 0.01646090534979424 * pow(-1 + sh, -5) *
2933 pow(sh, -6)*(-21. + 10. * sh +
2934 2. * (-847. + 972. * Lshb) * pow(sh, 2) +
2935 (11470. - 9648. * Lshb) * pow(sh, 3) +
2936 (-27667. + 19584. * Lshb -
2937 600. * pow(3.141592653589793, 2)) * pow(sh, 4) +
2938 8. * (3982. - 2565. * Lshb +
2939 150. * pow(3.141592653589793, 2)) * pow(sh, 5) +
2940 2. * (-9179. + 5256. * Lshb -
2941 300. * pow(3.141592653589793, 2)) * pow(sh, 6) -
2942 6. * (-607. + 348. * Lshb) * pow(sh, 7) -
2943 6. * (7. + 36. * Lshb) * pow(sh, 8) +
2944 6. * (5. + 18. * Lshb) * pow(sh, 9))) -
2945 0.09876543209876543 * pow(Lsh, 2) * pow(-1 + sh, -6) *
2946 pow(sh, -6)*(-6. - 14. * sh - 1215. * pow(sh, 2) +
2947 7695. * pow(sh, 3) - 19257. * pow(sh, 4) -
2948 60. * Lshb * pow(-1 + sh, 3) * pow(sh, 4) +
2949 25233. * pow(sh, 5) - 18375. * pow(sh, 6) +
2950 6983. * pow(sh, 7) - 1056. * pow(sh, 8) - 27. * pow(sh, 9) +
2952 0.04938271604938271 * pow(lz, 2) * pow(sh, -6)*
2953 (-3. - 25. * sh - 2580. * pow(sh, 2) + 73. * pow(sh, 4) +
2954 3. * pow(sh, 5) + 180. * pow(sh, 6) + 36. * pow(sh, 7) -
2955 738. * pow(sh, 8) + 612. * pow(sh, 9) - 180. * pow(sh, 10) +
2956 18. * pow(sh, 11)) -
2957 lz * (128. * lmu * pow(sh, -4) +
2958 5.925925925925926 * pow(Lsh, 2) * pow(-1 + sh, -3) *
2959 pow(sh, -2) + 7.111111111111111 * ash * sqrt4sh * sqrtsh *
2960 (-6. - 13. * sh + 20. * pow(sh, 2) - 8. * pow(sh, 3) +
2961 pow(sh, 4)) + 0.19753086419753085 * Lsh *
2962 pow(-1 + sh, -5) * pow(sh, -6)*
2963 (-3. - 10. * sh - 1240. * pow(sh, 2) +
2964 6380. * pow(sh, 3) - 12722. * pow(sh, 4) +
2965 12478. * pow(sh, 5) - 6074. * pow(sh, 6) +
2966 1158. * pow(sh, 7) + 3. * pow(sh, 8) + 3. * pow(sh, 9)) +
2967 0.0016460905349794238 * pow(-1 + sh, -4) * pow(sh, -6)*
2968 (105. + 55. * sh + 18530. * pow(sh, 2) -
2969 72330. * pow(sh, 3) + 91890. * pow(sh, 4) -
2970 26202. * pow(sh, 5) - 21156. * pow(sh, 6) +
2971 1154. * pow(sh, 7) + 54349. * pow(sh, 8) -
2972 76131. * pow(sh, 9) + 4716. * pow(sh, 10) +
2973 82440. * pow(sh, 11) - 83970. * pow(sh, 12) +
2974 36900. * pow(sh, 13) - 7740. * pow(sh, 14) +
2975 630. * pow(sh, 15))) -
2976 0.000027434842249657064 * pow(-1 + sh, -5) * pow(sh, -6)*
2977 (2775. - 5550. * sh - 9.841575e6 * pow(sh, 2) +
2978 4.85214e7 * pow(sh, 3) - 9.54242e7 * pow(sh, 4) +
2979 9.3237052e7 * pow(sh, 5) - 4.3743086e7 * pow(sh, 6) +
2980 4.46555e6 * pow(sh, 7) + 7.670645e6 * pow(sh, 8) -
2981 4.88458e6 * pow(sh, 9) - 1.3483707e7 * pow(sh, 10) +
2982 3.4549776e7 * pow(sh, 11) - 3.627225e7 * pow(sh, 12) +
2983 2.100825e7 * pow(sh, 13) - 6.9246e6 * pow(sh, 14) +
2984 1.21095e6 * pow(sh, 15) - 86850. * pow(sh, 16) +
2985 1200. * pow(3.141592653589793, 2)*
2986 (-9. - 30. * sh - 2694. * pow(sh, 2) +
2987 13896. * pow(sh, 3) - 27086. * pow(sh, 4) +
2988 25318. * pow(sh, 5) - 11444. * pow(sh, 6) +
2989 2810. * pow(sh, 7) - 1033. * pow(sh, 8) -
2990 2857. * pow(sh, 9) + 10122. * pow(sh, 10) -
2991 14418. * pow(sh, 11) + 11736. * pow(sh, 12) -
2992 5778. * pow(sh, 13) + 1692. * pow(sh, 14) -
2993 270. * pow(sh, 15) + 18. * pow(sh, 16)))) * pow(z, 4);
2995 res += pow(z, 4)*(-1424.5308091578474 -
2996 5728.7251525464435 * (1. - sh) -
2997 14040.07057463595 * pow(1. - sh, 2) -
2998 27178.850860660903 * pow(1. - sh, 3) +
2999 log(muh)*(-117.33333333333333 -
3000 597.3333333333334 * (1. - sh) -
3001 1749.3333333333333 * pow(1. - sh, 2) -
3002 3925.3333333333335 * pow(1. - sh, 3) -
3003 7509.333333333333 * pow(1. - sh, 4)) -
3004 46368.872894680535 * pow(1. - sh, 4) +
3005 log(z)*(34.25652635469211 +
3006 169.5621944588141 * (1. - sh) +
3007 865.728317571408 * pow(1. - sh, 2) +
3008 2770.286236159471 * pow(1. - sh, 3) +
3009 log(muh)*(-128. - 512. * (1. - sh) -
3010 1280. * pow(1. - sh, 2) -
3011 2560. * pow(1. - sh, 3) - 4480. * pow(1. - sh, 4)
3012 ) + 6293.656266280077 * pow(1. - sh, 4)) +
3013 (128.59259259259258 +
3014 497.77777777777777 * (1. - sh) +
3015 1274.962962962963 * pow(1. - sh, 2) +
3016 2588.641975308642 * pow(1. - sh, 3) +
3017 4550.617283950617 * pow(1. - sh, 4)) * pow(log(z), 2));
3022 res += 0.006449987402368355 * pow(3.141592653589793, 2) * pow(sh, -7)*
3023 (35. + 330. * sh + 1389. * pow(sh, 2) + 3212. * pow(sh, 3) +
3024 1389. * pow(sh, 4) + 330. * pow(sh, 5) + 35. * pow(sh, 6)) *
3027 res += (427.7872976387514 + 1711.1491905550056 * (1. - sh) +
3028 4470.377260324953 * pow(1. - sh, 2) +
3029 9518.267372462218 * pow(1. - sh, 3) +
3030 17894.4954414067 * pow(1. - sh, 4)) * pow(z, 4.5);
3035 res += (-1. * lmu * (352.7111111111111 * pow(sh, -5) -
3036 341.3333333333333 * Lsh * pow(sh, -5)) +
3037 29.917860700861013 * pow(-1 + sh, -3) * pow(sh, -3) -
3038 24.88888888888889 * Li3sh * pow(-1 + sh, -3) * pow(sh, -3) +
3039 8.296296296296296 * pow(Lsh, 3) * pow(-1 + sh, -3) *
3040 pow(sh, -3) + 2.3703703703703702 * cl2 * sh * sqrt4sh * sqrtsh *
3041 (44. + 108. * sh - 221. * pow(sh, 2) + 132. * pow(sh, 3) -
3042 33. * pow(sh, 4) + 3. * pow(sh, 5)) +
3043 0.23703703703703705 * ash * sh * sqrt4sh * sqrtsh *
3044 pow(-1 + sh, -3)*(-496. + 326. * sh + 6147. * pow(sh, 2) -
3045 16495. * pow(sh, 3) + 18680. * pow(sh, 4) -
3046 11251. * pow(sh, 5) + 3760. * pow(sh, 6) -
3047 658. * pow(sh, 7) + 47. * pow(sh, 8)) +
3048 1.1851851851851851 * Li2sh * pow(-1 + sh, -6) * pow(sh, -5)*
3049 (-72. + 429. * sh - 1073. * pow(sh, 2) + 1442. * pow(sh, 3) -
3050 1092. * pow(sh, 4) + 435. * pow(sh, 5) - 47. * pow(sh, 6) -
3051 6. * pow(sh, 7) - pow(sh, 8) + pow(sh, 9)) -
3052 7.111111111111111 * sh * pow(ash, 2) * pow(-1 + sh, -4)*
3053 (-24. + 162. * sh + 114. * pow(sh, 2) - 1809. * pow(sh, 3) +
3054 4224. * pow(sh, 4) - 4957. * pow(sh, 5) +
3055 3423. * pow(sh, 6) - 1445. * pow(sh, 7) +
3056 366. * pow(sh, 8) - 51. * pow(sh, 9) + 3. * pow(sh, 10)) -
3057 Lsh * (-24.88888888888889 * Li2sh * pow(-1 + sh, -3) *
3058 pow(sh, -3) - 2.3703703703703702 * ash * sh * sqrt4sh *
3059 sqrtsh * (44. + 108. * sh - 221. * pow(sh, 2) +
3060 132. * pow(sh, 3) - 33. * pow(sh, 4) + 3. * pow(sh, 5)) +
3061 0.0003292181069958848 * pow(-1 + sh, -7) * pow(sh, -8)*
3062 (486. - 594. * sh - 8634. * pow(sh, 2) -
3063 132632. * pow(sh, 3) + 806792. * pow(sh, 4) -
3064 1.649394e6 * pow(sh, 5) +
3065 63000. * pow(3.141592653589793, 2) * pow(-1 + sh, 4) *
3066 pow(sh, 5) + 1.299724e6 * pow(sh, 6) +
3067 289454. * pow(sh, 7) - 1.252893e6 * pow(sh, 8) +
3068 839053. * pow(sh, 9) -
3069 3600. * Lshb * (-1. + sh) * pow(sh, 3)*
3070 (-72. + 429. * sh - 1073. * pow(sh, 2) +
3071 1442. * pow(sh, 3) - 1092. * pow(sh, 4) +
3072 435. * pow(sh, 5) - 47. * pow(sh, 6) -
3073 6. * pow(sh, 7) - pow(sh, 8) + pow(sh, 9)) -
3074 206662. * pow(sh, 10) + 15750. * pow(sh, 11) +
3075 75. * pow(sh, 12) - 525. * pow(sh, 13))) -
3076 0.009876543209876543 * pow(Lsh, 2) * pow(-1 + sh, -7) *
3077 pow(sh, -8)*(36. + 126. * sh - 84. * pow(sh, 2) +
3078 37024. * pow(sh, 3) - 269605. * pow(sh, 4) +
3079 818292. * pow(sh, 5) -
3080 1260. * Lshb * pow(-1 + sh, 4) * pow(sh, 5) -
3081 1.360586e6 * pow(sh, 6) + 1.344224e6 * pow(sh, 7) -
3082 788604. * pow(sh, 8) + 252535. * pow(sh, 9) -
3083 33268. * pow(sh, 10) - 300. * pow(sh, 11) -
3084 120. * pow(sh, 12) + 60. * pow(sh, 13)) +
3085 0.0049382716049382715 * pow(lz, 2) * pow(sh, -8)*
3086 (18. + 189. * sh + 903. * pow(sh, 2) + 78416. * pow(sh, 3) -
3087 630. * pow(sh, 4) - 1383. * pow(sh, 5) -
3088 1839. * pow(sh, 6) - 18. * pow(sh, 7) - 3120. * pow(sh, 8) -
3089 8160. * pow(sh, 9) - 1440. * pow(sh, 10) +
3090 46560. * pow(sh, 11) - 51720. * pow(sh, 12) +
3091 23040. * pow(sh, 13) - 4680. * pow(sh, 14) +
3092 360. * pow(sh, 15)) -
3093 lz * (341.3333333333333 * lmu * pow(sh, -5) +
3094 12.444444444444445 * pow(Lsh, 2) * pow(-1 + sh, -3) *
3095 pow(sh, -3) - 4.7407407407407405 * ash * sh * sqrt4sh *
3096 sqrtsh * (44. + 108. * sh - 221. * pow(sh, 2) +
3097 132. * pow(sh, 3) - 33. * pow(sh, 4) + 3. * pow(sh, 5)) +
3098 0.019753086419753086 * Lsh * pow(-1 + sh, -6) * pow(sh, -8)*
3099 (18. + 81. * sh + 39. * pow(sh, 2) + 37525. * pow(sh, 3) -
3100 231198. * pow(sh, 4) + 582087. * pow(sh, 5) -
3101 774221. * pow(sh, 6) + 575781. * pow(sh, 7) -
3102 227535. * pow(sh, 8) + 38728. * pow(sh, 9) -
3103 450. * pow(sh, 10) + 90. * pow(sh, 11) + 15. * pow(sh, 12)
3104 ) - 0.000047031158142269256 * pow(-1 + sh, -5) *
3105 pow(sh, -8)*(1701. + 1323. * sh - 29274. * pow(sh, 2) -
3106 957341. * pow(sh, 3) + 5.622085e6 * pow(sh, 4) -
3107 1.3372289e7 * pow(sh, 5) + 1.7244458e7 * pow(sh, 6) -
3108 1.2810919e7 * pow(sh, 7) + 6.090011e6 * pow(sh, 8) -
3109 4.03983e6 * pow(sh, 9) + 4.371768e6 * pow(sh, 10) +
3110 4.551324e6 * pow(sh, 11) - 1.1897697e7 * pow(sh, 12) -
3111 1.323264e7 * pow(sh, 13) + 6.277908e7 * pow(sh, 14) -
3112 8.508528e7 * pow(sh, 15) + 6.212682e7 * pow(sh, 16) -
3113 2.701314e7 * pow(sh, 17) + 6.98418e6 * pow(sh, 18) -
3114 990360. * pow(sh, 19) + 59220. * pow(sh, 20))) -
3115 1.679684219366759e-7 * pow(-1 + sh, -6) * pow(sh, -8)*
3116 (-161406. + 318843. * sh + 1.910657e6 * pow(sh, 2) +
3117 3.755399935e9 * pow(sh, 3) -
3118 2.2256359629e10 * pow(sh, 4) +
3119 5.4715328248e10 * pow(sh, 5) -
3120 7.1399970054e10 * pow(sh, 6) +
3121 5.2256870574e10 * pow(sh, 7) -
3122 2.1211105332e10 * pow(sh, 8) +
3123 5.981785659e9 * pow(sh, 9) -
3124 2.944201467e9 * pow(sh, 10) -
3125 1.246759469e9 * pow(sh, 11) +
3126 2.466843039e9 * pow(sh, 12) +
3127 1.8171607842e10 * pow(sh, 13) -
3128 6.053910752e10 * pow(sh, 14) +
3129 8.996385496e10 * pow(sh, 15) -
3130 7.96292826e10 * pow(sh, 16) +
3131 4.509057616e10 * pow(sh, 17) -
3132 1.647058952e10 * pow(sh, 18) +
3133 3.74875284e9 * pow(sh, 19) - 4.8245988e8 * pow(sh, 20) +
3134 2.674812e7 * pow(sh, 21) -
3135 58800. * pow(3.141592653589793, 2)*
3136 (-18. - 81. * sh - 39. * pow(sh, 2) -
3137 28525. * pow(sh, 3) + 176088. * pow(sh, 4) -
3138 439209. * pow(sh, 5) + 573427. * pow(sh, 6) -
3139 414687. * pow(sh, 7) + 157597. * pow(sh, 8) -
3140 24058. * pow(sh, 9) + 3132. * pow(sh, 10) -
3141 3736. * pow(sh, 11) - 77479. * pow(sh, 12) +
3142 318588. * pow(sh, 13) - 608560. * pow(sh, 14) +
3143 702440. * pow(sh, 15) - 529920. * pow(sh, 16) +
3144 267160. * pow(sh, 17) - 89120. * pow(sh, 18) +
3145 18840. * pow(sh, 19) - 2280. * pow(sh, 20) +
3146 120. * pow(sh, 21)))) * pow(z, 5);
3148 res += pow(z, 5)*(-4202.986795121401 -
3149 21563.981801901395 * (1. - sh) -
3150 63115.765734768596 * pow(1. - sh, 2) -
3151 142119.7662710934 * pow(1. - sh, 3) +
3152 log(muh)*(-352.7111111111111 -
3153 2104.8888888888887 * (1. - sh) -
3154 7168. * pow(1. - sh, 2) -
3155 18432. * pow(1. - sh, 3) -
3156 39850.666666666664 * pow(1. - sh, 4)) -
3157 277830.4069154384 * pow(1. - sh, 4) +
3158 log(z)*(251.91238577751216 +
3159 956.915478778669 * (1. - sh) +
3160 4881.515489147685 * pow(1. - sh, 2) +
3161 16699.40413911222 * pow(1. - sh, 3) +
3162 log(muh)*(-341.3333333333333 -
3163 1706.6666666666667 * (1. - sh) -
3164 5120. * pow(1. - sh, 2) -
3165 11946.666666666666 * pow(1. - sh, 3) -
3166 23893.333333333332 * pow(1. - sh, 4)) +
3167 41147.84749220197 * pow(1. - sh, 4)) +
3168 (377.758024691358 + 1862.716049382716 * (1. - sh) +
3169 5810.814814814815 * pow(1. - sh, 2) +
3170 13779.160493827161 * pow(1. - sh, 3) +
3171 27666.172839506173 * pow(1. - sh, 4)) * pow(log(z), 2));
3176 res += 0.0004560597153189746 * pow(3.141592653589793, 2) * pow(sh, -9)*
3177 (315. + 3640. * sh + 19380. * pow(sh, 2) + 62088. * pow(sh, 3) +
3178 124834. * pow(sh, 4) + 62088. * pow(sh, 5) +
3179 19380. * pow(sh, 6) + 3640. * pow(sh, 7) + 315. * pow(sh, 8)) *
3182 res += (1330.8938148761156 + 6654.469074380578 * (1. - sh) +
3183 20761.943512067406 * pow(1. - sh, 2) +
3184 51372.50125421806 * pow(1. - sh, 3) +
3185 110145.72275758653 * pow(1. - sh, 4)) * pow(z, 5.5);
3190 res += (-1. * lmu * (1097.4814814814815 * pow(sh, -6) -
3191 995.5555555555555 * Lsh * pow(sh, -6)) +
3192 71.80286568206643 * pow(-1 + sh, -3) * pow(sh, -4) -
3193 59.733333333333334 * Li3sh * pow(-1 + sh, -3) * pow(sh, -4) +
3194 19.91111111111111 * pow(Lsh, 3) * pow(-1 + sh, -3) *
3195 pow(sh, -4) - 2.765432098765432 * cl2 * sh * sqrt4sh * sqrtsh *
3196 (20. - 194. * sh - 657. * pow(sh, 2) + 1581. * pow(sh, 3) -
3197 1239. * pow(sh, 4) + 462. * pow(sh, 5) - 84. * pow(sh, 6) +
3198 6. * pow(sh, 7)) - 0.02633744855967078 * ash * sh * sqrt4sh *
3199 sqrtsh * pow(-1 + sh, -3)*
3200 (-3230. + 42428. * sh + 5415. * pow(sh, 2) -
3201 547510. * pow(sh, 3) + 1.454143e6 * pow(sh, 4) -
3202 1.806873e6 * pow(sh, 5) + 1.280646e6 * pow(sh, 6) -
3203 546171. * pow(sh, 7) + 138714. * pow(sh, 8) -
3204 19329. * pow(sh, 9) + 1137. * pow(sh, 10)) +
3205 8.296296296296296 * pow(ash, 2) * pow(-1 + sh, -4) * pow(sh, 2)*
3206 (210. - 870. * sh - 1395. * pow(sh, 2) +
3207 13113. * pow(sh, 3) - 31009. * pow(sh, 4) +
3208 39370. * pow(sh, 5) - 30843. * pow(sh, 6) +
3209 15585. * pow(sh, 7) - 5091. * pow(sh, 8) +
3210 1038. * pow(sh, 9) - 120. * pow(sh, 10) + 6. * pow(sh, 11))\
3211 + 0.019753086419753086 * Li2sh * pow(-1 + sh, -7) *
3212 pow(sh, -6)*(12600. - 87696. * sh + 262296. * pow(sh, 2) -
3213 436590. * pow(sh, 3) + 435240. * pow(sh, 4) -
3214 257040. * pow(sh, 5) + 79436. * pow(sh, 6) -
3215 8456. * pow(sh, 7) - 360. * pow(sh, 9) + 45. * pow(sh, 11))\
3216 - Lsh * (-59.733333333333334 * Li2sh * pow(-1 + sh, -3) *
3217 pow(sh, -4) + 2.765432098765432 * ash * sh * sqrt4sh *
3218 sqrtsh * (20. - 194. * sh - 657. * pow(sh, 2) +
3219 1581. * pow(sh, 3) - 1239. * pow(sh, 4) +
3220 462. * pow(sh, 5) - 84. * pow(sh, 6) + 6. * pow(sh, 7)) +
3221 0.00006584362139917696 * pow(-1 + sh, -8) * pow(sh, -10)*
3222 (-1320. + 1416. * sh + 30540. * pow(sh, 2) +
3223 40260. * pow(sh, 3) + 6.70004e6 * pow(sh, 4) -
3224 5.379436e7 * pow(sh, 5) + 1.79724938e8 * pow(sh, 6) +
3225 756000. * pow(3.141592653589793, 2) * pow(-1 + sh, 5) *
3226 pow(sh, 6) - 3.3536944e8 * pow(sh, 7) +
3227 3.8325104e8 * pow(sh, 8) - 2.7317989e8 * pow(sh, 9) +
3228 1.17246435e8 * pow(sh, 10) -
3229 2.7037239e7 * pow(sh, 11) -
3230 300. * Lshb * (-1. + sh) * pow(sh, 4)*
3231 (12600. - 87696. * sh + 262296. * pow(sh, 2) -
3232 436590. * pow(sh, 3) + 435240. * pow(sh, 4) -
3233 257040. * pow(sh, 5) + 79436. * pow(sh, 6) -
3234 8456. * pow(sh, 7) - 360. * pow(sh, 9) +
3235 45. * pow(sh, 11)) + 2.36385e6 * pow(sh, 12) -
3236 31650. * pow(sh, 13) + 60375. * pow(sh, 14) -
3237 3780. * pow(sh, 15) - 1215. * pow(sh, 16))) -
3238 0.0019753086419753087 * pow(Lsh, 2) * pow(-1 + sh, -8) *
3239 pow(sh, -10)*(-120. - 552. * sh - 120. * pow(sh, 2) +
3240 3480. * pow(sh, 3) - 612380. * pow(sh, 4) +
3241 4.99354e6 * pow(sh, 5) - 1.7689517e7 * pow(sh, 6) -
3242 15120. * Lshb * pow(-1 + sh, 5) * pow(sh, 6) +
3243 3.5499895e7 * pow(sh, 7) - 4.421963e7 * pow(sh, 8) +
3244 3.502108e7 * pow(sh, 9) - 1.719707e7 * pow(sh, 10) +
3245 4.767574e6 * pow(sh, 11) - 567440. * pow(sh, 12) +
3246 1800. * pow(sh, 13) - 1800. * pow(sh, 14) -
3247 225. * pow(sh, 15) + 225. * pow(sh, 16)) -
3248 0.0019753086419753087 * pow(lz, 2) * pow(sh, -10)*
3249 (-30. - 378. * sh - 2214. * pow(sh, 2) - 7938. * pow(sh, 3) -
3250 638000. * pow(sh, 4) + 7560. * pow(sh, 5) +
3251 6273. * pow(sh, 6) + 13689. * pow(sh, 7) +
3252 10908. * pow(sh, 8) + 30. * pow(sh, 9) +
3253 20300. * pow(sh, 10) + 27300. * pow(sh, 11) +
3254 113400. * pow(sh, 12) + 49700. * pow(sh, 13) -
3255 880950. * pow(sh, 14) + 1.19385e6 * pow(sh, 15) -
3256 706650. * pow(sh, 16) + 216300. * pow(sh, 17) -
3257 33600. * pow(sh, 18) + 2100. * pow(sh, 19)) -
3258 lz * (995.5555555555555 * lmu * pow(sh, -6) +
3259 29.866666666666667 * pow(Lsh, 2) * pow(-1 + sh, -3) *
3260 pow(sh, -4) + 5.530864197530864 * ash * sh * sqrt4sh *
3261 sqrtsh * (20. - 194. * sh - 657. * pow(sh, 2) +
3262 1581. * pow(sh, 3) - 1239. * pow(sh, 4) +
3263 462. * pow(sh, 5) - 84. * pow(sh, 6) + 6. * pow(sh, 7)) +
3264 0.003950617283950617 * Lsh * pow(-1 + sh, -7) *
3265 pow(sh, -10)*(-60. - 336. * sh - 396. * pow(sh, 2) +
3266 1344. * pow(sh, 3) - 614396. * pow(sh, 4) +
3267 4.375784e6 * pow(sh, 5) - 1.3247019e7 * pow(sh, 6) +
3268 2.2150786e7 * pow(sh, 7) - 2.2142404e7 * pow(sh, 8) +
3269 1.3260636e7 * pow(sh, 9) - 4.429894e6 * pow(sh, 10) +
3270 647380. * pow(sh, 11) - 2940. * pow(sh, 12) -
3271 1740. * pow(sh, 13) + 585. * pow(sh, 14) +
3272 45. * pow(sh, 15)) +
3273 9.40623162845385e-6 * pow(-1 + sh, -6) * pow(sh, -10)*
3274 (4620. + 4284. * sh - 102942. * pow(sh, 2) -
3275 351078. * pow(sh, 3) - 4.7830454e7 * pow(sh, 4) +
3276 3.0838269e8 * pow(sh, 5) - 8.19066822e8 * pow(sh, 6) +
3277 1.178613196e9 * pow(sh, 7) -
3278 9.89156856e8 * pow(sh, 8) +
3279 4.74656612e8 * pow(sh, 9) - 9.292233e7 * pow(sh, 10) -
3280 5.9056617e7 * pow(sh, 11) +
3281 8.9385341e7 * pow(sh, 12) -
3282 1.01542571e8 * pow(sh, 13) +
3283 3.36952887e8 * pow(sh, 14) -
3284 4.6151941e8 * pow(sh, 15) -
3285 1.1490745e9 * pow(sh, 16) +
3286 5.1849126e9 * pow(sh, 17) -
3287 8.94761595e9 * pow(sh, 18) +
3288 9.1281316e9 * pow(sh, 19) -
3289 6.0729627e9 * pow(sh, 20) +
3290 2.7087795e9 * pow(sh, 21) - 8.039283e8 * pow(sh, 22) +
3291 1.522416e8 * pow(sh, 23) - 1.66278e7 * pow(sh, 24) +
3292 795900. * pow(sh, 25))) -
3293 7.465263197185596e-9 * pow(-1 + sh, -7) * pow(sh, -10)*
3294 (1.60524e6 - 3.372768e6 * sh - 2.1255318e7 * pow(sh, 2) +
3295 5.1364152e7 * pow(sh, 3) -
3296 2.17838821368e11 * pow(sh, 4) +
3297 1.503468880032e12 * pow(sh, 5) -
3298 4.435848139401e12 * pow(sh, 6) +
3299 7.250465882724e12 * pow(sh, 7) -
3300 7.100583095886e12 * pow(sh, 8) +
3301 4.194932382824e12 * pow(sh, 9) -
3302 1.400004505326e12 * pow(sh, 10) +
3303 1.30381507564e11 * pow(sh, 11) +
3304 2.09548894394e11 * pow(sh, 12) -
3305 3.42677229216e11 * pow(sh, 13) +
3306 7.84252942969e11 * pow(sh, 14) -
3307 7.91328639346e11 * pow(sh, 15) -
3308 3.49537340577e12 * pow(sh, 16) +
3309 1.51956042427e13 * pow(sh, 17) -
3310 2.91706462341e13 * pow(sh, 18) +
3311 3.468763704735e13 * pow(sh, 19) -
3312 2.790223133985e13 * pow(sh, 20) +
3313 1.564098014115e13 * pow(sh, 21) -
3314 6.1223554161e12 * pow(sh, 22) +
3315 1.63963647855e12 * pow(sh, 23) -
3316 2.859293178e11 * pow(sh, 24) +
3317 2.92020939e10 * pow(sh, 25) - 1.3226913e9 * pow(sh, 26) +
3318 88200. * pow(3.141592653589793, 2)*
3319 (-180. - 1008. * sh - 1188. * pow(sh, 2) +
3320 4032. * pow(sh, 3) - 1.452588e6 * pow(sh, 4) +
3321 1.0340232e7 * pow(sh, 5) - 3.1138149e7 * pow(sh, 6) +
3322 5.1560286e7 * pow(sh, 7) - 5.0820204e7 * pow(sh, 8) +
3323 2.9901036e7 * pow(sh, 9) - 9.847894e6 * pow(sh, 10) +
3324 1.626412e6 * pow(sh, 11) - 564648. * pow(sh, 12) +
3325 1.515872e6 * pow(sh, 13) - 1.613573e6 * pow(sh, 14) -
3326 1.0017483e7 * pow(sh, 15) + 5.154742e7 * pow(sh, 16) -
3327 1.211917e8 * pow(sh, 17) + 1.785651e8 * pow(sh, 18) -
3328 1.800631e8 * pow(sh, 19) + 1.286215e8 * pow(sh, 20) -
3329 6.57363e7 * pow(sh, 21) + 2.38644e7 * pow(sh, 22) -
3330 5.9997e6 * pow(sh, 23) + 991200. * pow(sh, 24) -
3331 96600. * pow(sh, 25) + 4200. * pow(sh, 26)))) * pow(z, 6);
3333 res += pow(z, 6)*(-13382.609103673856 -
3334 83909.57333198964 * (1. - sh) -
3335 285314.7964864959 * pow(1. - sh, 2) -
3336 733141.6984787509 * pow(1. - sh, 3) +
3337 log(muh)*(-1097.4814814814815 -
3338 7580.444444444444 * (1. - sh) -
3339 29518.222222222223 * pow(1. - sh, 2) -
3340 85684.14814814815 * pow(1. - sh, 3) -
3341 206727.11111111112 * pow(1. - sh, 4)) -
3342 1.6179936283094224e6 * pow(1. - sh, 4) +
3343 log(z)*(1155.1126760626482 +
3344 4506.419624502263 * (1. - sh) +
3345 25308.484542876675 * pow(1. - sh, 2) +
3346 94030.12199742932 * pow(1. - sh, 3) +
3347 log(muh)*(-995.5555555555555 -
3348 5973.333333333333 * (1. - sh) -
3349 20906.666666666668 * pow(1. - sh, 2) -
3350 55751.11111111111 * pow(1. - sh, 3) -
3351 125440. * pow(1. - sh, 4)) +
3352 250854.98047286444 * pow(1. - sh, 4)) +
3353 (1201.679012345679 + 7214.074074074074 * (1. - sh) +
3354 26538.192592592593 * pow(1. - sh, 2) +
3355 71813.39259259259 * pow(1. - sh, 3) +
3356 162002.66666666666 * pow(1. - sh, 4)) * pow(log(z), 2));
3361 res += 0.00014351529503044654 * pow(3.141592653589793, 2) * pow(sh, -11)*
3362 (693. + 9450. * sh + 60305. * pow(sh, 2) +
3363 238200. * pow(sh, 3) + 642234. * pow(sh, 4) +
3364 1.173308e6 * pow(sh, 5) + 642234. * pow(sh, 6) +
3365 238200. * pow(sh, 7) + 60305. * pow(sh, 8) +
3366 9450. * pow(sh, 9) + 693. * pow(sh, 10)) * pow(z, 6.5);
3368 res += (4355.652485049105 + 26133.914910294636 * (1. - sh) +
3369 94735.44154981805 * pow(1. - sh, 2) +
3370 266783.7147092577 * pow(1. - sh, 3) +
3371 641447.6079321424 * pow(1. - sh, 4)) * pow(z, 6.5);
3376 res += (-1. * lmu * (3527.9238095238097 * pow(sh, -7) -
3377 3072. * Lsh * pow(sh, -7)) +
3378 188.0551244054121 * pow(-1 + sh, -3) * pow(sh, -5) -
3379 156.44444444444446 * Li3sh * pow(-1 + sh, -3) * pow(sh, -5) +
3380 52.148148148148145 * pow(Lsh, 3) * pow(-1 + sh, -3) *
3381 pow(sh, -5) + 4.266666666666667 * cl2 * sqrt4sh * sqrtsh *
3382 pow(sh, 2)*(-204. + 654. * sh + 3621. * pow(sh, 2) -
3383 9390. * pow(sh, 3) + 8928. * pow(sh, 4) -
3384 4384. * pow(sh, 5) + 1190. * pow(sh, 6) -
3385 170. * pow(sh, 7) + 10. * pow(sh, 8)) +
3386 0.06772486772486773 * ash * sqrt4sh * sqrtsh * pow(-1 + sh, -3) *
3387 pow(sh, 2)*(23244. - 145995. * sh - 135165. * pow(sh, 2) +
3388 2.214474e6 * pow(sh, 3) - 5.812671e6 * pow(sh, 4) +
3389 7.761562e6 * pow(sh, 5) - 6.239858e6 * pow(sh, 6) +
3390 3.192422e6 * pow(sh, 7) - 1.048133e6 * pow(sh, 8) +
3391 214001. * pow(sh, 9) - 24740. * pow(sh, 10) +
3392 1237. * pow(sh, 11)) +
3393 0.7111111111111111 * Li2sh * pow(-1 + sh, -8) * pow(sh, -7)*
3394 (-1080. + 8600. * sh - 29988. * pow(sh, 2) +
3395 59752. * pow(sh, 3) - 74250. * pow(sh, 4) +
3396 58630. * pow(sh, 5) - 28380. * pow(sh, 6) +
3397 7472. * pow(sh, 7) - 698. * pow(sh, 8) - 10. * pow(sh, 11) +
3398 pow(sh, 12) + pow(sh, 13)) -
3399 12.8 * pow(ash, 2) * pow(-1 + sh, -4) * pow(sh, 2)*
3400 (60. - 1390. * sh + 3920. * pow(sh, 2) +
3401 11563. * pow(sh, 3) - 80966. * pow(sh, 4) +
3402 191832. * pow(sh, 5) - 258720. * pow(sh, 6) +
3403 224115. * pow(sh, 7) - 130416. * pow(sh, 8) +
3404 51582. * pow(sh, 9) - 13684. * pow(sh, 10) +
3405 2330. * pow(sh, 11) - 230. * pow(sh, 12) + 10. * pow(sh, 13))
3406 - Lsh * (-156.44444444444446 * Li2sh * pow(-1 + sh, -3) *
3407 pow(sh, -5) - 4.266666666666667 * ash * sqrt4sh * sqrtsh *
3408 pow(sh, 2)*(-204. + 654. * sh + 3621. * pow(sh, 2) -
3409 9390. * pow(sh, 3) + 8928. * pow(sh, 4) -
3410 4384. * pow(sh, 5) + 1190. * pow(sh, 6) -
3411 170. * pow(sh, 7) + 10. * pow(sh, 8)) +
3412 0.000013437473754934073 * pow(-1 + sh, -9) * pow(sh, -12)*
3413 (3900. - 3800. * sh - 116136. * pow(sh, 2) -
3414 213534. * pow(sh, 3) + 423060. * pow(sh, 4) -
3415 1.71541326e8 * pow(sh, 5) +
3416 1.547854104e9 * pow(sh, 6) -
3417 6.099085398e9 * pow(sh, 7) +
3418 9.702e6 * pow(3.141592653589793, 2) * pow(-1 + sh, 6) *
3419 pow(sh, 7) + 1.3863208764e10 * pow(sh, 8) -
3420 2.0080171644e10 * pow(sh, 9) +
3421 1.9223278801e10 * pow(sh, 10) -
3422 1.2149332644e10 * pow(sh, 11) +
3423 4.878367842e9 * pow(sh, 12) -
3424 1.128611791e9 * pow(sh, 13) -
3425 52920. * Lshb * (-1. + sh) * pow(sh, 5)*
3426 (-1080. + 8600. * sh - 29988. * pow(sh, 2) +
3427 59752. * pow(sh, 3) - 74250. * pow(sh, 4) +
3428 58630. * pow(sh, 5) - 28380. * pow(sh, 6) +
3429 7472. * pow(sh, 7) - 698. * pow(sh, 8) -
3430 10. * pow(sh, 11) + pow(sh, 12) + pow(sh, 13)) +
3431 1.16995017e8 * pow(sh, 14) - 1.22745e6 * pow(sh, 15) -
3432 53410. * pow(sh, 16) + 252693. * pow(sh, 17) -
3433 23814. * pow(sh, 18) - 3234. * pow(sh, 19))) -
3434 0.0028218694885361554 * pow(Lsh, 2) * pow(-1 + sh, -9) *
3435 pow(sh, -12)*(60. + 340. * sh + 312. * pow(sh, 2) -
3436 1872. * pow(sh, 3) - 4980. * pow(sh, 4) +
3437 1.468248e6 * pow(sh, 5) - 1.3307922e7 * pow(sh, 6) +
3438 5.3745563e7 * pow(sh, 7) -
3439 27720. * Lshb * pow(-1 + sh, 6) * pow(sh, 7) -
3440 1.2604899e8 * pow(sh, 8) + 1.89165441e8 * pow(sh, 9) -
3441 1.88411563e8 * pow(sh, 10) + 1.24499847e8 * pow(sh, 11) -
3442 5.2564902e7 * pow(sh, 12) + 1.2834946e7 * pow(sh, 13) -
3443 1.375368e6 * pow(sh, 14) + 1260. * pow(sh, 16) -
3444 1386. * pow(sh, 17) + 126. * pow(sh, 19)) +
3445 0.0028218694885361554 * pow(lz, 2) * pow(sh, -12)*
3446 (15. + 220. * sh + 1518. * pow(sh, 2) + 6534. * pow(sh, 3) +
3447 19503. * pow(sh, 4) + 1.511448e6 * pow(sh, 5) -
3448 20790. * pow(sh, 6) - 6546. * pow(sh, 7) -
3449 24741. * pow(sh, 8) - 29490. * pow(sh, 9) -
3450 19183. * pow(sh, 10) - 15. * pow(sh, 11) -
3451 42336. * pow(sh, 12) - 50652. * pow(sh, 13) -
3452 42336. * pow(sh, 14) - 522396. * pow(sh, 15) -
3453 514080. * pow(sh, 16) + 5.330178e6 * pow(sh, 17) -
3454 8.247204e6 * pow(sh, 18) + 6.008688e6 * pow(sh, 19) -
3455 2.443392e6 * pow(sh, 20) + 570780. * pow(sh, 21) -
3456 71820. * pow(sh, 22) + 3780. * pow(sh, 23)) -
3457 lz * (3072. * lmu * pow(sh, -7) +
3458 78.22222222222223 * pow(Lsh, 2) * pow(-1 + sh, -3) *
3459 pow(sh, -5) - 8.533333333333333 * ash * sqrt4sh * sqrtsh *
3460 pow(sh, 2)*(-204. + 654. * sh + 3621. * pow(sh, 2) -
3461 9390. * pow(sh, 3) + 8928. * pow(sh, 4) -
3462 4384. * pow(sh, 5) + 1190. * pow(sh, 6) -
3463 170. * pow(sh, 7) + 10. * pow(sh, 8)) +
3464 0.005643738977072311 * Lsh * pow(-1 + sh, -8) *
3465 pow(sh, -12)*(30. + 200. * sh + 356. * pow(sh, 2) -
3466 580. * pow(sh, 3) - 3070. * pow(sh, 4) +
3467 1.465988e6 * pow(sh, 5) - 1.1840494e7 * pow(sh, 6) +
3468 4.1785315e7 * pow(sh, 7) - 8.3994971e7 * pow(sh, 8) +
3469 1.0527886e8 * pow(sh, 9) - 8.4349693e7 * pow(sh, 10) +
3470 4.2263834e7 * pow(sh, 11) -
3471 1.2154288e7 * pow(sh, 12) + 1.55972e6 * pow(sh, 13) -
3472 2380. * pow(sh, 14) - 2065. * pow(sh, 15) -
3473 945. * pow(sh, 16) + 462. * pow(sh, 17) +
3474 21. * pow(sh, 18)) -
3475 4.071961743919416e-7 * pow(-1 + sh, -7) * pow(sh, -12)*
3476 (64350. + 66000. * sh - 1.848594e6 * pow(sh, 2) -
3477 7.286499e6 * pow(sh, 3) - 5.743914e6 * pow(sh, 4) -
3478 5.685999132e9 * pow(sh, 5) +
3479 4.1685525882e10 * pow(sh, 6) -
3480 1.30007412909e11 * pow(sh, 7) +
3481 2.26326493503e11 * pow(sh, 8) -
3482 2.39590844427e11 * pow(sh, 9) +
3483 1.56064921276e11 * pow(sh, 10) -
3484 5.9156118388e10 * pow(sh, 11) +
3485 1.1660658501e10 * pow(sh, 12) -
3486 5.482108499e9 * pow(sh, 13) +
3487 1.2697937753e10 * pow(sh, 14) -
3488 1.820498844e10 * pow(sh, 15) +
3489 4.184087512e9 * pow(sh, 16) +
3490 6.5568484406e10 * pow(sh, 17) -
3491 9.5193663855e10 * pow(sh, 18) -
3492 3.83610629052e11 * pow(sh, 19) +
3493 1.851933873636e12 * pow(sh, 20) -
3494 3.850977196836e12 * pow(sh, 21) +
3495 4.997388254166e12 * pow(sh, 22) -
3496 4.45400972478e12 * pow(sh, 23) +
3497 2.82695858676e12 * pow(sh, 24) -
3498 1.2912095196e12 * pow(sh, 25) +
3499 4.2141221892e11 * pow(sh, 26) -
3500 9.582123474e10 * pow(sh, 27) +
3501 1.44039357e10 * pow(sh, 28) -
3502 1.28490516e9 * pow(sh, 29) + 5.143446e7 * pow(sh, 30)))
3503 + 2.9379233361611947e-11 * pow(-1 + sh, -8) * pow(sh, -12)*
3504 (2.074545e8 - 4.76982e8 * sh -
3505 3.223561968e9 * pow(sh, 2) +
3506 1.2229671465e10 * pow(sh, 3) +
3507 5.5118379195e10 * pow(sh, 4) -
3508 1.50774673843134e14 * pow(sh, 5) +
3509 1.186047300987252e15 * pow(sh, 6) -
3510 4.079182881529008e15 * pow(sh, 7) +
3511 8.005291590723075e15 * pow(sh, 8) -
3512 9.81400671207825e15 * pow(sh, 9) +
3513 7.724627706932279e15 * pow(sh, 10) -
3514 3.853787364524662e15 * pow(sh, 11) +
3515 1.164725088632579e15 * pow(sh, 12) -
3516 2.634808520857e14 * pow(sh, 13) +
3517 2.5170701729432e14 * pow(sh, 14) -
3518 4.0013657985217e14 * pow(sh, 15) +
3519 1.47801666642455e14 * pow(sh, 16) +
3520 1.2877818678578e15 * pow(sh, 17) -
3521 2.076397177767913e15 * pow(sh, 18) -
3522 8.385630050032431e15 * pow(sh, 19) +
3523 4.538622316201646e16 * pow(sh, 20) -
3524 1.061838620189715e17 * pow(sh, 21) +
3525 1.575009771316346e17 * pow(sh, 22) -
3526 1.6342684496066314e17 * pow(sh, 23) +
3527 1.2333678370646925e17 * pow(sh, 24) -
3528 6.8698738259084936e16 * pow(sh, 25) +
3529 2.8233898574272776e16 * pow(sh, 26) -
3530 8.447633225048556e15 * pow(sh, 27) +
3531 1.786675943700648e15 * pow(sh, 28) -
3532 2.5272282124782e14 * pow(sh, 29) +
3533 2.141343983394e13 * pow(sh, 30) -
3534 8.2047522942e11 * pow(sh, 31) +
3535 3.20166e7 * pow(3.141592653589793, 2)*
3536 (-90. - 600. * sh - 1068. * pow(sh, 2) +
3537 1740. * pow(sh, 3) + 9210. * pow(sh, 4) -
3538 3.558804e6 * pow(sh, 5) + 2.8685982e7 * pow(sh, 6) -
3539 1.00848529e8 * pow(sh, 7) +
3540 2.01465905e8 * pow(sh, 8) - 2.5035526e8 * pow(sh, 9) +
3541 1.98432319e8 * pow(sh, 10) -
3542 9.8218082e7 * pow(sh, 11) +
3543 2.7878488e7 * pow(sh, 12) - 3.151676e6 * pow(sh, 13) -
3544 1.140968e6 * pow(sh, 14) + 600307. * pow(sh, 15) +
3545 5.766227e6 * pow(sh, 16) - 8.667286e6 * pow(sh, 17) -
3546 7.4435635e7 * pow(sh, 18) +
3547 4.29448764e8 * pow(sh, 19) -
3548 1.174957056e9 * pow(sh, 20) +
3549 2.075532984e9 * pow(sh, 21) -
3550 2.59118244e9 * pow(sh, 22) +
3551 2.377353888e9 * pow(sh, 23) -
3552 1.631212128e9 * pow(sh, 24) +
3553 8.40926772e8 * pow(sh, 25) -
3554 3.23868888e8 * pow(sh, 26) +
3555 9.1648368e7 * pow(sh, 27) -
3556 1.8464544e7 * pow(sh, 28) + 2.50236e6 * pow(sh, 29) -
3557 204120. * pow(sh, 30) + 7560. * pow(sh, 31)))) * pow(z, 7);
3559 res += pow(z, 7)*(-44491.606444877616 -
3560 330794.4598052619 * (1. - sh) -
3561 1.2804206866827896e6 * pow(1. - sh, 2) -
3562 3.6984879243897516e6 * pow(1. - sh, 3) +
3563 log(muh)*(-3527.9238095238097 -
3564 27767.466666666667 * (1. - sh) -
3565 121821.86666666667 * pow(1. - sh, 2) -
3566 394137.6 * pow(1. - sh, 3) -
3567 1.049856e6 * pow(1. - sh, 4)) -
3568 9.101629224853888e6 * pow(1. - sh, 4) +
3569 log(z)*(4793.772650900163 +
3570 20076.271464118738 * (1. - sh) +
3571 125799.80299624978 * pow(1. - sh, 2) +
3572 506715.246202913 * pow(1. - sh, 3) +
3573 log(muh)*(-3072. - 21504. * (1. - sh) -
3574 86016. * pow(1. - sh, 2) -
3575 258048. * pow(1. - sh, 3) -
3576 645120. * pow(1. - sh, 4)) +
3577 1.4556579496221067e6 * pow(1. - sh, 4)) +
3578 (4000.5164021164023 +
3579 28350.044444444444 * (1. - sh) +
3580 120184.65185185186 * pow(1. - sh, 2) +
3581 364920.47407407407 * pow(1. - sh, 3) +
3582 913431.7037037037 * pow(1. - sh, 4)) * pow(log(z), 2));
3587 res += 0.000024287203774383262 * pow(3.141592653589793, 2) *
3588 pow(sh, -13)*(3003. + 47124. * sh + 349398. * pow(sh, 2) +
3589 1.62554e6 * pow(sh, 3) + 5.297685e6 * pow(sh, 4) +
3590 1.265412e7 * pow(sh, 5) + 2.15477e7 * pow(sh, 6) +
3591 1.265412e7 * pow(sh, 7) + 5.297685e6 * pow(sh, 8) +
3592 1.62554e6 * pow(sh, 9) + 349398. * pow(sh, 10) +
3593 47124. * pow(sh, 11) + 3003. * pow(sh, 12)) * pow(z, 7.5);
3595 res += (14742.208410935435 + 103195.45887654803 * (1. - sh) +
3596 426049.8230760341 * pow(1. - sh, 2) +
3597 1.3444894070773118e6 * pow(1. - sh, 3) +
3598 3.579434527547286e6 * pow(1. - sh, 4)) * pow(z, 7.5);
3603 res += (-1. * lmu * (11644.08888888889 * pow(sh, -8) -
3604 9856. * Lsh * pow(sh, -8)) +
3605 523.8678465579336 * pow(-1 + sh, -3) * pow(sh, -6) -
3606 435.8095238095238 * Li3sh * pow(-1 + sh, -3) * pow(sh, -6) +
3607 145.26984126984127 * pow(Lsh, 3) * pow(-1 + sh, -3) *
3608 pow(sh, -6) - 7.822222222222222 * cl2 * sqrt4sh * sqrtsh *
3609 pow(sh, 2)*(-42. + 1287. * sh - 1614. * pow(sh, 2) -
3610 19017. * pow(sh, 3) + 50556. * pow(sh, 4) -
3611 55330. * pow(sh, 5) + 33315. * pow(sh, 6) -
3612 11945. * pow(sh, 7) + 2550. * pow(sh, 8) -
3613 300. * pow(sh, 9) + 15. * pow(sh, 10)) -
3614 0.0028218694885361554 * ash * sqrt4sh * sqrtsh *
3615 pow(-1 + sh, -3) * pow(sh, 2)*
3616 (228522. - 7.741593e6 * sh + 3.0723729e7 * pow(sh, 2) +
3617 6.0003996e7 * pow(sh, 3) - 5.8092981e8 * pow(sh, 4) +
3618 1.504380659e9 * pow(sh, 5) - 2.124842634e9 * pow(sh, 6) +
3619 1.888618168e9 * pow(sh, 7) - 1.114436468e9 * pow(sh, 8) +
3620 4.4384322e8 * pow(sh, 9) - 1.18087858e8 * pow(sh, 10) +
3621 2.0123511e7 * pow(sh, 11) - 1.986441e6 * pow(sh, 12) +
3622 86367. * pow(sh, 13)) +
3623 23.466666666666665 * pow(ash, 2) * pow(-1 + sh, -4) *
3624 pow(sh, 3)*(-700. + 7875. * sh - 14966. * pow(sh, 2) -
3625 79023. * pow(sh, 3) + 455850. * pow(sh, 4) -
3626 1.077135e6 * pow(sh, 5) + 1.523775e6 * pow(sh, 6) -
3627 1.432728e6 * pow(sh, 7) + 934836. * pow(sh, 8) -
3628 430150. * pow(sh, 9) + 139230. * pow(sh, 10) -
3629 31025. * pow(sh, 11) + 4530. * pow(sh, 12) -
3630 390. * pow(sh, 13) + 15. * pow(sh, 14)) +
3631 0.016931216931216932 * Li2sh * pow(-1 + sh, -9) * pow(sh, -8)*
3632 (145530. - 1.30482e6 * sh + 5.20032e6 * pow(sh, 2) -
3633 1.208298e7 * pow(sh, 3) + 1.801254e7 * pow(sh, 4) -
3634 1.78178e7 * pow(sh, 5) + 1.163162e7 * pow(sh, 6) -
3635 4.77477e6 * pow(sh, 7) + 1.082136e6 * pow(sh, 8) -
3636 94842. * pow(sh, 9) - 420. * pow(sh, 13) +
3637 70. * pow(sh, 14) + 35. * pow(sh, 15)) -
3638 Lsh * (-435.8095238095238 * Li2sh * pow(-1 + sh, -3) *
3639 pow(sh, -6) + 7.822222222222222 * ash * sqrt4sh * sqrtsh *
3640 pow(sh, 2)*(-42. + 1287. * sh - 1614. * pow(sh, 2) -
3641 19017. * pow(sh, 3) + 50556. * pow(sh, 4) -
3642 55330. * pow(sh, 5) + 33315. * pow(sh, 6) -
3643 11945. * pow(sh, 7) + 2550. * pow(sh, 8) -
3644 300. * pow(sh, 9) + 15. * pow(sh, 10)) +
3645 4.799097769619311e-7 * pow(-1 + sh, -10) * pow(sh, -14)*
3646 (-70875. + 64125. * sh + 2.661225e6 * pow(sh, 2) +
3647 6.471633e6 * pow(sh, 3) - 8.447481e6 * pow(sh, 4) -
3648 4.7108817e7 * pow(sh, 5) +
3649 2.1458196025e10 * pow(sh, 6) -
3650 2.14166617315e11 * pow(sh, 7) +
3651 9.57541423647e11 * pow(sh, 8) +
3652 7.56756e8 * pow(3.141592653589793, 2) *
3653 pow(-1 + sh, 7) * pow(sh, 8) -
3654 2.521983619869e12 * pow(sh, 9) +
3655 4.337037789987e12 * pow(sh, 10) -
3656 5.090519317485e12 * pow(sh, 11) +
3657 4.129949513865e12 * pow(sh, 12) -
3658 2.286556939359e12 * pow(sh, 13) +
3659 8.27081528763e11 * pow(sh, 14) -
3660 1.76924790709e11 * pow(sh, 15) -
3661 35280. * Lshb * (-1. + sh) * pow(sh, 6)*
3662 (145530. - 1.30482e6 * sh + 5.20032e6 * pow(sh, 2) -
3663 1.208298e7 * pow(sh, 3) + 1.801254e7 * pow(sh, 4) -
3664 1.78178e7 * pow(sh, 5) + 1.163162e7 * pow(sh, 6) -
3665 4.77477e6 * pow(sh, 7) + 1.082136e6 * pow(sh, 8) -
3666 94842. * pow(sh, 9) - 420. * pow(sh, 13) +
3667 70. * pow(sh, 14) + 35. * pow(sh, 15)) +
3668 1.717636396e10 * pow(sh, 16) -
3669 1.22598e7 * pow(sh, 17) - 4.2757596e7 * pow(sh, 18) +
3670 2.333772e6 * pow(sh, 19) + 6.325704e6 * pow(sh, 20) -
3671 688800. * pow(sh, 21) - 54600. * pow(sh, 22))) -
3672 0.00040312421264802217 * pow(Lsh, 2) * pow(-1 + sh, -10) *
3673 pow(sh, -14)*(-315. - 2115. * sh - 3255. * pow(sh, 2) +
3674 9681. * pow(sh, 3) + 37863. * pow(sh, 4) +
3675 21231. * pow(sh, 5) - 3.5869645e7 * pow(sh, 6) +
3676 3.5987455e8 * pow(sh, 7) - 1.631139771e9 * pow(sh, 8) -
3677 540540. * Lshb * pow(-1 + sh, 7) * pow(sh, 8) +
3678 4.372215162e9 * pow(sh, 9) -
3679 7.668316866e9 * pow(sh, 10) +
3680 9.19386972e9 * pow(sh, 11) - 7.62973092e9 * pow(sh, 12) +
3681 4.325428674e9 * pow(sh, 13) -
3682 1.601657043e9 * pow(sh, 14) +
3683 3.49216864e8 * pow(sh, 15) - 3.3959485e7 * pow(sh, 16) +
3684 8820. * pow(sh, 19) - 10290. * pow(sh, 20) +
3685 735. * pow(sh, 21) + 735. * pow(sh, 22)) -
3686 0.00010078105316200554 * pow(lz, 2) * pow(sh, -14)*
3687 (-315. - 5265. * sh - 41730. * pow(sh, 2) -
3688 208494. * pow(sh, 3) - 734877. * pow(sh, 4) -
3689 1.926639e6 * pow(sh, 5) - 1.4698138e8 * pow(sh, 6) +
3690 2.16216e6 * pow(sh, 7) + 48987. * pow(sh, 8) +
3691 1.669041e6 * pow(sh, 9) + 2.629962e6 * pow(sh, 10) +
3692 2.41431e6 * pow(sh, 11) + 1.436625e6 * pow(sh, 12) +
3693 315. * pow(sh, 13) + 3.764376e6 * pow(sh, 14) +
3694 4.230072e6 * pow(sh, 15) + 5.006232e6 * pow(sh, 16) -
3695 7.37352e6 * pow(sh, 17) + 1.0643094e8 * pow(sh, 18) +
3696 1.95359472e8 * pow(sh, 19) -
3697 1.425010356e9 * pow(sh, 20) +
3698 2.406212424e9 * pow(sh, 21) -
3699 2.05080876e9 * pow(sh, 22) + 1.03141962e9 * pow(sh, 23) -
3700 3.2026302e8 * pow(sh, 24) + 6.054048e7 * pow(sh, 25) -
3701 6.40332e6 * pow(sh, 26) + 291060. * pow(sh, 27)) -
3702 lz * (9856. * lmu * pow(sh, -8) +
3703 217.9047619047619 * pow(Lsh, 2) * pow(-1 + sh, -3) *
3704 pow(sh, -6) + 15.644444444444444 * ash * sqrt4sh * sqrtsh *
3705 pow(sh, 2)*(-42. + 1287. * sh - 1614. * pow(sh, 2) -
3706 19017. * pow(sh, 3) + 50556. * pow(sh, 4) -
3707 55330. * pow(sh, 5) + 33315. * pow(sh, 6) -
3708 11945. * pow(sh, 7) + 2550. * pow(sh, 8) -
3709 300. * pow(sh, 9) + 15. * pow(sh, 10)) +
3710 0.00040312421264802217 * Lsh * pow(-1 + sh, -9) *
3711 pow(sh, -14)*(-315. - 2430. * sh - 5685. * pow(sh, 2) +
3712 3996. * pow(sh, 3) + 41859. * pow(sh, 4) +
3713 63090. * pow(sh, 5) - 7.1605955e7 * pow(sh, 6) +
3714 6.4801935e8 * pow(sh, 7) -
3715 2.609788929e9 * pow(sh, 8) +
3716 6.121801554e9 * pow(sh, 9) -
3717 9.217982115e9 * pow(sh, 10) +
3718 9.24507918e9 * pow(sh, 11) -
3719 6.181354605e9 * pow(sh, 12) +
3720 2.66089059e9 * pow(sh, 13) -
3721 6.71469273e8 * pow(sh, 14) +
3722 7.6209266e7 * pow(sh, 15) + 37800. * pow(sh, 16) -
3723 47460. * pow(sh, 17) - 29232. * pow(sh, 18) -
3724 9828. * pow(sh, 19) + 6930. * pow(sh, 20) +
3725 210. * pow(sh, 21)) -
3726 1.6780062131536054e-9 * pow(-1 + sh, -8) * pow(sh, -14)*
3727 (-1.0135125e7 - 1.1100375e7 * sh +
3728 3.6848955e8 * pow(sh, 2) +
3729 1.673522994e9 * pow(sh, 3) +
3730 1.770566655e9 * pow(sh, 4) -
3731 4.868950515e9 * pow(sh, 5) +
3732 6.12638023426e12 * pow(sh, 6) -
3733 5.053531208726e13 * pow(sh, 7) +
3734 1.81942427597788e14 * pow(sh, 8) -
3735 3.7509726864038e14 * pow(sh, 9) +
3736 4.863909914149e14 * pow(sh, 10) -
3737 4.0837905886362e14 * pow(sh, 11) +
3738 2.185978661919e14 * pow(sh, 12) -
3739 6.9223437871788e13 * pow(sh, 13) +
3740 9.478097027448e12 * pow(sh, 14) +
3741 4.822522515092e12 * pow(sh, 15) -
3742 1.3702469432263e13 * pow(sh, 16) +
3743 2.0211432885167e13 * pow(sh, 17) -
3744 1.5340657706842e13 * pow(sh, 18) +
3745 3.8517131981282e13 * pow(sh, 19) -
3746 1.98231128437323e14 * pow(sh, 20) +
3747 2.48611036890867e14 * pow(sh, 21) +
3748 1.312843416208632e15 * pow(sh, 22) -
3749 6.858549767677608e15 * pow(sh, 23) +
3750 1.6510967898168744e16 * pow(sh, 24) -
3751 2.57336419311972e16 * pow(sh, 25) +
3752 2.851348331813982e16 * pow(sh, 26) -
3753 2.336070503975208e16 * pow(sh, 27) +
3754 1.439716999494774e16 * pow(sh, 28) -
3755 6.70336841266092e15 * pow(sh, 29) +
3756 2.34397396677444e15 * pow(sh, 30) -
3757 6.0527470649646e14 * pow(sh, 31) +
3758 1.1182510300962e14 * pow(sh, 32) -
3759 1.396347000048e13 * pow(sh, 33) +
3760 1.054350297e12 * pow(sh, 34) -
3761 3.631041414e10 * pow(sh, 35))) -
3762 4.656471897973153e-15 * pow(-1 + sh, -9) * pow(sh, -14)*
3763 (7.34806697625e11 - 1.86465740175e12 * sh -
3764 1.3539275231625e13 * pow(sh, 2) +
3765 6.7027116030732e13 * pow(sh, 3) +
3766 3.63276406618983e14 * pow(sh, 4) +
3767 7.707518005473e13 * pow(sh, 5) -
3768 2.670939112625576e18 * pow(sh, 6) +
3769 2.358705949072117e19 * pow(sh, 7) -
3770 9.25798205589377e19 * pow(sh, 8) +
3771 2.1183909722959626e20 * pow(sh, 9) -
3772 3.1161735718840415e20 * pow(sh, 10) +
3773 3.06304944728649e20 * pow(sh, 11) -
3774 2.023336897561395e20 * pow(sh, 12) +
3775 8.778312486562154e19 * pow(sh, 13) -
3776 2.3206885685742117e19 * pow(sh, 14) +
3777 1.3935768215051139e18 * pow(sh, 15) +
3778 5.883932779473471e18 * pow(sh, 16) -
3779 1.0545119614382879e19 * pow(sh, 17) +
3780 1.1056449646096937e19 * pow(sh, 18) -
3781 2.3434959475666616e19 * pow(sh, 19) +
3782 1.0419598368525897e20 * pow(sh, 20) -
3783 1.688591180885733e20 * pow(sh, 21) -
3784 5.308738122339353e20 * pow(sh, 22) +
3785 3.5306726064196524e21 * pow(sh, 23) -
3786 9.689077766441178e21 * pow(sh, 24) +
3787 1.7097533811227756e22 * pow(sh, 25) -
3788 2.1588630869532427e22 * pow(sh, 26) +
3789 2.038422378040867e22 * pow(sh, 27) -
3790 1.4690695476739364e22 * pow(sh, 28) +
3791 8.14476034656345e21 * pow(sh, 29) -
3792 3.4697517653715583e21 * pow(sh, 30) +
3793 1.1250567404926722e21 * pow(sh, 31) -
3794 2.7233981875541377e20 * pow(sh, 32) +
3795 4.759269487750347e19 * pow(sh, 33) -
3796 5.663805036845516e18 * pow(sh, 34) +
3797 4.101841359381566e17 * pow(sh, 35) -
3798 1.362249094884678e16 * pow(sh, 36) +
3799 4.32864432e10 * pow(3.141592653589793, 2)*
3800 (-315. - 2430. * sh - 5685. * pow(sh, 2) +
3801 3996. * pow(sh, 3) + 41859. * pow(sh, 4) +
3802 63090. * pow(sh, 5) - 5.9090375e7 * pow(sh, 6) +
3803 5.3345259e8 * pow(sh, 7) -
3804 2.141709439e9 * pow(sh, 8) +
3805 5.001401214e9 * pow(sh, 9) -
3806 7.486485645e9 * pow(sh, 10) +
3807 7.45371022e9 * pow(sh, 11) -
3808 4.941375175e9 * pow(sh, 12) +
3809 2.10797493e9 * pow(sh, 13) -
3810 5.28751827e8 * pow(sh, 14) +
3811 6.6874034e7 * pow(sh, 15) - 2.302929e7 * pow(sh, 16) +
3812 4.790626e7 * pow(sh, 17) -
3813 1.24781262e8 * pow(sh, 18) +
3814 4.35359256e8 * pow(sh, 19) -
3815 5.45760408e8 * pow(sh, 20) -
3816 3.99911175e9 * pow(sh, 21) +
3817 2.559251793e10 * pow(sh, 22) -
3818 7.87509723e10 * pow(sh, 23) +
3819 1.60158041736e11 * pow(sh, 24) -
3820 2.35861111332e11 * pow(sh, 25) +
3821 2.6202440418e11 * pow(sh, 26) -
3822 2.2419142284e11 * pow(sh, 27) +
3823 1.491808626e11 * pow(sh, 28) -
3824 7.7359627884e10 * pow(sh, 29) +
3825 3.1118686368e10 * pow(sh, 30) -
3826 9.59731542e9 * pow(sh, 31) +
3827 2.22259884e9 * pow(sh, 32) -
3828 3.733653e8 * pow(sh, 33) + 4.288284e7 * pow(sh, 34) -
3829 3.00762e6 * pow(sh, 35) + 97020. * pow(sh, 36)))) *
3832 res += pow(z, 8)*(-152175.3347907426 -
3833 1.3121004907612824e6 * (1. - sh) -
3834 5.692547653098881e6 * pow(1. - sh, 2) -
3835 1.8262414685026634e7 * pow(1. - sh, 3) +
3836 log(muh)*(-11644.08888888889 -
3837 103008.71111111112 * (1. - sh) -
3838 502963.2 * pow(1. - sh, 2) -
3839 1.794816e6 * pow(1. - sh, 3) -
3840 5.231424e6 * pow(1. - sh, 4)) -
3841 4.9606666429006e7 * pow(1. - sh, 4) +
3842 log(z)*(19181.53823244784 +
3843 86972.5498690609 * (1. - sh) +
3844 607120.085934602 * pow(1. - sh, 2) +
3845 2.6398902680343925e6 * pow(1. - sh, 3) +
3846 log(muh)*(-9856. - 78848. * (1. - sh) -
3847 354816. * pow(1. - sh, 2) -
3848 1.18272e6 * pow(1. - sh, 3) -
3849 3.25248e6 * pow(1. - sh, 4)) +
3850 8.124517405933031e6 * pow(1. - sh, 4)) +
3851 (13720.493827160493 +
3852 112232.75061728395 * (1. - sh) +
3853 538923.3777777777 * pow(1. - sh, 2) +
3854 1.8125184592592593e6 * pow(1. - sh, 3) +
3855 4.984340004232804e6 * pow(1. - sh, 4)) * pow(log(z), 2)
3861 res += 8.66719820968187e-6 * pow(3.141592653589793, 2) * pow(sh, -15)*
3862 (6435. + 114114. * sh + 962577. * pow(sh, 2) +
3863 5.137524e6 * pow(sh, 3) + 1.9445755e7 * pow(sh, 4) +
3864 5.526435e7 * pow(sh, 5) + 1.20557265e8 * pow(sh, 6) +
3865 1.9446652e8 * pow(sh, 7) + 1.20557265e8 * pow(sh, 8) +
3866 5.526435e7 * pow(sh, 9) + 1.9445755e7 * pow(sh, 10) +
3867 5.137524e6 * pow(sh, 11) + 962577. * pow(sh, 12) +
3868 114114. * pow(sh, 13) + 6435. * pow(sh, 14)) * pow(z, 8.5);
3870 res += (51106.322491242834 + 408850.5799299427 * (1. - sh) +
3871 1.8934892483005472e6 * pow(1. - sh, 2) +
3872 6.615713446491386e6 * pow(1. - sh, 3) +
3873 1.9308607466222685e7 * pow(1. - sh, 4)) * pow(z, 8.5);
3878 res += (-1. * lmu * (39253.5139329806 * pow(sh, -9) -
3879 32540.444444444445 * Lsh * pow(sh, -9)) +
3880 1527.947885793973 * pow(-1 + sh, -3) * pow(sh, -7) -
3881 1271.111111111111 * Li3sh * pow(-1 + sh, -3) * pow(sh, -7) +
3882 423.7037037037037 * pow(Lsh, 3) * pow(-1 + sh, -3) *
3883 pow(sh, -7) + 16.141093474426807 * cl2 * sqrt4sh * sqrtsh *
3884 pow(sh, 3)*(552. - 6504. * sh + 954. * pow(sh, 2) +
3885 96534. * pow(sh, 3) - 256765. * pow(sh, 4) +
3886 311884. * pow(sh, 5) - 219753. * pow(sh, 6) +
3887 97149. * pow(sh, 7) - 27453. * pow(sh, 8) +
3888 4830. * pow(sh, 9) - 483. * pow(sh, 10) + 21. * pow(sh, 11))\
3889 + 0.00062708210856359 * ash * sqrt4sh * sqrtsh *
3890 pow(-1 + sh, -3) * pow(sh, 3)*
3891 (-2.9602032e7 + 4.3918554e8 * sh -
3892 1.183665042e9 * pow(sh, 2) - 4.124801388e9 * pow(sh, 3) +
3893 2.9786749451e10 * pow(sh, 4) -
3894 7.6073083277e10 * pow(sh, 5) +
3895 1.12343648214e11 * pow(sh, 6) -
3896 1.08358670795e11 * pow(sh, 7) +
3897 7.1769426953e10 * pow(sh, 8) -
3898 3.3301199019e10 * pow(sh, 9) +
3899 1.0824913791e10 * pow(sh, 10) -
3900 2.416552923e9 * pow(sh, 11) +
3901 3.53029242e8 * pow(sh, 12) - 3.0393246e7 * pow(sh, 13) +
3902 1.168971e6 * pow(sh, 14)) -
3903 48.423280423280424 * pow(ash, 2) * pow(-1 + sh, -4) *
3904 pow(sh, 3)*(-112. + 5446. * sh - 40068. * pow(sh, 2) +
3905 45402. * pow(sh, 3) + 481386. * pow(sh, 4) -
3906 2.416545e6 * pow(sh, 5) + 5.68308e6 * pow(sh, 6) -
3907 8.356089e6 * pow(sh, 7) + 8.410689e6 * pow(sh, 8) -
3908 6.034798e6 * pow(sh, 9) + 3.142636e6 * pow(sh, 10) -
3909 1.191618e6 * pow(sh, 11) + 325686. * pow(sh, 12) -
3910 62481. * pow(sh, 13) + 7980. * pow(sh, 14) -
3911 609. * pow(sh, 15) + 21. * pow(sh, 16)) +
3912 0.02821869488536155 * Li2sh * pow(-1 + sh, -10) * pow(sh, -9)*
3913 (-288288. + 2.873871e6 * sh - 1.2889503e7 * pow(sh, 2) +
3914 3.423222e7 * pow(sh, 3) - 5.956146e7 * pow(sh, 4) +
3915 7.0825482e7 * pow(sh, 5) - 5.8121154e7 * pow(sh, 6) +
3916 3.2324292e7 * pow(sh, 7) - 1.153152e7 * pow(sh, 8) +
3917 2.327611e6 * pow(sh, 9) - 186043. * pow(sh, 10) -
3918 252. * pow(sh, 15) + 54. * pow(sh, 16) + 18. * pow(sh, 17))\
3919 - Lsh * (-1271.111111111111 * Li2sh * pow(-1 + sh, -3) *
3920 pow(sh, -7) - 16.141093474426807 * ash * sqrt4sh * sqrtsh *
3921 pow(sh, 3)*(552. - 6504. * sh + 954. * pow(sh, 2) +
3922 96534. * pow(sh, 3) - 256765. * pow(sh, 4) +
3923 311884. * pow(sh, 5) - 219753. * pow(sh, 6) +
3924 97149. * pow(sh, 7) - 27453. * pow(sh, 8) +
3925 4830. * pow(sh, 9) - 483. * pow(sh, 10) +
3926 21. * pow(sh, 11)) +
3927 9.331578996481995e-7 * pow(-1 + sh, -11) * pow(sh, -16)*
3928 (24990. - 21315. * sh - 1.15923e6 * pow(sh, 2) -
3929 3.54162e6 * pow(sh, 3) + 2.396436e6 * pow(sh, 4) +
3930 2.6605446e7 * pow(sh, 5) + 2.5110204e7 * pow(sh, 6) -
3931 4.6121470412e10 * pow(sh, 7) +
3932 5.06007096796e11 * pow(sh, 8) -
3933 2.5214921261e12 * pow(sh, 9) +
3934 1.135134e9 * pow(3.141592653589793, 2) *
3935 pow(-1 + sh, 8) * pow(sh, 9) +
3936 7.515033973326e12 * pow(sh, 10) -
3937 1.4885226217656e13 * pow(sh, 11) +
3938 2.0578551617862e13 * pow(sh, 12) -
3939 2.0264715219756e13 * pow(sh, 13) +
3940 1.4215757817264e13 * pow(sh, 14) -
3941 6.963287491704e12 * pow(sh, 15) +
3942 2.269595684207e12 * pow(sh, 16) -
3943 4.4385323008e11 * pow(sh, 17) -
3944 30240. * Lshb * (-1. + sh) * pow(sh, 7)*
3945 (-288288. + 2.873871e6 * sh -
3946 1.2889503e7 * pow(sh, 2) +
3947 3.423222e7 * pow(sh, 3) - 5.956146e7 * pow(sh, 4) +
3948 7.0825482e7 * pow(sh, 5) -
3949 5.8121154e7 * pow(sh, 6) +
3950 3.2324292e7 * pow(sh, 7) -
3951 1.153152e7 * pow(sh, 8) + 2.327611e6 * pow(sh, 9) -
3952 186043. * pow(sh, 10) - 252. * pow(sh, 15) +
3953 54. * pow(sh, 16) + 18. * pow(sh, 17)) +
3954 3.9661774136e10 * pow(sh, 18) +
3955 8.107722e7 * pow(sh, 19) - 2.4528798e7 * pow(sh, 20) -
3956 2.4018498e7 * pow(sh, 21) + 3.225852e6 * pow(sh, 22) +
3957 2.9835e6 * pow(sh, 23) - 343845. * pow(sh, 24) -
3958 18225. * pow(sh, 25))) -
3959 0.0005878894767783657 * pow(Lsh, 2) * pow(-1 + sh, -11) *
3960 pow(sh, -16)*(168. + 1302. * sh + 2760. * pow(sh, 2) -
3961 4440. * pow(sh, 3) - 27552. * pow(sh, 4) -
3962 26796. * pow(sh, 5) + 59136. * pow(sh, 6) +
3963 8.709724e7 * pow(sh, 7) - 9.61414058e8 * pow(sh, 8) +
3964 4.834002949e9 * pow(sh, 9) -
3965 1.08108e6 * Lshb * pow(-1 + sh, 8) * pow(sh, 9) -
3966 1.4570114508e10 * pow(sh, 10) +
3967 2.922310728e10 * pow(sh, 11) -
3968 4.0941204696e10 * pow(sh, 12) +
3969 4.087575366e10 * pow(sh, 13) -
3970 2.9075948544e10 * pow(sh, 14) +
3971 1.443451164e10 * pow(sh, 15) -
3972 4.759554808e9 * pow(sh, 16) +
3973 9.37068023e8 * pow(sh, 17) - 8.3312536e7 * pow(sh, 18) +
3974 6048. * pow(sh, 22) - 7344. * pow(sh, 23) +
3975 864. * pow(sh, 24) + 432. * pow(sh, 25)) -
3976 0.00029394473838918284 * pow(lz, 2) * pow(sh, -16)*
3977 (-84. - 1575. * sh - 14085. * pow(sh, 2) -
3978 79950. * pow(sh, 3) - 323154. * pow(sh, 4) -
3979 987129. * pow(sh, 5) - 2.347059e6 * pow(sh, 6) -
3980 1.78300736e8 * pow(sh, 7) + 2.7027e6 * pow(sh, 8) -
3981 473685. * pow(sh, 9) + 1.243629e6 * pow(sh, 10) +
3982 2.698722e6 * pow(sh, 11) + 2.985906e6 * pow(sh, 12) +
3983 2.433285e6 * pow(sh, 13) + 1.388619e6 * pow(sh, 14) +
3984 84. * pow(sh, 15) + 4.25568e6 * pow(sh, 16) +
3985 4.612608e6 * pow(sh, 17) + 5.189184e6 * pow(sh, 18) +
3986 7.605312e6 * pow(sh, 19) - 6.1281792e7 * pow(sh, 20) +
3987 2.75438592e8 * pow(sh, 21) + 8.71398528e8 * pow(sh, 22) -
3988 5.028236928e9 * pow(sh, 23) +
3989 8.99602704e9 * pow(sh, 24) -
3990 8.653636992e9 * pow(sh, 25) +
3991 5.137333344e9 * pow(sh, 26) -
3992 1.978438176e9 * pow(sh, 27) +
3993 4.97379168e8 * pow(sh, 28) - 7.8990912e7 * pow(sh, 29) +
3994 7.2072e6 * pow(sh, 30) - 288288. * pow(sh, 31)) -
3995 lz * (32540.444444444445 * lmu * pow(sh, -9) +
3996 635.5555555555555 * pow(Lsh, 2) * pow(-1 + sh, -3) *
3997 pow(sh, -7) - 32.282186948853614 * ash * sqrt4sh * sqrtsh *
3998 pow(sh, 3)*(552. - 6504. * sh + 954. * pow(sh, 2) +
3999 96534. * pow(sh, 3) - 256765. * pow(sh, 4) +
4000 311884. * pow(sh, 5) - 219753. * pow(sh, 6) +
4001 97149. * pow(sh, 7) - 27453. * pow(sh, 8) +
4002 4830. * pow(sh, 9) - 483. * pow(sh, 10) +
4003 21. * pow(sh, 11)) +
4004 0.0011757789535567313 * Lsh * pow(-1 + sh, -10) *
4005 pow(sh, -16)*(84. + 735. * sh + 2115. * pow(sh, 2) -
4006 105. * pow(sh, 3) - 13881. * pow(sh, 4) -
4007 27279. * pow(sh, 5) + 2289. * pow(sh, 6) +
4008 8.7045013e7 * pow(sh, 7) - 8.74124185e8 * pow(sh, 8) +
4009 3.955482141e9 * pow(sh, 9) -
4010 1.0600174539e10 * pow(sh, 10) +
4011 1.8626345409e10 * pow(sh, 11) -
4012 2.2429329447e10 * pow(sh, 12) +
4013 1.8752939835e10 * pow(sh, 13) -
4014 1.0757383293e10 * pow(sh, 14) +
4015 4.056843039e9 * pow(sh, 15) -
4016 9.10168681e8 * pow(sh, 16) +
4017 9.2636188e7 * pow(sh, 17) + 75600. * pow(sh, 18) +
4018 9072. * pow(sh, 19) - 22680. * pow(sh, 20) -
4019 9828. * pow(sh, 21) - 2268. * pow(sh, 22) +
4020 2484. * pow(sh, 23) + 54. * pow(sh, 24)) +
4021 3.2627898589097882e-9 * pow(-1 + sh, -9) * pow(sh, -16)*
4022 (-3.57357e6 - 4.099095e6 * sh +
4023 1.6114527e8 * pow(sh, 2) + 8.32841295e8 * pow(sh, 3) +
4024 1.161846972e9 * pow(sh, 4) -
4025 2.313726129e9 * pow(sh, 5) -
4026 9.380058402e9 * pow(sh, 6) +
4027 1.3161069788009e13 * pow(sh, 7) -
4028 1.20977154695826e14 * pow(sh, 8) +
4029 4.94447733407199e14 * pow(sh, 9) -
4030 1.180536055693566e15 * pow(sh, 10) +
4031 1.818479082421521e15 * pow(sh, 11) -
4032 1.879273746539124e15 * pow(sh, 12) +
4033 1.307997164466393e15 * pow(sh, 13) -
4034 5.94838522786482e14 * pow(sh, 14) +
4035 1.62091018204671e14 * pow(sh, 15) -
4036 2.194598322815e13 * pow(sh, 16) +
4037 9.520266140208e12 * pow(sh, 17) -
4038 3.1476144603672e13 * pow(sh, 18) +
4039 6.08108801358e13 * pow(sh, 19) -
4040 8.0522179628688e13 * pow(sh, 20) +
4041 7.1735262216768e13 * pow(sh, 21) +
4042 1.85972284144884e14 * pow(sh, 22) -
4043 1.24257330785952e15 * pow(sh, 23) +
4044 1.516348491395922e15 * pow(sh, 24) +
4045 9.623404446386016e15 * pow(sh, 25) -
4046 5.409785863372176e16 * pow(sh, 26) +
4047 1.4651551524606797e17 * pow(sh, 27) -
4048 2.639914855668742e17 * pow(sh, 28) +
4049 3.46785263569046e17 * pow(sh, 29) -
4050 3.4581010205224454e17 * pow(sh, 30) +
4051 2.67102371124581e17 * pow(sh, 31) -
4052 1.612804363486636e17 * pow(sh, 32) +
4053 7.625393556018638e16 * pow(sh, 33) -
4054 2.809203040360517e16 * pow(sh, 34) +
4055 7.968174811781184e15 * pow(sh, 35) -
4056 1.704011763662208e15 * pow(sh, 36) +
4057 2.65362784294608e14 * pow(sh, 37) -
4058 2.836102317048e13 * pow(sh, 38) +
4059 1.857724700832e12 * pow(sh, 39) -
4060 5.6166718608e10 * pow(sh, 40))) -
4061 1.3581376369088363e-14 * pow(-1 + sh, -10) * pow(sh, -16)*
4062 (-1.522039519e11 + 4.26276275425e11 * sh +
4063 3.315607505925e12 * pow(sh, 2) -
4064 1.9935815474575e13 * pow(sh, 3) -
4065 1.24680380503679e14 * pow(sh, 4) -
4066 9.4734104241777e13 * pow(sh, 5) +
4067 6.72219759402191e14 * pow(sh, 6) +
4068 2.6103747782620196e18 * pow(sh, 7) -
4069 2.5573475986293776e19 * pow(sh, 8) +
4070 1.1271511598880529e20 * pow(sh, 9) -
4071 2.943261559227364e20 * pow(sh, 10) +
4072 5.045168745362524e20 * pow(sh, 11) -
4073 5.941440025879828e20 * pow(sh, 12) +
4074 4.8854165231351536e20 * pow(sh, 13) -
4075 2.7897519097641733e20 * pow(sh, 14) +
4076 1.0745710257767067e20 * pow(sh, 15) -
4077 2.6239676246066704e19 * pow(sh, 16) +
4078 5.344996730119183e18 * pow(sh, 17) -
4079 8.145618788181729e18 * pow(sh, 18) +
4080 1.8132101139625396e19 * pow(sh, 19) -
4081 2.789986937416512e19 * pow(sh, 20) +
4082 2.916708794791665e19 * pow(sh, 21) +
4083 4.506612585214779e19 * pow(sh, 22) -
4084 3.9814881989842764e20 * pow(sh, 23) +
4085 7.61642854944141e20 * pow(sh, 24) +
4086 1.912231127650318e21 * pow(sh, 25) -
4087 1.5547276029752026e22 * pow(sh, 26) +
4088 4.875241535562659e22 * pow(sh, 27) -
4089 9.909337425361018e22 * pow(sh, 28) +
4090 1.4647795672750128e23 * pow(sh, 29) -
4091 1.651330016200192e23 * pow(sh, 30) +
4092 1.453924699389752e23 * pow(sh, 31) -
4093 1.011735575962528e23 * pow(sh, 32) +
4094 5.58883301567305e22 * pow(sh, 33) -
4095 2.447160559972393e22 * pow(sh, 34) +
4096 8.433457866753118e21 * pow(sh, 35) -
4097 2.2566304052798968e21 * pow(sh, 36) +
4098 4.585243554016694e20 * pow(sh, 37) -
4099 6.826335499960051e19 * pow(sh, 38) +
4100 7.011890262150629e18 * pow(sh, 39) -
4101 4.434709029248097e17 * pow(sh, 40) +
4102 1.2998152420597344e16 * pow(sh, 41) -
4103 4.32864432e10 * pow(3.141592653589793, 2)*
4104 (-84. - 735. * sh - 2115. * pow(sh, 2) +
4105 105. * pow(sh, 3) + 13881. * pow(sh, 4) +
4106 27279. * pow(sh, 5) - 2289. * pow(sh, 6) -
4107 7.2918901e7 * pow(sh, 7) + 7.30592797e8 * pow(sh, 8) -
4108 3.297157631e9 * pow(sh, 9) +
4109 8.805108353e9 * pow(sh, 10) -
4110 1.5403793139e10 * pow(sh, 11) +
4111 1.8449564077e10 * pow(sh, 12) -
4112 1.5329503273e10 * pow(sh, 13) +
4113 8.733035583e9 * pow(sh, 14) -
4114 3.270393893e9 * pow(sh, 15) +
4115 7.28627507e8 * pow(sh, 16) -
4116 6.571293e7 * pow(sh, 17) - 3.2675014e7 * pow(sh, 18) +
4117 7.242013e7 * pow(sh, 19) - 7.6735758e7 * pow(sh, 20) -
4118 2.40223726e8 * pow(sh, 21) +
4119 1.627847962e9 * pow(sh, 22) -
4120 2.273776974e9 * pow(sh, 23) -
4121 1.712105386e10 * pow(sh, 24) +
4122 1.18341944808e11 * pow(sh, 25) -
4123 3.9995878208e11 * pow(sh, 26) +
4124 9.10773328608e11 * pow(sh, 27) -
4125 1.531097366656e12 * pow(sh, 28) +
4126 1.98076916752e12 * pow(sh, 29) -
4127 2.016366553344e12 * pow(sh, 30) +
4128 1.634962897568e12 * pow(sh, 31) -
4129 1.06202452928e12 * pow(sh, 32) +
4130 5.53053223008e11 * pow(sh, 33) -
4131 2.30009473408e11 * pow(sh, 34) +
4132 7.5735678304e10 * pow(sh, 35) -
4133 1.945628256e10 * pow(sh, 36) +
4134 3.810741792e9 * pow(sh, 37) -
4135 5.48735616e8 * pow(sh, 38) +
4136 5.4678624e7 * pow(sh, 39) - 3.36336e6 * pow(sh, 40) +
4137 96096. * pow(sh, 41)))) * pow(z, 9);
4139 res += pow(z, 9)*(-531170.0836599696 -
4140 5.221397600319266e6 * (1. - sh) -
4141 2.50820565037545e7 * pow(1. - sh, 2) -
4142 8.849720248983555e7 * pow(1. - sh, 3) +
4143 log(muh)*(-39253.5139329806 -
4144 385822.06984126987 * (1. - sh) -
4145 2.075542349206349e6 * pow(1. - sh, 2) -
4146 8.098428613756614e6 * pow(1. - sh, 3) -
4147 2.5637579174603175e7 * pow(1. - sh, 4)) -
4148 2.6309546792497513e8 * pow(1. - sh, 4) +
4149 log(z)*(75507.81523340371 +
4150 370284.8608070748 * (1. - sh) +
4151 2.8644751808652747e6 * pow(1. - sh, 2) +
4152 1.3385139112281743e7 * pow(1. - sh, 3) +
4153 log(muh)*(-32540.444444444445 -
4154 292864. * (1. - sh) -
4155 1.46432e6 * pow(1. - sh, 2) -
4156 5.369173333333333e6 * pow(1. - sh, 3) -
4157 1.610752e7 * pow(1. - sh, 4)) +
4158 4.39338775383846e7 * pow(1. - sh, 4)) +
4159 (48060.242210464436 + 446131.6613756614 * (1. - sh) +
4160 2.3941590476190476e6 * pow(1. - sh, 2) +
4161 8.828220007054674e6 * pow(1. - sh, 3) +
4162 2.6455532740740743e7 * pow(1. - sh, 4)) *
4168 res += 4.0250146484590727e-7 * pow(3.141592653589793, 2) * pow(sh, -17)*
4169 (109395. + 2.16216e6 * sh + 2.0430696e7 * pow(sh, 2) +
4170 1.22866128e8 * pow(sh, 3) + 5.27900436e8 * pow(sh, 4) +
4171 1.72236568e9 * pow(sh, 5) + 4.41295356e9 * pow(sh, 6) +
4172 8.96714376e9 * pow(sh, 7) + 1.385377093e10 * pow(sh, 8) +
4173 8.96714376e9 * pow(sh, 9) + 4.41295356e9 * pow(sh, 10) +
4174 1.72236568e9 * pow(sh, 11) + 5.27900436e8 * pow(sh, 12) +
4175 1.22866128e8 * pow(sh, 13) + 2.0430696e7 * pow(sh, 14) +
4176 2.16216e6 * pow(sh, 15) + 109395. * pow(sh, 16)) * pow(z, 9.5);
4178 res += (180375.25585144534 + 1.6233773026630078e6 * (1. - sh) +
4179 8.333336820336774e6 * pow(1. - sh, 2) +
4180 3.192642028570582e7 * pow(1. - sh, 3) +
4181 1.0132579997454941e8 * pow(1. - sh, 4)) * pow(z, 9.5);
4186 res += (-1. * lmu * (134616.0761904762 * pow(sh, -10) -
4187 109824. * Lsh * pow(sh, -10)) +
4188 4617.798054844007 * pow(-1 + sh, -3) * pow(sh, -8) -
4189 3841.58024691358 * Li3sh * pow(-1 + sh, -3) * pow(sh, -8) +
4190 1280.5267489711935 * pow(Lsh, 3) * pow(-1 + sh, -3) *
4191 pow(sh, -8) - 36.317460317460316 * cl2 * sqrt4sh * sqrtsh *
4192 pow(sh, 3)*(72. - 4364. * sh + 28750. * pow(sh, 2) +
4193 23168. * pow(sh, 3) - 476659. * pow(sh, 4) +
4194 1.255353e6 * pow(sh, 5) - 1.649661e6 * pow(sh, 6) +
4195 1.316274e6 * pow(sh, 7) - 686175. * pow(sh, 8) +
4196 239883. * pow(sh, 9) - 55951. * pow(sh, 10) +
4197 8372. * pow(sh, 11) - 728. * pow(sh, 12) + 28. * pow(sh, 13))
4198 - 0.005643738977072311 * ash * sqrt4sh * sqrtsh *
4199 pow(-1 + sh, -3) * pow(sh, 3)*
4200 (-1.003446e6 + 6.399262e7 * sh -
4201 5.88476654e8 * pow(sh, 2) + 1.052685443e9 * pow(sh, 3) +
4202 6.558619027e9 * pow(sh, 4) -
4203 3.9015984112e10 * pow(sh, 5) +
4204 9.831642973e10 * pow(sh, 6) -
4205 1.5050270954e11 * pow(sh, 7) +
4206 1.55340604554e11 * pow(sh, 8) -
4207 1.13220265218e11 * pow(sh, 9) +
4208 5.9521934912e10 * pow(sh, 10) -
4209 2.2692174847e10 * pow(sh, 11) +
4210 6.219574195e9 * pow(sh, 12) -
4211 1.194652762e9 * pow(sh, 13) + 1.5263422e8 * pow(sh, 14) -
4212 1.1648401e7 * pow(sh, 15) + 401669. * pow(sh, 16)) +
4213 108.95238095238095 * pow(ash, 2) * pow(-1 + sh, -4) *
4214 pow(sh, 4)*(1638. - 34902. * sh + 187614. * pow(sh, 2) -
4215 73296. * pow(sh, 3) - 2.717613e6 * pow(sh, 4) +
4216 1.2278721e7 * pow(sh, 5) - 2.8715127e7 * pow(sh, 6) +
4217 4.3584072e7 * pow(sh, 7) - 4.6464222e7 * pow(sh, 8) +
4218 3.6125334e7 * pow(sh, 9) - 2.0866662e7 * pow(sh, 10) +
4219 9.012312e6 * pow(sh, 11) - 2.899701e6 * pow(sh, 12) +
4220 684873. * pow(sh, 13) - 115255. * pow(sh, 14) +
4221 13076. * pow(sh, 15) - 896. * pow(sh, 16) + 28. * pow(sh, 17)
4222 ) + 0.007054673721340388 * Li2sh * pow(-1 + sh, -11) *
4223 pow(sh, -10)*(3.891888e6 - 4.2698656e7 * sh +
4224 2.12860648e8 * pow(sh, 2) - 6.36197562e8 * pow(sh, 3) +
4225 1.265813472e9 * pow(sh, 4) - 1.758542016e9 * pow(sh, 5) +
4226 1.737624672e9 * pow(sh, 6) - 1.217575632e9 * pow(sh, 7) +
4227 5.89815408e8 * pow(sh, 8) - 1.86234048e8 * pow(sh, 9) +
4228 3.3715396e7 * pow(sh, 10) - 2.523664e6 * pow(sh, 11) -
4229 1008. * pow(sh, 17) + 252. * pow(sh, 18) + 63. * pow(sh, 19))
4230 - Lsh * (-3841.58024691358 * Li2sh * pow(-1 + sh, -3) *
4231 pow(sh, -8) + 36.317460317460316 * ash * sqrt4sh * sqrtsh *
4232 pow(sh, 3)*(72. - 4364. * sh + 28750. * pow(sh, 2) +
4233 23168. * pow(sh, 3) - 476659. * pow(sh, 4) +
4234 1.255353e6 * pow(sh, 5) - 1.649661e6 * pow(sh, 6) +
4235 1.316274e6 * pow(sh, 7) - 686175. * pow(sh, 8) +
4236 239883. * pow(sh, 9) - 55951. * pow(sh, 10) +
4237 8372. * pow(sh, 11) - 728. * pow(sh, 12) +
4238 28. * pow(sh, 13)) +
4239 6.22105266432133e-8 * pow(-1 + sh, -12) * pow(sh, -18)*
4240 (-268128. + 217952. * sh + 1.5087492e7 * pow(sh, 2) +
4241 5.5721052e7 * pow(sh, 3) - 6.214128e6 * pow(sh, 4) -
4242 4.11907536e8 * pow(sh, 5) -
4243 6.39424296e8 * pow(sh, 6) +
4244 9.51140904e8 * pow(sh, 7) +
4245 2.797636941792e12 * pow(sh, 8) -
4246 3.3528009426912e13 * pow(sh, 9) +
4247 1.8405229744577e14 * pow(sh, 10) +
4248 5.1459408e10 * pow(3.141592653589793, 2) *
4249 pow(-1 + sh, 9) * pow(sh, 10) -
4250 6.11232376478328e14 * pow(sh, 11) +
4251 1.367454094140762e15 * pow(sh, 12) -
4252 2.171310899030364e15 * pow(sh, 13) +
4253 2.509372042935534e15 * pow(sh, 14) -
4254 2.126984603171022e15 * pow(sh, 15) +
4255 1.312509475395804e15 * pow(sh, 16) -
4256 5.75196391765836e14 * pow(sh, 17) +
4257 1.70043077466117e14 * pow(sh, 18) -
4258 3.0485673234577e13 * pow(sh, 19) -
4259 113400. * Lshb * (-1. + sh) * pow(sh, 8)*
4260 (3.891888e6 - 4.2698656e7 * sh +
4261 2.12860648e8 * pow(sh, 2) -
4262 6.36197562e8 * pow(sh, 3) +
4263 1.265813472e9 * pow(sh, 4) -
4264 1.758542016e9 * pow(sh, 5) +
4265 1.737624672e9 * pow(sh, 6) -
4266 1.217575632e9 * pow(sh, 7) +
4267 5.89815408e8 * pow(sh, 8) -
4268 1.86234048e8 * pow(sh, 9) +
4269 3.3715396e7 * pow(sh, 10) -
4270 2.523664e6 * pow(sh, 11) - 1008. * pow(sh, 17) +
4271 252. * pow(sh, 18) + 63. * pow(sh, 19)) +
4272 2.505502772136e12 * pow(sh, 20) +
4273 3.810922668e9 * pow(sh, 21) +
4274 9.25325478e8 * pow(sh, 22) -
4275 6.18090102e8 * pow(sh, 23) -
4276 3.71771532e8 * pow(sh, 24) + 7.86267e7 * pow(sh, 25) +
4277 4.1725125e7 * pow(sh, 26) - 4.8951e6 * pow(sh, 27) -
4278 187425. * pow(sh, 28))) -
4279 0.00007838526357044875 * pow(Lsh, 2) * pow(-1 + sh, -12) *
4280 pow(sh, -18)*(-1008. - 8848. * sh - 23688. * pow(sh, 2) +
4281 17352. * pow(sh, 3) + 212592. * pow(sh, 4) +
4282 311184. * pow(sh, 5) - 400176. * pow(sh, 6) -
4283 1.492176e6 * pow(sh, 7) - 2.343842424e9 * pow(sh, 8) +
4284 2.82465414e10 * pow(sh, 9) -
4285 1.56076722399e11 * pow(sh, 10) -
4286 2.450448e7 * Lshb * pow(-1 + sh, 9) * pow(sh, 10) +
4287 5.22413429229e11 * pow(sh, 11) -
4288 1.178919945486e12 * pow(sh, 12) +
4289 1.88907310284e12 * pow(sh, 13) -
4290 2.203548370128e12 * pow(sh, 14) +
4291 1.885005025848e12 * pow(sh, 15) -
4292 1.173372074448e12 * pow(sh, 16) +
4293 5.18137391592e11 * pow(sh, 17) -
4294 1.53980764326e11 * pow(sh, 18) +
4295 2.7628439638e10 * pow(sh, 19) -
4296 2.260854288e9 * pow(sh, 20) + 45360. * pow(sh, 25) -
4297 56700. * pow(sh, 26) + 8505. * pow(sh, 27) +
4298 2835. * pow(sh, 28)) -
4299 0.00007838526357044875 * pow(lz, 2) * pow(sh, -18)*
4300 (-252. - 5236. * sh - 52122. * pow(sh, 2) -
4301 330990. * pow(sh, 3) - 1.50586e6 * pow(sh, 4) -
4302 5.221788e6 * pow(sh, 5) - 1.4306214e7 * pow(sh, 6) -
4303 3.1500898e7 * pow(sh, 7) - 2.399232528e9 * pow(sh, 8) +
4304 3.675672e7 * pow(sh, 9) - 1.1865305e7 * pow(sh, 10) +
4305 7.819275e6 * pow(sh, 11) + 2.920308e7 * pow(sh, 12) +
4306 3.763735e7 * pow(sh, 13) + 3.5618775e7 * pow(sh, 14) +
4307 2.7424047e7 * pow(sh, 15) + 1.5311086e7 * pow(sh, 16) +
4308 252. * pow(sh, 17) + 5.397678e7 * pow(sh, 18) +
4309 5.722002e7 * pow(sh, 19) + 6.231654e7 * pow(sh, 20) +
4310 6.996132e7 * pow(sh, 21) + 2.7057888e8 * pow(sh, 22) -
4311 3.18810492e9 * pow(sh, 23) + 7.5729654e9 * pow(sh, 24) +
4312 4.194296964e10 * pow(sh, 25) -
4313 2.0547002619e11 * pow(sh, 26) +
4314 3.8218629735e11 * pow(sh, 27) -
4315 4.0430890071e11 * pow(sh, 28) +
4316 2.7416659842e11 * pow(sh, 29) -
4317 1.2521582073e11 * pow(sh, 30) +
4318 3.911266359e10 * pow(sh, 31) -
4319 8.26458633e9 * pow(sh, 32) + 1.13189076e9 * pow(sh, 33) -
4320 9.081072e7 * pow(sh, 34) + 3.24324e6 * pow(sh, 35)) -
4321 lz * (109824. * lmu * pow(sh, -10) +
4322 1920.79012345679 * pow(Lsh, 2) * pow(-1 + sh, -3) *
4323 pow(sh, -8) + 72.63492063492063 * ash * sqrt4sh * sqrtsh *
4324 pow(sh, 3)*(72. - 4364. * sh + 28750. * pow(sh, 2) +
4325 23168. * pow(sh, 3) - 476659. * pow(sh, 4) +
4326 1.255353e6 * pow(sh, 5) - 1.649661e6 * pow(sh, 6) +
4327 1.316274e6 * pow(sh, 7) - 686175. * pow(sh, 8) +
4328 239883. * pow(sh, 9) - 55951. * pow(sh, 10) +
4329 8372. * pow(sh, 11) - 728. * pow(sh, 12) +
4330 28. * pow(sh, 13)) +
4331 0.0001567705271408975 * Lsh * pow(-1 + sh, -11) *
4332 pow(sh, -18)*(-504. - 4928. * sh - 16772. * pow(sh, 2) -
4333 8096. * pow(sh, 3) + 98200. * pow(sh, 4) +
4334 253792. * pow(sh, 5) + 53704. * pow(sh, 6) -
4335 692384. * pow(sh, 7) - 2.34466232e9 * pow(sh, 8) +
4336 2.5895265088e10 * pow(sh, 9) -
4337 1.30081243153e11 * pow(sh, 10) +
4338 3.91980399854e11 * pow(sh, 11) -
4339 7.87080144424e11 * pow(sh, 12) +
4340 1.105852892984e12 * pow(sh, 13) -
4341 1.109622439576e12 * pow(sh, 14) +
4342 7.95460832936e11 * pow(sh, 15) -
4343 3.99495040168e11 * pow(sh, 16) +
4344 1.33989748088e11 * pow(sh, 17) -
4345 2.7053178394e10 * pow(sh, 18) +
4346 2.494260648e9 * pow(sh, 19) + 864864. * pow(sh, 20) +
4347 769608. * pow(sh, 21) - 13608. * pow(sh, 22) -
4348 222264. * pow(sh, 23) - 69120. * pow(sh, 24) -
4349 9000. * pow(sh, 25) + 19215. * pow(sh, 26) +
4350 315. * pow(sh, 27)) +
4351 2.5590508697331672e-11 * pow(-1 + sh, -10) * pow(sh, -18)*
4352 (3.25909584e8 + 3.86898512e8 * sh -
4353 1.7890959086e10 * pow(sh, 2) -
4354 1.0389775539e11 * pow(sh, 3) -
4355 1.8235127911e11 * pow(sh, 4) +
4356 2.39868807178e11 * pow(sh, 5) +
4357 1.439309125254e12 * pow(sh, 6) +
4358 1.482637674518e12 * pow(sh, 7) -
4359 6.794302525868858e15 * pow(sh, 8) +
4360 6.901772268373983e16 * pow(sh, 9) -
4361 3.158878275063268e17 * pow(sh, 10) +
4362 8.578247944222483e17 * pow(sh, 11) -
4363 1.5323085794390446e18 * pow(sh, 12) +
4364 1.8841304231440591e18 * pow(sh, 13) -
4365 1.6184014070975416e18 * pow(sh, 14) +
4366 9.617660963186045e17 * pow(sh, 15) -
4367 3.8015225593619155e17 * pow(sh, 16) +
4368 9.092524285521288e16 * pow(sh, 17) -
4369 9.590076560304648e15 * pow(sh, 18) -
4370 4.556468500926695e15 * pow(sh, 19) +
4371 1.7572146654042552e16 * pow(sh, 20) -
4372 3.9615878221273704e16 * pow(sh, 21) +
4373 5.638676809242492e16 * pow(sh, 22) -
4374 4.406121643569693e16 * pow(sh, 23) +
4375 3.2318412281568824e16 * pow(sh, 24) -
4376 4.259492022110168e17 * pow(sh, 25) +
4377 2.1201430825186813e18 * pow(sh, 26) -
4378 2.4923037486423537e18 * pow(sh, 27) -
4379 1.7838580318483542e19 * pow(sh, 28) +
4380 1.0701104247446746e20 * pow(sh, 29) -
4381 3.1929857696090705e20 * pow(sh, 30) +
4382 6.474487284281146e20 * pow(sh, 31) -
4383 9.760291519177703e20 * pow(sh, 32) +
4384 1.1392452492311442e21 * pow(sh, 33) -
4385 1.0520263768751638e21 * pow(sh, 34) +
4386 7.775724274489533e20 * pow(sh, 35) -
4387 4.624318941820542e20 * pow(sh, 36) +
4388 2.213656928475123e20 * pow(sh, 37) -
4389 8.494727327998501e19 * pow(sh, 38) +
4390 2.5900145606169035e19 * pow(sh, 39) -
4391 6.181939718636802e18 * pow(sh, 40) +
4392 1.1286121824774998e18 * pow(sh, 41) -
4393 1.5195925385464826e17 * pow(sh, 42) +
4394 1.420123992415128e16 * pow(sh, 43) -
4395 8.2168151561496e14 * pow(sh, 44) +
4396 2.214605244852e13 * pow(sh, 45))) -
4397 4.177278391107532e-18 * pow(-1 + sh, -11) * pow(sh, -18)*
4398 (3.16414088198208e14 - 9.74539823921664e14 * sh -
4399 8.087033764735926e15 * pow(sh, 2) +
4400 5.684694862939435e16 * pow(sh, 3) +
4401 4.0319420043801754e17 * pow(sh, 4) +
4402 4.882604186351892e17 * pow(sh, 5) -
4403 1.9325128069129482e18 * pow(sh, 6) -
4404 6.292691026868566e18 * pow(sh, 7) -
4405 2.4287339249894217e22 * pow(sh, 8) +
4406 2.612488077307783e23 * pow(sh, 9) -
4407 1.276336786564487e24 * pow(sh, 10) +
4408 3.7407720820676834e24 * pow(sh, 11) -
4409 7.312076867879097e24 * pow(sh, 12) +
4410 1.0021549960712915e25 * pow(sh, 13) -
4411 9.851452036875564e24 * pow(sh, 14) +
4412 6.978601203194415e24 * pow(sh, 15) -
4413 3.521127295712042e24 * pow(sh, 16) +
4414 1.224012074235529e24 * pow(sh, 17) -
4415 2.690422527577656e23 * pow(sh, 18) +
4416 6.92847965398663e21 * pow(sh, 19) +
4417 1.0487308963151102e23 * pow(sh, 20) -
4418 2.7045906540143195e23 * pow(sh, 21) +
4419 4.4609820546041925e23 * pow(sh, 22) -
4420 4.3341866014896366e23 * pow(sh, 23) +
4421 3.0284128932535685e23 * pow(sh, 24) -
4422 2.8502065219122033e24 * pow(sh, 25) +
4423 1.6744213714265902e25 * pow(sh, 26) -
4424 3.4546092096924117e25 * pow(sh, 27) -
4425 6.6958214497636675e25 * pow(sh, 28) +
4426 6.719624033223081e26 * pow(sh, 29) -
4427 2.3634332855143025e27 * pow(sh, 30) +
4428 5.406861258181357e27 * pow(sh, 31) -
4429 9.10419041097793e27 * pow(sh, 32) +
4430 1.1865952156817548e28 * pow(sh, 33) -
4431 1.2283961231268607e28 * pow(sh, 34) +
4432 1.0245250602869844e28 * pow(sh, 35) -
4433 6.935097181125567e27 * pow(sh, 36) +
4434 3.819594455600949e27 * pow(sh, 37) -
4435 1.7090015768187537e27 * pow(sh, 38) +
4436 6.177750381749111e26 * pow(sh, 39) -
4437 1.7862394817725466e26 * pow(sh, 40) +
4438 4.0667139734465495e25 * pow(sh, 41) -
4439 7.117942930727479e24 * pow(sh, 42) +
4440 9.229373088621336e23 * pow(sh, 43) -
4441 8.339200912205272e22 * pow(sh, 44) +
4442 4.6814118287209613e21 * pow(sh, 45) -
4443 1.2279937968685318e20 * pow(sh, 46) +
4444 6.2548910424e12 * pow(3.141592653589793, 2)*
4445 (-1512. - 14784. * sh - 50316. * pow(sh, 2) -
4446 24288. * pow(sh, 3) + 294600. * pow(sh, 4) +
4447 761376. * pow(sh, 5) + 161112. * pow(sh, 6) -
4448 2.077152e6 * pow(sh, 7) - 5.96371776e9 * pow(sh, 8) +
4449 6.5734816224e10 * pow(sh, 9) -
4450 3.29439299647e11 * pow(sh, 10) +
4451 9.89864674466e11 * pow(sh, 11) -
4452 1.980655533976e12 * pow(sh, 12) +
4453 2.77128605852e12 * pow(sh, 13) -
4454 2.767417361128e12 * pow(sh, 14) +
4455 1.973276912096e12 * pow(sh, 15) -
4456 9.85353747448e11 * pow(sh, 16) +
4457 3.28634709368e11 * pow(sh, 17) -
4458 6.613484449e10 * pow(sh, 18) +
4459 6.792399064e9 * pow(sh, 19) -
4460 3.03847082e9 * pow(sh, 20) +
4461 7.83321994e9 * pow(sh, 21) -
4462 1.36334135e10 * pow(sh, 22) +
4463 2.612527414e10 * pow(sh, 23) -
4464 1.2043334018e11 * pow(sh, 24) +
4465 5.11152445252e11 * pow(sh, 25) -
4466 7.00772118731e11 * pow(sh, 26) -
4467 5.082162807869e12 * pow(sh, 27) +
4468 3.7496386828908e13 * pow(sh, 28) -
4469 1.3686666259266e14 * pow(sh, 29) +
4470 3.4205865415722e14 * pow(sh, 30) -
4471 6.4118721989322e14 * pow(sh, 31) +
4472 9.3979081997658e14 * pow(sh, 32) -
4473 1.10202690808782e15 * pow(sh, 33) +
4474 1.04788397968254e15 * pow(sh, 34) -
4475 8.141721939351e14 * pow(sh, 35) +
4476 5.1869434238622e14 * pow(sh, 36) -
4477 2.7096438795426e14 * pow(sh, 37) +
4478 1.1569642322238e14 * pow(sh, 38) -
4479 4.010774591538e13 * pow(sh, 39) +
4480 1.116674149734e13 * pow(sh, 40) -
4481 2.4564705165e12 * pow(sh, 41) +
4482 4.1666228604e11 * pow(sh, 42) -
4483 5.249021778e10 * pow(sh, 43) +
4484 4.61837376e9 * pow(sh, 44) -
4485 2.5297272e8 * pow(sh, 45) + 6.48648e6 * pow(sh, 46)))) *
4488 res += pow(z, 10)*(-1.8827679840714834e6 -
4489 2.081717517867475e7 * (1. - sh) -
4490 1.0963380905010003e8 * pow(1. - sh, 2) -
4491 4.219709458628686e8 * pow(1. - sh, 3) +
4492 log(muh)*(-134616.0761904762 -
4493 1.455984761904762e6 * (1. - sh) -
4494 8.55703619047619e6 * pow(1. - sh, 2) -
4495 3.624158476190476e7 * pow(1. - sh, 3) -
4496 1.2382547047619048e8 * pow(1. - sh, 4)) -
4497 1.3631906422358756e9 * pow(1. - sh, 4) +
4498 log(z)*(294861.71513582155 +
4499 1.5577045832132744e6 * (1. - sh) +
4500 1.3273873434568334e7 * pow(1. - sh, 2) +
4501 6.637266945174166e7 * pow(1. - sh, 3) +
4502 log(muh)*(-109824. - 1.09824e6 * (1. - sh) -
4503 6.04032e6 * pow(1. - sh, 2) -
4504 2.416128e7 * pow(1. - sh, 3) -
4505 7.852416e7 * pow(1. - sh, 4)) +
4506 2.3142944320276135e8 * pow(1. - sh, 4)) +
4507 (171022.4 + 1.777716148148148e6 * (1. - sh) +
4508 1.0547910696296297e7 * pow(1. - sh, 2) +
4509 4.229024415943562e7 * pow(1. - sh, 3) +
4510 1.3716763538765433e8 * pow(1. - sh, 4)) *