2478 double lmu = log(muh);
2479 double Lsh = log(sh);
2480 double sh2 = pow(sh, 2);
2481 double sh3 = pow(sh, 3);
2482 double lz2 = pow(lz, 2);
2483 double lz3 = pow(lz, 3);
2484 double sqrtz = pow(z, 0.5);
2485 double z2 = pow(z, 2);
2486 double z3 = pow(z, 3);
2492 res = -4.61812 - 0.0493827 * Lsh - 0.06584362139917696 * lmu * Lsh +
2493 2.814815 * lz + lmu * (-1.953360768175583 +
2494 1.1851851851851851 * lz) + 4.47441 * sh + 0.74074 * lz * sh +
2495 14.4621 * z - 0.592593 * Lsh * z + 4.796605 * lz * z +
2496 71.3855 * sh * z - 2.66667 * Lsh * sh * z + 4.238385 * lz * sh * z -
2497 0.592595 * Lsh * lz * sh * z - 1.0534979423868314 * pow(lmu, 2) -
2498 0.2962975 * z * lz2 + 3.134725 * sh * z * lz2 -
2499 0.098765375 * z * lz3 - 0.098765375 * sh * z * lz3 +
2500 0.0493826875 * sh * z * pow(lz, 4) + 37.1282 * sh2 +
2501 0.0164609 * Lsh * sh2 + 11.03105 * lz * sh2 +
2502 212.74 * z * sh2 - 5.33333 * Lsh * z * sh2 -
2503 10.96075 * lz * z * sh2 - 1.185185 * Lsh * lz * z * sh2 +
2504 1.3333325 * lz2 * sh2 +
2505 14.2931 * z * lz2 * sh2 +
2506 0.098765375 * lz3 * sh2 -
2507 0.29629625 * z * lz3 * sh2 +
2508 0.148148125 * z * pow(lz, 4) * sh2 + 79.7475 * sh3 +
2509 0.0219479 * Lsh * sh3 + 28.65855 * lz * sh3 +
2510 425.579 * z * sh3 - 8.2963 * Lsh * z * sh3 -
2511 34.4008 * lz * z * sh3 - 1.77778 * Lsh * lz * z * sh3 +
2512 2.814825 * lz2 * sh3 +
2513 32.33925 * z * lz2 * sh3 +
2514 0.29629625 * lz3 * sh3 -
2515 0.6913575 * z * lz3 * sh3 +
2516 0.29629625 * z * pow(lz, 4) * sh3 +
2517 lmu * sh3*(0.00020902736952119668 -
2518 0.0037624926513815404 * pow(z, -3)) -
2519 0.0759415 * sh3 * pow(z, -3) +
2520 lmu * sh2*(0.0014109347442680777 -
2521 0.025396825396825397 * pow(z, -2)) -
2522 0.403158 * sh2 * pow(z, -2) -
2523 0.00480894 * sh3 * pow(z, -2) +
2524 lmu * sh * (0.01316872427983539 -
2525 0.23703703703703705 * pow(z, -1)) -
2526 2.48507 * sh * pow(z, -1) - 0.0613169 * sh2 * pow(z, -1) -
2527 1.81002 * sh3 * pow(z, -1) -
2528 0.4597295 * lz * sh3 * pow(z, -1) -
2529 0.04938275 * lz2 * sh3 * pow(z, -1) +
2530 sh3*(7.79821 * pow(z, -0.5) - 319.726 * sqrtz) -
2531 23.3946 * sh * sqrtz - 140.368 * sh2 * sqrtz +
2532 3.8991 * pow(z, 1.5) - 16.0864 * z2 +
2533 4.95977 * Lsh * z2 + 27.12195 * lz * z2 -
2534 0.592595 * Lsh * lz * z2 - 18.1301 * sh * z2 +
2535 18.6539 * Lsh * sh * z2 + 74.798 * lz * sh * z2 -
2536 2.37037 * Lsh * lz * sh * z2 -
2537 3.85185 * lz2 * z2 -
2538 0.5925925 * Lsh * lz2 * z2 -
2539 12.2963 * sh * lz2 * z2 -
2540 2.37037 * Lsh * sh * lz2 * z2 -
2541 0.4938275 * lz3 * z2 -
2542 1.3827125 * sh * lz3 * z2 -
2543 44.6829 * sh2 * z2 +
2544 40.786 * Lsh * sh2 * z2 +
2545 136.0075 * lz * sh2 * z2 -
2546 7.1111 * Lsh * lz * sh2 * z2 -
2547 29.77775 * lz2 * sh2 * z2 -
2548 5.333325 * Lsh * lz2 * sh2 * z2 -
2549 2.6666625 * lz3 * sh2 * z2 -
2550 87.8946 * sh3 * z2 +
2551 70.2698 * Lsh * sh3 * z2 +
2552 208.806 * lz * sh3 * z2 -
2553 15.80245 * Lsh * lz * sh3 * z2 -
2554 56.79 * lz2 * sh3 * z2 -
2555 9.481475 * Lsh * lz2 * sh3 * z2 -
2556 4.345675 * lz3 * sh3 * z2 - 14.73 * z3 -
2557 9.20287 * Lsh * z3 - 14.28805 * lz * z3 -
2558 0.52675 * Lsh * lz * z3 - 72.89 * sh * z3 -
2559 41.6104 * Lsh * sh * z3 - 34.0675 * lz * sh * z3 -
2560 1.185185 * Lsh * lz * sh * z3 +
2561 5.037025 * lz2 * z3 +
2562 0.7901225 * Lsh * lz2 * z3 +
2563 15.901225 * sh * lz2 * z3 +
2564 3.55555 * Lsh * sh * lz2 * z3 -
2565 137.203 * sh2 * z3 -
2566 111.356 * Lsh * sh2 * z3 -
2567 49.7185 * lz * sh2 * z3 +
2568 42.22225 * lz2 * sh2 * z3 +
2569 9.481475 * Lsh * lz2 * sh2 * z3 -
2570 279.268 * sh3 * z3 -
2571 231.893 * Lsh * sh3 * z3 -
2572 73.4265 * lz * sh3 * z3 +
2573 5.92595 * Lsh * lz * sh3 * z3 +
2574 82.81475 * lz2 * sh3 * z3 +
2575 19.753075 * Lsh * lz2 * sh3 * z3;
2577 double Lshb = log(1. - sh);
2578 double sqrtsh = sqrt(sh);
2579 double sqrt4sh = sqrt(4. - sh);
2580 double Li2sh = real(
Li2(sh));
2581 double Li3sh = real(
Li3(sh));
2582 double Li3shb = real(
Li3(1. - sh));
2583 double Li4sh = real(
Li4(sh));
2584 double ash = asin(sqrtsh / 2.);
2585 double cl2 =
Cl2(2. * ash);
2586 double cl3 =
Cl3(2. * ash);
2589 res += (11.39728026699467 - 1.0534979423868314 * pow(lmu, 2) -
2590 0.03292181069958848 * Li2sh * (-8. + 17. * sh) *
2591 pow(-1 + sh, -1) - 0.03292181069958848 * cl2 * (2. + sh) *
2592 sqrt4sh * sqrtsh * pow(sh, -2) -
2593 lmu * (-1.1193415637860082 * Lsh +
2594 0.13168724279835392 * ash * (2. + sh) * sqrt4sh * sqrtsh *
2595 pow(sh, -2) + 0.13168724279835392 * (-2. + 29. * sh) *
2596 pow(sh, -1)) + 0.14814814814814814 * pow(Lsh, 2) *
2597 pow(-1 + sh, -3) * pow(sh, 2) -
2598 Lsh * (0.03292181069958848 * ash * (2. + sh) * sqrt4sh *
2599 sqrtsh * pow(sh, -2) +
2600 0.01646090534979424 * pow(-1 + sh, -2)*
2601 (-168. + 307. * sh - 121. * pow(sh, 2) +
2602 2. * Lshb * (8. - 25. * sh + 17. * pow(sh, 2)))) -
2603 0.03292181069958848 * ash * sqrt4sh * sqrtsh * pow(-1 + sh, -3) *
2604 pow(sh, -2)*(-4. + 3. * sh + 18. * pow(sh, 2) -
2605 16. * pow(sh, 3) + 5. * pow(sh, 4)) -
2606 0.09876543209876543 * pow(ash, 2) * pow(-1 + sh, -4) *
2607 pow(sh, -1)*(12. - 33. * sh + 18. * pow(sh, 2) -
2608 4. * pow(sh, 4) + pow(sh, 5)) +
2609 0.001828989483310471 * pow(-1 + sh, -4) * pow(sh, -1)*
2610 (144. + 2. * sh * (-2068. +
2611 75. * pow(3.141592653589793, 2)) +
2612 (15077. - 618. * pow(3.141592653589793, 2)) *
2614 (-214. + 9. * pow(3.141592653589793, 2)) * pow(sh, 3) +
2615 (14249. - 615. * pow(3.141592653589793, 2)) *
2616 pow(sh, 4) + 2. * (-1753. +
2617 78. * pow(3.141592653589793, 2)) * pow(sh, 5))) *
2620 res += pow(z, 0)*(5.605483340182181 -
2621 1.9301201034078073 * (1. - sh) -
2622 1.1708211541924651 * pow(1. - sh, 2) -
2623 0.8303057591438543 * pow(1. - sh, 3) +
2624 log(muh)*(-3.9138369114536733 -
2625 1.136144688522592 * (1. - sh) -
2626 0.5573265536192934 * pow(1. - sh, 2) -
2627 0.37357014711105585 * pow(1. - sh, 3) -
2628 0.27973271946953165 * pow(1. - sh, 4)) -
2629 0.6410250317311895 * pow(1. - sh, 4) -
2630 1.0534979423868314 * pow(log(muh), 2));
2642 res += (-0.09876543209876543 * sh * pow(Lsh, 4) * pow(-1 + sh, -3) +
2643 0.5925925925925926 * Li2sh * (1. + 2. * sh) *
2644 pow(-1 + sh, -2) - 0.19753086419753085 * (-3. + sh) *
2645 pow(Lsh, 3) * pow(-1 + sh, -2) -
2646 7.111111111111111 * lmu * pow(sh, -1) -
2647 lz * (0.5925925925925926 * Lsh *
2648 (2. + 2. * Lshb * (-1. + sh) + sh) * pow(-1 + sh, -2) +
2649 1.1851851851851851 * Li2sh * pow(-1 + sh, -1) -
2650 0.5925925925925926 * pow(Lsh, 2) * pow(-1 + sh, -1) -
2651 0.19753086419753085 *
2652 (144. + sh * (-135. + pow(3.141592653589793, 2))) *
2653 pow(-1 + sh, -1) * pow(sh, -1)) +
2654 0.7123300166871669 * pow(-1 + sh, -3)*
2655 (-3. + 3. * sh + 2. * pow(sh, 2)) -
2656 0.5925925925925926 * ash * sqrt4sh * sqrtsh * pow(-1 + sh, -3) *
2657 pow(sh, -1)*(2. - 11. * sh + 3. * pow(sh, 2)) +
2658 0.5925925925925926 * Li3sh * pow(-1 + sh, -3)*
2659 (-3. - 3. * sh + 4. * pow(sh, 2)) -
2660 pow(Lsh, 2)*(1.1851851851851851 * Li2sh * sh *
2662 0.09876543209876543 * pow(-1 + sh, -3)*
2663 (12. + 3. * Lshb * (-7. + 5. * sh) +
2664 2. * sh * (-6. + pow(3.141592653589793, 2)) +
2666 3.5555555555555554 * pow(ash, 2) * pow(-1 + sh, -4)*
2667 (-1. + 8. * sh - 5. * pow(sh, 2) + pow(sh, 3)) -
2668 Lsh * (5.698640133497335 * sh * pow(-1 + sh, -3) -
2669 4.7407407407407405 * Li3sh * sh * pow(-1 + sh, -3) +
2670 0.5925925925925926 * Li2sh * pow(-1 + sh, -3)*
2671 (-5. + sh + 2. * pow(sh, 2)) -
2672 0.09876543209876543 * pow(-1 + sh, -3) * pow(sh, -1)*
2673 (-288. + sh * (906. - 6. * Lshb -
2674 7. * pow(3.141592653589793, 2)) +
2675 (-903. - 6. * Lshb + 5. * pow(3.141592653589793, 2)) *
2676 pow(sh, 2) + 3. * (95. + 4. * Lshb) * pow(sh, 3))) -
2677 0.019753086419753086 * pow(-1 + sh, -4) * pow(sh, -1)*
2678 (-4. * (-1. + sh) * pow(3.141592653589793, 4) * pow(sh, 2) +
2679 15. * (-1. + sh)*(-20. + 48. * sh +
2680 3. * (-13. + 8. * Li4sh) * pow(sh, 2) + 11. * pow(sh, 3))\
2681 + 5. * pow(3.141592653589793, 2)*
2682 (-3. + 8. * sh - 10. * pow(sh, 2) + 2. * pow(sh, 4)))) *
2685 res += (-0.6372090343808106 - 30.841194974865694 * (1. - sh) -
2686 45.8723194923496 * pow(1. - sh, 2) -
2687 55.64666652217136 * pow(1. - sh, 3) +
2688 log(z)*(-29.333333333333332 -
2689 29.037037037037038 * (1. - sh) -
2690 28.921810699588477 * pow(1. - sh, 2) -
2691 28.854320987654322 * pow(1. - sh, 3) -
2692 28.8079012345679 * pow(1. - sh, 4)) +
2693 log(muh)*(-7.111111111111111 -
2694 7.111111111111111 * (1. - sh) -
2695 7.111111111111111 * pow(1. - sh, 2) -
2696 7.111111111111111 * pow(1. - sh, 3) -
2697 7.111111111111111 * pow(1. - sh, 4)) -
2698 62.91148562558235 * pow(1. - sh, 4)) * pow(z, 1);
2703 res += 0.3950617283950617 * (-2. + sh) * pow(3.141592653589793, 2) *
2704 pow(sh, -1) * pow(z, 1.5);
2706 res += (-3.8991029732698697 - 7.798205946539739 * (1. - sh) -
2707 7.798205946539739 * pow(1. - sh, 2) -
2708 7.798205946539739 * pow(1. - sh, 3) -
2709 7.798205946539739 * pow(1. - sh, 4)) * pow(z, 1.5);
2714 res += (-0.5925925925925926 * Li3shb * (6. + sh) * pow(-1 + sh, -3) -
2715 0.5925925925925926 * Li3sh * (-11. + 3. * sh) *
2716 pow(-1 + sh, -3) - 0.09876543209876543 * (2. + sh) *
2717 pow(Lsh, 3) * pow(-1 + sh, -3) -
2718 lmu * (3.5555555555555554 * pow(sh, -2) +
2719 7.111111111111111 * Lsh * pow(sh, -2)) +
2720 0.5925925925925926 * Li2sh * pow(-1 + sh, -3) * pow(sh, -2)*
2721 (-3. + 10. * sh - 9. * pow(sh, 2) + pow(sh, 3)) +
2722 0.5925925925925926 * cl2 * sqrt4sh * sqrtsh * pow(-1 + sh, -3) *
2723 pow(sh, -1)*(2. + 8. * pow(sh, 2) - 5. * pow(sh, 3) +
2724 pow(sh, 4)) - 0.5925925925925926 * ash * sqrt4sh * sqrtsh *
2725 pow(-1 + sh, -3) * pow(sh, -1)*
2726 (-2. - 4. * sh + 24. * pow(sh, 2) - 15. * pow(sh, 3) +
2727 3. * pow(sh, 4)) - 0.7123300166871669 * pow(-1 + sh, -4) *
2728 pow(sh, -2)*(64. - 256. * sh + 375. * pow(sh, 2) -
2729 236. * pow(sh, 3) + 59. * pow(sh, 4)) +
2730 0.14814814814814814 * pow(Lsh, 2) * pow(-1 + sh, -4) *
2731 pow(sh, -2)*(185. - 751. * sh +
2732 3. * (379. + 2. * Lshb) * pow(sh, 2) -
2733 (777. + 8. * Lshb) * pow(sh, 3) +
2734 2. * (100. + Lshb) * pow(sh, 4)) -
2735 Lsh * (-0.5925925925925926 * Li2sh * (-7. + 2. * sh) *
2737 3.5555555555555554 * pow(ash, 2) * pow(-1 + sh, -4)*
2738 (-1. - 3. * sh + pow(sh, 2)) -
2739 0.04938271604938271 * pow(-1 + sh, -3) * pow(sh, -2)*
2740 (-423. + 1236. * sh -
2741 (1179. + 46. * pow(3.141592653589793, 2)) *
2743 (225. + pow(3.141592653589793, 2)) * pow(sh, 3) +
2744 12. * Lshb * (-3. + 10. * sh - 9. * pow(sh, 2) +
2746 0.5925925925925926 * ash * sqrt4sh * sqrtsh *
2747 pow(-1 + sh, -3) * pow(sh, -1)*
2748 (2. + 8. * pow(sh, 2) - 5. * pow(sh, 3) + pow(sh, 4))) -
2749 1.7777777777777777 * pow(ash, 2) * pow(-1 + sh, -4)*
2750 (4. - 2. * sh - 24. * pow(sh, 2) + 23. * pow(sh, 3) -
2751 8. * pow(sh, 4) + pow(sh, 5)) -
2752 pow(lz, 2)*(-0.2962962962962963 * Lsh * (6. + sh) *
2754 0.14814814814814814 * pow(-1 + sh, -2) * pow(sh, -2)*
2755 (188. - 378. * sh + 169. * pow(sh, 2) + 12. * pow(sh, 3) -
2756 6. * pow(sh, 4) + pow(sh, 5))) -
2757 0.024691358024691357 * pow(-1 + sh, -4) * pow(sh, -2)*
2758 (88. - 448. * sh - 3. * (-415. + 48. * cl3) * pow(sh, 2) -
2759 2. * (917. + 216. * cl3) * pow(sh, 3) +
2760 (1453. + 144. * cl3) * pow(sh, 4) - 693. * pow(sh, 5) +
2761 210. * pow(sh, 6) - 21. * pow(sh, 7) -
2762 2. * pow(3.141592653589793, 2)*
2763 (-603. + 2441. * sh - 3677. * pow(sh, 2) +
2764 2457. * pow(sh, 3) - 662. * pow(sh, 4) +
2765 52. * pow(sh, 5) - 16. * pow(sh, 6) + 2. * pow(sh, 7))) -
2766 lz * (0.5925925925925926 * Li2sh * (6. + sh) *
2768 0.2962962962962963 * pow(Lsh, 2) * pow(-1 + sh, -2) -
2769 7.111111111111111 * lmu * pow(sh, -2) -
2770 7.111111111111111 * pow(ash, 2) * pow(-1 + sh, -4)*
2771 (-1. - 3. * sh + pow(sh, 2)) +
2772 0.2962962962962963 * Lsh * pow(-1 + sh, -3) * pow(sh, -2)*
2773 (-187. + 566. * sh +
2774 3. * (-191. + 4. * Lshb) * pow(sh, 2) +
2775 2. * (96. + Lshb) * pow(sh, 3)) -
2776 1.1851851851851851 * ash * sqrt4sh * sqrtsh *
2777 pow(-1 + sh, -3) * pow(sh, -1)*
2778 (2. + 8. * pow(sh, 2) - 5. * pow(sh, 3) + pow(sh, 4)) +
2779 0.04938271604938271 * pow(-1 + sh, -4) * pow(sh, -2)*
2780 (432. - 1728. * sh +
2781 (2511. + 8. * pow(3.141592653589793, 2)) *
2783 (735. + 11. * pow(3.141592653589793, 2)) * pow(sh, 3)\
2784 + (99. + 2. * pow(3.141592653589793, 2)) *
2785 pow(sh, 4) + 225. * pow(sh, 5) - 78. * pow(sh, 6) +
2786 9. * pow(sh, 7)))) * pow(z, 2);
2788 res += pow(z, 2)*(-336.0382986459167 -
2789 695.8939576145283 * (1. - sh) -
2790 1040.3241576582361 * pow(1. - sh, 2) -
2791 1365.7546667417848 * pow(1. - sh, 3) -
2792 1671.7824637501126 * pow(1. - sh, 4) +
2793 log(muh)*(-3.5555555555555554 +
2794 7.111111111111111 * pow(1. - sh, 2) +
2795 16.59259259259259 * pow(1. - sh, 3) +
2796 27.85185185185185 * pow(1. - sh, 4)) +
2797 log(z)*(-21.887216272530303 +
2798 9.991951790611958 * (1. - sh) +
2799 71.95973563934483 * pow(1. - sh, 2) +
2800 151.94549827445954 * pow(1. - sh, 3) +
2801 245.72125954180052 * pow(1. - sh, 4) +
2802 log(muh)*(7.111111111111111 +
2803 14.222222222222221 * (1. - sh) +
2804 21.333333333333332 * pow(1. - sh, 2) +
2805 28.444444444444443 * pow(1. - sh, 3) +
2806 35.55555555555556 * pow(1. - sh, 4))) +
2807 (27.65432098765432 + 55.67901234567901 * (1. - sh) +
2808 83.6 * pow(1. - sh, 2) +
2809 111.39753086419753 * pow(1. - sh, 3) +
2810 139.20987654320987 * pow(1. - sh, 4)) * pow(log(z), 2));
2815 res += -0.05267489711934156 * pow(3.141592653589793, 2) * pow(sh, -3)*
2816 (3. + 14. * sh + 3. * pow(sh, 2)) * pow(z, 2.5);
2818 res += (-10.397607928719653 - 20.795215857439306 * (1. - sh) -
2819 32.75246497546691 * pow(1. - sh, 2) -
2820 46.26935528280246 * pow(1. - sh, 3) -
2821 61.34588677944595 * pow(1. - sh, 4)) * pow(z, 2.5);
2826 res += (-0.7901234567901234 * Li3shb * (2. + sh) * pow(-1 + sh, -4) -
2827 lmu * (-6.320987654320987 * pow(sh, -3) +
2828 9.481481481481481 * Lsh * pow(sh, -3)) +
2829 1.1851851851851851 * Li3sh * pow(-1 + sh, -3) * pow(sh, -1) -
2830 0.3950617283950617 * pow(Lsh, 3) * pow(-1 + sh, -3) *
2831 pow(sh, -1) + 0.4748866777914446 * pow(-1 + sh, -4) *
2832 pow(sh, -1)*(3. + sh + 2. * pow(sh, 2)) -
2833 0.06584362139917696 * Li2sh * pow(-1 + sh, -4) * pow(sh, -3)*
2834 (-36. + 144. * sh - 234. * pow(sh, 2) + 208. * pow(sh, 3) -
2835 73. * pow(sh, 4) + 9. * pow(sh, 5)) -
2836 0.01646090534979424 * pow(Lsh, 2) * pow(-1 + sh, -5) *
2837 pow(sh, -4)*(-12. + 3682. * sh - 18347. * pow(sh, 2) +
2838 (36889. + 36. * Lshb) * pow(sh, 3) -
2839 2. * (18563. + 36. * Lshb) * pow(sh, 4) +
2840 4. * (4679. + 9. * Lshb) * pow(sh, 5) - 3748. * pow(sh, 6) -
2842 0.3950617283950617 * pow(ash, 2) * pow(-1 + sh, -4)*
2843 (16. - 22. * sh - 171. * pow(sh, 2) + 336. * pow(sh, 3) -
2844 300. * pow(sh, 4) + 141. * pow(sh, 5) - 33. * pow(sh, 6) +
2845 3. * pow(sh, 7)) - pow(lz, 2)*
2846 (-0.3950617283950617 * Lsh * (2. + sh) * pow(-1 + sh, -4) +
2847 0.01646090534979424 * pow(-1 + sh, -3) * pow(sh, -4)*
2848 (-3. + 3685. * sh - 11046. * pow(sh, 2) +
2849 11061. * pow(sh, 3) - 3604. * pow(sh, 4) -
2850 132. * pow(sh, 5) + 321. * pow(sh, 6) -
2851 378. * pow(sh, 7) + 222. * pow(sh, 8) -
2852 60. * pow(sh, 9) + 6. * pow(sh, 10))) -
2853 0.13168724279835392 * ash * pow(-1 + sh, -3)*
2854 (-8. - 21. * sh - 17. * pow(sh, 2) + 23. * pow(sh, 3) -
2855 8. * pow(sh, 4) + pow(sh, 5)) * pow(-((-4 + sh) * sh), 0.5)\
2856 - Lsh * (-2.3703703703703702 * (2. + sh) * pow(ash, 2) *
2858 1.1851851851851851 * Li2sh * pow(-1 + sh, -3) *
2859 pow(sh, -1) + 0.0054869684499314125 * pow(-1 + sh, -4) *
2860 pow(sh, -4)*(30. + (12452. - 432. * Lshb) * sh +
2861 2. * (-25315. + 864. * Lshb) * pow(sh, 2) +
2862 36. * (-78. * Lshb +
2863 5. * (424. + pow(3.141592653589793, 2))) *
2864 pow(sh, 3) + (-50317. + 2496. * Lshb -
2865 36. * pow(3.141592653589793, 2)) * pow(sh, 4) +
2866 (12496. - 876. * Lshb +
2867 72. * pow(3.141592653589793, 2)) * pow(sh, 5) +
2868 27. * (3. + 4. * Lshb) * pow(sh, 6)) +
2869 0.3950617283950617 * ash * pow(-1 + sh, -3)*
2870 (4. + 15. * sh - 29. * pow(sh, 2) + 23. * pow(sh, 3) -
2871 8. * pow(sh, 4) + pow(sh, 5)) *
2872 pow(-((-4 + sh) * sh), 0.5)) -
2873 lz * (0.7901234567901234 * Li2sh * (2. + sh) *
2875 4.7407407407407405 * (2. + sh) * pow(ash, 2) *
2877 9.481481481481481 * lmu * pow(sh, -3) -
2878 0.5925925925925926 * pow(Lsh, 2) * pow(-1 + sh, -3) *
2879 pow(sh, -1) + 0.03292181069958848 * Lsh *
2880 pow(-1 + sh, -4) * pow(sh, -4)*
2881 (-6. + 3682. * sh - 14713. * pow(sh, 2) +
2882 22170. * pow(sh, 3) +
2883 6. * (-2491. + 8. * Lshb) * pow(sh, 4) +
2884 6. * (631. + 4. * Lshb) * pow(sh, 5) - 9. * pow(sh, 6)) +
2885 0.0054869684499314125 * pow(-1 + sh, -3) * pow(sh, -4)*
2886 (15. + 12482. * sh - 37473. * pow(sh, 2) +
2887 37278. * pow(sh, 3) - 12305. * pow(sh, 4) +
2888 498. * pow(sh, 5) - 1215. * pow(sh, 6) +
2889 606. * pow(sh, 7) - 84. * pow(sh, 8) - 24. * pow(sh, 9) +
2891 0.7901234567901234 * ash * pow(-1 + sh, -3)*
2892 (4. + 15. * sh - 29. * pow(sh, 2) + 23. * pow(sh, 3) -
2893 8. * pow(sh, 4) + pow(sh, 5)) *
2894 pow(-((-4 + sh) * sh), 0.5)) -
2895 0.0027434842249657062 * pow(-1 + sh, -4) * pow(sh, -4)*
2896 (19. + 31332. * sh - 125971. * pow(sh, 2) +
2897 189649. * pow(sh, 3) - 94. * pow(sh, 11) +
2898 6. * pow(3.141592653589793, 2)*
2899 (-6. + 3838. * sh - 15355. * pow(sh, 2) +
2900 23258. * pow(sh, 3) - 15790. * pow(sh, 4) +
2901 4134. * pow(sh, 5) - 299. * pow(sh, 6) +
2902 440. * pow(sh, 7) - 400. * pow(sh, 8) +
2903 188. * pow(sh, 9) - 44. * pow(sh, 10) + 4. * pow(sh, 11))\
2904 + pow(sh, 4)*(-125525. + 1728. * cl3 -
2905 576. * cl2 * pow(-((-4 + sh) * sh), 0.5)) +
2906 12. * pow(sh, 5)*(2563. + 72. * cl3 -
2907 132. * cl2 * pow(-((-4 + sh) * sh), 0.5)) +
2908 2. * pow(sh, 10)*(451. +
2909 72. * cl2 * pow(-((-4 + sh) * sh), 0.5)) -
2910 4. * pow(sh, 9)*(791. +
2911 324. * cl2 * pow(-((-4 + sh) * sh), 0.5)) +
2913 4464. * cl2 * pow(-((-4 + sh) * sh), 0.5)) +
2915 6336. * cl2 * pow(-((-4 + sh) * sh), 0.5)) -
2917 7488. * cl2 * pow(-((-4 + sh) * sh), 0.5)))) * pow(z, 3);
2919 res += pow(z, 3)*(-691.9182670329545 -
2920 2017.613655726338 * (1. - sh) -
2921 3881.550454689942 * pow(1. - sh, 2) -
2922 6209.7474017909 * pow(1. - sh, 3) -
2923 8925.087722919276 * pow(1. - sh, 4) +
2924 log(muh)*(6.320987654320987 +
2925 28.444444444444443 * (1. - sh) +
2926 71.11111111111111 * pow(1. - sh, 2) +
2927 137.4814814814815 * pow(1. - sh, 3) +
2928 229.92592592592592 * pow(1. - sh, 4)) +
2929 log(z)*(65.87937511351353 +
2930 322.6369826409675 * (1. - sh) +
2931 831.0074785910778 * pow(1. - sh, 2) +
2932 1624.0261584651576 * pow(1. - sh, 3) +
2933 2734.9197954441984 * pow(1. - sh, 4) +
2934 log(muh)*(9.481481481481481 +
2935 28.444444444444443 * (1. - sh) +
2936 56.888888888888886 * pow(1. - sh, 2) +
2937 94.81481481481481 * pow(1. - sh, 3) +
2938 142.22222222222223 * pow(1. - sh, 4))) +
2939 (59.7037037037037 + 181.14567901234568 * (1. - sh) +
2940 362.4493827160494 * pow(1. - sh, 2) +
2941 603.5673133450912 * pow(1. - sh, 3) +
2942 905.1428571428571 * pow(1. - sh, 4)) * pow(log(z), 2));
2947 res += -0.004514991181657848 * pow(3.141592653589793, 2) * pow(sh, -5)*
2948 (15. + 108. * sh + 314. * pow(sh, 2) + 108. * pow(sh, 3) +
2949 15. * pow(sh, 4)) * pow(z, 3.5);
2951 res += (-24.954259028927165 - 74.8627770867815 * (1. - sh) -
2952 157.21183188224114 * pow(1. - sh, 2) -
2953 279.48770112398427 * pow(1. - sh, 3) -
2954 449.8450801732495 * pow(1. - sh, 4)) * pow(z, 3.5);
2959 res += (-1. * lmu * (-19.555555555555557 * pow(sh, -4) +
2960 21.333333333333332 * Lsh * pow(sh, -4)) -
2961 2.374433388957223 * pow(-1 + sh, -3) * pow(sh, -2) +
2962 1.9753086419753085 * Li3sh * pow(-1 + sh, -3) * pow(sh, -2) -
2963 0.6584362139917695 * pow(Lsh, 3) * pow(-1 + sh, -3) *
2964 pow(sh, -2) - 0.09876543209876543 * Li2sh *
2965 pow(-1 + sh, -5) * pow(sh, -4)*
2966 (54. - 268. * sh + 544. * pow(sh, 2) - 570. * pow(sh, 3) +
2967 292. * pow(sh, 4) - 58. * pow(sh, 5) - 6. * pow(sh, 6) +
2968 3. * pow(sh, 7)) - 1.7777777777777777 * sh * pow(ash, 2) *
2969 pow(-1 + sh, -4)*(-24. - 3. * sh + 174. * pow(sh, 2) -
2970 390. * pow(sh, 3) + 411. * pow(sh, 4) - 241. * pow(sh, 5) +
2971 80. * pow(sh, 6) - 14. * pow(sh, 7) + pow(sh, 8)) -
2972 0.01646090534979424 * pow(Lsh, 2) * pow(-1 + sh, -6) *
2973 pow(sh, -6)*(6. + 14. * sh - 9747. * pow(sh, 2) +
2974 58077. * pow(sh, 3) -
2975 3. * (48391. + 20. * Lshb) * pow(sh, 4) +
2976 3. * (64669. + 60. * Lshb) * pow(sh, 5) -
2977 15. * (9737. + 12. * Lshb) * pow(sh, 6) +
2978 (58789. + 60. * Lshb) * pow(sh, 7) - 9906. * pow(sh, 8) +
2979 27. * pow(sh, 9) - 9. * pow(sh, 10)) +
2980 0.00823045267489712 * pow(lz, 2) * pow(sh, -6)*
2981 (-3. - 25. * sh + 19344. * pow(sh, 2) + 73. * pow(sh, 4) +
2982 3. * pow(sh, 5) + 180. * pow(sh, 6) + 36. * pow(sh, 7) -
2983 738. * pow(sh, 8) + 612. * pow(sh, 9) - 180. * pow(sh, 10) +
2984 18. * pow(sh, 11)) +
2985 0.09876543209876543 * ash * pow(-1 + sh, -3)*
2986 (-6. - 41. * sh - 239. * pow(sh, 2) + 671. * pow(sh, 3) -
2987 680. * pow(sh, 4) + 329. * pow(sh, 5) - 77. * pow(sh, 6) +
2988 7. * pow(sh, 7)) * pow(-((-4 + sh) * sh), 0.5) -
2989 lz * (-21.333333333333332 * lmu * pow(sh, -4) -
2990 0.9876543209876543 * pow(Lsh, 2) * pow(-1 + sh, -3) *
2991 pow(sh, -2) - 0.03292181069958848 * Lsh *
2992 pow(-1 + sh, -5) * pow(sh, -6)*
2993 (-3. - 10. * sh + 9722. * pow(sh, 2) -
2994 48430. * pow(sh, 3) + 96898. * pow(sh, 4) -
2995 97142. * pow(sh, 5) + 48736. * pow(sh, 6) -
2996 9804. * pow(sh, 7) + 3. * pow(sh, 8) + 3. * pow(sh, 9)) -
2997 0.00027434842249657066 * pow(-1 + sh, -4) * pow(sh, -6)*
2998 (105. + 55. * sh + 963530. * pow(sh, 2) -
2999 3.85233e6 * pow(sh, 3) + 5.76189e6 * pow(sh, 4) -
3000 3.806202e6 * pow(sh, 5) + 923844. * pow(sh, 6) +
3001 1154. * pow(sh, 7) + 54349. * pow(sh, 8) -
3002 76131. * pow(sh, 9) + 4716. * pow(sh, 10) +
3003 82440. * pow(sh, 11) - 83970. * pow(sh, 12) +
3004 36900. * pow(sh, 13) - 7740. * pow(sh, 14) +
3005 630. * pow(sh, 15)) -
3006 1.1851851851851851 * ash *
3007 (-6. - 13. * sh + 20. * pow(sh, 2) - 8. * pow(sh, 3) +
3008 pow(sh, 4)) * pow(-((-4 + sh) * sh), 0.5)) -
3009 Lsh * (1.9753086419753085 * Li2sh * pow(-1 + sh, -3) *
3010 pow(sh, -2) + 0.0027434842249657062 * pow(-1 + sh, -5) *
3011 pow(sh, -6)*(-21. + 10. * sh +
3012 2. * (-48097. + 972. * Lshb) * pow(sh, 2) +
3013 (483970. - 9648. * Lshb) * pow(sh, 3) +
3014 (-972667. + 19584. * Lshb -
3015 600. * pow(3.141592653589793, 2)) * pow(sh, 4) +
3016 8. * (122107. - 2565. * Lshb +
3017 150. * pow(3.141592653589793, 2)) * pow(sh, 5) +
3018 2. * (-245429. + 5256. * Lshb -
3019 300. * pow(3.141592653589793, 2)) * pow(sh, 6) +
3020 (98142. - 2088. * Lshb) * pow(sh, 7) -
3021 6. * (7. + 36. * Lshb) * pow(sh, 8) +
3022 6. * (5. + 18. * Lshb) * pow(sh, 9)) -
3023 0.5925925925925926 * ash *
3024 (-6. - 13. * sh + 20. * pow(sh, 2) - 8. * pow(sh, 3) +
3025 pow(sh, 4)) * pow(-((-4 + sh) * sh), 0.5)) -
3026 4.572473708276177e-6 * pow(-1 + sh, -5) * pow(sh, -6)*
3027 (-1200. * pow(3.141592653589793, 2)*
3028 (-9. - 30. * sh + 30192. * pow(sh, 2) -
3029 150534. * pow(sh, 3) + 301774. * pow(sh, 4) -
3030 303542. * pow(sh, 5) + 152986. * pow(sh, 6) -
3031 30076. * pow(sh, 7) - 1033. * pow(sh, 8) -
3032 2857. * pow(sh, 9) + 10122. * pow(sh, 10) -
3033 14418. * pow(sh, 11) + 11736. * pow(sh, 12) -
3034 5778. * pow(sh, 13) + 1692. * pow(sh, 14) -
3035 270. * pow(sh, 15) + 18. * pow(sh, 16)) +
3036 (-1. + sh)*(2775. - 2775. * sh +
3037 1.836255e7 * pow(sh, 2) - 7.415055e7 * pow(sh, 3) +
3038 1.1249425e8 * pow(sh, 4) - 7.6337698e7 * pow(sh, 5) +
3039 86850. * pow(sh, 15) +
3040 pow(sh, 8)*(4.883011e6 -
3041 4.6656e6 * cl2 * pow(-((-4 + sh) * sh), 0.5)) -
3043 (1249. + 144. * cl2 * pow(-((-4 + sh) * sh), 0.5)) +
3045 (6445. + 1728. * cl2 * pow(-((-4 + sh) * sh), 0.5)) -
3046 4050. * pow(sh, 12)*
3047 (3755. + 1856. * cl2 * pow(-((-4 + sh) * sh), 0.5)) +
3048 4500. * pow(sh, 11)*
3049 (4681. + 4176. * cl2 * pow(-((-4 + sh) * sh), 0.5)) +
3052 64800. * cl2 * pow(-((-4 + sh) * sh), 0.5)) -
3053 2. * pow(sh, 7)*(1.393817e6 +
3054 712800. * cl2 * pow(-((-4 + sh) * sh), 0.5)) -
3056 (374591. + 716400. * cl2 * pow(-((-4 + sh) * sh), 0.5))
3058 (-523. + 6.1344e6 * cl2 * pow(-((-4 + sh) * sh), 0.5)))
3061 res += pow(z, 4)*(-1672.4687901566804 -
3062 6425.515582048289 * (1. - sh) -
3063 15411.783043279445 * pow(1. - sh, 2) -
3064 29568.2502704632 * pow(1. - sh, 3) -
3065 49534.42028818623 * pow(1. - sh, 4) +
3066 log(muh)*(19.555555555555557 +
3067 99.55555555555556 * (1. - sh) +
3068 291.55555555555554 * pow(1. - sh, 2) +
3069 654.2222222222222 * pow(1. - sh, 3) +
3070 1251.5555555555557 * pow(1. - sh, 4)) +
3071 log(z)*(253.5498382001439 +
3072 1369.6655601827902 * (1. - sh) +
3073 4072.304539664025 * pow(1. - sh, 2) +
3074 9174.433775454903 * pow(1. - sh, 3) +
3075 17618.760992657026 * pow(1. - sh, 4) +
3076 log(muh)*(21.333333333333332 +
3077 85.33333333333333 * (1. - sh) +
3078 213.33333333333334 * pow(1. - sh, 2) +
3079 426.6666666666667 * pow(1. - sh, 3) +
3080 746.6666666666666 * pow(1. - sh, 4))) +
3081 (159.01234567901236 + 638.8148148148148 * (1. - sh) +
3082 1591.9506172839506 * pow(1. - sh, 2) +
3083 3177.448559670782 * pow(1. - sh, 3) +
3084 5557.119341563786 * pow(1. - sh, 4)) * pow(log(z), 2));
3089 res += -0.0010749979003947259 * pow(3.141592653589793, 2) * pow(sh, -7)*
3090 (35. + 330. * sh + 1389. * pow(sh, 2) + 3212. * pow(sh, 3) +
3091 1389. * pow(sh, 4) + 330. * pow(sh, 5) + 35. * pow(sh, 6)) *
3094 res += (-71.2978829397919 - 285.1915317591676 * (1. - sh) -
3095 745.0628767208254 * pow(1. - sh, 2) -
3096 1586.3778954103698 * pow(1. - sh, 3) -
3097 2982.415906901117 * pow(1. - sh, 4)) * pow(z, 4.5);
3102 res += (-1. * lmu * (-58.785185185185185 * pow(sh, -5) +
3103 56.888888888888886 * Lsh * pow(sh, -5)) -
3104 4.986310116810168 * pow(-1 + sh, -3) * pow(sh, -3) +
3105 4.148148148148148 * Li3sh * pow(-1 + sh, -3) * pow(sh, -3) -
3106 1.382716049382716 * pow(Lsh, 3) * pow(-1 + sh, -3) *
3107 pow(sh, -3) - 0.3950617283950617 * cl2 * sh * sqrt4sh * sqrtsh *
3108 (44. + 108. * sh - 221. * pow(sh, 2) + 132. * pow(sh, 3) -
3109 33. * pow(sh, 4) + 3. * pow(sh, 5)) -
3110 0.03950617283950617 * ash * sh * sqrt4sh * sqrtsh *
3111 pow(-1 + sh, -3)*(-496. + 326. * sh + 6147. * pow(sh, 2) -
3112 16495. * pow(sh, 3) + 18680. * pow(sh, 4) -
3113 11251. * pow(sh, 5) + 3760. * pow(sh, 6) -
3114 658. * pow(sh, 7) + 47. * pow(sh, 8)) -
3115 0.19753086419753085 * Li2sh * pow(-1 + sh, -6) * pow(sh, -5)*
3116 (-72. + 429. * sh - 1073. * pow(sh, 2) + 1442. * pow(sh, 3) -
3117 1092. * pow(sh, 4) + 435. * pow(sh, 5) - 47. * pow(sh, 6) -
3118 6. * pow(sh, 7) - pow(sh, 8) + pow(sh, 9)) +
3119 1.1851851851851851 * sh * pow(ash, 2) * pow(-1 + sh, -4)*
3120 (-24. + 162. * sh + 114. * pow(sh, 2) - 1809. * pow(sh, 3) +
3121 4224. * pow(sh, 4) - 4957. * pow(sh, 5) +
3122 3423. * pow(sh, 6) - 1445. * pow(sh, 7) +
3123 366. * pow(sh, 8) - 51. * pow(sh, 9) + 3. * pow(sh, 10)) -
3124 Lsh * (4.148148148148148 * Li2sh * pow(-1 + sh, -3) *
3125 pow(sh, -3) + 0.3950617283950617 * ash * sh * sqrt4sh *
3126 sqrtsh * (44. + 108. * sh - 221. * pow(sh, 2) +
3127 132. * pow(sh, 3) - 33. * pow(sh, 4) + 3. * pow(sh, 5)) +
3128 0.00005486968449931413 * pow(-1 + sh, -6) * pow(sh, -8)*
3129 (486. - 108. * sh - 8742. * pow(sh, 2) +
3130 (1.6072738e7 - 259200. * Lshb) * pow(sh, 3) +
3131 6. * (-1.6103209e7 + 257400. * Lshb) * pow(sh, 4) +
3132 24. * (1.0092821e7 - 160950. * Lshb +
3133 2625. * pow(3.141592653589793, 2)) * pow(sh, 5) +
3134 4. * (-8.0991623e7 + 1.2978e6 * Lshb -
3135 47250. * pow(3.141592653589793, 2)) * pow(sh, 6) +
3136 6. * (4.0636147e7 - 655200. * Lshb +
3137 31500. * pow(3.141592653589793, 2)) * pow(sh, 7) +
3138 3. * (-3.2644121e7 + 522000. * Lshb -
3139 21000. * pow(3.141592653589793, 2)) * pow(sh, 8) -
3140 2. * (-8.202737e6 + 84600. * Lshb) * pow(sh, 9) -
3141 900. * (17. + 24. * Lshb) * pow(sh, 10) -
3142 450. * (-1. + 8. * Lshb) * pow(sh, 11) +
3143 75. * (7. + 48. * Lshb) * pow(sh, 12))) -
3144 0.0016460905349794238 * pow(Lsh, 2) * pow(-1 + sh, -7) *
3145 pow(sh, -8)*(-36. - 126. * sh + 84. * pow(sh, 2) +
3146 299648. * pow(sh, 3) - 2.087099e6 * pow(sh, 4) +
3147 60. * (104197. + 21. * Lshb) * pow(sh, 5) -
3148 2. * (5.211467e6 + 2520. * Lshb) * pow(sh, 6) +
3149 56. * (186416. + 135. * Lshb) * pow(sh, 7) -
3150 12. * (523459. + 420. * Lshb) * pow(sh, 8) +
3151 (2.104169e6 + 1260. * Lshb) * pow(sh, 9) -
3152 303404. * pow(sh, 10) + 300. * pow(sh, 11) +
3153 120. * pow(sh, 12) - 60. * pow(sh, 13)) +
3154 0.0008230452674897119 * pow(lz, 2) * pow(sh, -8)*
3155 (-18. - 189. * sh - 903. * pow(sh, 2) +
3156 594928. * pow(sh, 3) + 630. * pow(sh, 4) +
3157 1383. * pow(sh, 5) + 1839. * pow(sh, 6) + 18. * pow(sh, 7) +
3158 3120. * pow(sh, 8) + 8160. * pow(sh, 9) +
3159 1440. * pow(sh, 10) - 46560. * pow(sh, 11) +
3160 51720. * pow(sh, 12) - 23040. * pow(sh, 13) +
3161 4680. * pow(sh, 14) - 360. * pow(sh, 15)) -
3162 lz * (-56.888888888888886 * lmu * pow(sh, -5) -
3163 2.074074074074074 * pow(Lsh, 2) * pow(-1 + sh, -3) *
3164 pow(sh, -3) + 0.7901234567901234 * ash * sh * sqrt4sh *
3165 sqrtsh * (44. + 108. * sh - 221. * pow(sh, 2) +
3166 132. * pow(sh, 3) - 33. * pow(sh, 4) + 3. * pow(sh, 5)) -
3167 0.0032921810699588477 * Lsh * pow(-1 + sh, -6) *
3168 pow(sh, -8)*(18. + 81. * sh + 39. * pow(sh, 2) -
3169 299147. * pow(sh, 3) + 1.788834e6 * pow(sh, 4) -
3170 4.467993e6 * pow(sh, 5) + 5.959219e6 * pow(sh, 6) -
3171 4.474299e6 * pow(sh, 7) + 1.792497e6 * pow(sh, 8) -
3172 297944. * pow(sh, 9) - 450. * pow(sh, 10) +
3173 90. * pow(sh, 11) + 15. * pow(sh, 12)) +
3174 7.838526357044876e-6 * pow(-1 + sh, -5) * pow(sh, -8)*
3175 (1701. + 1323. * sh - 29274. * pow(sh, 2) +
3176 1.12541443e8 * pow(sh, 3) -
3177 5.61871835e8 * pow(sh, 4) +
3178 1.121615551e9 * pow(sh, 5) -
3179 1.117743382e9 * pow(sh, 6) +
3180 5.54683001e8 * pow(sh, 7) -
3181 1.07408773e8 * pow(sh, 8) - 4.03983e6 * pow(sh, 9) +
3182 4.371768e6 * pow(sh, 10) + 4.551324e6 * pow(sh, 11) -
3183 1.1897697e7 * pow(sh, 12) - 1.323264e7 * pow(sh, 13) +
3184 6.277908e7 * pow(sh, 14) - 8.508528e7 * pow(sh, 15) +
3185 6.212682e7 * pow(sh, 16) - 2.701314e7 * pow(sh, 17) +
3186 6.98418e6 * pow(sh, 18) - 990360. * pow(sh, 19) +
3187 59220. * pow(sh, 20))) -
3188 2.7994736989445985e-8 * pow(-1 + sh, -6) * pow(sh, -8)*
3189 (161406. - 318843. * sh - 1.910657e6 * pow(sh, 2) +
3190 2.6587841e7 * pow(sh, 3) - 4.35567027e8 * pow(sh, 4) +
3191 2.014488392e9 * pow(sh, 5) - 4.239785466e9 * pow(sh, 6) +
3192 4.472946066e9 * pow(sh, 7) - 1.480821324e9 * pow(sh, 8) -
3193 2.199797883e9 * pow(sh, 9) +
3194 2.944201467e9 * pow(sh, 10) +
3195 1.246759469e9 * pow(sh, 11) -
3196 2.466843039e9 * pow(sh, 12) -
3197 1.8171607842e10 * pow(sh, 13) +
3198 6.053910752e10 * pow(sh, 14) -
3199 8.996385496e10 * pow(sh, 15) +
3200 7.96292826e10 * pow(sh, 16) -
3201 4.509057616e10 * pow(sh, 17) +
3202 1.647058952e10 * pow(sh, 18) -
3203 3.74875284e9 * pow(sh, 19) + 4.8245988e8 * pow(sh, 20) -
3204 2.674812e7 * pow(sh, 21) +
3205 58800. * pow(3.141592653589793, 2)*
3206 (-18. - 81. * sh - 39. * pow(sh, 2) +
3207 308147. * pow(sh, 3) - 1.843944e6 * pow(sh, 4) +
3208 4.610871e6 * pow(sh, 5) - 6.160013e6 * pow(sh, 6) +
3209 4.635393e6 * pow(sh, 7) - 1.862435e6 * pow(sh, 8) +
3210 312614. * pow(sh, 9) + 3132. * pow(sh, 10) -
3211 3736. * pow(sh, 11) - 77479. * pow(sh, 12) +
3212 318588. * pow(sh, 13) - 608560. * pow(sh, 14) +
3213 702440. * pow(sh, 15) - 529920. * pow(sh, 16) +
3214 267160. * pow(sh, 17) - 89120. * pow(sh, 18) +
3215 18840. * pow(sh, 19) - 2280. * pow(sh, 20) +
3216 120. * pow(sh, 21)))) * pow(z, 5);
3218 res += pow(z, 5)*(-4875.039604802492 -
3219 23394.026843086856 * (1. - sh) -
3220 67666.42652379921 * pow(1. - sh, 2) -
3221 152266.3652960926 * pow(1. - sh, 3) -
3222 292874.8144929396 * pow(1. - sh, 4) +
3223 log(muh)*(58.785185185185185 +
3224 350.81481481481484 * (1. - sh) +
3225 1194.6666666666667 * pow(1. - sh, 2) +
3226 3072. * pow(1. - sh, 3) +
3227 6641.777777777777 * pow(1. - sh, 4)) +
3228 log(z)*(847.677812246958 +
3229 5397.215321438122 * (1. - sh) +
3230 18627.480751808722 * pow(1. - sh, 2) +
3231 48121.18079162945 * pow(1. - sh, 3) +
3232 104649.22541796633 * pow(1. - sh, 4) +
3233 log(muh)*(56.888888888888886 +
3234 284.44444444444446 * (1. - sh) +
3235 853.3333333333334 * pow(1. - sh, 2) +
3236 1991.111111111111 * pow(1. - sh, 3) +
3237 3982.222222222222 * pow(1. - sh, 4))) +
3238 (491.23292181069957 +
3239 2460.5102880658437 * (1. - sh) +
3240 7344.419753086419 * pow(1. - sh, 2) +
3241 17100.213991769546 * pow(1. - sh, 3) +
3242 34182.45267489712 * pow(1. - sh, 4)) * pow(log(z), 2));
3247 res += -0.00007600995255316243 * pow(3.141592653589793, 2) * pow(sh, -9)*
3248 (315. + 3640. * sh + 19380. * pow(sh, 2) + 62088. * pow(sh, 3) +
3249 124834. * pow(sh, 4) + 62088. * pow(sh, 5) +
3250 19380. * pow(sh, 6) + 3640. * pow(sh, 7) + 315. * pow(sh, 8)) *
3253 res += (-221.8156358126859 - 1109.0781790634296 * (1. - sh) -
3254 3460.3239186779006 * pow(1. - sh, 2) -
3255 8562.083542369677 * pow(1. - sh, 3) -
3256 18357.620459597754 * pow(1. - sh, 4)) * pow(z, 5.5);
3261 res += (-1. * lmu * (-182.91358024691357 * pow(sh, -6) +
3262 165.92592592592592 * Lsh * pow(sh, -6)) -
3263 11.967144280344405 * pow(-1 + sh, -3) * pow(sh, -4) +
3264 9.955555555555556 * Li3sh * pow(-1 + sh, -3) * pow(sh, -4) -
3265 3.3185185185185184 * pow(Lsh, 3) * pow(-1 + sh, -3) *
3266 pow(sh, -4) + 0.4609053497942387 * cl2 * sh * sqrt4sh * sqrtsh *
3267 (20. - 194. * sh - 657. * pow(sh, 2) + 1581. * pow(sh, 3) -
3268 1239. * pow(sh, 4) + 462. * pow(sh, 5) - 84. * pow(sh, 6) +
3269 6. * pow(sh, 7)) + 0.0043895747599451305 * ash * sh *
3270 sqrt4sh * sqrtsh * pow(-1 + sh, -3)*
3271 (-3230. + 42428. * sh + 5415. * pow(sh, 2) -
3272 547510. * pow(sh, 3) + 1.454143e6 * pow(sh, 4) -
3273 1.806873e6 * pow(sh, 5) + 1.280646e6 * pow(sh, 6) -
3274 546171. * pow(sh, 7) + 138714. * pow(sh, 8) -
3275 19329. * pow(sh, 9) + 1137. * pow(sh, 10)) -
3276 1.382716049382716 * pow(ash, 2) * pow(-1 + sh, -4) * pow(sh, 2)*
3277 (210. - 870. * sh - 1395. * pow(sh, 2) +
3278 13113. * pow(sh, 3) - 31009. * pow(sh, 4) +
3279 39370. * pow(sh, 5) - 30843. * pow(sh, 6) +
3280 15585. * pow(sh, 7) - 5091. * pow(sh, 8) +
3281 1038. * pow(sh, 9) - 120. * pow(sh, 10) + 6. * pow(sh, 11))\
3282 - 0.0032921810699588477 * Li2sh * pow(-1 + sh, -7) *
3283 pow(sh, -6)*(12600. - 87696. * sh + 262296. * pow(sh, 2) -
3284 436590. * pow(sh, 3) + 435240. * pow(sh, 4) -
3285 257040. * pow(sh, 5) + 79436. * pow(sh, 6) -
3286 8456. * pow(sh, 7) - 360. * pow(sh, 9) + 45. * pow(sh, 11))\
3287 - Lsh * (9.955555555555556 * Li2sh * pow(-1 + sh, -3) *
3288 pow(sh, -4) - 0.4609053497942387 * ash * sh * sqrt4sh *
3289 sqrtsh * (20. - 194. * sh - 657. * pow(sh, 2) +
3290 1581. * pow(sh, 3) - 1239. * pow(sh, 4) +
3291 462. * pow(sh, 5) - 84. * pow(sh, 6) + 6. * pow(sh, 7)) +
3292 0.000010973936899862826 * pow(-1 + sh, -7) * pow(sh, -10)*
3293 (-1320. + 96. * sh + 30636. * pow(sh, 2) +
3294 70896. * pow(sh, 3) +
3295 8. * (-3.4452983e7 + 472500. * Lshb) * pow(sh, 4) -
3296 112. * (-1.7229823e7 + 234900. * Lshb) * pow(sh, 5) +
3297 6. * (-9.66264881e8 + 1.31148e7 * Lshb -
3298 126000. * pow(3.141592653589793, 2)) * pow(sh, 6) +
3299 (9.681150074e9 - 1.30977e8 * Lshb +
3300 3.024e6 * pow(3.141592653589793, 2)) * pow(sh, 7) +
3301 2. * (-4.851617443e9 + 6.5286e7 * Lshb -
3302 2.268e6 * pow(3.141592653589793, 2)) * pow(sh, 8) +
3303 24. * (2.43237251e8 - 3.213e6 * Lshb +
3304 126000. * pow(3.141592653589793, 2)) * pow(sh, 9) +
3305 3. * (-6.50704647e8 + 7.9436e6 * Lshb -
3306 252000. * pow(3.141592653589793, 2)) * pow(sh, 10)\
3307 - 60. * (-4.666787e6 + 42280. * Lshb) * pow(sh, 11) -
3308 23730. * pow(sh, 12) -
3309 60. * (923. + 1800. * Lshb) * pow(sh, 13) +
3310 4995. * pow(sh, 14) +
3311 135. * (9. + 100. * Lshb) * pow(sh, 15))) -
3312 0.0003292181069958848 * pow(Lsh, 2) * pow(-1 + sh, -8) *
3313 pow(sh, -10)*(120. + 552. * sh + 120. * pow(sh, 2) -
3314 3480. * pow(sh, 3) - 4.90282e6 * pow(sh, 4) +
3315 3.912806e7 * pow(sh, 5) -
3316 (1.36736083e8 + 15120. * Lshb) * pow(sh, 6) +
3317 5. * (5.4670261e7 + 15120. * Lshb) * pow(sh, 7) -
3318 70. * (4.883491e6 + 2160. * Lshb) * pow(sh, 8) +
3319 40. * (6.845753e6 + 3780. * Lshb) * pow(sh, 9) -
3320 10. * (1.3722853e7 + 7560. * Lshb) * pow(sh, 10) +
3321 2. * (1.9677013e7 + 7560. * Lshb) * pow(sh, 11) -
3322 4.94776e6 * pow(sh, 12) - 1800. * pow(sh, 13) +
3323 1800. * pow(sh, 14) + 225. * pow(sh, 15) - 225. * pow(sh, 16)
3324 ) + 0.0003292181069958848 * pow(lz, 2) * pow(sh, -10)*
3325 (-30. - 378. * sh - 2214. * pow(sh, 2) - 7938. * pow(sh, 3) +
3326 4.8772e6 * pow(sh, 4) + 7560. * pow(sh, 5) +
3327 6273. * pow(sh, 6) + 13689. * pow(sh, 7) +
3328 10908. * pow(sh, 8) + 30. * pow(sh, 9) +
3329 20300. * pow(sh, 10) + 27300. * pow(sh, 11) +
3330 113400. * pow(sh, 12) + 49700. * pow(sh, 13) -
3331 880950. * pow(sh, 14) + 1.19385e6 * pow(sh, 15) -
3332 706650. * pow(sh, 16) + 216300. * pow(sh, 17) -
3333 33600. * pow(sh, 18) + 2100. * pow(sh, 19)) -
3334 lz * (-165.92592592592592 * lmu * pow(sh, -6) -
3335 4.977777777777778 * pow(Lsh, 2) * pow(-1 + sh, -3) *
3336 pow(sh, -4) - 0.9218106995884774 * ash * sh * sqrt4sh *
3337 sqrtsh * (20. - 194. * sh - 657. * pow(sh, 2) +
3338 1581. * pow(sh, 3) - 1239. * pow(sh, 4) +
3339 462. * pow(sh, 5) - 84. * pow(sh, 6) + 6. * pow(sh, 7)) -
3340 0.0006584362139917696 * Lsh * pow(-1 + sh, -7) *
3341 pow(sh, -10)*(-60. - 336. * sh - 396. * pow(sh, 2) +
3342 1344. * pow(sh, 3) + 4.900804e6 * pow(sh, 4) -
3343 3.4230616e7 * pow(sh, 5) + 1.02572181e8 * pow(sh, 6) -
3344 1.70881214e8 * pow(sh, 7) +
3345 1.70889596e8 * pow(sh, 8) -
3346 1.02558564e8 * pow(sh, 9) +
3347 3.4176506e7 * pow(sh, 10) - 4.86782e6 * pow(sh, 11) -
3348 2940. * pow(sh, 12) - 1740. * pow(sh, 13) +
3349 585. * pow(sh, 14) + 45. * pow(sh, 15)) +
3350 1.567705271408975e-6 * pow(-1 + sh, -6) * pow(sh, -10)*
3351 (-4620. - 4284. * sh + 102942. * pow(sh, 2) +
3352 351078. * pow(sh, 3) - 1.928933146e9 * pow(sh, 4) +
3353 1.155219891e10 * pow(sh, 5) -
3354 2.8832387178e10 * pow(sh, 6) +
3355 3.8356658804e10 * pow(sh, 7) -
3356 2.8662297144e10 * pow(sh, 8) +
3357 1.1385924988e10 * pow(sh, 9) -
3358 1.88384127e9 * pow(sh, 10) +
3359 5.9056617e7 * pow(sh, 11) -
3360 8.9385341e7 * pow(sh, 12) +
3361 1.01542571e8 * pow(sh, 13) -
3362 3.36952887e8 * pow(sh, 14) +
3363 4.6151941e8 * pow(sh, 15) +
3364 1.1490745e9 * pow(sh, 16) -
3365 5.1849126e9 * pow(sh, 17) +
3366 8.94761595e9 * pow(sh, 18) -
3367 9.1281316e9 * pow(sh, 19) +
3368 6.0729627e9 * pow(sh, 20) -
3369 2.7087795e9 * pow(sh, 21) + 8.039283e8 * pow(sh, 22) -
3370 1.522416e8 * pow(sh, 23) + 1.66278e7 * pow(sh, 24) -
3371 795900. * pow(sh, 25))) -
3372 1.244210532864266e-9 * pow(-1 + sh, -7) * pow(sh, -10)*
3373 (-1.60524e6 + 3.372768e6 * sh + 2.1255318e7 * pow(sh, 2) -
3374 5.1364152e7 * pow(sh, 3) +
3375 3.79257522168e11 * pow(sh, 4) -
3376 2.633399785632e12 * pow(sh, 5) +
3377 7.825640856201e12 * pow(sh, 6) -
3378 1.2900120410724e13 * pow(sh, 7) +
3379 1.2750237623886e13 * pow(sh, 8) -
3380 7.584725099624e12 * pow(sh, 9) +
3381 2.529935410926e12 * pow(sh, 10) -
3382 2.91800208364e11 * pow(sh, 11) -
3383 2.09548894394e11 * pow(sh, 12) +
3384 3.42677229216e11 * pow(sh, 13) -
3385 7.84252942969e11 * pow(sh, 14) +
3386 7.91328639346e11 * pow(sh, 15) +
3387 3.49537340577e12 * pow(sh, 16) -
3388 1.51956042427e13 * pow(sh, 17) +
3389 2.91706462341e13 * pow(sh, 18) -
3390 3.468763704735e13 * pow(sh, 19) +
3391 2.790223133985e13 * pow(sh, 20) -
3392 1.564098014115e13 * pow(sh, 21) +
3393 6.1223554161e12 * pow(sh, 22) -
3394 1.63963647855e12 * pow(sh, 23) +
3395 2.859293178e11 * pow(sh, 24) -
3396 2.92020939e10 * pow(sh, 25) + 1.3226913e9 * pow(sh, 26) -
3397 88200. * pow(3.141592653589793, 2)*
3398 (-180. - 1008. * sh - 1188. * pow(sh, 2) +
3399 4032. * pow(sh, 3) + 1.5093012e7 * pow(sh, 4) -
3400 1.05478968e8 * pow(sh, 5) +
3401 3.16319451e8 * pow(sh, 6) -
3402 5.27535714e8 * pow(sh, 7) +
3403 5.28275796e8 * pow(sh, 8) -
3404 3.17556564e8 * pow(sh, 9) +
3405 1.05971306e8 * pow(sh, 10) -
3406 1.4919188e7 * pow(sh, 11) - 564648. * pow(sh, 12) +
3407 1.515872e6 * pow(sh, 13) - 1.613573e6 * pow(sh, 14) -
3408 1.0017483e7 * pow(sh, 15) + 5.154742e7 * pow(sh, 16) -
3409 1.211917e8 * pow(sh, 17) + 1.785651e8 * pow(sh, 18) -
3410 1.800631e8 * pow(sh, 19) + 1.286215e8 * pow(sh, 20) -
3411 6.57363e7 * pow(sh, 21) + 2.38644e7 * pow(sh, 22) -
3412 5.9997e6 * pow(sh, 23) + 991200. * pow(sh, 24) -
3413 96600. * pow(sh, 25) + 4200. * pow(sh, 26)))) * pow(z, 6);
3415 res += pow(z, 6)*(-15489.003566799385 -
3416 89232.71889975584 * (1. - sh) -
3417 302596.6493265383 * pow(1. - sh, 2) -
3418 781979.762945469 * pow(1. - sh, 3) -
3419 1.6992401790855133e6 * pow(1. - sh, 4) +
3420 log(muh)*(182.91358024691357 +
3421 1263.4074074074074 * (1. - sh) +
3422 4919.7037037037035 * pow(1. - sh, 2) +
3423 14280.691358024691 * pow(1. - sh, 3) +
3424 34454.51851851852 * pow(1. - sh, 4)) +
3425 log(z)*(2906.463936705608 +
3426 21474.23376628666 * (1. - sh) +
3427 84464.70442803907 * pow(1. - sh, 2) +
3428 246235.5920127741 * pow(1. - sh, 3) +
3429 598321.9180693374 * pow(1. - sh, 4) +
3430 log(muh)*(165.92592592592592 +
3431 995.5555555555555 * (1. - sh) +
3432 3484.4444444444443 * pow(1. - sh, 2) +
3433 9291.851851851852 * pow(1. - sh, 3) +
3434 20906.666666666668 * pow(1. - sh, 4))) +
3435 (1615.4238683127571 + 9691.876543209877 * (1. - sh) +
3436 33706.74567901235 * pow(1. - sh, 2) +
3437 89710.5086419753 * pow(1. - sh, 3) +
3438 201778.22222222222 * pow(1. - sh, 4)) * pow(log(z), 2));
3443 res += -0.000023919215838407757 * pow(3.141592653589793, 2) *
3444 pow(sh, -11)*(693. + 9450. * sh + 60305. * pow(sh, 2) +
3445 238200. * pow(sh, 3) + 642234. * pow(sh, 4) +
3446 1.173308e6 * pow(sh, 5) + 642234. * pow(sh, 6) +
3447 238200. * pow(sh, 7) + 60305. * pow(sh, 8) +
3448 9450. * pow(sh, 9) + 693. * pow(sh, 10)) * pow(z, 6.5);
3450 res += (-725.9420808415176 - 4355.652485049105 * (1. - sh) -
3451 15789.240258303007 * pow(1. - sh, 2) -
3452 44463.95245154296 * pow(1. - sh, 3) -
3453 106907.93465535706 * pow(1. - sh, 4)) * pow(z, 6.5);
3458 res += (-1. * lmu * (-587.9873015873015 * pow(sh, -7) +
3459 512. * Lsh * pow(sh, -7)) -
3460 31.342520734235343 * pow(-1 + sh, -3) * pow(sh, -5) +
3461 26.074074074074073 * Li3sh * pow(-1 + sh, -3) * pow(sh, -5) -
3462 8.691358024691358 * pow(Lsh, 3) * pow(-1 + sh, -3) *
3463 pow(sh, -5) - 0.7111111111111111 * cl2 * sqrt4sh * sqrtsh *
3464 pow(sh, 2)*(-204. + 654. * sh + 3621. * pow(sh, 2) -
3465 9390. * pow(sh, 3) + 8928. * pow(sh, 4) -
3466 4384. * pow(sh, 5) + 1190. * pow(sh, 6) -
3467 170. * pow(sh, 7) + 10. * pow(sh, 8)) -
3468 0.011287477954144622 * ash * sqrt4sh * sqrtsh *
3469 pow(-1 + sh, -3) * pow(sh, 2)*
3470 (23244. - 145995. * sh - 135165. * pow(sh, 2) +
3471 2.214474e6 * pow(sh, 3) - 5.812671e6 * pow(sh, 4) +
3472 7.761562e6 * pow(sh, 5) - 6.239858e6 * pow(sh, 6) +
3473 3.192422e6 * pow(sh, 7) - 1.048133e6 * pow(sh, 8) +
3474 214001. * pow(sh, 9) - 24740. * pow(sh, 10) +
3475 1237. * pow(sh, 11)) -
3476 0.11851851851851852 * Li2sh * pow(-1 + sh, -8) * pow(sh, -7)*
3477 (-1080. + 8600. * sh - 29988. * pow(sh, 2) +
3478 59752. * pow(sh, 3) - 74250. * pow(sh, 4) +
3479 58630. * pow(sh, 5) - 28380. * pow(sh, 6) +
3480 7472. * pow(sh, 7) - 698. * pow(sh, 8) - 10. * pow(sh, 11) +
3481 pow(sh, 12) + pow(sh, 13)) +
3482 2.1333333333333333 * pow(ash, 2) * pow(-1 + sh, -4) *
3483 pow(sh, 2)*(60. - 1390. * sh + 3920. * pow(sh, 2) +
3484 11563. * pow(sh, 3) - 80966. * pow(sh, 4) +
3485 191832. * pow(sh, 5) - 258720. * pow(sh, 6) +
3486 224115. * pow(sh, 7) - 130416. * pow(sh, 8) +
3487 51582. * pow(sh, 9) - 13684. * pow(sh, 10) +
3488 2330. * pow(sh, 11) - 230. * pow(sh, 12) + 10. * pow(sh, 13))
3489 - Lsh * (26.074074074074073 * Li2sh * pow(-1 + sh, -3) *
3490 pow(sh, -5) + 0.7111111111111111 * ash * sqrt4sh * sqrtsh *
3491 pow(sh, 2)*(-204. + 654. * sh + 3621. * pow(sh, 2) -
3492 9390. * pow(sh, 3) + 8928. * pow(sh, 4) -
3493 4384. * pow(sh, 5) + 1190. * pow(sh, 6) -
3494 170. * pow(sh, 7) + 10. * pow(sh, 8)) +
3495 2.239578959155679e-6 * pow(-1 + sh, -8) * pow(sh, -12)*
3496 (3900. + 100. * sh - 116036. * pow(sh, 2) -
3497 329570. * pow(sh, 3) + 93490. * pow(sh, 4) +
3498 (4.740596068e9 - 5.71536e7 * Lshb) * pow(sh, 5) +
3499 4. * (-9.479986241e9 + 1.13778e8 * Lshb) * pow(sh, 6) +
3500 2. * (6.6407275091e10 - 7.9348248e8 * Lshb +
3501 4.851e6 * pow(3.141592653589793, 2)) * pow(sh, 7) +
3502 10. * (-2.6593392899e10 + 3.16207584e8 * Lshb -
3503 4.851e6 * pow(3.141592653589793, 2)) * pow(sh, 8) +
3504 10. * (3.3290343127e10 - 3.92931e8 * Lshb +
3505 9.702e6 * pow(3.141592653589793, 2)) * pow(sh, 9) +
3506 (-2.66790821833e11 + 3.1026996e9 * Lshb -
3507 9.702e7 * pow(3.141592653589793, 2)) * pow(sh, 10)\
3508 + (1.33671533459e11 - 1.5018696e9 * Lshb +
3509 4.851e7 * pow(3.141592653589793, 2)) * pow(sh, 11)\
3510 + (-3.8283679243e10 + 3.9541824e8 * Lshb -
3511 9.702e6 * pow(3.141592653589793, 2)) * pow(sh, 12)\
3512 + (4.796104102e9 - 3.693816e7 * Lshb) *
3513 pow(sh, 13) + 1.055215e6 * pow(sh, 14) -
3514 172235. * pow(sh, 15) -
3515 735. * (307. + 720. * Lshb) * pow(sh, 16) +
3516 1176. * (23. + 45. * Lshb) * pow(sh, 17) +
3517 294. * (11. + 180. * Lshb) * pow(sh, 18))) -
3518 0.00047031158142269255 * pow(Lsh, 2) * pow(-1 + sh, -9) *
3519 pow(sh, -12)*(-60. - 340. * sh - 312. * pow(sh, 2) +
3520 1872. * pow(sh, 3) + 4980. * pow(sh, 4) +
3521 1.1658504e7 * pow(sh, 5) - 1.04832846e8 * pow(sh, 6) +
3522 11. * (3.8074319e7 + 2520. * Lshb) * pow(sh, 7) -
3523 6. * (1.62766363e8 + 27720. * Lshb) * pow(sh, 8) +
3524 3. * (4.88268437e8 + 138600. * Lshb) * pow(sh, 9) -
3525 (1.465559189e9 + 554400. * Lshb) * pow(sh, 10) +
3526 (9.78147321e8 + 415800. * Lshb) * pow(sh, 11) -
3527 30. * (1.3999939e7 + 5544. * Lshb) * pow(sh, 12) +
3528 2. * (5.2652911e7 + 13860. * Lshb) * pow(sh, 13) -
3529 1.1751384e7 * pow(sh, 14) - 1260. * pow(sh, 16) +
3530 1386. * pow(sh, 17) - 126. * pow(sh, 19)) -
3531 0.00047031158142269255 * pow(lz, 2) * pow(sh, -12)*
3532 (15. + 220. * sh + 1518. * pow(sh, 2) + 6534. * pow(sh, 3) +
3533 19503. * pow(sh, 4) - 1.1615304e7 * pow(sh, 5) -
3534 20790. * pow(sh, 6) - 6546. * pow(sh, 7) -
3535 24741. * pow(sh, 8) - 29490. * pow(sh, 9) -
3536 19183. * pow(sh, 10) - 15. * pow(sh, 11) -
3537 42336. * pow(sh, 12) - 50652. * pow(sh, 13) -
3538 42336. * pow(sh, 14) - 522396. * pow(sh, 15) -
3539 514080. * pow(sh, 16) + 5.330178e6 * pow(sh, 17) -
3540 8.247204e6 * pow(sh, 18) + 6.008688e6 * pow(sh, 19) -
3541 2.443392e6 * pow(sh, 20) + 570780. * pow(sh, 21) -
3542 71820. * pow(sh, 22) + 3780. * pow(sh, 23)) -
3543 lz * (-512. * lmu * pow(sh, -7) -
3544 13.037037037037036 * pow(Lsh, 2) * pow(-1 + sh, -3) *
3545 pow(sh, -5) + 1.4222222222222223 * ash * sqrt4sh * sqrtsh *
3546 pow(sh, 2)*(-204. + 654. * sh + 3621. * pow(sh, 2) -
3547 9390. * pow(sh, 3) + 8928. * pow(sh, 4) -
3548 4384. * pow(sh, 5) + 1190. * pow(sh, 6) -
3549 170. * pow(sh, 7) + 10. * pow(sh, 8)) -
3550 0.0009406231628453851 * Lsh * pow(-1 + sh, -8) *
3551 pow(sh, -12)*(30. + 200. * sh + 356. * pow(sh, 2) -
3552 580. * pow(sh, 3) - 3070. * pow(sh, 4) -
3553 1.1660764e7 * pow(sh, 5) + 9.3173522e7 * pow(sh, 6) -
3554 3.25763741e8 * pow(sh, 7) +
3555 6.51103141e8 * pow(sh, 8) - 8.1359378e8 * pow(sh, 9) +
3556 6.50748419e8 * pow(sh, 10) -
3557 3.25285222e8 * pow(sh, 11) +
3558 9.2859728e7 * pow(sh, 12) -
3559 1.1567032e7 * pow(sh, 13) - 2380. * pow(sh, 14) -
3560 2065. * pow(sh, 15) - 945. * pow(sh, 16) +
3561 462. * pow(sh, 17) + 21. * pow(sh, 18)) +
3562 6.78660290653236e-8 * pow(-1 + sh, -7) * pow(sh, -12)*
3563 (64350. + 66000. * sh - 1.848594e6 * pow(sh, 2) -
3564 7.286499e6 * pow(sh, 3) - 5.743914e6 * pow(sh, 4) +
3565 1.564114497e11 * pow(sh, 5) -
3566 1.092996615942e12 * pow(sh, 6) +
3567 3.274039012563e12 * pow(sh, 7) -
3568 5.447084215617e12 * pow(sh, 8) +
3569 5.433819864693e12 * pow(sh, 9) -
3570 3.247981504196e12 * pow(sh, 10) +
3571 1.075526023436e12 * pow(sh, 11) -
3572 1.50436790331e11 * pow(sh, 12) -
3573 5.482108499e9 * pow(sh, 13) +
3574 1.2697937753e10 * pow(sh, 14) -
3575 1.820498844e10 * pow(sh, 15) +
3576 4.184087512e9 * pow(sh, 16) +
3577 6.5568484406e10 * pow(sh, 17) -
3578 9.5193663855e10 * pow(sh, 18) -
3579 3.83610629052e11 * pow(sh, 19) +
3580 1.851933873636e12 * pow(sh, 20) -
3581 3.850977196836e12 * pow(sh, 21) +
3582 4.997388254166e12 * pow(sh, 22) -
3583 4.45400972478e12 * pow(sh, 23) +
3584 2.82695858676e12 * pow(sh, 24) -
3585 1.2912095196e12 * pow(sh, 25) +
3586 4.2141221892e11 * pow(sh, 26) -
3587 9.582123474e10 * pow(sh, 27) +
3588 1.44039357e10 * pow(sh, 28) -
3589 1.28490516e9 * pow(sh, 29) + 5.143446e7 * pow(sh, 30)))
3590 - 4.896538893601991e-12 * pow(-1 + sh, -8) * pow(sh, -12)*
3591 (2.074545e8 - 4.76982e8 * sh -
3592 3.223561968e9 * pow(sh, 2) +
3593 1.2229671465e10 * pow(sh, 3) +
3594 5.5118379195e10 * pow(sh, 4) -
3595 5.39918795738718e14 * pow(sh, 5) +
3596 4.299200276151924e15 * pow(sh, 6) -
3597 1.497521829460536e16 * pow(sh, 7) +
3598 2.979736241687578e16 * pow(sh, 8) -
3599 3.705409524476913e16 * pow(sh, 9) +
3600 2.9516698533084984e16 * pow(sh, 10) -
3601 1.4749822777601014e16 * pow(sh, 11) +
3602 4.277878063797251e15 * pow(sh, 12) -
3603 6.52624973981284e14 * pow(sh, 13) +
3604 2.5170701729432e14 * pow(sh, 14) -
3605 4.0013657985217e14 * pow(sh, 15) +
3606 1.47801666642455e14 * pow(sh, 16) +
3607 1.2877818678578e15 * pow(sh, 17) -
3608 2.076397177767913e15 * pow(sh, 18) -
3609 8.385630050032431e15 * pow(sh, 19) +
3610 4.538622316201646e16 * pow(sh, 20) -
3611 1.061838620189715e17 * pow(sh, 21) +
3612 1.575009771316346e17 * pow(sh, 22) -
3613 1.6342684496066314e17 * pow(sh, 23) +
3614 1.2333678370646925e17 * pow(sh, 24) -
3615 6.8698738259084936e16 * pow(sh, 25) +
3616 2.8233898574272776e16 * pow(sh, 26) -
3617 8.447633225048556e15 * pow(sh, 27) +
3618 1.786675943700648e15 * pow(sh, 28) -
3619 2.5272282124782e14 * pow(sh, 29) +
3620 2.141343983394e13 * pow(sh, 30) -
3621 8.2047522942e11 * pow(sh, 31) +
3622 3.20166e7 * pow(3.141592653589793, 2)*
3623 (-90. - 600. * sh - 1068. * pow(sh, 2) +
3624 1740. * pow(sh, 3) + 9210. * pow(sh, 4) +
3625 3.5821452e7 * pow(sh, 5) - 2.86356066e8 * pow(sh, 6) +
3626 1.001798639e9 * pow(sh, 7) -
3627 2.003828431e9 * pow(sh, 8) +
3628 2.50626266e9 * pow(sh, 9) -
3629 2.006862017e9 * pow(sh, 10) +
3630 1.004429086e9 * pow(sh, 11) -
3631 2.8716356e8 * pow(sh, 12) + 3.622858e7 * pow(sh, 13) -
3632 1.140968e6 * pow(sh, 14) + 600307. * pow(sh, 15) +
3633 5.766227e6 * pow(sh, 16) - 8.667286e6 * pow(sh, 17) -
3634 7.4435635e7 * pow(sh, 18) +
3635 4.29448764e8 * pow(sh, 19) -
3636 1.174957056e9 * pow(sh, 20) +
3637 2.075532984e9 * pow(sh, 21) -
3638 2.59118244e9 * pow(sh, 22) +
3639 2.377353888e9 * pow(sh, 23) -
3640 1.631212128e9 * pow(sh, 24) +
3641 8.40926772e8 * pow(sh, 25) -
3642 3.23868888e8 * pow(sh, 26) +
3643 9.1648368e7 * pow(sh, 27) -
3644 1.8464544e7 * pow(sh, 28) + 2.50236e6 * pow(sh, 29) -
3645 204120. * pow(sh, 30) + 7560. * pow(sh, 31)))) * pow(z, 7);
3647 res += pow(z, 7)*(-51610.88930322329 -
3648 347049.77916686254 * (1. - sh) -
3649 1.3506484596566306e6 * pow(1. - sh, 2) -
3650 3.942197589163436e6 * pow(1. - sh, 3) -
3651 9.550309146362534e6 * pow(1. - sh, 4) +
3652 log(muh)*(587.9873015873015 +
3653 4627.9111111111115 * (1. - sh) +
3654 20303.644444444446 * pow(1. - sh, 2) +
3655 65689.6 * pow(1. - sh, 3) +
3656 174976. * pow(1. - sh, 4)) +
3657 log(z)*(10201.948065363957 +
3658 86007.6529570384 * (1. - sh) +
3659 379663.80341596896 * pow(1. - sh, 2) +
3660 1.2326804892307315e6 * pow(1. - sh, 3) +
3661 3.309516784057691e6 * pow(1. - sh, 4) +
3662 log(muh)*(512. + 3584. * (1. - sh) +
3663 14336. * pow(1. - sh, 2) +
3664 43008. * pow(1. - sh, 3) +
3665 107520. * pow(1. - sh, 4))) +
3666 (5506.910758377425 + 38490.637037037035 * (1. - sh) +
3667 152831.80246913582 * pow(1. - sh, 2) +
3668 457767.65432098764 * pow(1. - sh, 3) +
3669 1.1442307160493827e6 * pow(1. - sh, 4)) *
3675 res += -4.047867295730543e-6 * pow(3.141592653589793, 2) * pow(sh, -13)*
3676 (3003. + 47124. * sh + 349398. * pow(sh, 2) +
3677 1.62554e6 * pow(sh, 3) + 5.297685e6 * pow(sh, 4) +
3678 1.265412e7 * pow(sh, 5) + 2.15477e7 * pow(sh, 6) +
3679 1.265412e7 * pow(sh, 7) + 5.297685e6 * pow(sh, 8) +
3680 1.62554e6 * pow(sh, 9) + 349398. * pow(sh, 10) +
3681 47124. * pow(sh, 11) + 3003. * pow(sh, 12)) * pow(z, 7.5);
3683 res += (-2457.0347351559058 - 17199.24314609134 * (1. - sh) -
3684 71008.30384600567 * pow(1. - sh, 2) -
3685 224081.5678462186 * pow(1. - sh, 3) -
3686 596572.421257881 * pow(1. - sh, 4)) * pow(z, 7.5);
3691 res += (-1. * lmu * (-1940.6814814814816 * pow(sh, -8) +
3692 1642.6666666666667 * Lsh * pow(sh, -8)) -
3693 87.3113077596556 * pow(-1 + sh, -3) * pow(sh, -6) +
3694 72.63492063492063 * Li3sh * pow(-1 + sh, -3) * pow(sh, -6) -
3695 24.211640211640212 * pow(Lsh, 3) * pow(-1 + sh, -3) *
3696 pow(sh, -6) + 1.3037037037037038 * cl2 * sqrt4sh * sqrtsh *
3697 pow(sh, 2)*(-42. + 1287. * sh - 1614. * pow(sh, 2) -
3698 19017. * pow(sh, 3) + 50556. * pow(sh, 4) -
3699 55330. * pow(sh, 5) + 33315. * pow(sh, 6) -
3700 11945. * pow(sh, 7) + 2550. * pow(sh, 8) -
3701 300. * pow(sh, 9) + 15. * pow(sh, 10)) +
3702 0.00047031158142269255 * ash * sqrt4sh * sqrtsh *
3703 pow(-1 + sh, -3) * pow(sh, 2)*
3704 (228522. - 7.741593e6 * sh + 3.0723729e7 * pow(sh, 2) +
3705 6.0003996e7 * pow(sh, 3) - 5.8092981e8 * pow(sh, 4) +
3706 1.504380659e9 * pow(sh, 5) - 2.124842634e9 * pow(sh, 6) +
3707 1.888618168e9 * pow(sh, 7) - 1.114436468e9 * pow(sh, 8) +
3708 4.4384322e8 * pow(sh, 9) - 1.18087858e8 * pow(sh, 10) +
3709 2.0123511e7 * pow(sh, 11) - 1.986441e6 * pow(sh, 12) +
3710 86367. * pow(sh, 13)) -
3711 3.911111111111111 * pow(ash, 2) * pow(-1 + sh, -4) * pow(sh, 3)*
3712 (-700. + 7875. * sh - 14966. * pow(sh, 2) -
3713 79023. * pow(sh, 3) + 455850. * pow(sh, 4) -
3714 1.077135e6 * pow(sh, 5) + 1.523775e6 * pow(sh, 6) -
3715 1.432728e6 * pow(sh, 7) + 934836. * pow(sh, 8) -
3716 430150. * pow(sh, 9) + 139230. * pow(sh, 10) -
3717 31025. * pow(sh, 11) + 4530. * pow(sh, 12) -
3718 390. * pow(sh, 13) + 15. * pow(sh, 14)) -
3719 0.0028218694885361554 * Li2sh * pow(-1 + sh, -9) * pow(sh, -8)*
3720 (145530. - 1.30482e6 * sh + 5.20032e6 * pow(sh, 2) -
3721 1.208298e7 * pow(sh, 3) + 1.801254e7 * pow(sh, 4) -
3722 1.78178e7 * pow(sh, 5) + 1.163162e7 * pow(sh, 6) -
3723 4.77477e6 * pow(sh, 7) + 1.082136e6 * pow(sh, 8) -
3724 94842. * pow(sh, 9) - 420. * pow(sh, 13) +
3725 70. * pow(sh, 14) + 35. * pow(sh, 15)) -
3726 Lsh * (72.63492063492063 * Li2sh * pow(-1 + sh, -3) *
3727 pow(sh, -6) - 1.3037037037037038 * ash * sqrt4sh * sqrtsh *
3728 pow(sh, 2)*(-42. + 1287. * sh - 1614. * pow(sh, 2) -
3729 19017. * pow(sh, 3) + 50556. * pow(sh, 4) -
3730 55330. * pow(sh, 5) + 33315. * pow(sh, 6) -
3731 11945. * pow(sh, 7) + 2550. * pow(sh, 8) -
3732 300. * pow(sh, 9) + 15. * pow(sh, 10)) +
3733 7.998496282698853e-8 * pow(-1 + sh, -9) * pow(sh, -14)*
3734 (-70875. - 6750. * sh + 2.654475e6 * pow(sh, 2) +
3735 9.126108e6 * pow(sh, 3) + 678627. * pow(sh, 4) -
3736 4.643019e7 * pow(sh, 5) +
3737 25. * (-1.8964236617e10 + 2.05371936e8 * Lshb) *
3739 (-2.3705023777e10 + 2.5574472e8 * Lshb) * pow(sh, 7)\
3740 + 3. * (-5.691283317731e12 + 6.11557632e10 * Lshb -
3741 2.52252e8 * pow(3.141592653589793, 2)) * pow(sh, 8)\
3742 + 6. * (6.644381363023e12 - 7.10479224e10 * Lshb +
3743 7.56756e8 * pow(3.141592653589793, 2)) * pow(sh, 9)\
3744 + 225. * (-2.66023942651e11 +
3745 2.824366272e9 * Lshb -
3746 5.04504e7 * pow(3.141592653589793, 2)) * pow(sh, 10)
3747 + 120. * (4.99371243863e11 - 5.2384332e9 * Lshb +
3748 1.26126e8 * pow(3.141592653589793, 2)) * pow(sh, 11)
3749 + 225. * (-1.77796507943e11 +
3750 1.823838016e9 * Lshb -
3751 5.04504e7 * pow(3.141592653589793, 2)) * pow(sh, 12)
3752 + 2106. * (8.153537761e9 - 7.99876e7 * Lshb +
3753 2.156e6 * pow(3.141592653589793, 2)) * pow(sh, 13)\
3754 + (-4.299863603271e12 + 3.817775808e10 * Lshb -
3755 7.56756e8 * pow(3.141592653589793, 2)) * pow(sh, 14)
3756 - 140. * (-3.417060133e9 + 2.3900184e7 * Lshb) *
3757 pow(sh, 15) + 4.710132e7 * pow(sh, 16) +
3758 3.484152e7 * pow(sh, 17) - 7.916076e6 * pow(sh, 18) -
3759 6048. * (923. + 2450. * Lshb) * pow(sh, 19) +
3760 12600. * (59. + 196. * Lshb) * pow(sh, 20) +
3761 4200. * (13. + 294. * Lshb) * pow(sh, 21))) +
3762 0.00006718736877467037 * pow(Lsh, 2) * pow(-1 + sh, -10) *
3763 pow(sh, -14)*(-315. - 2115. * sh - 3255. * pow(sh, 2) +
3764 9681. * pow(sh, 3) + 37863. * pow(sh, 4) +
3765 21231. * pow(sh, 5) + 2.843603e8 * pow(sh, 6) -
3766 2.8424249e9 * pow(sh, 7) +
3767 6. * (2.129867959e9 + 90090. * Lshb) * pow(sh, 8) -
3768 42. * (8.10842339e8 + 90090. * Lshb) * pow(sh, 9) +
3769 84. * (7.09285376e8 + 135135. * Lshb) * pow(sh, 10) -
3770 1260. * (5.6749267e7 + 15015. * Lshb) * pow(sh, 11) +
3771 630. * (9.4632631e7 + 30030. * Lshb) * pow(sh, 12) -
3772 378. * (9.0217367e7 + 30030. * Lshb) * pow(sh, 13) +
3773 42. * (3.04968821e8 + 90090. * Lshb) * pow(sh, 14) -
3774 2. * (1.426541293e9 + 270270. * Lshb) * pow(sh, 15) +
3775 2.8627046e8 * pow(sh, 16) + 8820. * pow(sh, 19) -
3776 10290. * pow(sh, 20) + 735. * pow(sh, 21) +
3777 735. * pow(sh, 22)) +
3778 0.00001679684219366759 * pow(lz, 2) * pow(sh, -14)*
3779 (-315. - 5265. * sh - 41730. * pow(sh, 2) -
3780 208494. * pow(sh, 3) - 734877. * pow(sh, 4) -
3781 1.926639e6 * pow(sh, 5) + 1.1339384e9 * pow(sh, 6) +
3782 2.16216e6 * pow(sh, 7) + 48987. * pow(sh, 8) +
3783 1.669041e6 * pow(sh, 9) + 2.629962e6 * pow(sh, 10) +
3784 2.41431e6 * pow(sh, 11) + 1.436625e6 * pow(sh, 12) +
3785 315. * pow(sh, 13) + 3.764376e6 * pow(sh, 14) +
3786 4.230072e6 * pow(sh, 15) + 5.006232e6 * pow(sh, 16) -
3787 7.37352e6 * pow(sh, 17) + 1.0643094e8 * pow(sh, 18) +
3788 1.95359472e8 * pow(sh, 19) -
3789 1.425010356e9 * pow(sh, 20) +
3790 2.406212424e9 * pow(sh, 21) -
3791 2.05080876e9 * pow(sh, 22) + 1.03141962e9 * pow(sh, 23) -
3792 3.2026302e8 * pow(sh, 24) + 6.054048e7 * pow(sh, 25) -
3793 6.40332e6 * pow(sh, 26) + 291060. * pow(sh, 27)) -
3794 lz * (-1642.6666666666667 * lmu * pow(sh, -8) -
3795 36.317460317460316 * pow(Lsh, 2) * pow(-1 + sh, -3) *
3796 pow(sh, -6) - 2.6074074074074076 * ash * sqrt4sh * sqrtsh *
3797 pow(sh, 2)*(-42. + 1287. * sh - 1614. * pow(sh, 2) -
3798 19017. * pow(sh, 3) + 50556. * pow(sh, 4) -
3799 55330. * pow(sh, 5) + 33315. * pow(sh, 6) -
3800 11945. * pow(sh, 7) + 2550. * pow(sh, 8) -
3801 300. * pow(sh, 9) + 15. * pow(sh, 10)) -
3802 0.00006718736877467037 * Lsh * pow(-1 + sh, -9) *
3803 pow(sh, -14)*(-315. - 2430. * sh - 5685. * pow(sh, 2) +
3804 3996. * pow(sh, 3) + 41859. * pow(sh, 4) +
3805 63090. * pow(sh, 5) + 5.68853935e8 * pow(sh, 6) -
3806 5.11611966e9 * pow(sh, 7) +
3807 2.0446767111e10 * pow(sh, 8) -
3808 4.7676829206e10 * pow(sh, 9) +
3809 7.1479964025e10 * pow(sh, 10) -
3810 7.145286696e10 * pow(sh, 11) +
3811 4.7617276155e10 * pow(sh, 12) -
3812 2.039566545e10 * pow(sh, 13) +
3813 5.092669737e9 * pow(sh, 14) -
3814 5.64250624e8 * pow(sh, 15) + 37800. * pow(sh, 16) -
3815 47460. * pow(sh, 17) - 29232. * pow(sh, 18) -
3816 9828. * pow(sh, 19) + 6930. * pow(sh, 20) +
3817 210. * pow(sh, 21)) +
3818 2.796677021922676e-10 * pow(-1 + sh, -8) * pow(sh, -14)*
3819 (-1.0135125e7 - 1.1100375e7 * sh +
3820 3.6848955e8 * pow(sh, 2) +
3821 1.673522994e9 * pow(sh, 3) +
3822 1.770566655e9 * pow(sh, 4) -
3823 4.868950515e9 * pow(sh, 5) -
3824 1.355916766061e14 * pow(sh, 6) +
3825 1.08320914263562e15 * pow(sh, 7) -
3826 3.786163163932292e15 * pow(sh, 8) +
3827 7.56111391441978e15 * pow(sh, 9) -
3828 9.4338729874103e15 * pow(sh, 10) +
3829 7.52783212419654e15 * pow(sh, 11) -
3830 3.74950772533818e15 * pow(sh, 12) +
3831 1.064521016851092e15 * pow(sh, 13) -
3832 1.32239959812912e14 * pow(sh, 14) +
3833 4.822522515092e12 * pow(sh, 15) -
3834 1.3702469432263e13 * pow(sh, 16) +
3835 2.0211432885167e13 * pow(sh, 17) -
3836 1.5340657706842e13 * pow(sh, 18) +
3837 3.8517131981282e13 * pow(sh, 19) -
3838 1.98231128437323e14 * pow(sh, 20) +
3839 2.48611036890867e14 * pow(sh, 21) +
3840 1.312843416208632e15 * pow(sh, 22) -
3841 6.858549767677608e15 * pow(sh, 23) +
3842 1.6510967898168744e16 * pow(sh, 24) -
3843 2.57336419311972e16 * pow(sh, 25) +
3844 2.851348331813982e16 * pow(sh, 26) -
3845 2.336070503975208e16 * pow(sh, 27) +
3846 1.439716999494774e16 * pow(sh, 28) -
3847 6.70336841266092e15 * pow(sh, 29) +
3848 2.34397396677444e15 * pow(sh, 30) -
3849 6.0527470649646e14 * pow(sh, 31) +
3850 1.1182510300962e14 * pow(sh, 32) -
3851 1.396347000048e13 * pow(sh, 33) +
3852 1.054350297e12 * pow(sh, 34) -
3853 3.631041414e10 * pow(sh, 35))) +
3854 7.760786496621921e-16 * pow(-1 + sh, -9) * pow(sh, -14)*
3855 (7.34806697625e11 - 1.86465740175e12 * sh -
3856 1.3539275231625e13 * pow(sh, 2) +
3857 6.7027116030732e13 * pow(sh, 3) +
3858 3.63276406618983e14 * pow(sh, 4) +
3859 7.707518005473e13 * pow(sh, 5) -
3860 1.517387130121548e19 * pow(sh, 6) +
3861 1.361134491880303e20 * pow(sh, 7) -
3862 5.426853793481742e20 * pow(sh, 8) +
3863 1.262085401071148e21 * pow(sh, 9) -
3864 1.886986812950732e21 * pow(sh, 10) +
3865 1.8816744004909767e21 * pow(sh, 11) -
3866 1.2525799935976913e21 * pow(sh, 12) +
3867 5.37888683654858e20 * pow(sh, 13) -
3868 1.3573327538305124e20 * pow(sh, 14) +
3869 1.3896509010095016e19 * pow(sh, 15) +
3870 5.883932779473471e18 * pow(sh, 16) -
3871 1.0545119614382879e19 * pow(sh, 17) +
3872 1.1056449646096937e19 * pow(sh, 18) -
3873 2.3434959475666616e19 * pow(sh, 19) +
3874 1.0419598368525897e20 * pow(sh, 20) -
3875 1.688591180885733e20 * pow(sh, 21) -
3876 5.308738122339353e20 * pow(sh, 22) +
3877 3.5306726064196524e21 * pow(sh, 23) -
3878 9.689077766441178e21 * pow(sh, 24) +
3879 1.7097533811227756e22 * pow(sh, 25) -
3880 2.1588630869532427e22 * pow(sh, 26) +
3881 2.038422378040867e22 * pow(sh, 27) -
3882 1.4690695476739364e22 * pow(sh, 28) +
3883 8.14476034656345e21 * pow(sh, 29) -
3884 3.4697517653715583e21 * pow(sh, 30) +
3885 1.1250567404926722e21 * pow(sh, 31) -
3886 2.7233981875541377e20 * pow(sh, 32) +
3887 4.759269487750347e19 * pow(sh, 33) -
3888 5.663805036845516e18 * pow(sh, 34) +
3889 4.101841359381566e17 * pow(sh, 35) -
3890 1.362249094884678e16 * pow(sh, 36) +
3891 4.32864432e10 * pow(3.141592653589793, 2)*
3892 (-315. - 2430. * sh - 5685. * pow(sh, 2) +
3893 3996. * pow(sh, 3) + 41859. * pow(sh, 4) +
3894 63090. * pow(sh, 5) + 5.81369515e8 * pow(sh, 6) -
3895 5.23068642e9 * pow(sh, 7) +
3896 2.0914846601e10 * pow(sh, 8) -
3897 4.8797229546e10 * pow(sh, 9) +
3898 7.3211460495e10 * pow(sh, 10) -
3899 7.324423592e10 * pow(sh, 11) +
3900 4.8857255585e10 * pow(sh, 12) -
3901 2.094858111e10 * pow(sh, 13) +
3902 5.235387183e9 * pow(sh, 14) -
3903 5.73585856e8 * pow(sh, 15) -
3904 2.302929e7 * pow(sh, 16) + 4.790626e7 * pow(sh, 17) -
3905 1.24781262e8 * pow(sh, 18) +
3906 4.35359256e8 * pow(sh, 19) -
3907 5.45760408e8 * pow(sh, 20) -
3908 3.99911175e9 * pow(sh, 21) +
3909 2.559251793e10 * pow(sh, 22) -
3910 7.87509723e10 * pow(sh, 23) +
3911 1.60158041736e11 * pow(sh, 24) -
3912 2.35861111332e11 * pow(sh, 25) +
3913 2.6202440418e11 * pow(sh, 26) -
3914 2.2419142284e11 * pow(sh, 27) +
3915 1.491808626e11 * pow(sh, 28) -
3916 7.7359627884e10 * pow(sh, 29) +
3917 3.1118686368e10 * pow(sh, 30) -
3918 9.59731542e9 * pow(sh, 31) +
3919 2.22259884e9 * pow(sh, 32) -
3920 3.733653e8 * pow(sh, 33) + 4.288284e7 * pow(sh, 34) -
3921 3.00762e6 * pow(sh, 35) + 97020. * pow(sh, 36)))) *
3924 res += pow(z, 8)*(-177282.7451111828 -
3925 1.362845028834542e6 * (1. - sh) -
3926 5.9880687949730605e6 * pow(1. - sh, 2) -
3927 1.9481491807919256e7 * pow(1. - sh, 3) -
3928 5.2053685182060525e7 * pow(1. - sh, 4) +
3929 log(muh)*(1940.6814814814816 +
3930 17168.118518518517 * (1. - sh) +
3931 83827.2 * pow(1. - sh, 2) +
3932 299136. * pow(1. - sh, 3) +
3933 871904. * pow(1. - sh, 4)) +
3934 log(z)*(36437.04027695534 +
3935 345607.09636221116 * (1. - sh) +
3936 1.691397924301905e6 * pow(1. - sh, 2) +
3937 6.051670084075397e6 * pow(1. - sh, 3) +
3938 1.778888064751204e7 * pow(1. - sh, 4) +
3939 log(muh)*(1642.6666666666667 +
3940 13141.333333333334 * (1. - sh) +
3941 59136. * pow(1. - sh, 2) +
3942 197120. * pow(1. - sh, 3) +
3943 542080. * pow(1. - sh, 4))) +
3944 (19228.65843621399 + 153417.80082304525 * (1. - sh) +
3945 684734.1037037037 * pow(1. - sh, 2) +
3946 2.279762479012346e6 * pow(1. - sh, 3) +
3947 6.2693611104056435e6 * pow(1. - sh, 4)) *
3953 res += -1.4445330349469783e-6 * pow(3.141592653589793, 2) * pow(sh, -15)*
3954 (6435. + 114114. * sh + 962577. * pow(sh, 2) +
3955 5.137524e6 * pow(sh, 3) + 1.9445755e7 * pow(sh, 4) +
3956 5.526435e7 * pow(sh, 5) + 1.20557265e8 * pow(sh, 6) +
3957 1.9446652e8 * pow(sh, 7) + 1.20557265e8 * pow(sh, 8) +
3958 5.526435e7 * pow(sh, 9) + 1.9445755e7 * pow(sh, 10) +
3959 5.137524e6 * pow(sh, 11) + 962577. * pow(sh, 12) +
3960 114114. * pow(sh, 13) + 6435. * pow(sh, 14)) * pow(z, 8.5);
3962 res += (-8517.720415207139 - 68141.76332165711 * (1. - sh) -
3963 315581.54138342454 * pow(1. - sh, 2) -
3964 1.1026189077485641e6 * pow(1. - sh, 3) -
3965 3.2181012443704475e6 * pow(1. - sh, 4)) * pow(z, 8.5);
3970 res += (-1. * lmu * (-6542.252322163433 * pow(sh, -9) +
3971 5423.407407407408 * Lsh * pow(sh, -9)) -
3972 254.65798096566218 * pow(-1 + sh, -3) * pow(sh, -7) +
3973 211.85185185185185 * Li3sh * pow(-1 + sh, -3) * pow(sh, -7) -
3974 70.61728395061728 * pow(Lsh, 3) * pow(-1 + sh, -3) *
3975 pow(sh, -7) - 2.6901822457378013 * cl2 * sqrt4sh * sqrtsh *
3976 pow(sh, 3)*(552. - 6504. * sh + 954. * pow(sh, 2) +
3977 96534. * pow(sh, 3) - 256765. * pow(sh, 4) +
3978 311884. * pow(sh, 5) - 219753. * pow(sh, 6) +
3979 97149. * pow(sh, 7) - 27453. * pow(sh, 8) +
3980 4830. * pow(sh, 9) - 483. * pow(sh, 10) + 21. * pow(sh, 11))\
3981 - 0.00010451368476059834 * ash * sqrt4sh * sqrtsh *
3982 pow(-1 + sh, -3) * pow(sh, 3)*
3983 (-2.9602032e7 + 4.3918554e8 * sh -
3984 1.183665042e9 * pow(sh, 2) - 4.124801388e9 * pow(sh, 3) +
3985 2.9786749451e10 * pow(sh, 4) -
3986 7.6073083277e10 * pow(sh, 5) +
3987 1.12343648214e11 * pow(sh, 6) -
3988 1.08358670795e11 * pow(sh, 7) +
3989 7.1769426953e10 * pow(sh, 8) -
3990 3.3301199019e10 * pow(sh, 9) +
3991 1.0824913791e10 * pow(sh, 10) -
3992 2.416552923e9 * pow(sh, 11) +
3993 3.53029242e8 * pow(sh, 12) - 3.0393246e7 * pow(sh, 13) +
3994 1.168971e6 * pow(sh, 14)) +
3995 8.070546737213403 * pow(ash, 2) * pow(-1 + sh, -4) * pow(sh, 3)*
3996 (-112. + 5446. * sh - 40068. * pow(sh, 2) +
3997 45402. * pow(sh, 3) + 481386. * pow(sh, 4) -
3998 2.416545e6 * pow(sh, 5) + 5.68308e6 * pow(sh, 6) -
3999 8.356089e6 * pow(sh, 7) + 8.410689e6 * pow(sh, 8) -
4000 6.034798e6 * pow(sh, 9) + 3.142636e6 * pow(sh, 10) -
4001 1.191618e6 * pow(sh, 11) + 325686. * pow(sh, 12) -
4002 62481. * pow(sh, 13) + 7980. * pow(sh, 14) -
4003 609. * pow(sh, 15) + 21. * pow(sh, 16)) -
4004 0.004703115814226925 * Li2sh * pow(-1 + sh, -10) * pow(sh, -9)*
4005 (-288288. + 2.873871e6 * sh - 1.2889503e7 * pow(sh, 2) +
4006 3.423222e7 * pow(sh, 3) - 5.956146e7 * pow(sh, 4) +
4007 7.0825482e7 * pow(sh, 5) - 5.8121154e7 * pow(sh, 6) +
4008 3.2324292e7 * pow(sh, 7) - 1.153152e7 * pow(sh, 8) +
4009 2.327611e6 * pow(sh, 9) - 186043. * pow(sh, 10) -
4010 252. * pow(sh, 15) + 54. * pow(sh, 16) + 18. * pow(sh, 17))\
4011 - Lsh * (211.85185185185185 * Li2sh * pow(-1 + sh, -3) *
4012 pow(sh, -7) + 2.6901822457378013 * ash * sqrt4sh * sqrtsh *
4013 pow(sh, 3)*(552. - 6504. * sh + 954. * pow(sh, 2) +
4014 96534. * pow(sh, 3) - 256765. * pow(sh, 4) +
4015 311884. * pow(sh, 5) - 219753. * pow(sh, 6) +
4016 97149. * pow(sh, 7) - 27453. * pow(sh, 8) +
4017 4830. * pow(sh, 9) - 483. * pow(sh, 10) +
4018 21. * pow(sh, 11)) +
4019 1.5552631660803325e-7 * pow(-1 + sh, -10) * pow(sh, -16)*
4020 (24990. + 3675. * sh - 1.155555e6 * pow(sh, 2) -
4021 4.697175e6 * pow(sh, 3) - 2.300739e6 * pow(sh, 4) +
4022 2.4304707e7 * pow(sh, 5) + 4.9414911e7 * pow(sh, 6) +
4023 (8.82608902003e11 - 8.71782912e9 * Lshb) *
4024 pow(sh, 7) + (-8.826874533745e12 +
4025 8.690585904e10 * Lshb) * pow(sh, 8) +
4026 15. * (2.648605733525e12 - 2.5985238048e10 * Lshb +
4027 7.56756e7 * pow(3.141592653589793, 2)) * pow(sh, 9)\
4028 + 3. * (-3.5329412670653e13 +
4029 3.450607776e11 * Lshb -
4030 2.648646e9 * pow(3.141592653589793, 2)) *
4032 105. * (1.767535730921e12 - 1.715370048e10 * Lshb +
4033 2.270268e8 * pow(3.141592653589793, 2)) *
4035 21. * (-1.0613371381061e13 +
4036 1.0198869408e11 * Lshb -
4037 1.89189e9 * pow(3.141592653589793, 2)) * pow(sh, 12)
4038 + 189. * (9.83624804999e11 - 9.29938464e9 * Lshb +
4039 2.1021e8 * pow(3.141592653589793, 2)) * pow(sh, 13)\
4040 + 15. * (-7.089591334283e12 +
4041 6.5165772672e10 * Lshb -
4042 1.5891876e9 * pow(3.141592653589793, 2)) *
4044 33. * (1.209854560067e12 - 1.05670656e10 * Lshb +
4045 2.40786e8 * pow(3.141592653589793, 2)) * pow(sh, 15)
4046 + 2. * (-4.441328248151e12 +
4047 3.519347832e10 * Lshb -
4048 5.67567e8 * pow(3.141592653589793, 2)) * pow(sh, 16)
4049 + (8.88980806162e11 - 5.62594032e9 * Lshb) *
4050 pow(sh, 17) - 3.8377206e7 * pow(sh, 18) +
4051 4.2700014e7 * pow(sh, 19) +
4052 1.8171216e7 * pow(sh, 20) - 5.847282e6 * pow(sh, 21) -
4053 1890. * (1387. + 4032. * Lshb) * pow(sh, 22) +
4054 2430. * (149. + 672. * Lshb) * pow(sh, 23) +
4055 1215. * (15. + 448. * Lshb) * pow(sh, 24))) -
4056 0.00009798157946306095 * pow(Lsh, 2) * pow(-1 + sh, -11) *
4057 pow(sh, -16)*(-168. - 1302. * sh - 2760. * pow(sh, 2) +
4058 4440. * pow(sh, 3) + 27552. * pow(sh, 4) +
4059 26796. * pow(sh, 5) - 59136. * pow(sh, 6) +
4060 6.91904744e8 * pow(sh, 7) - 7.607607766e9 * pow(sh, 8) +
4061 (3.8011106171e10 + 1.08108e6 * Lshb) * pow(sh, 9) -
4062 36. * (3.165700357e9 + 240240. * Lshb) * pow(sh, 10) +
4063 720. * (3.16454927e8 + 42042. * Lshb) * pow(sh, 11) -
4064 504. * (6.32852603e8 + 120120. * Lshb) * pow(sh, 12) +
4065 84. * (3.797894797e9 + 900900. * Lshb) * pow(sh, 13) -
4066 192. * (1.187472428e9 + 315315. * Lshb) * pow(sh, 14) +
4067 120. * (9.50840131e8 + 252252. * Lshb) * pow(sh, 15) -
4068 104. * (3.66207253e8 + 83160. * Lshb) * pow(sh, 16) +
4069 (7.631953801e9 + 1.08108e6 * Lshb) * pow(sh, 17) -
4070 6.95689448e8 * pow(sh, 18) - 6048. * pow(sh, 22) +
4071 7344. * pow(sh, 23) - 864. * pow(sh, 24) - 432. * pow(sh, 25)
4072 ) + 0.00004899078973153047 * pow(lz, 2) * pow(sh, -16)*
4073 (-84. - 1575. * sh - 14085. * pow(sh, 2) -
4074 79950. * pow(sh, 3) - 323154. * pow(sh, 4) -
4075 987129. * pow(sh, 5) - 2.347059e6 * pow(sh, 6) +
4076 1.379703232e9 * pow(sh, 7) + 2.7027e6 * pow(sh, 8) -
4077 473685. * pow(sh, 9) + 1.243629e6 * pow(sh, 10) +
4078 2.698722e6 * pow(sh, 11) + 2.985906e6 * pow(sh, 12) +
4079 2.433285e6 * pow(sh, 13) + 1.388619e6 * pow(sh, 14) +
4080 84. * pow(sh, 15) + 4.25568e6 * pow(sh, 16) +
4081 4.612608e6 * pow(sh, 17) + 5.189184e6 * pow(sh, 18) +
4082 7.605312e6 * pow(sh, 19) - 6.1281792e7 * pow(sh, 20) +
4083 2.75438592e8 * pow(sh, 21) + 8.71398528e8 * pow(sh, 22) -
4084 5.028236928e9 * pow(sh, 23) +
4085 8.99602704e9 * pow(sh, 24) -
4086 8.653636992e9 * pow(sh, 25) +
4087 5.137333344e9 * pow(sh, 26) -
4088 1.978438176e9 * pow(sh, 27) +
4089 4.97379168e8 * pow(sh, 28) - 7.8990912e7 * pow(sh, 29) +
4090 7.2072e6 * pow(sh, 30) - 288288. * pow(sh, 31)) -
4091 lz * (-5423.407407407408 * lmu * pow(sh, -9) -
4092 105.92592592592592 * pow(Lsh, 2) * pow(-1 + sh, -3) *
4093 pow(sh, -7) + 5.380364491475603 * ash * sqrt4sh * sqrtsh *
4094 pow(sh, 3)*(552. - 6504. * sh + 954. * pow(sh, 2) +
4095 96534. * pow(sh, 3) - 256765. * pow(sh, 4) +
4096 311884. * pow(sh, 5) - 219753. * pow(sh, 6) +
4097 97149. * pow(sh, 7) - 27453. * pow(sh, 8) +
4098 4830. * pow(sh, 9) - 483. * pow(sh, 10) +
4099 21. * pow(sh, 11)) +
4100 0.0001959631589261219 * Lsh * pow(-1 + sh, -10) *
4101 pow(sh, -16)*(-84. - 735. * sh - 2115. * pow(sh, 2) +
4102 105. * pow(sh, 3) + 13881. * pow(sh, 4) +
4103 27279. * pow(sh, 5) - 2289. * pow(sh, 6) +
4104 6.91956971e8 * pow(sh, 7) -
4105 6.915895655e9 * pow(sh, 8) +
4106 3.1099607139e10 * pow(sh, 9) -
4107 8.2880063541e10 * pow(sh, 10) +
4108 1.44964071231e11 * pow(sh, 11) -
4109 1.73879170521e11 * pow(sh, 12) +
4110 1.44837476805e11 * pow(sh, 13) -
4111 8.2722854787e10 * pow(sh, 14) +
4112 3.0998246241e10 * pow(sh, 15) -
4113 6.879851159e9 * pow(sh, 16) +
4114 6.86365796e8 * pow(sh, 17) - 75600. * pow(sh, 18) -
4115 9072. * pow(sh, 19) + 22680. * pow(sh, 20) +
4116 9828. * pow(sh, 21) + 2268. * pow(sh, 22) -
4117 2484. * pow(sh, 23) - 54. * pow(sh, 24)) +
4118 5.43798309818298e-10 * pow(-1 + sh, -9) * pow(sh, -16)*
4119 (3.57357e6 + 4.099095e6 * sh -
4120 1.6114527e8 * pow(sh, 2) - 8.32841295e8 * pow(sh, 3) -
4121 1.161846972e9 * pow(sh, 4) +
4122 2.313726129e9 * pow(sh, 5) +
4123 9.380058402e9 * pow(sh, 6) +
4124 2.52441684058135e14 * pow(sh, 7) -
4125 2.26944762991947e15 * pow(sh, 8) +
4126 9.067251405053984e15 * pow(sh, 9) -
4127 2.113009526738253e16 * pow(sh, 10) +
4128 3.1647467902192624e16 * pow(sh, 11) -
4129 3.158667323807502e16 * pow(sh, 12) +
4130 2.1002634158609704e16 * pow(sh, 13) -
4131 8.966860615674702e15 * pow(sh, 14) +
4132 2.228333766410625e15 * pow(sh, 15) -
4133 2.43656770617994e14 * pow(sh, 16) -
4134 9.520266140208e12 * pow(sh, 17) +
4135 3.1476144603672e13 * pow(sh, 18) -
4136 6.08108801358e13 * pow(sh, 19) +
4137 8.0522179628688e13 * pow(sh, 20) -
4138 7.1735262216768e13 * pow(sh, 21) -
4139 1.85972284144884e14 * pow(sh, 22) +
4140 1.24257330785952e15 * pow(sh, 23) -
4141 1.516348491395922e15 * pow(sh, 24) -
4142 9.623404446386016e15 * pow(sh, 25) +
4143 5.409785863372176e16 * pow(sh, 26) -
4144 1.4651551524606797e17 * pow(sh, 27) +
4145 2.639914855668742e17 * pow(sh, 28) -
4146 3.46785263569046e17 * pow(sh, 29) +
4147 3.4581010205224454e17 * pow(sh, 30) -
4148 2.67102371124581e17 * pow(sh, 31) +
4149 1.612804363486636e17 * pow(sh, 32) -
4150 7.625393556018638e16 * pow(sh, 33) +
4151 2.809203040360517e16 * pow(sh, 34) -
4152 7.968174811781184e15 * pow(sh, 35) +
4153 1.704011763662208e15 * pow(sh, 36) -
4154 2.65362784294608e14 * pow(sh, 37) +
4155 2.836102317048e13 * pow(sh, 38) -
4156 1.857724700832e12 * pow(sh, 39) +
4157 5.6166718608e10 * pow(sh, 40))) -
4158 2.263562728181394e-15 * pow(-1 + sh, -10) * pow(sh, -16)*
4159 (1.522039519e11 - 4.26276275425e11 * sh -
4160 3.315607505925e12 * pow(sh, 2) +
4161 1.9935815474575e13 * pow(sh, 3) +
4162 1.24680380503679e14 * pow(sh, 4) +
4163 9.4734104241777e13 * pow(sh, 5) -
4164 6.72219759402191e14 * pow(sh, 6) -
4165 2.148676883608435e19 * pow(sh, 7) +
4166 2.1433741656451708e20 * pow(sh, 8) -
4167 9.621528485908101e20 * pow(sh, 9) +
4168 2.559493442861416e21 * pow(sh, 10) -
4169 4.468559626678942e21 * pow(sh, 11) +
4170 5.35099530515921e21 * pow(sh, 12) -
4171 4.4525844044562044e21 * pow(sh, 13) +
4172 2.544142477915097e21 * pow(sh, 14) -
4173 9.568948351796755e20 * pow(sh, 15) +
4174 2.1500361682429e20 * pow(sh, 16) -
4175 2.422139078794151e19 * pow(sh, 17) +
4176 8.145618788181729e18 * pow(sh, 18) -
4177 1.8132101139625396e19 * pow(sh, 19) +
4178 2.789986937416512e19 * pow(sh, 20) -
4179 2.916708794791665e19 * pow(sh, 21) -
4180 4.506612585214779e19 * pow(sh, 22) +
4181 3.9814881989842764e20 * pow(sh, 23) -
4182 7.61642854944141e20 * pow(sh, 24) -
4183 1.912231127650318e21 * pow(sh, 25) +
4184 1.5547276029752026e22 * pow(sh, 26) -
4185 4.875241535562659e22 * pow(sh, 27) +
4186 9.909337425361018e22 * pow(sh, 28) -
4187 1.4647795672750128e23 * pow(sh, 29) +
4188 1.651330016200192e23 * pow(sh, 30) -
4189 1.453924699389752e23 * pow(sh, 31) +
4190 1.011735575962528e23 * pow(sh, 32) -
4191 5.58883301567305e22 * pow(sh, 33) +
4192 2.447160559972393e22 * pow(sh, 34) -
4193 8.433457866753118e21 * pow(sh, 35) +
4194 2.2566304052798968e21 * pow(sh, 36) -
4195 4.585243554016694e20 * pow(sh, 37) +
4196 6.826335499960051e19 * pow(sh, 38) -
4197 7.011890262150629e18 * pow(sh, 39) +
4198 4.434709029248097e17 * pow(sh, 40) -
4199 1.2998152420597344e16 * pow(sh, 41) +
4200 4.32864432e10 * pow(3.141592653589793, 2)*
4201 (-84. - 735. * sh - 2115. * pow(sh, 2) +
4202 105. * pow(sh, 3) + 13881. * pow(sh, 4) +
4203 27279. * pow(sh, 5) - 2289. * pow(sh, 6) +
4204 7.06083083e8 * pow(sh, 7) -
4205 7.059427043e9 * pow(sh, 8) +
4206 3.1757931649e10 * pow(sh, 9) -
4207 8.4675129727e10 * pow(sh, 10) +
4208 1.48186623501e11 * pow(sh, 11) -
4209 1.77858935891e11 * pow(sh, 12) +
4210 1.48260913367e11 * pow(sh, 13) -
4211 8.4747202497e10 * pow(sh, 14) +
4212 3.1784695387e10 * pow(sh, 15) -
4213 7.061392333e9 * pow(sh, 16) +
4214 7.13289054e8 * pow(sh, 17) -
4215 3.2675014e7 * pow(sh, 18) + 7.242013e7 * pow(sh, 19) -
4216 7.6735758e7 * pow(sh, 20) -
4217 2.40223726e8 * pow(sh, 21) +
4218 1.627847962e9 * pow(sh, 22) -
4219 2.273776974e9 * pow(sh, 23) -
4220 1.712105386e10 * pow(sh, 24) +
4221 1.18341944808e11 * pow(sh, 25) -
4222 3.9995878208e11 * pow(sh, 26) +
4223 9.10773328608e11 * pow(sh, 27) -
4224 1.531097366656e12 * pow(sh, 28) +
4225 1.98076916752e12 * pow(sh, 29) -
4226 2.016366553344e12 * pow(sh, 30) +
4227 1.634962897568e12 * pow(sh, 31) -
4228 1.06202452928e12 * pow(sh, 32) +
4229 5.53053223008e11 * pow(sh, 33) -
4230 2.30009473408e11 * pow(sh, 34) +
4231 7.5735678304e10 * pow(sh, 35) -
4232 1.945628256e10 * pow(sh, 36) +
4233 3.810741792e9 * pow(sh, 37) -
4234 5.48735616e8 * pow(sh, 38) +
4235 5.4678624e7 * pow(sh, 39) - 3.36336e6 * pow(sh, 40) +
4236 96096. * pow(sh, 41)))) * pow(z, 9);
4238 res += pow(z, 9)*(-622069.3836274818 -
4239 5.380712316126128e6 * (1. - sh) -
4240 2.6348101173329648e7 * pow(1. - sh, 2) -
4241 9.453816935977913e7 * pow(1. - sh, 3) -
4242 2.761541131481308e8 * pow(1. - sh, 4) +
4243 log(muh)*(6542.252322163433 +
4244 64303.678306878304 * (1. - sh) +
4245 345923.7248677249 * pow(1. - sh, 2) +
4246 1.349738102292769e6 * pow(1. - sh, 3) +
4247 4.272929862433863e6 * pow(1. - sh, 4)) +
4248 log(z)*(131849.69275238464 +
4249 1.3908505037480777e6 * (1. - sh) +
4250 7.472361309110008e6 * pow(1. - sh, 2) +
4251 2.920814956913386e7 * pow(1. - sh, 3) +
4252 9.329175253624554e7 * pow(1. - sh, 4) +
4253 log(muh)*(5423.407407407408 +
4254 48810.666666666664 * (1. - sh) +
4255 244053.33333333334 * pow(1. - sh, 2) +
4256 894862.2222222222 * pow(1. - sh, 3) +
4257 2.6845866666666665e6 * pow(1. - sh, 4))) +
4258 (68317.8044287674 + 612595.3262786596 * (1. - sh) +
4259 3.035726507936508e6 * pow(1. - sh, 2) +
4260 1.1122724390358612e7 * pow(1. - sh, 3) +
4261 3.3373027717813052e7 * pow(1. - sh, 4)) *
4267 res += -6.708357747431788e-8 * pow(3.141592653589793, 2) * pow(sh, -17)*
4268 (109395. + 2.16216e6 * sh + 2.0430696e7 * pow(sh, 2) +
4269 1.22866128e8 * pow(sh, 3) + 5.27900436e8 * pow(sh, 4) +
4270 1.72236568e9 * pow(sh, 5) + 4.41295356e9 * pow(sh, 6) +
4271 8.96714376e9 * pow(sh, 7) + 1.385377093e10 * pow(sh, 8) +
4272 8.96714376e9 * pow(sh, 9) + 4.41295356e9 * pow(sh, 10) +
4273 1.72236568e9 * pow(sh, 11) + 5.27900436e8 * pow(sh, 12) +
4274 1.22866128e8 * pow(sh, 13) + 2.0430696e7 * pow(sh, 14) +
4275 2.16216e6 * pow(sh, 15) + 109395. * pow(sh, 16)) * pow(z, 9.5);
4277 res += (-30062.542641907552 - 270562.883777168 * (1. - sh) -
4278 1.388889470056129e6 * pow(1. - sh, 2) -
4279 5.321070047617637e6 * pow(1. - sh, 3) -
4280 1.6887633329091568e7 * pow(1. - sh, 4)) * pow(z, 9.5);
4285 res += (-1. * lmu * (-22436.0126984127 * pow(sh, -10) +
4286 18304. * Lsh * pow(sh, -10)) -
4287 769.6330091406679 * pow(-1 + sh, -3) * pow(sh, -8) +
4288 640.2633744855967 * Li3sh * pow(-1 + sh, -3) * pow(sh, -8) -
4289 213.42112482853224 * pow(Lsh, 3) * pow(-1 + sh, -3) *
4290 pow(sh, -8) + 6.052910052910053 * cl2 * sqrt4sh * sqrtsh *
4291 pow(sh, 3)*(72. - 4364. * sh + 28750. * pow(sh, 2) +
4292 23168. * pow(sh, 3) - 476659. * pow(sh, 4) +
4293 1.255353e6 * pow(sh, 5) - 1.649661e6 * pow(sh, 6) +
4294 1.316274e6 * pow(sh, 7) - 686175. * pow(sh, 8) +
4295 239883. * pow(sh, 9) - 55951. * pow(sh, 10) +
4296 8372. * pow(sh, 11) - 728. * pow(sh, 12) + 28. * pow(sh, 13))
4297 + 0.0009406231628453851 * ash * sqrt4sh * sqrtsh *
4298 pow(-1 + sh, -3) * pow(sh, 3)*
4299 (-1.003446e6 + 6.399262e7 * sh -
4300 5.88476654e8 * pow(sh, 2) + 1.052685443e9 * pow(sh, 3) +
4301 6.558619027e9 * pow(sh, 4) -
4302 3.9015984112e10 * pow(sh, 5) +
4303 9.831642973e10 * pow(sh, 6) -
4304 1.5050270954e11 * pow(sh, 7) +
4305 1.55340604554e11 * pow(sh, 8) -
4306 1.13220265218e11 * pow(sh, 9) +
4307 5.9521934912e10 * pow(sh, 10) -
4308 2.2692174847e10 * pow(sh, 11) +
4309 6.219574195e9 * pow(sh, 12) -
4310 1.194652762e9 * pow(sh, 13) + 1.5263422e8 * pow(sh, 14) -
4311 1.1648401e7 * pow(sh, 15) + 401669. * pow(sh, 16)) -
4312 18.158730158730158 * pow(ash, 2) * pow(-1 + sh, -4) *
4313 pow(sh, 4)*(1638. - 34902. * sh + 187614. * pow(sh, 2) -
4314 73296. * pow(sh, 3) - 2.717613e6 * pow(sh, 4) +
4315 1.2278721e7 * pow(sh, 5) - 2.8715127e7 * pow(sh, 6) +
4316 4.3584072e7 * pow(sh, 7) - 4.6464222e7 * pow(sh, 8) +
4317 3.6125334e7 * pow(sh, 9) - 2.0866662e7 * pow(sh, 10) +
4318 9.012312e6 * pow(sh, 11) - 2.899701e6 * pow(sh, 12) +
4319 684873. * pow(sh, 13) - 115255. * pow(sh, 14) +
4320 13076. * pow(sh, 15) - 896. * pow(sh, 16) + 28. * pow(sh, 17)
4321 ) - 0.0011757789535567313 * Li2sh * pow(-1 + sh, -11) *
4322 pow(sh, -10)*(3.891888e6 - 4.2698656e7 * sh +
4323 2.12860648e8 * pow(sh, 2) - 6.36197562e8 * pow(sh, 3) +
4324 1.265813472e9 * pow(sh, 4) - 1.758542016e9 * pow(sh, 5) +
4325 1.737624672e9 * pow(sh, 6) - 1.217575632e9 * pow(sh, 7) +
4326 5.89815408e8 * pow(sh, 8) - 1.86234048e8 * pow(sh, 9) +
4327 3.3715396e7 * pow(sh, 10) - 2.523664e6 * pow(sh, 11) -
4328 1008. * pow(sh, 17) + 252. * pow(sh, 18) + 63. * pow(sh, 19))
4329 - Lsh * (640.2633744855967 * Li2sh * pow(-1 + sh, -3) *
4330 pow(sh, -8) - 6.052910052910053 * ash * sqrt4sh * sqrtsh *
4331 pow(sh, 3)*(72. - 4364. * sh + 28750. * pow(sh, 2) +
4332 23168. * pow(sh, 3) - 476659. * pow(sh, 4) +
4333 1.255353e6 * pow(sh, 5) - 1.649661e6 * pow(sh, 6) +
4334 1.316274e6 * pow(sh, 7) - 686175. * pow(sh, 8) +
4335 239883. * pow(sh, 9) - 55951. * pow(sh, 10) +
4336 8372. * pow(sh, 11) - 728. * pow(sh, 12) +
4337 28. * pow(sh, 13)) +
4338 1.0368421107202216e-8 * pow(-1 + sh, -11) * pow(sh, -18)*
4339 (-268128. - 50176. * sh + 1.5037316e7 * pow(sh, 2) +
4340 7.0758368e7 * pow(sh, 3) + 6.454424e7 * pow(sh, 4) -
4341 3.47363296e8 * pow(sh, 5) -
4342 9.86787592e8 * pow(sh, 6) - 3.5646688e7 * pow(sh, 7) +
4343 448. * (-1.08098395051e11 + 9.8513415e8 * Lshb) *
4344 pow(sh, 8) - 2464. *
4345 (-2.16214325051e11 + 1.9651086e9 * Lshb) * pow(sh, 9)
4346 + (-2.664090635973398e15 +
4347 2.41383974832e13 * Lshb -
4348 5.1459408e10 * pow(3.141592653589793, 2)) *
4350 (7.994327088697714e15 - 7.21448035308e13 * Lshb +
4351 4.11675264e11 * pow(3.141592653589793, 2)) *
4353 28. * (-5.71247555169563e14 +
4354 5.1265445616e12 * Lshb -
4355 5.1459408e10 * pow(3.141592653589793, 2)) *
4357 112. * (2.00040160003213e14 -
4358 1.7805237912e12 * Lshb +
4359 2.5729704e10 * pow(3.141592653589793, 2)) *
4361 14. * (-1.601332890109447e15 +
4362 1.40747598432e13 * Lshb -
4363 2.5729704e11 * pow(3.141592653589793, 2)) *
4365 112. * (1.43081208030667e14 -
4366 1.2327953274e12 * Lshb +
4367 2.5729704e10 * pow(3.141592653589793, 2)) *
4369 28. * (-2.86396712598419e14 +
4370 2.3887524024e12 * Lshb -
4371 5.1459408e10 * pow(3.141592653589793, 2)) *
4373 32. * (8.3604554894621e13 - 6.599669076e11 * Lshb +
4374 1.2864852e10 * pow(3.141592653589793, 2)) *
4376 (-5.35506196250843e14 + 3.8233259064e12 * Lshb -
4377 5.1459408e10 * pow(3.141592653589793, 2)) *
4379 252. * (-1.93318721627e11 + 1.1356488e9 * Lshb) *
4380 pow(sh, 19) - 3.861655812e9 * pow(sh, 20) -
4381 5.0733144e7 * pow(sh, 21) +
4382 8.74592334e8 * pow(sh, 22) +
4383 2.56502232e8 * pow(sh, 23) -
4384 1.152693e8 * pow(sh, 24) -
4385 1800. * (20357. + 63504. * Lshb) * pow(sh, 25) +
4386 11025. * (461. + 2592. * Lshb) * pow(sh, 26) +
4387 11025. * (17. + 648. * Lshb) * pow(sh, 27))) -
4388 0.000013064210595074793 * pow(Lsh, 2) * pow(-1 + sh, -12) *
4389 pow(sh, -18)*(1008. + 8848. * sh + 23688. * pow(sh, 2) -
4390 17352. * pow(sh, 3) - 212592. * pow(sh, 4) -
4391 311184. * pow(sh, 5) + 400176. * pow(sh, 6) +
4392 1.492176e6 * pow(sh, 7) - 1.8665467128e10 * pow(sh, 8) +
4393 2.23865173224e11 * pow(sh, 9) -
4394 3. * (4.10179236011e11 + 8.16816e6 * Lshb) * pow(sh, 10) +
4395 1053. * (3.893290287e9 + 209440. * Lshb) * pow(sh, 11) -
4396 234. * (3.9404650781e10 + 3.76992e6 * Lshb) *
4397 pow(sh, 12) + 1848. *
4398 (7.981764103e9 + 1.11384e6 * Lshb) * pow(sh, 13) -
4399 5040. * (3.414494773e9 + 612612. * Lshb) * pow(sh, 14) +
4400 3528. * (4.182077137e9 + 875160. * Lshb) * pow(sh, 15) -
4401 1344. * (6.864759043e9 + 1.53153e6 * Lshb) * pow(sh, 16) +
4402 1224. * (3.352868227e9 + 720720. * Lshb) * pow(sh, 17) -
4403 18. * (6.8479648117e10 + 1.225224e7 * Lshb) *
4404 pow(sh, 18) + (2.24483274986e11 + 2.450448e7 * Lshb) *
4405 pow(sh, 19) - 1.8748455264e10 * pow(sh, 20) -
4406 45360. * pow(sh, 25) + 56700. * pow(sh, 26) -
4407 8505. * pow(sh, 27) - 2835. * pow(sh, 28)) +
4408 0.000013064210595074793 * pow(lz, 2) * pow(sh, -18)*
4409 (-252. - 5236. * sh - 52122. * pow(sh, 2) -
4410 330990. * pow(sh, 3) - 1.50586e6 * pow(sh, 4) -
4411 5.221788e6 * pow(sh, 5) - 1.4306214e7 * pow(sh, 6) -
4412 3.1500898e7 * pow(sh, 7) + 1.8610077024e10 * pow(sh, 8) +
4413 3.675672e7 * pow(sh, 9) - 1.1865305e7 * pow(sh, 10) +
4414 7.819275e6 * pow(sh, 11) + 2.920308e7 * pow(sh, 12) +
4415 3.763735e7 * pow(sh, 13) + 3.5618775e7 * pow(sh, 14) +
4416 2.7424047e7 * pow(sh, 15) + 1.5311086e7 * pow(sh, 16) +
4417 252. * pow(sh, 17) + 5.397678e7 * pow(sh, 18) +
4418 5.722002e7 * pow(sh, 19) + 6.231654e7 * pow(sh, 20) +
4419 6.996132e7 * pow(sh, 21) + 2.7057888e8 * pow(sh, 22) -
4420 3.18810492e9 * pow(sh, 23) + 7.5729654e9 * pow(sh, 24) +
4421 4.194296964e10 * pow(sh, 25) -
4422 2.0547002619e11 * pow(sh, 26) +
4423 3.8218629735e11 * pow(sh, 27) -
4424 4.0430890071e11 * pow(sh, 28) +
4425 2.7416659842e11 * pow(sh, 29) -
4426 1.2521582073e11 * pow(sh, 30) +
4427 3.911266359e10 * pow(sh, 31) -
4428 8.26458633e9 * pow(sh, 32) + 1.13189076e9 * pow(sh, 33) -
4429 9.081072e7 * pow(sh, 34) + 3.24324e6 * pow(sh, 35)) -
4430 lz * (-18304. * lmu * pow(sh, -10) -
4431 320.13168724279836 * pow(Lsh, 2) * pow(-1 + sh, -3) *
4432 pow(sh, -8) - 12.105820105820106 * ash * sqrt4sh * sqrtsh *
4433 pow(sh, 3)*(72. - 4364. * sh + 28750. * pow(sh, 2) +
4434 23168. * pow(sh, 3) - 476659. * pow(sh, 4) +
4435 1.255353e6 * pow(sh, 5) - 1.649661e6 * pow(sh, 6) +
4436 1.316274e6 * pow(sh, 7) - 686175. * pow(sh, 8) +
4437 239883. * pow(sh, 9) - 55951. * pow(sh, 10) +
4438 8372. * pow(sh, 11) - 728. * pow(sh, 12) +
4439 28. * pow(sh, 13)) -
4440 0.000026128421190149586 * Lsh * pow(-1 + sh, -11) *
4441 pow(sh, -18)*(-504. - 4928. * sh - 16772. * pow(sh, 2) -
4442 8096. * pow(sh, 3) + 98200. * pow(sh, 4) +
4443 253792. * pow(sh, 5) + 53704. * pow(sh, 6) -
4444 692384. * pow(sh, 7) + 1.8664647232e10 * pow(sh, 8) -
4445 2.05207139984e11 * pow(sh, 9) +
4446 1.025430782207e12 * pow(sh, 10) -
4447 3.074555676226e12 * pow(sh, 11) +
4448 6.145992007736e12 * pow(sh, 12) -
4449 8.60044812004e12 * pow(sh, 13) +
4450 8.596678573448e12 * pow(sh, 14) -
4451 6.137611319224e12 * pow(sh, 15) +
4452 3.067041035912e12 * pow(sh, 16) -
4453 1.021522277272e12 * pow(sh, 17) +
4454 2.04049226678e11 * pow(sh, 18) -
4455 1.8515048904e10 * pow(sh, 19) +
4456 864864. * pow(sh, 20) + 769608. * pow(sh, 21) -
4457 13608. * pow(sh, 22) - 222264. * pow(sh, 23) -
4458 69120. * pow(sh, 24) - 9000. * pow(sh, 25) +
4459 19215. * pow(sh, 26) + 315. * pow(sh, 27)) -
4460 4.265084782888612e-12 * pow(-1 + sh, -10) * pow(sh, -18)*
4461 (3.25909584e8 + 3.86898512e8 * sh -
4462 1.7890959086e10 * pow(sh, 2) -
4463 1.0389775539e11 * pow(sh, 3) -
4464 1.8235127911e11 * pow(sh, 4) +
4465 2.39868807178e11 * pow(sh, 5) +
4466 1.439309125254e12 * pow(sh, 6) +
4467 1.482637674518e12 * pow(sh, 7) +
4468 1.1773533109183245e17 * pow(sh, 8) -
4469 1.1762786134932733e18 * pow(sh, 9) +
4470 5.287945685290232e18 * pow(sh, 10) -
4471 1.408573123970191e19 * pow(sh, 11) +
4472 2.461891448027823e19 * pow(sh, 12) -
4473 2.949733724851667e19 * pow(sh, 13) +
4474 2.4532821652619735e19 * pow(sh, 14) -
4475 1.3981789937805554e19 * pow(sh, 15) +
4476 5.223681256860368e18 * pow(sh, 16) -
4477 1.1543710933218002e18 * pow(sh, 17) +
4478 1.1493955705739667e17 * pow(sh, 18) -
4479 4.556468500926695e15 * pow(sh, 19) +
4480 1.7572146654042552e16 * pow(sh, 20) -
4481 3.9615878221273704e16 * pow(sh, 21) +
4482 5.638676809242492e16 * pow(sh, 22) -
4483 4.406121643569693e16 * pow(sh, 23) +
4484 3.2318412281568824e16 * pow(sh, 24) -
4485 4.259492022110168e17 * pow(sh, 25) +
4486 2.1201430825186813e18 * pow(sh, 26) -
4487 2.4923037486423537e18 * pow(sh, 27) -
4488 1.7838580318483542e19 * pow(sh, 28) +
4489 1.0701104247446746e20 * pow(sh, 29) -
4490 3.1929857696090705e20 * pow(sh, 30) +
4491 6.474487284281146e20 * pow(sh, 31) -
4492 9.760291519177703e20 * pow(sh, 32) +
4493 1.1392452492311442e21 * pow(sh, 33) -
4494 1.0520263768751638e21 * pow(sh, 34) +
4495 7.775724274489533e20 * pow(sh, 35) -
4496 4.624318941820542e20 * pow(sh, 36) +
4497 2.213656928475123e20 * pow(sh, 37) -
4498 8.494727327998501e19 * pow(sh, 38) +
4499 2.5900145606169035e19 * pow(sh, 39) -
4500 6.181939718636802e18 * pow(sh, 40) +
4501 1.1286121824774998e18 * pow(sh, 41) -
4502 1.5195925385464826e17 * pow(sh, 42) +
4503 1.420123992415128e16 * pow(sh, 43) -
4504 8.2168151561496e14 * pow(sh, 44) +
4505 2.214605244852e13 * pow(sh, 45))) -
4506 6.962130651845887e-19 * pow(-1 + sh, -11) * pow(sh, -18)*
4507 (-3.16414088198208e14 + 9.74539823921664e14 * sh +
4508 8.087033764735926e15 * pow(sh, 2) -
4509 5.684694862939435e16 * pow(sh, 3) -
4510 4.0319420043801754e17 * pow(sh, 4) -
4511 4.882604186351892e17 * pow(sh, 5) +
4512 1.9325128069129482e18 * pow(sh, 6) +
4513 6.292691026868566e18 * pow(sh, 7) +
4514 2.792081386627729e23 * pow(sh, 8) -
4515 3.0653776012724437e24 * pow(sh, 9) +
4516 1.5296980754272814e25 * pow(sh, 10) -
4517 4.580270398519267e25 * pow(sh, 11) +
4518 9.143594067412906e25 * pow(sh, 12) -
4519 1.2779495928946286e26 * pow(sh, 13) +
4520 1.2762486136562552e26 * pow(sh, 14) -
4521 9.110246500944438e25 * pow(sh, 15) +
4522 4.558305919883703e25 * pow(sh, 16) -
4523 1.5244656041943856e25 * pow(sh, 17) +
4524 3.073171046299431e24 * pow(sh, 18) -
4525 2.6184927906686532e23 * pow(sh, 19) -
4526 1.0487308963151102e23 * pow(sh, 20) +
4527 2.7045906540143195e23 * pow(sh, 21) -
4528 4.4609820546041925e23 * pow(sh, 22) +
4529 4.3341866014896366e23 * pow(sh, 23) -
4530 3.0284128932535685e23 * pow(sh, 24) +
4531 2.8502065219122033e24 * pow(sh, 25) -
4532 1.6744213714265902e25 * pow(sh, 26) +
4533 3.4546092096924117e25 * pow(sh, 27) +
4534 6.6958214497636675e25 * pow(sh, 28) -
4535 6.719624033223081e26 * pow(sh, 29) +
4536 2.3634332855143025e27 * pow(sh, 30) -
4537 5.406861258181357e27 * pow(sh, 31) +
4538 9.10419041097793e27 * pow(sh, 32) -
4539 1.1865952156817548e28 * pow(sh, 33) +
4540 1.2283961231268607e28 * pow(sh, 34) -
4541 1.0245250602869844e28 * pow(sh, 35) +
4542 6.935097181125567e27 * pow(sh, 36) -
4543 3.819594455600949e27 * pow(sh, 37) +
4544 1.7090015768187537e27 * pow(sh, 38) -
4545 6.177750381749111e26 * pow(sh, 39) +
4546 1.7862394817725466e26 * pow(sh, 40) -
4547 4.0667139734465495e25 * pow(sh, 41) +
4548 7.117942930727479e24 * pow(sh, 42) -
4549 9.229373088621336e23 * pow(sh, 43) +
4550 8.339200912205272e22 * pow(sh, 44) -
4551 4.6814118287209613e21 * pow(sh, 45) +
4552 1.2279937968685318e20 * pow(sh, 46) -
4553 6.2548910424e12 * pow(3.141592653589793, 2)*
4554 (-1512. - 14784. * sh - 50316. * pow(sh, 2) -
4555 24288. * pow(sh, 3) + 294600. * pow(sh, 4) +
4556 761376. * pow(sh, 5) + 161112. * pow(sh, 6) -
4557 2.077152e6 * pow(sh, 7) +
4558 5.7064210896e10 * pow(sh, 8) -
4559 6.27572398992e11 * pow(sh, 9) +
4560 3.137096776433e12 * pow(sh, 10) -
4561 9.409743553774e12 * pow(sh, 11) +
4562 1.8818560922504e13 * pow(sh, 12) -
4563 2.6347616980552e13 * pow(sh, 13) +
4564 2.6351485677944e13 * pow(sh, 14) -
4565 1.8825939544384e13 * pow(sh, 15) +
4566 9.414254480792e12 * pow(sh, 16) -
4567 3.137901366712e12 * pow(sh, 17) +
4568 6.27172370726e11 * pow(sh, 18) -
4569 5.6235529592e10 * pow(sh, 19) -
4570 3.03847082e9 * pow(sh, 20) +
4571 7.83321994e9 * pow(sh, 21) -
4572 1.36334135e10 * pow(sh, 22) +
4573 2.612527414e10 * pow(sh, 23) -
4574 1.2043334018e11 * pow(sh, 24) +
4575 5.11152445252e11 * pow(sh, 25) -
4576 7.00772118731e11 * pow(sh, 26) -
4577 5.082162807869e12 * pow(sh, 27) +
4578 3.7496386828908e13 * pow(sh, 28) -
4579 1.3686666259266e14 * pow(sh, 29) +
4580 3.4205865415722e14 * pow(sh, 30) -
4581 6.4118721989322e14 * pow(sh, 31) +
4582 9.3979081997658e14 * pow(sh, 32) -
4583 1.10202690808782e15 * pow(sh, 33) +
4584 1.04788397968254e15 * pow(sh, 34) -
4585 8.141721939351e14 * pow(sh, 35) +
4586 5.1869434238622e14 * pow(sh, 36) -
4587 2.7096438795426e14 * pow(sh, 37) +
4588 1.1569642322238e14 * pow(sh, 38) -
4589 4.010774591538e13 * pow(sh, 39) +
4590 1.116674149734e13 * pow(sh, 40) -
4591 2.4564705165e12 * pow(sh, 41) +
4592 4.1666228604e11 * pow(sh, 42) -
4593 5.249021778e10 * pow(sh, 43) +
4594 4.61837376e9 * pow(sh, 44) -
4595 2.5297272e8 * pow(sh, 45) + 6.48648e6 * pow(sh, 46)))) *
4598 res += pow(z, 10)*(-2.2176369034656202e6 -
4599 2.1313657032968003e7 * (1. - sh) -
4600 1.1510510548217365e8 * pow(1. - sh, 2) -
4601 4.5152247350887704e8 * pow(1. - sh, 3) -
4602 1.4313252283712986e9 * pow(1. - sh, 4) +
4603 log(muh)*(22436.0126984127 +
4604 242664.12698412698 * (1. - sh) +
4605 1.4261726984126985e6 * pow(1. - sh, 2) +
4606 6.040264126984127e6 * pow(1. - sh, 3) +
4607 2.0637578412698414e7 * pow(1. - sh, 4)) +
4608 log(z)*(481985.82617224846 +
4609 5.600617111968864e6 * (1. - sh) +
4610 3.2763678189123973e7 * pow(1. - sh, 2) +
4611 1.3890575176721412e8 * pow(1. - sh, 3) +
4612 4.790156868849013e8 * pow(1. - sh, 4) +
4613 log(muh)*(18304. + 183040. * (1. - sh) +
4614 1.00672e6 * pow(1. - sh, 2) +
4615 4.02688e6 * pow(1. - sh, 3) +
4616 1.308736e7 * pow(1. - sh, 4))) +
4617 (245966.3111111111 +
4618 2.4484144197530863e6 * (1. - sh) +
4619 1.3337867328395061e7 * pow(1. - sh, 2) +
4620 5.3335035751205176e7 * pow(1. - sh, 3) +
4621 1.733848092131687e8 * pow(1. - sh, 4)) * pow(log(z), 2)