4462 const double s_sqrt = sqrt(s);
4463 const double spin_deg_factor = 1.0;
4466 double tmax = mandelstam_t[0];
4467 double tmin = mandelstam_t[1];
4473 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4476 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4479 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
4480 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
4482 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
4483 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4487 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
4488 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4490 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
4492 (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
4493 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
4494 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4497 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4498 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s))) +
4502 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
4504 (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
4505 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
4506 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4509 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4510 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
4511 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
4514 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
4515 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
4517 (1. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s))))) /
4518 (1. * pow(
a1_mass, 2) - 1. * tmin) +
4520 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
4523 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
4524 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s -
4530 0.5 *
eta2 * pow(m_rho, 2) + 0.5 *
eta1 * s - 1. *
eta2 * s) *
4534 (12. + 1. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
4535 delta * (-8. + 8. *
C4 * pow(m_rho, 2))) +
4536 (-4. - 3. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
4537 delta * (8. + 8. *
C4 * pow(m_rho, 2))) *
4553 0.5 * pow(m_rho, 2) + 1. * s)) *
4555 (0.25 * (-2. + 1. *
delta) *
4556 (-8. *
C4 * pow(m_rho, 4) +
4558 (-2. - 3. *
delta) * s +
4559 pow(m_rho, 2) * (2. + 1. *
delta + 16. *
C4 * s)) *
4563 (32 * pow(
C4, 2) * pow(m_rho, 8) + 2 * pow(
delta, 2) * pow(s, 2) +
4564 8 *
C4 * pow(m_rho, 6) * (-6 +
delta - 8 *
C4 * s) +
4565 2 *
delta * pow(m_rho, 2) * s * (-6 +
delta - 8 *
C4 * s) +
4566 pow(m_rho, 4) * (12 - pow(
delta, 2) + 8 *
C4 * (6 +
delta) * s +
4567 32 * pow(
C4, 2) * pow(s, 2))) *
4572 (0.75 * pow(m_rho, 4) - 0.125 *
delta * pow(m_rho, 4) -
4573 1. *
C4 * pow(m_rho, 6) +
4575 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4577 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) -
4578 0.25 * pow(m_rho, 2) * s - 0.375 *
delta * pow(m_rho, 2) * s +
4579 2. *
C4 * pow(m_rho, 4) * s + 0.25 *
delta * pow(s, 2) -
4580 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
4581 eta1 * (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
4583 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) +
4585 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4586 0.25 *
delta * pow(s, 2) +
4587 1. *
C4 * pow(m_rho, 2) * pow(s, 2))) *
4590 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
4593 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
4594 2. * pow(m_rho, 4) +
4596 pow(m_rho, 2) - 1. * s) -
4597 3. * pow(m_rho, 2) * s + pow(s, 2)) +
4600 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
4602 pow(m_rho, 2) - 1. * s) +
4603 pow(m_rho, 2) * s + pow(s, 2)) +
4607 4. * pow(
pion_mass, 2) * pow(m_rho, 2) +
4608 2. * pow(m_rho, 4) - 2. * pow(s, 2) +
4610 2. * pow(m_rho, 2) + 2. * s))) *
4612 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4615 2. * pow(m_rho, 4) +
4617 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4622 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
4623 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4625 4. * pow(m_rho, 2) + 4. * s)) +
4629 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
4630 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4632 4. * pow(m_rho, 2) + 4. * s))) *
4634 (0.125 * (-2. + 1. *
delta) *
4635 (2. + 1. *
delta - 8. *
C4 * pow(m_rho, 2)) * pow(tmin, 2)) /
4637 0.5 * pow(1. *
eta1 - 1. *
eta2, 2) *
4638 (-0.5 + 1. *
C4 * pow(m_rho, 2)) * pow(tmin, 2) -
4640 (0.5 - 0.125 * pow(
delta, 2) - 2. *
C4 * pow(m_rho, 2) +
4641 1. *
C4 *
delta * pow(m_rho, 2)) *
4644 0.0625 * pow(1. *
eta1 - 1. *
eta2, 4) *
4645 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 0.5 * s) *
4647 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
4649 1. * pow(m_rho, 2) - 1. * s) +
4651 1. * pow(m_rho, 2) + s)) *
4653 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) * pow(tmin, 3) -
4654 0.020833333333333332 * pow(1. *
eta1 - 1. *
eta2, 4) * pow(tmin, 3) +
4655 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4659 2. * pow(m_rho, 4) +
4661 4. * pow(m_rho, 2) * s - 4. * pow(s, 2)) +
4665 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
4666 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4668 4. * pow(m_rho, 2) + 4. * s)) +
4672 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
4673 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4675 4. * pow(m_rho, 2) + 4. * s))) *
4676 (-1. * pow(m_rho, 2) + s + tmin) -
4677 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
4678 (
eta2 * (-1. * pow(
a1_mass, 2) - 1. * pow(m_rho, 2) - 1. * s) +
4679 eta1 * (pow(
a1_mass, 2) - 1. * pow(m_rho, 2) + s)) *
4680 pow(-1. * pow(m_rho, 2) + s + tmin, 2) +
4681 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) *
4682 pow(-1. * pow(m_rho, 2) + s + tmin, 3) +
4684 (-1. + 2. *
C4 * pow(m_rho, 2)) *
4686 1. * pow(m_rho, 2) + s + tmin,
4689 (
eta2 * (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
4690 0.5 *
C4 * pow(m_rho, 6) +
4692 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4694 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4695 0.125 * pow(m_rho, 2) * s -
4696 0.1875 *
delta * pow(m_rho, 2) * s +
4697 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
4698 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
4699 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
4701 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4703 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4704 0.125 *
delta * pow(s, 2) +
4705 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
4707 1. * s + 1. * tmin)) /
4710 (
eta2 * (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
4711 0.5 *
C4 * pow(m_rho, 6) +
4713 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4715 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4716 0.125 * pow(m_rho, 2) * s -
4717 0.1875 *
delta * pow(m_rho, 2) * s +
4718 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
4719 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
4720 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
4722 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4724 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4725 0.125 *
delta * pow(s, 2) +
4726 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
4735 (-1. * pow(
a1_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
4737 (2. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
4738 pow(
pion_mass, 4) * (-1.5 * pow(m_rho, 2) + 1. * s) +
4740 (-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 2. * s)) +
4747 (1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
4748 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) +
4749 1. * pow(m_rho, 2) * s))) *
4758 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
4759 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
4761 (2. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) - 1.5 * s) +
4763 3. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4764 0.5 * pow(m_rho, 4) * s - 2. * pow(
pion_mass, 2) * pow(s, 2) -
4765 1. * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
4769 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
4770 1. * pow(m_rho, 2) * s +
4771 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 1. * s))) +
4774 4. * pow(
pion_mass, 4) * pow(m_rho, 2) +
4775 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
4780 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
4781 0.5 * pow(m_rho, 4) * s + 3. * pow(
pion_mass, 2) * pow(s, 2) +
4782 1. * pow(m_rho, 2) * pow(s, 2) - 0.5 * pow(s, 3) +
4785 (1. * pow(m_rho, 2) - 1. * s) * s +
4786 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 3. * s)))) *
4792 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
4793 2. * pow(m_rho, 2) * s + pow(s, 2) +
4795 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) +
4796 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4802 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
4804 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
4805 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
4806 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4809 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4810 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s)) +
4813 1. * pow(m_rho, 4) +
4815 6. * pow(m_rho, 2) - 6. * s) +
4816 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
4817 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s))) +
4821 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4824 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
4825 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
4827 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
4828 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4832 2. * pow(m_rho, 4) - 8. * pow(
pion_mass, 2) * s -
4833 4. * pow(m_rho, 2) * s + 4. * pow(s, 2) +
4838 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
4840 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
4841 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
4842 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4845 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4846 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
4847 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
4850 1. * pow(m_rho, 4) +
4852 6. * pow(m_rho, 2) - 6. * s) +
4853 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
4854 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s)) +
4857 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
4858 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
4860 (pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)))) *
4869 28. * pow(
pion_mass, 6) * pow(m_rho, 2) +
4870 20. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4871 10. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) -
4873 16. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
4874 12. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
4875 4. * pow(m_rho, 6) * s - 4. * pow(
pion_mass, 4) * pow(s, 2) -
4876 6. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
4878 4. * pow(m_rho, 2) * pow(s, 3) - 2. * pow(s, 4) +
4880 28. * pow(m_rho, 2) + 28. * s) +
4882 (84. * pow(
pion_mass, 4) + 24. * pow(m_rho, 4) +
4883 pow(
pion_mass, 2) * (84. * pow(m_rho, 2) - 72. * s) -
4884 36. * pow(m_rho, 2) * s + 12. * pow(s, 2)) +
4888 4. * pow(m_rho, 2) - 4. * s) +
4889 (4. * pow(m_rho, 2) - 4. * s) * s +
4890 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 8. * s)) +
4892 (-56. * pow(
pion_mass, 6) - 10. * pow(m_rho, 6) +
4893 18. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) -
4895 pow(
pion_mass, 4) * (-84. * pow(m_rho, 2) + 60. * s) +
4897 (-48. * pow(m_rho, 4) + 60. * pow(m_rho, 2) * s -
4898 12. * pow(s, 2)))) +
4902 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
4903 7. * pow(
pion_mass, 4) * pow(m_rho, 4) -
4904 2. * pow(
pion_mass, 2) * pow(m_rho, 6) +
4906 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
4908 11. * pow(
pion_mass, 4) * pow(m_rho, 2) * s +
4909 6. * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
4910 1. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) -
4911 1. * pow(m_rho, 4) * pow(s, 2) -
4913 1. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
4916 1. * pow(m_rho, 4) +
4917 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
4918 1. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4920 2. * pow(m_rho, 2) + 2. * s)) +
4922 (-42. * pow(
pion_mass, 4) - 9. * pow(m_rho, 4) +
4923 21. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
4924 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
4926 (28. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
4927 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) -
4928 9. * pow(m_rho, 4) * s + 6. * pow(m_rho, 2) * pow(s, 2) +
4931 (18. * pow(m_rho, 4) - 36. * pow(m_rho, 2) * s +
4936 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
4937 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4938 6. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) +
4940 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
4942 1. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
4943 16. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
4944 7. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) +
4945 14. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
4946 9. * pow(m_rho, 4) * pow(s, 2) -
4948 5. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
4951 3. * pow(m_rho, 4) +
4952 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
4953 5. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4955 2. * pow(m_rho, 2) + 2. * s)) +
4957 (-42. * pow(
pion_mass, 4) - 3. * pow(m_rho, 4) +
4958 9. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
4959 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
4961 (28. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
4962 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) +
4963 9. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) +
4966 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
4967 6. * pow(s, 2))))) *
4973 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
4974 2. * pow(m_rho, 2) * s + pow(s, 2) +
4976 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)) +
4977 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
4982 (4. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) -
4983 2. * pow(
pion_mass, 2) * s + 2. * pow(m_rho, 2) * s) +
4985 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
4986 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4991 pow(m_rho, 2) * (pow(m_rho, 2) - 1. * s) * s +
4992 pow(
pion_mass, 4) * (-3. * pow(m_rho, 2) + s) +
4994 3. * pow(m_rho, 2) + 3. * s) +
4996 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
4997 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
4998 4. * pow(m_rho, 2) * s + 2. * pow(s, 2))) +
5001 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5002 pow(
pion_mass, 4) * (3. * pow(m_rho, 2) + s) +
5004 3. * pow(m_rho, 2) + 3. * s) +
5006 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
5007 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
5008 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)))) *
5009 log(abs(-1. * pow(
a1_mass, 2) + tmin)) -
5013 0.5 * pow(
pion_mass, 4) * pow(m_rho, 2) +
5015 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5017 (0.5 * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s)) +
5020 (1. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
5021 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
5022 0.5 * pow(m_rho, 2) * s) +
5024 (-2. * pow(
pion_mass, 4) - 0.5 * pow(m_rho, 2) * s +
5025 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
5026 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5031 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2)) +
5032 pow(
pion_mass, 4) * (0.5 * pow(m_rho, 2) - 1. * s) +
5033 pow(
pion_mass, 2) * (-0.5 * pow(m_rho, 2) + 0.5 * s) * s +
5036 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5037 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s)) +
5040 pow(
pion_mass, 2) * (1.5 * pow(m_rho, 2) - 2. * s) * s +
5041 (-0.5 * pow(m_rho, 2) + 0.5 * s) * pow(s, 2) +
5042 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1.5 * s) +
5044 (-2. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 2) * s +
5045 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
5046 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5048 1. * pow(m_rho, 2) + 1. * s) -
5049 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5051 0.5 * pow(m_rho, 2) + 0.5 * s) *
5055 (8. * pow(
pion_mass, 2) + 4. * pow(m_rho, 2) - 4. * s) +
5057 (-12. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) +
5058 4. * pow(m_rho, 2) * s +
5059 pow(
pion_mass, 2) * (-12. * pow(m_rho, 2) + 8. * s)) +
5062 4. * pow(
pion_mass, 4) * pow(m_rho, 2) -
5063 2. * pow(m_rho, 6) + 4. * pow(m_rho, 4) * s -
5064 2. * pow(m_rho, 2) * pow(s, 2) +
5066 (-8. * pow(m_rho, 4) + 8. * pow(m_rho, 2) * s)) +
5068 (8. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
5069 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 4. * s) -
5070 2. * pow(m_rho, 4) * s - 2. * pow(m_rho, 2) * pow(s, 2) +
5073 (8. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s -
5078 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
5080 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
5081 pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
5083 (6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5084 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) +
5085 pow(m_rho, 2) * s) +
5087 (-4. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) -
5088 5. * pow(m_rho, 4) * s + 4. * pow(m_rho, 2) * pow(s, 2) -
5090 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
5092 (2. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s +
5096 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
5098 (-2. * pow(m_rho, 2) + 2. * s) +
5100 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
5102 (3. * pow(m_rho, 4) - 5. * pow(m_rho, 2) * s) +
5104 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
5105 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
5106 3. * pow(m_rho, 2) * s) +
5108 (-4. * pow(
pion_mass, 6) + pow(m_rho, 4) * s -
5110 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
5112 (-2. * pow(m_rho, 4) + 4. * pow(m_rho, 2) * s +
5113 2. * pow(s, 2))))) *
5114 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5117 0.25 * pow(m_rho, 4) - 1. * pow(
pion_mass, 2) * s -
5118 0.5 * pow(m_rho, 2) * s + 0.25 * pow(s, 2) +
5120 (-2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 1. * s)) -
5124 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5126 (pow(
pion_mass, 2) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
5127 (-0.25 + 0.125 *
delta) * (pow(m_rho, 2) + s)) +
5129 (-1. *
C4 * pow(m_rho, 6) +
5131 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) +
5132 0.25 *
delta * pow(s, 2) +
5134 (-0.25 - 0.375 *
delta - 1. *
C4 * s) +
5135 pow(m_rho, 4) * (0.75 - 0.125 *
delta + 2. *
C4 * s))) +
5138 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
5140 (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
5142 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
5143 0.25 *
delta * pow(s, 2) +
5144 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5146 (pow(
pion_mass, 4) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
5147 s * ((0.25 - 0.125 *
delta) * pow(m_rho, 2) +
5148 (-0.25 + 0.125 *
delta) * s) +
5150 (2. *
C4 * pow(m_rho, 4) + (0.5 + 0.25 *
delta) * s +
5151 pow(m_rho, 2) * (-1. - 2. *
C4 * s))))) *
5152 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5155 log(abs(-1. * pow(
pion_mass, 2) + tmin)) +
5158 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
5160 5.999999999999999 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 2) *
5162 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 2) * pow(s, 2)) *
5163 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5164 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
5166 (0. +
eta2 * pow(
pion_mass, 2) * (4. * pow(m_rho, 2) - 4. * s) +
5167 eta1 * (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s +
5168 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 4. * s))) *
5169 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5171 (2. * (-2. + 1. *
delta) *
5172 (0. + (-0.25 + 0.125 *
delta) * pow(m_rho, 2) * s +
5174 pow(m_rho, 2) * (1. + 2. *
C4 * s))) *
5175 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5178 (
eta1 * (pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1. * s) +
5180 (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5181 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
5183 (-1. * pow(m_rho, 4) + 2.5 * pow(m_rho, 2) * s -
5185 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
5187 eta2 * (0.5 * pow(m_rho, 6) +
5188 pow(
pion_mass, 4) * (1. * pow(m_rho, 2) - 1. * s) -
5189 1.5 * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) -
5192 (1.5 * pow(m_rho, 4) - 3. * pow(m_rho, 2) * s +
5195 (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
5197 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
5198 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5201 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
5202 2. * pow(m_rho, 2) * s + pow(s, 2) +
5204 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
5206 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s +
5211 (-1. * pow(m_rho, 2) + 1. * s) +
5213 (-0.5 * pow(m_rho, 4) +
5214 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5215 0.5 * pow(m_rho, 2) * s) +
5217 (0.5 * pow(m_rho, 4) - 0.5 * pow(m_rho, 2) * s +
5218 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
5219 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
5223 (
eta1 * (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5224 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
5225 eta2 * (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
5227 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
5228 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
5230 1. * pow(m_rho, 2) + 1. * s) -
5233 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
5235 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 4) * s +
5237 (
C4 * (32. - 16. *
delta) * pow(m_rho, 6) +
5240 (-8. + 24. *
delta - 10. * pow(
delta, 2) + 32. *
C4 * s -
5242 pow(m_rho, 4) * (-16. + 8. *
delta - 64. *
C4 * s +
5244 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
5245 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
5246 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5254 (-4. * pow(
pion_mass, 4) + 4. * pow(m_rho, 4) +
5255 2. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s) +
5257 (12. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) -
5258 8. * pow(
pion_mass, 2) * s - 4. * pow(m_rho, 2) * s +
5262 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
5264 (6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) - 3. * s) -
5266 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
5267 1. * pow(m_rho, 4) * s + pow(m_rho, 2) * pow(s, 2) +
5270 1. * pow(m_rho, 2) + s) +
5272 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5273 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
5274 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s))) +
5277 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
5279 (6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) - 3. * s) -
5284 pow(m_rho, 2) + s) +
5286 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5287 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
5288 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s)))) *
5292 2. * pow(
a1_mass, 2) * (-1. * pow(m_rho, 2) + s + tmin) -
5294 (-1. * pow(m_rho, 2) + s + tmin) +
5295 pow(-1. * pow(m_rho, 2) + s + tmin, 2))) -
5298 (0.5 - 1. *
C4 * pow(m_rho, 2)) +
5300 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
5302 (pow(
pion_mass, 2) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
5303 (0.25 - 0.125 *
delta) * (pow(m_rho, 2) + s)) +
5305 (1. *
C4 * pow(m_rho, 6) +
5307 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
5308 0.25 *
delta * pow(s, 2) +
5309 pow(m_rho, 4) * (-0.75 + 0.125 *
delta - 2. *
C4 * s) +
5311 (0.25 + 0.375 *
delta + 1. *
C4 * s))) +
5313 (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
5315 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5317 (-0.5 * pow(m_rho, 4) + 1. *
C4 * pow(m_rho, 6) +
5319 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) +
5320 0.25 *
delta * pow(s, 2) -
5321 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5323 (pow(
pion_mass, 4) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
5324 s * ((-0.25 + 0.125 *
delta) * pow(m_rho, 2) +
5325 (0.25 - 0.125 *
delta) * s) +
5327 (-2. *
C4 * pow(m_rho, 4) +
5328 (-0.5 - 0.25 *
delta) * s +
5329 pow(m_rho, 2) * (1. + 2. *
C4 * s))))) *
5332 1. * pow(m_rho, 2) + s + tmin,
5340 0.5 * pow(
pion_mass, 6) * pow(m_rho, 2) +
5343 (1. * pow(
pion_mass, 2) + 1.5 * pow(m_rho, 2) -
5346 0.5 * pow(m_rho, 2) - 1. * s)) +
5348 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5350 (-1. * pow(
pion_mass, 2) - 0.5 * pow(m_rho, 2) + 1. * s)) +
5354 1. * pow(
pion_mass, 2) * pow(m_rho, 4) -
5358 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
5359 0.5 * pow(m_rho, 2) * s) +
5361 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5362 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 0.5 * s) +
5363 0.5 * pow(m_rho, 2) * s) +
5367 0.5 * pow(m_rho, 2) * s +
5368 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1.5 * s)))) *
5371 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5372 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
5373 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
5374 2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
5376 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
5382 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5385 4. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
5386 1.5 * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
5388 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
5389 0.5 * pow(m_rho, 4) * pow(s, 2) -
5391 1. * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
5395 0.5 * pow(
pion_mass, 2) * s + 0.5 * pow(s, 2)) +
5398 (1. * pow(m_rho, 2) - 0.5 * s) * s +
5399 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 2.5 * s)) +
5402 pow(
pion_mass, 4) * (2. * pow(m_rho, 2) - 1.5 * s) -
5403 0.5 * pow(m_rho, 4) * s + 0.5 * pow(s, 3) +
5405 (1. * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s -
5411 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 0.5 * s) +
5413 pow(
a1_mass, 4) * (-0.5 * pow(m_rho, 4) +
5414 (-0.5 * pow(
pion_mass, 2) + 0.5 * s) * s) +
5416 (-0.5 * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
5419 (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
5422 (1. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 4) +
5423 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) -
5424 1. * pow(m_rho, 2) * s + 0.5 * pow(s, 2) +
5426 1. * pow(m_rho, 2) + 1.5 * s)) +
5429 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 0.5 * s) +
5431 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
5433 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
5434 0.5 * pow(s, 2))))) *
5437 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5438 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
5439 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
5440 2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
5442 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
5446 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
5447 2. * pow(m_rho, 2) * s + pow(s, 2) +
5449 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
5450 (0.0625 * pow(
eta1 - 1. *
eta2, 2) *
5454 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
5457 0.5 * pow(m_rho, 2) + 0.5 * s) +
5459 (10. * pow(
pion_mass, 6) - 2.5 * pow(m_rho, 6) +
5460 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) +
5461 6. * pow(m_rho, 4) * s - 4.5 * pow(m_rho, 2) * pow(s, 2) +
5465 (1. * pow(m_rho, 2) - 1. * s) * s +
5466 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
5468 (1. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
5469 pow(
pion_mass, 4) * (2.5 * pow(m_rho, 2) - 0.5 * s) -
5470 1. * pow(m_rho, 4) * s + 0.5 * pow(m_rho, 2) * pow(s, 2) +
5472 (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)) +
5475 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
5476 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
5477 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
5479 6. * pow(m_rho, 2) + 6. * s) +
5481 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
5484 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
5485 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
5486 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
5488 (-5. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) -
5489 3.5 * pow(m_rho, 6) * s +
5490 4.5 * pow(m_rho, 4) * pow(s, 2) -
5491 2.5 * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
5492 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
5494 (-2. * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s +
5497 (3. * pow(m_rho, 6) - 8. * pow(m_rho, 4) * s +
5498 7. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)))) +
5502 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
5505 0.5 * pow(m_rho, 2) + 0.5 * s) +
5507 (-10. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) +
5508 5. * pow(m_rho, 2) * s - 1. * pow(s, 2) +
5509 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
5511 (10. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
5512 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) -
5513 3. * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) +
5516 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s)) +
5519 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
5520 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
5521 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
5523 6. * pow(m_rho, 2) + 6. * s) +
5525 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
5528 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
5529 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
5530 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
5533 pow(
pion_mass, 6) * (2.5 * pow(m_rho, 2) - 0.5 * s) +
5535 (4. * pow(m_rho, 4) - 6. * pow(m_rho, 2) * s) +
5537 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
5540 (1.5 * pow(m_rho, 6) - 6. * pow(m_rho, 4) * s +
5541 4.5 * pow(m_rho, 2) * pow(s, 2))) +
5544 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
5546 (-8. * pow(m_rho, 4) + 13. * pow(m_rho, 2) * s +
5549 (-1. * pow(m_rho, 6) + 6. * pow(m_rho, 4) * s -
5550 3. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)) +
5551 s * (0.5 * pow(m_rho, 6) - 0.5 * pow(m_rho, 4) * s -
5552 0.5 * pow(m_rho, 2) * pow(s, 2) +
5553 0.5 * pow(s, 3)))) +
5558 1. * pow(m_rho, 2) - 1. * s) +
5560 (-10. * pow(
pion_mass, 2) - 5. * pow(m_rho, 2) + 5. * s) +
5562 (20. * pow(
pion_mass, 4) + 6. * pow(m_rho, 4) +
5563 pow(
pion_mass, 2) * (20. * pow(m_rho, 2) - 16. * s) -
5564 8. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5566 (-20. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
5567 6. * pow(m_rho, 4) * s - 2. * pow(s, 3) +
5568 pow(
pion_mass, 4) * (-30. * pow(m_rho, 2) + 18. * s) +
5570 (-18. * pow(m_rho, 4) + 18. * pow(m_rho, 2) * s)) +
5572 (-2. * pow(
pion_mass, 8) - 1. * pow(m_rho, 8) +
5573 3. * pow(m_rho, 6) * s - 3. * pow(m_rho, 4) * pow(s, 2) +
5574 1. * pow(m_rho, 2) * pow(s, 3) +
5575 pow(
pion_mass, 6) * (-5. * pow(m_rho, 2) + 1. * s) +
5577 (-10. * pow(m_rho, 4) + 10. * pow(m_rho, 2) * s) +
5579 (-6. * pow(m_rho, 6) + 12. * pow(m_rho, 4) * s -
5580 6. * pow(m_rho, 2) * pow(s, 2))) +
5582 (10. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) +
5583 pow(
pion_mass, 6) * (20. * pow(m_rho, 2) - 8. * s) -
5584 2. * pow(m_rho, 6) * s + 2. * pow(m_rho, 2) * pow(s, 3) -
5587 (22. * pow(m_rho, 4) - 20. * pow(m_rho, 2) * s -
5590 (8. * pow(m_rho, 6) - 12. * pow(m_rho, 4) * s +
5594 1. * pow(m_rho, 6) +
5595 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 12. * s) +
5597 12. * pow(m_rho, 2) - 12. * s) -
5598 3. * pow(m_rho, 4) * s + 3. * pow(m_rho, 2) * pow(s, 2) -
5601 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
5604 (-24. * pow(
pion_mass, 4) - 6. * pow(m_rho, 4) +
5605 12. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
5607 (-24. * pow(m_rho, 2) + 24. * s))))) *
5610 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5611 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
5612 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
5613 2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
5615 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
5619 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
5620 2. * pow(m_rho, 2) * s + pow(s, 2) +
5622 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)))) /
5626 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5629 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5632 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
5633 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
5635 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
5636 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
5640 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
5641 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5643 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
5645 (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
5646 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
5647 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5650 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5651 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s))) +
5655 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
5657 (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
5658 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
5659 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5662 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5663 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
5664 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
5667 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
5668 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
5670 (1. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s))))) /
5671 (1. * pow(
a1_mass, 2) - 1. * tmax) +
5673 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
5676 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
5677 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s -
5683 0.5 *
eta2 * pow(m_rho, 2) + 0.5 *
eta1 * s - 1. *
eta2 * s) *
5687 (12. + 1. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
5688 delta * (-8. + 8. *
C4 * pow(m_rho, 2))) +
5689 (-4. - 3. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
5690 delta * (8. + 8. *
C4 * pow(m_rho, 2))) *
5705 0.5 * pow(m_rho, 2) + 1. * s)) *
5707 (0.25 * (-2. + 1. *
delta) *
5708 (-8. *
C4 * pow(m_rho, 4) +
5710 (-2. - 3. *
delta) * s +
5711 pow(m_rho, 2) * (2. + 1. *
delta + 16. *
C4 * s)) *
5715 (32 * pow(
C4, 2) * pow(m_rho, 8) + 2 * pow(
delta, 2) * pow(s, 2) +
5716 8 *
C4 * pow(m_rho, 6) * (-6 +
delta - 8 *
C4 * s) +
5717 2 *
delta * pow(m_rho, 2) * s * (-6 +
delta - 8 *
C4 * s) +
5718 pow(m_rho, 4) * (12 - pow(
delta, 2) + 8 *
C4 * (6 +
delta) * s +
5719 32 * pow(
C4, 2) * pow(s, 2))) *
5724 (0.75 * pow(m_rho, 4) - 0.125 *
delta * pow(m_rho, 4) -
5725 1. *
C4 * pow(m_rho, 6) +
5727 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5729 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) -
5730 0.25 * pow(m_rho, 2) * s - 0.375 *
delta * pow(m_rho, 2) * s +
5731 2. *
C4 * pow(m_rho, 4) * s + 0.25 *
delta * pow(s, 2) -
5732 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5733 eta1 * (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
5735 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) +
5737 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5738 0.25 *
delta * pow(s, 2) +
5739 1. *
C4 * pow(m_rho, 2) * pow(s, 2))) *
5742 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
5746 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + 2. * pow(m_rho, 4) +
5748 (2. * pow(
pion_mass, 2) + pow(m_rho, 2) - 1. * s) -
5749 3. * pow(m_rho, 2) * s + pow(s, 2)) +
5752 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
5754 pow(m_rho, 2) - 1. * s) +
5755 pow(m_rho, 2) * s + pow(s, 2)) +
5759 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + 2. * pow(m_rho, 4) -
5762 2. * pow(m_rho, 2) + 2. * s))) *
5764 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5767 2. * pow(m_rho, 4) +
5769 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
5774 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
5775 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5777 4. * pow(m_rho, 2) + 4. * s)) +
5781 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
5782 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5784 4. * pow(m_rho, 2) + 4. * s))) *
5786 (0.125 * (-2. + 1. *
delta) *
5787 (2. + 1. *
delta - 8. *
C4 * pow(m_rho, 2)) * pow(tmax, 2)) /
5789 0.5 * pow(1. *
eta1 - 1. *
eta2, 2) *
5790 (-0.5 + 1. *
C4 * pow(m_rho, 2)) * pow(tmax, 2) -
5792 (0.5 - 0.125 * pow(
delta, 2) - 2. *
C4 * pow(m_rho, 2) +
5793 1. *
C4 *
delta * pow(m_rho, 2)) *
5796 0.0625 * pow(1. *
eta1 - 1. *
eta2, 4) *
5797 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 0.5 * s) *
5799 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
5801 1. * pow(m_rho, 2) - 1. * s) +
5803 1. * pow(m_rho, 2) + s)) *
5805 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) * pow(tmax, 3) -
5806 0.020833333333333332 * pow(1. *
eta1 - 1. *
eta2, 4) * pow(tmax, 3) +
5807 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5811 2. * pow(m_rho, 4) +
5813 4. * pow(m_rho, 2) * s - 4. * pow(s, 2)) +
5817 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5818 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5820 4. * pow(m_rho, 2) + 4. * s)) +
5824 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5825 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5827 4. * pow(m_rho, 2) + 4. * s))) *
5828 (-1. * pow(m_rho, 2) + s + tmax) -
5829 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
5830 (
eta2 * (-1. * pow(
a1_mass, 2) - 1. * pow(m_rho, 2) - 1. * s) +
5831 eta1 * (pow(
a1_mass, 2) - 1. * pow(m_rho, 2) + s)) *
5832 pow(-1. * pow(m_rho, 2) + s + tmax, 2) +
5833 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) *
5834 pow(-1. * pow(m_rho, 2) + s + tmax, 3) +
5836 (-1. + 2. *
C4 * pow(m_rho, 2)) *
5842 (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
5843 0.5 *
C4 * pow(m_rho, 6) +
5845 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5847 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5848 0.125 * pow(m_rho, 2) * s - 0.1875 *
delta * pow(m_rho, 2) * s +
5849 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
5850 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5851 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
5853 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5855 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5856 0.125 *
delta * pow(s, 2) +
5857 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
5859 1. * s + 1. * tmax)) /
5863 (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
5864 0.5 *
C4 * pow(m_rho, 6) +
5866 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5868 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5869 0.125 * pow(m_rho, 2) * s - 0.1875 *
delta * pow(m_rho, 2) * s +
5870 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
5871 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5872 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
5874 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5876 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5877 0.125 *
delta * pow(s, 2) +
5878 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
5887 (-1. * pow(
a1_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5889 (2. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5890 pow(
pion_mass, 4) * (-1.5 * pow(m_rho, 2) + 1. * s) +
5892 (-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 2. * s)) +
5899 (1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5900 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) +
5901 1. * pow(m_rho, 2) * s))) *
5910 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
5912 (2. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) - 1.5 * s) +
5914 3. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5915 0.5 * pow(m_rho, 4) * s - 2. * pow(
pion_mass, 2) * pow(s, 2) -
5916 1. * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
5920 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5921 1. * pow(m_rho, 2) * s +
5922 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 1. * s))) +
5925 4. * pow(
pion_mass, 4) * pow(m_rho, 2) +
5926 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
5931 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
5932 0.5 * pow(m_rho, 4) * s + 3. * pow(
pion_mass, 2) * pow(s, 2) +
5933 1. * pow(m_rho, 2) * pow(s, 2) - 0.5 * pow(s, 3) +
5936 (1. * pow(m_rho, 2) - 1. * s) * s +
5937 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 3. * s)))) *
5943 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
5944 2. * pow(m_rho, 2) * s + pow(s, 2) +
5946 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) +
5947 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5953 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
5955 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
5956 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
5957 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5960 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5961 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s)) +
5964 1. * pow(m_rho, 4) +
5966 6. * pow(m_rho, 2) - 6. * s) +
5967 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
5968 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s))) +
5972 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5975 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
5976 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
5978 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
5979 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
5983 2. * pow(m_rho, 4) - 8. * pow(
pion_mass, 2) * s -
5984 4. * pow(m_rho, 2) * s + 4. * pow(s, 2) +
5989 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
5991 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
5992 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
5993 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5996 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5997 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
5998 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
6001 1. * pow(m_rho, 4) +
6003 6. * pow(m_rho, 2) - 6. * s) +
6004 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
6005 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s)) +
6008 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
6009 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
6011 (pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)))) *
6020 28. * pow(
pion_mass, 6) * pow(m_rho, 2) +
6021 20. * pow(
pion_mass, 4) * pow(m_rho, 4) +
6022 10. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) -
6024 16. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
6025 12. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
6026 4. * pow(m_rho, 6) * s - 4. * pow(
pion_mass, 4) * pow(s, 2) -
6027 6. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
6029 4. * pow(m_rho, 2) * pow(s, 3) - 2. * pow(s, 4) +
6031 (-56. * pow(
pion_mass, 2) - 28. * pow(m_rho, 2) + 28. * s) +
6033 (84. * pow(
pion_mass, 4) + 24. * pow(m_rho, 4) +
6034 pow(
pion_mass, 2) * (84. * pow(m_rho, 2) - 72. * s) -
6035 36. * pow(m_rho, 2) * s + 12. * pow(s, 2)) +
6039 4. * pow(m_rho, 2) - 4. * s) +
6040 (4. * pow(m_rho, 2) - 4. * s) * s +
6041 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 8. * s)) +
6043 (-56. * pow(
pion_mass, 6) - 10. * pow(m_rho, 6) +
6044 18. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) -
6046 pow(
pion_mass, 4) * (-84. * pow(m_rho, 2) + 60. * s) +
6048 (-48. * pow(m_rho, 4) + 60. * pow(m_rho, 2) * s -
6049 12. * pow(s, 2)))) +
6053 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
6054 7. * pow(
pion_mass, 4) * pow(m_rho, 4) -
6055 2. * pow(
pion_mass, 2) * pow(m_rho, 6) +
6057 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
6059 11. * pow(
pion_mass, 4) * pow(m_rho, 2) * s +
6060 6. * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
6061 1. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) -
6062 1. * pow(m_rho, 4) * pow(s, 2) -
6064 1. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
6067 1. * pow(m_rho, 4) +
6068 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
6069 1. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
6071 2. * pow(m_rho, 2) + 2. * s)) +
6073 (-42. * pow(
pion_mass, 4) - 9. * pow(m_rho, 4) +
6074 21. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
6075 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
6077 (28. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
6078 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) -
6079 9. * pow(m_rho, 4) * s + 6. * pow(m_rho, 2) * pow(s, 2) +
6082 (18. * pow(m_rho, 4) - 36. * pow(m_rho, 2) * s +
6087 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
6088 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
6089 6. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) +
6091 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
6093 1. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
6094 16. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
6095 7. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) +
6096 14. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
6097 9. * pow(m_rho, 4) * pow(s, 2) -
6099 5. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
6102 3. * pow(m_rho, 4) +
6103 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
6104 5. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
6106 2. * pow(m_rho, 2) + 2. * s)) +
6108 (-42. * pow(
pion_mass, 4) - 3. * pow(m_rho, 4) +
6109 9. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
6110 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
6112 (28. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
6113 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) +
6114 9. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) +
6117 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
6118 6. * pow(s, 2))))) *
6124 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
6125 2. * pow(m_rho, 2) * s + pow(s, 2) +
6127 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)) +
6128 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
6133 (4. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) -
6134 2. * pow(
pion_mass, 2) * s + 2. * pow(m_rho, 2) * s) +
6136 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
6137 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
6142 pow(m_rho, 2) * (pow(m_rho, 2) - 1. * s) * s +
6143 pow(
pion_mass, 4) * (-3. * pow(m_rho, 2) + s) +
6145 (-6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) + 3. * s) +
6147 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
6148 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
6149 4. * pow(m_rho, 2) * s + 2. * pow(s, 2))) +
6152 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
6153 pow(
pion_mass, 4) * (3. * pow(m_rho, 2) + s) +
6155 (-6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) + 3. * s) +
6157 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
6158 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
6159 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)))) *
6160 log(abs(-1. * pow(
a1_mass, 2) + tmax)) -
6164 0.5 * pow(
pion_mass, 4) * pow(m_rho, 2) +
6166 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
6168 (0.5 * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s)) +
6171 (1. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
6172 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
6173 0.5 * pow(m_rho, 2) * s) +
6175 (-2. * pow(
pion_mass, 4) - 0.5 * pow(m_rho, 2) * s +
6176 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
6177 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6182 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2)) +
6183 pow(
pion_mass, 4) * (0.5 * pow(m_rho, 2) - 1. * s) +
6184 pow(
pion_mass, 2) * (-0.5 * pow(m_rho, 2) + 0.5 * s) * s +
6187 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6188 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s)) +
6191 pow(
pion_mass, 2) * (1.5 * pow(m_rho, 2) - 2. * s) * s +
6192 (-0.5 * pow(m_rho, 2) + 0.5 * s) * pow(s, 2) +
6193 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1.5 * s) +
6195 (-2. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 2) * s +
6196 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
6197 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6199 1. * pow(m_rho, 2) + 1. * s) -
6200 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
6206 (8. * pow(
pion_mass, 2) + 4. * pow(m_rho, 2) - 4. * s) +
6208 (-12. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) +
6209 4. * pow(m_rho, 2) * s +
6210 pow(
pion_mass, 2) * (-12. * pow(m_rho, 2) + 8. * s)) +
6213 4. * pow(
pion_mass, 4) * pow(m_rho, 2) -
6214 2. * pow(m_rho, 6) + 4. * pow(m_rho, 4) * s -
6215 2. * pow(m_rho, 2) * pow(s, 2) +
6217 (-8. * pow(m_rho, 4) + 8. * pow(m_rho, 2) * s)) +
6219 (8. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
6220 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 4. * s) -
6221 2. * pow(m_rho, 4) * s - 2. * pow(m_rho, 2) * pow(s, 2) +
6224 (8. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s -
6229 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
6231 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
6232 pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
6234 (6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6235 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) +
6236 pow(m_rho, 2) * s) +
6238 (-4. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) -
6239 5. * pow(m_rho, 4) * s + 4. * pow(m_rho, 2) * pow(s, 2) -
6241 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
6243 (2. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s +
6247 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
6249 (-2. * pow(m_rho, 2) + 2. * s) +
6251 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
6253 (3. * pow(m_rho, 4) - 5. * pow(m_rho, 2) * s) +
6255 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
6256 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
6257 3. * pow(m_rho, 2) * s) +
6259 (-4. * pow(
pion_mass, 6) + pow(m_rho, 4) * s -
6261 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
6263 (-2. * pow(m_rho, 4) + 4. * pow(m_rho, 2) * s +
6264 2. * pow(s, 2))))) *
6265 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6268 0.25 * pow(m_rho, 4) - 1. * pow(
pion_mass, 2) * s -
6269 0.5 * pow(m_rho, 2) * s + 0.25 * pow(s, 2) +
6271 (-2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 1. * s)) -
6275 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
6277 (pow(
pion_mass, 2) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
6278 (-0.25 + 0.125 *
delta) * (pow(m_rho, 2) + s)) +
6280 (-1. *
C4 * pow(m_rho, 6) +
6282 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) +
6283 0.25 *
delta * pow(s, 2) +
6284 pow(m_rho, 2) * s * (-0.25 - 0.375 *
delta - 1. *
C4 * s) +
6285 pow(m_rho, 4) * (0.75 - 0.125 *
delta + 2. *
C4 * s))) +
6288 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
6290 (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
6292 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
6293 0.25 *
delta * pow(s, 2) +
6294 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
6296 (pow(
pion_mass, 4) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
6297 s * ((0.25 - 0.125 *
delta) * pow(m_rho, 2) +
6298 (-0.25 + 0.125 *
delta) * s) +
6300 (2. *
C4 * pow(m_rho, 4) + (0.5 + 0.25 *
delta) * s +
6301 pow(m_rho, 2) * (-1. - 2. *
C4 * s))))) *
6302 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6305 log(abs(-1. * pow(
pion_mass, 2) + tmax)) +
6308 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
6310 5.999999999999999 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 2) *
6312 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 2) * pow(s, 2)) *
6313 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6314 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
6316 (0. +
eta2 * pow(
pion_mass, 2) * (4. * pow(m_rho, 2) - 4. * s) +
6317 eta1 * (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s +
6318 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 4. * s))) *
6319 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6321 (2. * (-2. + 1. *
delta) *
6322 (0. + (-0.25 + 0.125 *
delta) * pow(m_rho, 2) * s +
6324 pow(m_rho, 2) * (1. + 2. *
C4 * s))) *
6325 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6328 (
eta1 * (pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1. * s) +
6330 (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6331 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
6333 (-1. * pow(m_rho, 4) + 2.5 * pow(m_rho, 2) * s -
6335 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
6338 (0.5 * pow(m_rho, 6) +
6339 pow(
pion_mass, 4) * (1. * pow(m_rho, 2) - 1. * s) -
6340 1.5 * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) -
6342 pow(
pion_mass, 2) * (1.5 * pow(m_rho, 4) -
6343 3. * pow(m_rho, 2) * s + 1.5 * pow(s, 2)) +
6345 (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
6347 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
6348 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6351 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
6352 2. * pow(m_rho, 2) * s + pow(s, 2) +
6354 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
6356 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax)) +
6360 (-1. * pow(m_rho, 2) + 1. * s) +
6362 (-0.5 * pow(m_rho, 4) +
6363 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6364 0.5 * pow(m_rho, 2) * s) +
6366 (0.5 * pow(m_rho, 4) - 0.5 * pow(m_rho, 2) * s +
6367 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
6368 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
6372 (
eta1 * (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6373 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
6374 eta2 * (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
6376 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
6377 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
6379 1. * pow(m_rho, 2) + 1. * s) -
6382 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
6384 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 4) * s +
6386 (
C4 * (32. - 16. *
delta) * pow(m_rho, 6) +
6389 (-8. + 24. *
delta - 10. * pow(
delta, 2) + 32. *
C4 * s -
6391 pow(m_rho, 4) * (-16. + 8. *
delta - 64. *
C4 * s +
6393 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
6394 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
6395 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
6403 (-4. * pow(
pion_mass, 4) + 4. * pow(m_rho, 4) +
6404 2. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s) +
6406 (12. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) -
6407 8. * pow(
pion_mass, 2) * s - 4. * pow(m_rho, 2) * s +
6411 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
6413 (6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) - 3. * s) -
6415 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
6416 1. * pow(m_rho, 4) * s + pow(m_rho, 2) * pow(s, 2) +
6419 1. * pow(m_rho, 2) + s) +
6421 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6422 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
6423 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s))) +
6426 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
6428 (6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) - 3. * s) -
6433 pow(m_rho, 2) + s) +
6435 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6436 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
6437 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s)))) *
6441 2. * pow(
a1_mass, 2) * (-1. * pow(m_rho, 2) + s + tmax) -
6442 4. * pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + s + tmax) +
6443 pow(-1. * pow(m_rho, 2) + s + tmax, 2))) -
6446 (0.5 - 1. *
C4 * pow(m_rho, 2)) +
6448 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
6450 (pow(
pion_mass, 2) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
6451 (0.25 - 0.125 *
delta) * (pow(m_rho, 2) + s)) +
6453 (1. *
C4 * pow(m_rho, 6) +
6455 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
6456 0.25 *
delta * pow(s, 2) +
6457 pow(m_rho, 4) * (-0.75 + 0.125 *
delta - 2. *
C4 * s) +
6459 (0.25 + 0.375 *
delta + 1. *
C4 * s))) +
6462 (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
6464 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
6466 (-0.5 * pow(m_rho, 4) + 1. *
C4 * pow(m_rho, 6) +
6468 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) +
6469 0.25 *
delta * pow(s, 2) -
6470 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
6472 (pow(
pion_mass, 4) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
6473 s * ((-0.25 + 0.125 *
delta) * pow(m_rho, 2) +
6474 (0.25 - 0.125 *
delta) * s) +
6476 (-2. *
C4 * pow(m_rho, 4) + (-0.5 - 0.25 *
delta) * s +
6477 pow(m_rho, 2) * (1. + 2. *
C4 * s))))) *
6480 1. * pow(m_rho, 2) + s + tmax,
6488 0.5 * pow(
pion_mass, 6) * pow(m_rho, 2) +
6491 (1. * pow(
pion_mass, 2) + 1.5 * pow(m_rho, 2) -
6494 0.5 * pow(m_rho, 2) - 1. * s)) +
6496 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
6498 (-1. * pow(
pion_mass, 2) - 0.5 * pow(m_rho, 2) + 1. * s)) +
6501 1. * pow(
pion_mass, 2) * pow(m_rho, 4) -
6505 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
6506 0.5 * pow(m_rho, 2) * s) +
6508 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6509 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 0.5 * s) +
6510 0.5 * pow(m_rho, 2) * s) +
6514 0.5 * pow(m_rho, 2) * s +
6515 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1.5 * s)))) *
6518 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
6519 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
6520 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
6521 2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
6523 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
6529 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
6532 4. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
6533 1.5 * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
6535 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
6536 0.5 * pow(m_rho, 4) * pow(s, 2) -
6538 1. * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
6542 0.5 * pow(
pion_mass, 2) * s + 0.5 * pow(s, 2)) +
6545 (1. * pow(m_rho, 2) - 0.5 * s) * s +
6546 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 2.5 * s)) +
6549 pow(
pion_mass, 4) * (2. * pow(m_rho, 2) - 1.5 * s) -
6550 0.5 * pow(m_rho, 4) * s + 0.5 * pow(s, 3) +
6552 (1. * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s -
6558 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 0.5 * s) +
6560 pow(
a1_mass, 4) * (-0.5 * pow(m_rho, 4) +
6561 (-0.5 * pow(
pion_mass, 2) + 0.5 * s) * s) +
6562 pow(
pion_mass, 4) * (-0.5 * pow(m_rho, 4) +
6563 2. * pow(m_rho, 2) * s - 1.5 * pow(s, 2)) +
6565 (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
6568 (1. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 4) +
6569 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) -
6570 1. * pow(m_rho, 2) * s + 0.5 * pow(s, 2) +
6572 1. * pow(m_rho, 2) + 1.5 * s)) +
6575 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 0.5 * s) +
6577 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
6579 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
6580 0.5 * pow(s, 2))))) *
6583 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
6584 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
6585 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
6586 2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
6588 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
6592 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
6593 2. * pow(m_rho, 2) * s + pow(s, 2) +
6595 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
6596 (0.0625 * pow(
eta1 - 1. *
eta2, 2) *
6600 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
6603 0.5 * pow(m_rho, 2) + 0.5 * s) +
6605 (10. * pow(
pion_mass, 6) - 2.5 * pow(m_rho, 6) +
6606 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) +
6607 6. * pow(m_rho, 4) * s - 4.5 * pow(m_rho, 2) * pow(s, 2) +
6611 (1. * pow(m_rho, 2) - 1. * s) * s +
6612 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
6614 (1. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
6615 pow(
pion_mass, 4) * (2.5 * pow(m_rho, 2) - 0.5 * s) -
6616 1. * pow(m_rho, 4) * s + 0.5 * pow(m_rho, 2) * pow(s, 2) +
6618 (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)) +
6621 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
6622 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
6623 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
6625 6. * pow(m_rho, 2) + 6. * s) +
6627 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
6630 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
6631 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
6632 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
6634 (-5. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) -
6635 3.5 * pow(m_rho, 6) * s + 4.5 * pow(m_rho, 4) * pow(s, 2) -
6636 2.5 * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
6637 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
6639 (-2. * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s +
6642 (3. * pow(m_rho, 6) - 8. * pow(m_rho, 4) * s +
6643 7. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)))) +
6647 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
6650 0.5 * pow(m_rho, 2) + 0.5 * s) +
6652 (-10. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) +
6653 5. * pow(m_rho, 2) * s - 1. * pow(s, 2) +
6654 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
6656 (10. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
6657 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) -
6658 3. * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) +
6661 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s)) +
6664 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
6665 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
6666 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
6668 6. * pow(m_rho, 2) + 6. * s) +
6670 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
6673 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
6674 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
6675 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
6678 pow(
pion_mass, 6) * (2.5 * pow(m_rho, 2) - 0.5 * s) +
6680 (4. * pow(m_rho, 4) - 6. * pow(m_rho, 2) * s) +
6682 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
6685 (1.5 * pow(m_rho, 6) - 6. * pow(m_rho, 4) * s +
6686 4.5 * pow(m_rho, 2) * pow(s, 2))) +
6689 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
6691 (-8. * pow(m_rho, 4) + 13. * pow(m_rho, 2) * s +
6694 (-1. * pow(m_rho, 6) + 6. * pow(m_rho, 4) * s -
6695 3. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)) +
6696 s * (0.5 * pow(m_rho, 6) - 0.5 * pow(m_rho, 4) * s -
6697 0.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3)))) +
6702 1. * pow(m_rho, 2) - 1. * s) +
6704 (-10. * pow(
pion_mass, 2) - 5. * pow(m_rho, 2) + 5. * s) +
6706 (20. * pow(
pion_mass, 4) + 6. * pow(m_rho, 4) +
6707 pow(
pion_mass, 2) * (20. * pow(m_rho, 2) - 16. * s) -
6708 8. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
6710 (-20. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
6711 6. * pow(m_rho, 4) * s - 2. * pow(s, 3) +
6712 pow(
pion_mass, 4) * (-30. * pow(m_rho, 2) + 18. * s) +
6714 (-18. * pow(m_rho, 4) + 18. * pow(m_rho, 2) * s)) +
6716 (-2. * pow(
pion_mass, 8) - 1. * pow(m_rho, 8) +
6717 3. * pow(m_rho, 6) * s - 3. * pow(m_rho, 4) * pow(s, 2) +
6718 1. * pow(m_rho, 2) * pow(s, 3) +
6719 pow(
pion_mass, 6) * (-5. * pow(m_rho, 2) + 1. * s) +
6721 (-10. * pow(m_rho, 4) + 10. * pow(m_rho, 2) * s) +
6723 (-6. * pow(m_rho, 6) + 12. * pow(m_rho, 4) * s -
6724 6. * pow(m_rho, 2) * pow(s, 2))) +
6726 (10. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) +
6727 pow(
pion_mass, 6) * (20. * pow(m_rho, 2) - 8. * s) -
6728 2. * pow(m_rho, 6) * s + 2. * pow(m_rho, 2) * pow(s, 3) -
6731 (22. * pow(m_rho, 4) - 20. * pow(m_rho, 2) * s -
6734 (8. * pow(m_rho, 6) - 12. * pow(m_rho, 4) * s +
6738 1. * pow(m_rho, 6) +
6739 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 12. * s) +
6741 12. * pow(m_rho, 2) - 12. * s) -
6742 3. * pow(m_rho, 4) * s + 3. * pow(m_rho, 2) * pow(s, 2) -
6745 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
6748 (-24. * pow(
pion_mass, 4) - 6. * pow(m_rho, 4) +
6749 12. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
6751 (-24. * pow(m_rho, 2) + 24. * s))))) *
6754 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
6755 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
6756 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
6757 2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
6759 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
6763 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
6764 2. * pow(m_rho, 2) * s + pow(s, 2) +
6766 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)))) /