#include #include #include #include "IRoneloopintegrals.h" #include "amplitudes.iso.h" #include "amplitudes.em.h" #include "amplitudes.IR.h" typedef std::complex dcomplex; const double pi=4.*atan(1.); const double pi16 = 1./(16.*pi*pi); inline double prop(double msq, double qsq){ return 1./(qsq-msq); } namespace inputpars{ extern double mpp,mpo,mkp,mko,me,fpip,fkp,mu,ZE,mud,spe,cpe,e2, L1r,L2r,L3r,L4r,L5r,L6r,L7r,L8r,L9r,L10r, K1r,K2r,K3r,K4r,K5r,K6r,K7r,K8r,K9r,K10r,K11r,K12r,K13r,K14r, G8,G8p,G27,GE, D1r,D2r,D3r,D4r,D5r,D6r,D7r,D13r,D15r,D26r,D27r,D28r,D29r,D30r,D31r, N1r,N2r,N3r,N4r,N5r,N6r,N7r,N8r,N9r,N10r,N11r,N12r,N13r,N14r,N15r, Z1r,Z2r,Z3r,Z4r,Z5r,Z6r,Z7r,Z8r,Z9r,Z10r,Z11r,Z12r,Z13r,Z14r,Z15r ; } namespace photoncutoffs{ double eg3pi,egfpi,egcoulomb; } dcomplex Ap4IKLooo(double s1, double s2,double s3){ using namespace inputpars; double mpp2 = mpp*mpp, mpo2 = mpo*mpo, mkp2 = mkp*mkp, mko2 = mko*mko, me2 = me*me, mu2 = mu*mu; double mp2 = mpo2,mk2=mko2; dcomplex K3Pp4Ix6333=dcomplex(0.); const double mg2 = 0.; K3Pp4Ix6333 = + B1b(mg2,mp2,mp2,mu2)*G8 * ( - 3*mk2 ); K3Pp4Ix6333 += + B1b(mg2,mp2,mp2,mu2)*G27 * ( 3*mk2 ); K3Pp4Ix6333 += + B1b(mg2,mk2,mk2,mu2)*G8 * ( - mk2 ); K3Pp4Ix6333 += + B1b(mg2,mk2,mk2,mu2)*G27 * ( mk2 ); K3Pp4Ix6333 += + dqsqBb(mg2,mp2,mp2,mu2)*G8 * ( 6*mp2*mk2 ); K3Pp4Ix6333 += + dqsqBb(mg2,mp2,mp2,mu2)*G27 * ( - 6*mp2*mk2 ); K3Pp4Ix6333 += + dqsqBb(mg2,mk2,mk2,mu2)*G8 * ( 2*pow(mk2,2) ); K3Pp4Ix6333 += + dqsqBb(mg2,mk2,mk2,mu2)*G27 * ( - 2*pow(mk2,2) ); K3Pp4Ix6333 += + dqsqB1b(mg2,mp2,mp2,mu2)*G8 * ( - 6*mp2*mk2 ); K3Pp4Ix6333 += + dqsqB1b(mg2,mp2,mp2,mu2)*G27 * ( 6*mp2*mk2 ); K3Pp4Ix6333 += + dqsqB1b(mg2,mk2,mk2,mu2)*G8 * ( - 2*pow(mk2,2) ); K3Pp4Ix6333 += + dqsqB1b(mg2,mk2,mk2,mu2)*G27 * ( 2*pow(mk2,2) ) ; return K3Pp4Ix6333*e2/(pow(fpip,3)*fkp); } dcomplex Ap4IKLpmo(double s1, double s2,double s3){ using namespace inputpars; double mko2 = mko*mko, mkp2 = mkp*mkp,mpp2 = mpp*mpp,mpo2=mpo*mpo, me2 = me*me,mu2 = mu*mu; double mp2 = (2.*mpp2+mpo2)/3.,mk2=mko2; dcomplex K3Pp4Ix10113 = dcomplex(0.); const double mg2 = 0.; double propkp = prop(mk2,mp2); K3Pp4Ix10113 = + Ab(mp2,mu2)*G8 * ( 1./6. - 3./2.*mp2*s3*pow(propkp,2) - mp2* propkp + 2*pow(mp2,2)*pow(propkp,2) ); K3Pp4Ix10113 += + Ab(mp2,mu2)*G27 * ( 1./3. + 8./3.*mp2*propkp - s3*propkp ); K3Pp4Ix10113 += + Ab(mk2,mu2)*G8 * ( - 1 + 3./2.*mp2*s3*pow( propkp,2) + 3*mp2*propkp - 2*pow(mp2,2)*pow(propkp,2) - 3./2.* s3*propkp ); K3Pp4Ix10113 += + Ab(mk2,mu2)*G27 * ( - 2./3. + 5./3.*mp2* propkp - s3*propkp ); K3Pp4Ix10113 += + Bb(mp2,mp2,s3,mu2)*G8 * ( - 5./2.*mp2*s3* propkp - 3./2.*mp2 + 2*pow(mp2,2)*propkp + 7./3.*s3 ); K3Pp4Ix10113 += + Bb(mp2,mp2,s3,mu2)*G27 * ( 5./2.*mp2*s3*propkp - mp2 - 2*pow(mp2,2)*propkp - 5./6.*mk2 + 1./6.*s3 ); K3Pp4Ix10113 += + Bb(mp2,mk2,s3,mu2)*G8 * ( 5./2.*mp2*s3*propkp + 2*mp2 - 2*pow(mp2,2)*propkp + 3./2.*mk2 - 5./2.*s3 ); K3Pp4Ix10113 += + Bb(mp2,mk2,s3,mu2)*G27 * ( 5./3.*mp2*s3*propkp + 4./3.*mp2 - 4./3.*pow(mp2,2)*propkp + mk2 - 5./3.*s3 ); K3Pp4Ix10113 += + Bb(mg2,mp2,mp2,mu2)*G8 * ( - 3./2.*mp2 - 3./2. *mk2 + 29./6.*s3 ); K3Pp4Ix10113 += + Bb(mg2,mp2,mp2,mu2)*G27 * ( 8./3.*mp2*s3* propkp - 9*mp2 + 4*pow(mp2,2)*s3*pow(propkp,2) + 20./3.*pow( mp2,2)*propkp - 16./3.*pow(mp2,3)*pow(propkp,2) - 13./3.*mk2 + 11*s3 ); K3Pp4Ix10113 += + Bb(mg2,mk2,mp2,mu2)*G8 * ( 3./2.*mp2*s3*propkp + 2*mp2 - 2*pow(mp2,2)*propkp + 3./2.*mk2 - 3./2.*s3 ); K3Pp4Ix10113 += + Bb(mg2,mk2,mp2,mu2)*G27 * ( mp2*s3*propkp + 4./ 3.*mp2 - 4./3.*pow(mp2,2)*propkp + mk2 - s3 ); K3Pp4Ix10113 += + Bb(mg2,mk2,mk2,mu2)*G27 * ( 4*mp2*s3*propkp + 4./3.*mp2 - 4*pow(mp2,2)*s3*pow(propkp,2) - 20./3.*pow(mp2,2)* propkp + 16./3.*pow(mp2,3)*pow(propkp,2) ); K3Pp4Ix10113 += + B1b(mp2,mk2,s3,mu2)*G8 * ( - mk2 ); K3Pp4Ix10113 += + B1b(mp2,mk2,s3,mu2)*G27 * ( - 2./3.*mk2 ); K3Pp4Ix10113 += + B1b(mg2,mp2,mp2,mu2)*G8 * ( 7./2.*mp2 + 3*pow( mp2,2)*s3*pow(propkp,2) + 3*pow(mp2,2)*propkp - 4*pow(mp2,3)* pow(propkp,2) - 1./2.*mk2 - 37./6.*s3 ); K3Pp4Ix10113 += + B1b(mg2,mp2,mp2,mu2)*G27 * ( - 71./6.*mp2*s3* propkp + 31./3.*mp2 - 2*pow(mp2,2)*s3*pow(propkp,2) + 16./3.* pow(mp2,2)*propkp + 8./3.*pow(mp2,3)*pow(propkp,2) + 11./2.* mk2 - 73./6.*s3 ); K3Pp4Ix10113 += + B1b(mg2,mk2,mp2,mu2)*G8 * ( 3*mp2*s3*propkp + 3*mp2 - 3*pow(mp2,2)*s3*pow(propkp,2) - 7*pow(mp2,2)*propkp + 4*pow(mp2,3)*pow(propkp,2) ); K3Pp4Ix10113 += + B1b(mg2,mk2,mp2,mu2)*G27 * ( 2*mp2*s3*propkp + 2*mp2 - 2*pow(mp2,2)*s3*pow(propkp,2) - 14./3.*pow(mp2,2)* propkp + 8./3.*pow(mp2,3)*pow(propkp,2) ); K3Pp4Ix10113 += + B1b(mg2,mk2,mk2,mu2)*G8 * ( mp2 - s3 ); K3Pp4Ix10113 += + B1b(mg2,mk2,mk2,mu2)*G27 * ( - 13./2.*mp2*s3* propkp - 2./3.*mp2 + 4*pow(mp2,2)*s3*pow(propkp,2) + 10*pow( mp2,2)*propkp - 16./3.*pow(mp2,3)*pow(propkp,2) + 5./6.*mk2 - 3./2.*s3 ); K3Pp4Ix10113 += + Cb(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( 6*mp2*s3 - 4*pow(mp2,2) - 2*pow(s3,2) ); K3Pp4Ix10113 += + Cb(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( - 10./ 3.*mp2*mk2 + 22./3.*mp2*s3 - 5*mp2*pow(s3,2)*propkp + 50./3.* pow(mp2,2)*s3*propkp - 8./3.*pow(mp2,2) - 40./3.*pow(mp2,3)* propkp + 5./3.*mk2*s3 - 3*pow(s3,2) ); K3Pp4Ix10113 += + Cb(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( - 9./2. *mp2*mk2 + 7./2.*mk2*s3 - 3./2.*pow(mk2,2) ); K3Pp4Ix10113 += + Cb(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( - 3* mp2*mk2 + 7./3.*mk2*s3 - pow(mk2,2) ); K3Pp4Ix10113 += + C11b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( - 29./ 2.*mp2*s3 - 7./2.*mp2*pow(s3,2)*propkp + 13*pow(mp2,2)*s3* propkp + 53./3.*pow(mp2,2) - 12*pow(mp2,3)*propkp + 7./3.*pow( s3,2) ); K3Pp4Ix10113 += + C11b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( 10./3. *mp2*mk2 + 1./3.*mp2*s3 + 11./6.*mp2*pow(s3,2)*propkp - 44./3. *pow(mp2,2)*s3*propkp - 6*pow(mp2,2) + 56./3.*pow(mp2,3)* propkp - 5./6.*mk2*s3 + 1./6.*pow(s3,2) ); K3Pp4Ix10113 += + C11b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( 3./2.* mp2*mk2 + 13*mp2*s3 + 7./2.*mp2*pow(s3,2)*propkp - 13*pow( mp2,2)*s3*propkp - 12*pow(mp2,2) + 12*pow(mp2,3)*propkp + mk2* s3 - 1./2.*pow(mk2,2) - 7./2.*pow(s3,2) ); K3Pp4Ix10113 += + C11b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( mp2* mk2 + 26./3.*mp2*s3 + 7./3.*mp2*pow(s3,2)*propkp - 26./3.*pow( mp2,2)*s3*propkp - 8*pow(mp2,2) + 8*pow(mp2,3)*propkp + 2./3.* mk2*s3 - 1./3.*pow(mk2,2) - 7./3.*pow(s3,2) ); K3Pp4Ix10113 += + C12b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( 5./2.* mp2*s3 + 7./2.*mp2*pow(s3,2)*propkp - 6*pow(mp2,2)*s3*propkp - 4*pow(mp2,2) ); K3Pp4Ix10113 += + C12b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( 5./3. *mp2*s3 + 7./3.*mp2*pow(s3,2)*propkp - 4*pow(mp2,2)*s3*propkp - 8./3.*pow(mp2,2) ); K3Pp4Ix10113 += + C12b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( - mp2 *mk2 - 6*mp2*s3 - 7./2.*mp2*pow(s3,2)*propkp + 6*pow(mp2,2)*s3 *propkp - 5./2.*mk2*s3 + pow(mk2,2) + 7./2.*pow(s3,2) ); K3Pp4Ix10113 += + C12b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( - 2./ 3.*mp2*mk2 - 4*mp2*s3 - 7./3.*mp2*pow(s3,2)*propkp + 4*pow( mp2,2)*s3*propkp - 5./3.*mk2*s3 + 2./3.*pow(mk2,2) + 7./3.* pow(s3,2) ); K3Pp4Ix10113 += + C21b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( 16./3. *mp2*s3 + 1./2.*mp2*pow(s3,2)*propkp - 5*pow(mp2,2)*s3*propkp - 35./3.*pow(mp2,2) + 8*pow(mp2,3)*propkp - 2./3.*pow(s3,2) ) ; K3Pp4Ix10113 += + C21b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( - 16. /3.*mp2*s3 - 1./2.*mp2*pow(s3,2)*propkp + 10./3.*pow(mp2,2)*s3 *propkp + 10*pow(mp2,2) - 8*pow(mp2,3)*propkp + 2./3.*pow( s3,2) ); K3Pp4Ix10113 += + C21b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( - mp2 *mk2 - 5*mp2*s3 - 1./2.*mp2*pow(s3,2)*propkp + 5*pow(mp2,2)*s3 *propkp + 8*pow(mp2,2) - 8*pow(mp2,3)*propkp + 1./2.*mk2*s3 + 1./2.*pow(s3,2) ); K3Pp4Ix10113 += + C21b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( - 2./ 3.*mp2*mk2 - 10./3.*mp2*s3 - 1./3.*mp2*pow(s3,2)*propkp + 10./ 3.*pow(mp2,2)*s3*propkp + 16./3.*pow(mp2,2) - 16./3.*pow( mp2,3)*propkp + 1./3.*mk2*s3 + 1./3.*pow(s3,2) ); K3Pp4Ix10113 += + C22b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( - mp2 *s3 - mp2*pow(s3,2)*propkp ); K3Pp4Ix10113 += + C22b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( - 2./ 3.*mp2*s3 - 2./3.*mp2*pow(s3,2)*propkp ); K3Pp4Ix10113 += + C22b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( mp2* pow(s3,2)*propkp + mk2*s3 - pow(s3,2) ); K3Pp4Ix10113 += + C22b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( 2./3. *mp2*pow(s3,2)*propkp + 2./3.*mk2*s3 - 2./3.*pow(s3,2) ); K3Pp4Ix10113 += + C23b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( 1./2.* mp2*s3 + 1./2.*mp2*pow(s3,2)*propkp + 2*pow(mp2,2)*s3*propkp + 2*pow(mp2,2) ); K3Pp4Ix10113 += + C23b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( 1./3. *mp2*s3 + 1./3.*mp2*pow(s3,2)*propkp + 4./3.*pow(mp2,2)*s3* propkp + 4./3.*pow(mp2,2) ); K3Pp4Ix10113 += + C23b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( 2*mp2* mk2 + 2*mp2*s3 - 1./2.*mp2*pow(s3,2)*propkp - 2*pow(mp2,2)*s3* propkp - 3./2.*mk2*s3 + 1./2.*pow(s3,2) ); K3Pp4Ix10113 += + C23b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( 4./3. *mp2*mk2 + 4./3.*mp2*s3 - 1./3.*mp2*pow(s3,2)*propkp - 4./3.* pow(mp2,2)*s3*propkp - mk2*s3 + 1./3.*pow(s3,2) ); K3Pp4Ix10113 += + C24b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( - 3* mp2*s3*propkp - 35./3.*mp2 + 8*pow(mp2,2)*propkp + 8./3.*s3 ); K3Pp4Ix10113 += + C24b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( 4./3. *mp2*s3*propkp + 10*mp2 - 8*pow(mp2,2)*propkp - 8./3.*s3 ); K3Pp4Ix10113 += + C24b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( 3*mp2* s3*propkp + 8*mp2 - 8*pow(mp2,2)*propkp + mk2 - 3*s3 ); K3Pp4Ix10113 += + C24b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( 2*mp2 *s3*propkp + 16./3.*mp2 - 16./3.*pow(mp2,2)*propkp + 2./3.*mk2 - 2*s3 ); K3Pp4Ix10113 += + dqsqBb(mg2,mp2,mp2,mu2)*G8 * ( 10*mp2*s3 - 10* pow(mp2,2) ); K3Pp4Ix10113 += + dqsqBb(mg2,mp2,mp2,mu2)*G27 * ( - 25./3.*mp2* mk2 + 15*mp2*s3 + 25*pow(mp2,2)*s3*propkp - 20./3.*pow(mp2,2) - 100./3.*pow(mp2,3)*propkp ); K3Pp4Ix10113 += + dqsqBb(mg2,mk2,mk2,mu2)*G8 * ( - 2*mp2*mk2 + 2*mk2*s3 ); K3Pp4Ix10113 += + dqsqBb(mg2,mk2,mk2,mu2)*G27 * ( - 4./3.*mp2* mk2 - 5*mp2*s3 + 5*pow(mp2,2)*s3*propkp + 20./3.*pow(mp2,2) - 20./3.*pow(mp2,3)*propkp + 3*mk2*s3 - 5./3.*pow(mk2,2) ); K3Pp4Ix10113 += + dqsqB1b(mg2,mp2,mp2,mu2)*G8 * ( - 10*mp2*s3 + 10*pow(mp2,2) ); K3Pp4Ix10113 += + dqsqB1b(mg2,mp2,mp2,mu2)*G27 * ( 25./3.*mp2* mk2 - 15*mp2*s3 - 25*pow(mp2,2)*s3*propkp + 20./3.*pow(mp2,2) + 100./3.*pow(mp2,3)*propkp ); K3Pp4Ix10113 += + dqsqB1b(mg2,mk2,mk2,mu2)*G8 * ( 2*mp2*mk2 - 2* mk2*s3 ); K3Pp4Ix10113 += + dqsqB1b(mg2,mk2,mk2,mu2)*G27 * ( 4./3.*mp2*mk2 + 5*mp2*s3 - 5*pow(mp2,2)*s3*propkp - 20./3.*pow(mp2,2) + 20./ 3.*pow(mp2,3)*propkp - 3*mk2*s3 + 5./3.*pow(mk2,2) ); return K3Pp4Ix10113*e2/(pow(fpip,3)*fkp); } dcomplex Ap4IKSpmo(double s1, double s2,double s3){ using namespace inputpars; double mko2 = mko*mko, mkp2 = mkp*mkp,mpp2 = mpp*mpp,mpo2=mpo*mpo, me2 = me*me,mu2 = mu*mu; double mp2 = (2.*mpp2+mpo2)/3.,mk2=mko2; dcomplex K3Pp4Ix10123 = dcomplex(0.); const double mg2 = 0.; double propkp = prop(mk2,mp2); K3Pp4Ix10123 = + Ab(mp2,mu2)*G8 * ( - 1./2.*mp2*s1*pow(propkp,2) + 1./2.*mp2* s2*pow(propkp,2) ); K3Pp4Ix10123 += + Ab(mp2,mu2)*G27 * ( - 1./3.*s1*propkp + 1./3. *s2*propkp ); K3Pp4Ix10123 += + Ab(mk2,mu2)*G8 * ( 1./2.*mp2*s1*pow(propkp,2) - 1./2.*mp2*s2*pow(propkp,2) - 1./2.*s1*propkp + 1./2.*s2* propkp ); K3Pp4Ix10123 += + Ab(mk2,mu2)*G27 * ( - 1./3.*s1*propkp + 1./3. *s2*propkp ); K3Pp4Ix10123 += + Bb(mg2,mp2,mp2,mu2)*G8 * ( 1./2.*s1 - 1./2.*s2 ); K3Pp4Ix10123 += + Bb(mg2,mp2,mp2,mu2)*G27 * ( - 4./3.*mp2*s1* propkp + 4./3.*mp2*s2*propkp + 4./3.*pow(mp2,2)*s1*pow( propkp,2) - 4./3.*pow(mp2,2)*s2*pow(propkp,2) + 1./3.*s1 - 1./ 3.*s2 ); K3Pp4Ix10123 += + Bb(mg2,mk2,mp2,mu2)*G8 * ( 1./2.*mp2*s1*propkp - 1./2.*mp2*s2*propkp - 1./2.*s1 + 1./2.*s2 ); K3Pp4Ix10123 += + Bb(mg2,mk2,mp2,mu2)*G27 * ( 1./3.*mp2*s1* propkp - 1./3.*mp2*s2*propkp - 1./3.*s1 + 1./3.*s2 ); K3Pp4Ix10123 += + Bb(mg2,mk2,mk2,mu2)*G27 * ( 4./3.*mp2*s1* propkp - 4./3.*mp2*s2*propkp - 4./3.*pow(mp2,2)*s1*pow( propkp,2) + 4./3.*pow(mp2,2)*s2*pow(propkp,2) ); K3Pp4Ix10123 += + B1b(mg2,mp2,mp2,mu2)*G8 * ( pow(mp2,2)*s1*pow( propkp,2) - pow(mp2,2)*s2*pow(propkp,2) - 1./2.*s1 + 1./2.*s2 ); K3Pp4Ix10123 += + B1b(mg2,mp2,mp2,mu2)*G27 * ( - 17./6.*mp2*s1* propkp + 17./6.*mp2*s2*propkp - 2./3.*pow(mp2,2)*s1*pow( propkp,2) + 2./3.*pow(mp2,2)*s2*pow(propkp,2) + 73./6.*s1 - 73. /6.*s2 ); K3Pp4Ix10123 += + B1b(mg2,mk2,mp2,mu2)*G8 * ( mp2*s1*propkp - mp2*s2*propkp - pow(mp2,2)*s1*pow(propkp,2) + pow(mp2,2)*s2* pow(propkp,2) ); K3Pp4Ix10123 += + B1b(mg2,mk2,mp2,mu2)*G27 * ( 2./3.*mp2*s1* propkp - 2./3.*mp2*s2*propkp - 2./3.*pow(mp2,2)*s1*pow( propkp,2) + 2./3.*pow(mp2,2)*s2*pow(propkp,2) ); K3Pp4Ix10123 += + B1b(mg2,mk2,mk2,mu2)*G27 * ( - 13./6.*mp2*s1* propkp + 13./6.*mp2*s2*propkp + 4./3.*pow(mp2,2)*s1*pow( propkp,2) - 4./3.*pow(mp2,2)*s2*pow(propkp,2) + 5./2.*s1 - 5./ 2.*s2 ); K3Pp4Ix10123 += + Cb(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( - mp2* s1*s3*propkp - mp2*s1 - mp2*pow(s1,2)*propkp + mp2*s2*s3* propkp + mp2*s2 + mp2*pow(s2,2)*propkp + 4*pow(mp2,2)*s1* propkp - 4*pow(mp2,2)*s2*propkp ); K3Pp4Ix10123 += + Cb(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( - 2./3. *mp2*s1*s3*propkp + 6*mp2*s1 + mp2*pow(s1,2)*propkp + 2./3.* mp2*s2*s3*propkp - 6*mp2*s2 - mp2*pow(s2,2)*propkp - 2./3.* pow(mp2,2)*s1*propkp + 2./3.*pow(mp2,2)*s2*propkp + 5*mk2*s1 - 5*mk2*s2 - 5*pow(s1,2) + 5*pow(s2,2) ); K3Pp4Ix10123 += + Cb(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( mp2*s1* s3*propkp + 4*mp2*s1 + mp2*pow(s1,2)*propkp - mp2*s2*s3*propkp - 4*mp2*s2 - mp2*pow(s2,2)*propkp - 4*pow(mp2,2)*s1*propkp + 4*pow(mp2,2)*s2*propkp + 3./2.*mk2*s1 - 3./2.*mk2*s2 - s1*s3 - pow(s1,2) + s2*s3 + pow(s2,2) ); K3Pp4Ix10123 += + Cb(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( 2./3.* mp2*s1*s3*propkp + 8./3.*mp2*s1 + 2./3.*mp2*pow(s1,2)*propkp - 2./3.*mp2*s2*s3*propkp - 8./3.*mp2*s2 - 2./3.*mp2*pow(s2,2) *propkp - 8./3.*pow(mp2,2)*s1*propkp + 8./3.*pow(mp2,2)*s2* propkp + mk2*s1 - mk2*s2 - 2./3.*s1*s3 - 2./3.*pow(s1,2) + 2./ 3.*s2*s3 + 2./3.*pow(s2,2) ); K3Pp4Ix10123 += + C11b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( mp2*s1 *s3*propkp + 3./2.*mp2*s1 + 3./2.*mp2*pow(s1,2)*propkp - mp2* s2*s3*propkp - 3./2.*mp2*s2 - 3./2.*mp2*pow(s2,2)*propkp + mp2 *s3 - 5*pow(mp2,2)*s1*propkp + 5*pow(mp2,2)*s2*propkp - 2*pow( mp2,2) ); K3Pp4Ix10123 += + C11b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( - 10 *mp2*mk2 + 2./3.*mp2*s1*s3*propkp - 4*mp2*s1 - 3./2.*mp2*pow( s1,2)*propkp - 2./3.*mp2*s2*s3*propkp + 4*mp2*s2 + 3./2.*mp2* pow(s2,2)*propkp - 8./3.*mp2*s3 + 12*pow(mp2,2) - 15./2.*mk2* s1 + 15./2.*mk2*s2 + 5*mk2*s3 + 15./2.*pow(s1,2) - 15./2.*pow( s2,2) - 5./3.*pow(s3,2) ); K3Pp4Ix10123 += + C11b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( - mp2 *s1*s3*propkp - 5*mp2*s1 - 3./2.*mp2*pow(s1,2)*propkp + mp2*s2 *s3*propkp + 5*mp2*s2 + 3./2.*mp2*pow(s2,2)*propkp + 5*pow( mp2,2)*s1*propkp - 5*pow(mp2,2)*s2*propkp - 2*mk2*s1 + 2*mk2* s2 + s1*s3 + 3./2.*pow(s1,2) - s2*s3 - 3./2.*pow(s2,2) ); K3Pp4Ix10123 += + C11b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( - 2./ 3.*mp2*s1*s3*propkp - 10./3.*mp2*s1 - mp2*pow(s1,2)*propkp + 2. /3.*mp2*s2*s3*propkp + 10./3.*mp2*s2 + mp2*pow(s2,2)*propkp + 10./3.*pow(mp2,2)*s1*propkp - 10./3.*pow(mp2,2)*s2*propkp - 4./ 3.*mk2*s1 + 4./3.*mk2*s2 + 2./3.*s1*s3 + pow(s1,2) - 2./3.*s2* s3 - pow(s2,2) ); K3Pp4Ix10123 += + C12b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( 1./2.* mp2*s1*s3*propkp - 1./2.*mp2*s2*s3*propkp - 2*mp2*s3 + 4*pow( mp2,2) ); K3Pp4Ix10123 += + C12b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( 20* mp2*mk2 + 1./3.*mp2*s1*s3*propkp - 1./3.*mp2*s2*s3*propkp + 16. /3.*mp2*s3 - 24*pow(mp2,2) - 10*mk2*s3 + 10./3.*pow(s3,2) ); K3Pp4Ix10123 += + C12b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( - 1./ 2.*mp2*s1*s3*propkp + 1./2.*mp2*s2*s3*propkp + 1./2.*s1*s3 - 1. /2.*s2*s3 ); K3Pp4Ix10123 += + C12b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( - 1./ 3.*mp2*s1*s3*propkp + 1./3.*mp2*s2*s3*propkp + 1./3.*s1*s3 - 1. /3.*s2*s3 ); K3Pp4Ix10123 += + C21b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( - 1./ 2.*mp2*s1 - 1./2.*mp2*pow(s1,2)*propkp + 1./2.*mp2*s2 + 1./2.* mp2*pow(s2,2)*propkp - 1./2.*mp2*s3 + pow(mp2,2)*s1*propkp - pow(mp2,2)*s2*propkp + 2*pow(mp2,2) ); K3Pp4Ix10123 += + C21b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( 10* mp2*mk2 - 2*mp2*s1 + 1./2.*mp2*pow(s1,2)*propkp + 2*mp2*s2 - 1. /2.*mp2*pow(s2,2)*propkp - 1./3.*mp2*s3 + 2./3.*pow(mp2,2)*s1* propkp - 2./3.*pow(mp2,2)*s2*propkp - 12*pow(mp2,2) + 5./2.* mk2*s1 - 5./2.*mk2*s2 - 5./2.*mk2*s3 - 5./2.*pow(s1,2) + 5./2. *pow(s2,2) + 5./6.*pow(s3,2) ); K3Pp4Ix10123 += + C21b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( mp2*s1 + 1./2.*mp2*pow(s1,2)*propkp - mp2*s2 - 1./2.*mp2*pow(s2,2)* propkp - pow(mp2,2)*s1*propkp + pow(mp2,2)*s2*propkp + 1./2.* mk2*s1 - 1./2.*mk2*s2 - 1./2.*pow(s1,2) + 1./2.*pow(s2,2) ); K3Pp4Ix10123 += + C21b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( 2./3. *mp2*s1 + 1./3.*mp2*pow(s1,2)*propkp - 2./3.*mp2*s2 - 1./3.* mp2*pow(s2,2)*propkp - 2./3.*pow(mp2,2)*s1*propkp + 2./3.*pow( mp2,2)*s2*propkp + 1./3.*mk2*s1 - 1./3.*mk2*s2 - 1./3.*pow( s1,2) + 1./3.*pow(s2,2) ); K3Pp4Ix10123 += + C23b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( - 1./ 2.*mp2*s1*s3*propkp + 1./2.*mp2*s2*s3*propkp + mp2*s3 - 4*pow( mp2,2) ); K3Pp4Ix10123 += + C23b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( - 20 *mp2*mk2 - 1./3.*mp2*s1*s3*propkp + 1./3.*mp2*s2*s3*propkp + 2. /3.*mp2*s3 + 24*pow(mp2,2) + 5*mk2*s3 - 5./3.*pow(s3,2) ); K3Pp4Ix10123 += + C23b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( 1./2.* mp2*s1*s3*propkp - 1./2.*mp2*s2*s3*propkp - 1./2.*s1*s3 + 1./2. *s2*s3 ); K3Pp4Ix10123 += + C23b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( 1./3. *mp2*s1*s3*propkp - 1./3.*mp2*s2*s3*propkp - 1./3.*s1*s3 + 1./ 3.*s2*s3 ); K3Pp4Ix10123 += + C24b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( - mp2 *s1*propkp + mp2*s2*propkp ); K3Pp4Ix10123 += + C24b(mg2,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( 8./3. *mp2*s1*propkp - 8./3.*mp2*s2*propkp - 10*s1 + 10*s2 ); K3Pp4Ix10123 += + C24b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( mp2*s1 *propkp - mp2*s2*propkp - s1 + s2 ); K3Pp4Ix10123 += + C24b(mg2,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( 2./3. *mp2*s1*propkp - 2./3.*mp2*s2*propkp - 2./3.*s1 + 2./3.*s2 ); K3Pp4Ix10123 += + dqsqBb(mg2,mp2,mp2,mu2)*G27 * ( - 25*mp2*s1 + 25*mp2*s2 + 25./3.*pow(mp2,2)*s1*propkp - 25./3.*pow(mp2,2) *s2*propkp ); K3Pp4Ix10123 += + dqsqBb(mg2,mk2,mk2,mu2)*G27 * ( - 5./3.*mp2* s1 + 5./3.*mp2*s2 + 5./3.*pow(mp2,2)*s1*propkp - 5./3.*pow( mp2,2)*s2*propkp - 5*mk2*s1 + 5*mk2*s2 ); K3Pp4Ix10123 += + dqsqB1b(mg2,mp2,mp2,mu2)*G27 * ( 25*mp2*s1 - 25*mp2*s2 - 25./3.*pow(mp2,2)*s1*propkp + 25./3.*pow(mp2,2)*s2 *propkp ); K3Pp4Ix10123 += + dqsqB1b(mg2,mk2,mk2,mu2)*G27 * ( 5./3.*mp2*s1 - 5./3.*mp2*s2 - 5./3.*pow(mp2,2)*s1*propkp + 5./3.*pow( mp2,2)*s2*propkp + 5*mk2*s1 - 5*mk2*s2 ); return K3Pp4Ix10123*e2/(pow(fpip,3)*fkp); } dcomplex Ap4EKSpmo(double s1, double s2,double s3){ using namespace inputpars; double mko2 = mko*mko, mkp2 = mkp*mkp,mpp2 = mpp*mpp,mpo2=mpo*mpo, me2 = me*me,mu2 = mu*mu; double mp2 = (2.*mpp2+mpo2)/3.,mk2=mko2; dcomplex K3Pp4Ex10123 = dcomplex(0.); double propkp = prop(mk2,mp2); double props3 = 1./s3; K3Pp4Ex10123 = + G8*props3 * ( 4./3.*mk2*pi16*s1 - 4./3.*mk2*pi16*s2 - 4./3.* N14r*s1*s3 + 4./3.*N14r*s2*s3 - 2./3.*N15r*s1*s3 + 2./3.*N15r* s2*s3 - 2./9.*pi16*s1*s3 + 2./9.*pi16*s2*s3 ); K3Pp4Ex10123 += + G27*props3 * ( - 20./9.*mp2*pi16*s1 + 20./9.* mp2*pi16*s2 - 2./9.*mk2*pi16*s1 + 2./9.*mk2*pi16*s2 + 2./3.* D15r*s1*s3 - 2./3.*D15r*s2*s3 + 16./27.*pi16*s1*s3 - 16./27.* pi16*s2*s3 ); K3Pp4Ex10123 += + Ab(mp2,mu2)*G27*props3 * ( - 10./9.*mp2*s1* propkp + 10./9.*mp2*s2*propkp + 10./3.*s1 - 10./3.*s2 ); K3Pp4Ex10123 += + Ab(mk2,mu2)*G8*props3 * ( - 4./3.*s1 + 4./3.* s2 ); K3Pp4Ex10123 += + Ab(mk2,mu2)*G27*props3 * ( 10./9.*mp2*s1* propkp - 10./9.*mp2*s2*propkp + 2./9.*s1 - 2./9.*s2 ); K3Pp4Ex10123 += + Bb(mp2,mp2,s3,mu2)*G27*props3 * ( - 5./18.* mp2*s1*s3*propkp - 10./3.*mp2*s1 + 5./18.*mp2*s2*s3*propkp + 10./3.*mp2*s2 + 10./9.*pow(mp2,2)*s1*propkp - 10./9.*pow( mp2,2)*s2*propkp + 5./6.*s1*s3 - 5./6.*s2*s3 ); K3Pp4Ex10123 += + Bb(mk2,mk2,s3,mu2)*G8*props3 * ( 4./3.*mk2*s1 - 4./3.*mk2*s2 - 1./3.*s1*s3 + 1./3.*s2*s3 ); K3Pp4Ex10123 += + Bb(mk2,mk2,s3,mu2)*G27*props3 * ( 5./18.*mp2* s1*s3*propkp + 10./9.*mp2*s1 - 5./18.*mp2*s2*s3*propkp - 10./9. *mp2*s2 - 10./9.*pow(mp2,2)*s1*propkp + 10./9.*pow(mp2,2)*s2* propkp - 2./9.*mk2*s1 + 2./9.*mk2*s2 + 1./18.*s1*s3 - 1./18.* s2*s3 ); return K3Pp4Ex10123*e2/(pow(fpip,3)*fkp); } dcomplex Ap4IKPoop(double s1, double s2,double s3){ using namespace inputpars; double mko2 = mko*mko, mkp2 = mkp*mkp,mpp2 = mpp*mpp,mpo2=mpo*mpo, me2 = me*me,mu2 = mu*mu; double mp2 = (mpp2+2.*mpo2)/3.,mk2=mkp2; double propkp = prop(mk2,mp2); dcomplex K3Pp4Ix10331=dcomplex(0.); const double mg2 = 0.; K3Pp4Ix10331 = + Ab(mp2,mu2)*G8 * ( 1./3. - 3./2.*mp2*s3*pow(propkp,2) + 2*pow( mp2,2)*pow(propkp,2) ); K3Pp4Ix10331 += + Ab(mp2,mu2)*G27 * ( - 16./9. + 5./2.*mp2* propkp - s3*propkp ); K3Pp4Ix10331 += + Ab(mk2,mu2)*G8 * ( 1./6. + 3./2.*mp2*s3*pow( propkp,2) + 2*mp2*propkp - 2*pow(mp2,2)*pow(propkp,2) - 3./2.* s3*propkp ); K3Pp4Ix10331 += + Ab(mk2,mu2)*G27 * ( 1./9. + 11./6.*mp2*propkp - s3*propkp ); K3Pp4Ix10331 += + Bb(mp2,mp2,s3,mu2)*G8 * ( 10./3.*mp2*s3*propkp - 2*pow(mp2,2)*propkp ); K3Pp4Ix10331 += + Bb(mp2,mp2,s3,mu2)*G27 * ( 20./9.*mp2*s3* propkp - 4./3.*pow(mp2,2)*propkp ); K3Pp4Ix10331 += + Bb(mp2,mk2,s3,mu2)*G8 * ( - 5./2.*mp2*s3* propkp - 5./6.*mp2 + 2*pow(mp2,2)*propkp - 2*mk2 + 7./3.*s3 ); K3Pp4Ix10331 += + Bb(mp2,mk2,s3,mu2)*G27 * ( 5./2.*mp2*s3*propkp + 95./18.*mp2 - 2*pow(mp2,2)*propkp - 4./3.*mk2 - 1./9.*s3 ); K3Pp4Ix10331 += + Bb(mk2,mk2,s3,mu2)*G8 * ( - 5./6.*mp2*s3* propkp + 5./6.*s3 ); K3Pp4Ix10331 += + Bb(mk2,mk2,s3,mu2)*G27 * ( - 5./9.*mp2*s3* propkp + 5./9.*s3 ); K3Pp4Ix10331 += + Bb(mg2,mp2,mp2,mu2)*G8 * ( 6*mp2 + 10./3.*mk2 - 20./3.*s3 ); K3Pp4Ix10331 += + Bb(mg2,mp2,mp2,mu2)*G27 * ( - 2./3.*mp2*s3* propkp + 8./3.*mp2 + 4*pow(mp2,2)*s3*pow(propkp,2) + 20./3.* pow(mp2,2)*propkp - 16./3.*pow(mp2,3)*pow(propkp,2) + 20./9.* mk2 - 70./9.*s3 ); K3Pp4Ix10331 += + Bb(mg2,mp2,mk2,mu2)*G8 * ( 2*mp2*s3*propkp - 2. /3.*mp2 - 2*pow(mp2,2)*propkp ); K3Pp4Ix10331 += + Bb(mg2,mp2,mk2,mu2)*G27 * ( 4./3.*mp2*s3* propkp - 4./9.*mp2 - 4./3.*pow(mp2,2)*propkp ); K3Pp4Ix10331 += + Bb(mg2,mk2,mp2,mu2)*G8 * ( - 1./2.*mp2*s3* propkp - 1./6.*mk2 + 1./2.*s3 ); K3Pp4Ix10331 += + Bb(mg2,mk2,mp2,mu2)*G27 * ( - 1./3.*mp2*s3* propkp - 1./9.*mk2 + 1./3.*s3 ); K3Pp4Ix10331 += + Bb(mg2,mk2,mk2,mu2)*G8 * ( 17./6.*mp2 + 5./6.* mk2 - 19./6.*s3 ); K3Pp4Ix10331 += + Bb(mg2,mk2,mk2,mu2)*G27 * ( 22./3.*mp2*s3* propkp + 59./9.*mp2 - 4*pow(mp2,2)*s3*pow(propkp,2) - 20./3.* pow(mp2,2)*propkp + 16./3.*pow(mp2,3)*pow(propkp,2) + 5./9.* mk2 - 79./9.*s3 ); K3Pp4Ix10331 += + B1b(mp2,mk2,s3,mu2)*G8 * ( - 2./3.*mp2 + 8./3. *mk2 ); K3Pp4Ix10331 += + B1b(mp2,mk2,s3,mu2)*G27 * ( - 34./9.*mp2 + 16. /9.*mk2 - 5./3.*s3 ); K3Pp4Ix10331 += + B1b(mg2,mp2,mp2,mu2)*G8 * ( - 25./3.*mp2 - pow(mp2,2)*s3*pow(propkp,2) - 1./3.*pow(mp2,2)*propkp - 2./3.* mk2 + 17./3.*s3 ); K3Pp4Ix10331 += + B1b(mg2,mp2,mp2,mu2)*G27 * ( - 23./3.*mp2*s3* propkp - 128./9.*mp2 - 14./3.*pow(mp2,2)*s3*pow(propkp,2) + 28. /9.*pow(mp2,2)*propkp + 16./3.*pow(mp2,3)*pow(propkp,2) - 34./ 9.*mk2 + 139./9.*s3 ); K3Pp4Ix10331 += + B1b(mg2,mp2,mk2,mu2)*G8 * ( - 4*mp2*s3*propkp + 4./3.*mp2 + 4*pow(mp2,2)*s3*pow(propkp,2) + 8./3.*pow( mp2,2)*propkp - 4*pow(mp2,3)*pow(propkp,2) ); K3Pp4Ix10331 += + B1b(mg2,mp2,mk2,mu2)*G27 * ( - 8./3.*mp2*s3* propkp + 8./9.*mp2 + 8./3.*pow(mp2,2)*s3*pow(propkp,2) + 16./9. *pow(mp2,2)*propkp - 8./3.*pow(mp2,3)*pow(propkp,2) ); K3Pp4Ix10331 += + B1b(mg2,mk2,mp2,mu2)*G8 * ( - mp2*s3*propkp - 1./3.*mp2 + pow(mp2,2)*s3*pow(propkp,2) + 1./3.*pow(mp2,2)* propkp ); K3Pp4Ix10331 += + B1b(mg2,mk2,mp2,mu2)*G27 * ( - 2./3.*mp2*s3* propkp - 2./9.*mp2 + 2./3.*pow(mp2,2)*s3*pow(propkp,2) + 2./9. *pow(mp2,2)*propkp ); K3Pp4Ix10331 += + B1b(mg2,mk2,mk2,mu2)*G8 * ( 8*mp2*s3*propkp - 1./2.*mp2 - 4*pow(mp2,2)*s3*pow(propkp,2) - 20./3.*pow(mp2,2)* propkp + 4*pow(mp2,3)*pow(propkp,2) - 5./2.*mk2 - 5./6.*s3 ); K3Pp4Ix10331 += + B1b(mg2,mk2,mk2,mu2)*G27 * ( - 16./3.*mp2*s3* propkp - 20./3.*mp2 + 4./3.*pow(mp2,2)*s3*pow(propkp,2) + 50./ 9.*pow(mp2,2)*propkp - 8./3.*pow(mp2,3)*pow(propkp,2) - 10./3. *mk2 + 70./9.*s3 ); K3Pp4Ix10331 += + Cb(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( 2*mp2* mk2 - 14./3.*mp2*s3 + 16./3.*pow(mp2,2) ); K3Pp4Ix10331 += + Cb(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( 4./3.* mp2*mk2 - 28./9.*mp2*s3 + 32./9.*pow(mp2,2) ); K3Pp4Ix10331 += + Cb(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( 2*mp2* mk2 - 4*mp2*s3 + 2*pow(mp2,2) - 2*mk2*s3 + 2*pow(s3,2) ); K3Pp4Ix10331 += + Cb(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( 29./3.* mp2*mk2 - 58./3.*mp2*s3 - 5*mp2*pow(s3,2)*propkp + 50./3.*pow( mp2,2)*s3*propkp + 44./3.*pow(mp2,2) - 40./3.*pow(mp2,3)* propkp - 8*mk2*s3 + 5./3.*pow(mk2,2) + 19./3.*pow(s3,2) ); K3Pp4Ix10331 += + Cb(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G8 * ( 1./2.* mp2*mk2 - 7./6.*mk2*s3 + 5./6.*pow(mk2,2) ); K3Pp4Ix10331 += + Cb(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G27 * ( 1./3.* mp2*mk2 - 7./9.*mk2*s3 + 5./9.*pow(mk2,2) ); K3Pp4Ix10331 += + C11b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( - 8./ 3.*mp2*mk2 + 22./3.*mp2*s3 + 14./3.*mp2*pow(s3,2)*propkp - 50./ 3.*pow(mp2,2)*s3*propkp - 14*pow(mp2,2) + 44./3.*pow(mp2,3)* propkp ); K3Pp4Ix10331 += + C11b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( - 16. /9.*mp2*mk2 + 44./9.*mp2*s3 + 28./9.*mp2*pow(s3,2)*propkp - 100./9.*pow(mp2,2)*s3*propkp - 28./3.*pow(mp2,2) + 88./9.*pow( mp2,3)*propkp ); K3Pp4Ix10331 += + C11b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( 29./3. *mp2*mk2 - 38./3.*mp2*s3 - 14./3.*mp2*pow(s3,2)*propkp + 43./3. *pow(mp2,2)*s3*propkp + 32./3.*pow(mp2,2) - 12*pow(mp2,3)* propkp - 7*mk2*s3 + 10./3.*pow(mk2,2) + 11./3.*pow(s3,2) ); K3Pp4Ix10331 += + C11b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( - 19. /18.*mp2*mk2 + 74./9.*mp2*s3 + 19./18.*mp2*pow(s3,2)*propkp - 124./9.*pow(mp2,2)*s3*propkp - 116./9.*pow(mp2,2) + 56./3.* pow(mp2,3)*propkp + 2*mk2*s3 - 5./18.*pow(mk2,2) - 31./18.* pow(s3,2) ); K3Pp4Ix10331 += + C11b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G8 * ( 5./6.* mp2*mk2 - 7./3.*mp2*s3 + 7./3.*pow(mp2,2)*s3*propkp + 8./3.* pow(mp2,2) - 8./3.*pow(mp2,3)*propkp + 1./2.*mk2*s3 - 1./2.* pow(mk2,2) ); K3Pp4Ix10331 += + C11b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G27 * ( 5./9. *mp2*mk2 - 14./9.*mp2*s3 + 14./9.*pow(mp2,2)*s3*propkp + 16./9. *pow(mp2,2) - 16./9.*pow(mp2,3)*propkp + 1./3.*mk2*s3 - 1./3.* pow(mk2,2) ); K3Pp4Ix10331 += + C12b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( - 14./ 3.*mp2*mk2 + 4*mp2*s3 - 14./3.*mp2*pow(s3,2)*propkp + 22./3.* pow(mp2,2)*s3*propkp - 8./3.*pow(mp2,2) ); K3Pp4Ix10331 += + C12b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( - 28. /9.*mp2*mk2 + 8./3.*mp2*s3 - 28./9.*mp2*pow(s3,2)*propkp + 44./ 9.*pow(mp2,2)*s3*propkp - 16./9.*pow(mp2,2) ); K3Pp4Ix10331 += + C12b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( 3./2.* mp2*mk2 + 7*mp2*s3 + 35./6.*mp2*pow(s3,2)*propkp - 26./3.*pow( mp2,2)*s3*propkp + 7./6.*pow(mp2,2) + 2./3.*mk2*s3 + 10./3.* pow(mk2,2) - 14./3.*pow(s3,2) ); K3Pp4Ix10331 += + C12b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( - 22. /3.*mp2*mk2 + 14./3.*mp2*s3 + 35./9.*mp2*pow(s3,2)*propkp - 52. /9.*pow(mp2,2)*s3*propkp - 23./9.*pow(mp2,2) + 19./9.*mk2*s3 + 5./9.*pow(mk2,2) + 2./9.*pow(s3,2) ); K3Pp4Ix10331 += + C12b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G8 * ( 7./6.* mp2*mk2 - 4./3.*mp2*s3 - 7./6.*mp2*pow(s3,2)*propkp + 4./3.* pow(mp2,2)*s3*propkp - 7./3.*mk2*s3 + 1./6.*pow(mk2,2) + 7./6. *pow(s3,2) ); K3Pp4Ix10331 += + C12b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G27 * ( 7./9. *mp2*mk2 - 8./9.*mp2*s3 - 7./9.*mp2*pow(s3,2)*propkp + 8./9.* pow(mp2,2)*s3*propkp - 14./9.*mk2*s3 + 1./9.*pow(mk2,2) + 7./9. *pow(s3,2) ); K3Pp4Ix10331 += + C21b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( 2./3.* mp2*mk2 - 4./3.*mp2*s3 - 2./3.*mp2*pow(s3,2)*propkp + 20./3.* pow(mp2,2)*s3*propkp + 6*pow(mp2,2) - 32./3.*pow(mp2,3)*propkp ); K3Pp4Ix10331 += + C21b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( 4./9. *mp2*mk2 - 8./9.*mp2*s3 - 4./9.*mp2*pow(s3,2)*propkp + 40./9.* pow(mp2,2)*s3*propkp + 4*pow(mp2,2) - 64./9.*pow(mp2,3)*propkp ); K3Pp4Ix10331 += + C21b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( - 6* mp2*mk2 + 20./3.*mp2*s3 + 2./3.*mp2*pow(s3,2)*propkp - 17./3.* pow(mp2,2)*s3*propkp - 29./3.*pow(mp2,2) + 8*pow(mp2,3)*propkp + 4./3.*mk2*s3 - 2./3.*pow(mk2,2) - 2./3.*pow(s3,2) ); K3Pp4Ix10331 += + C21b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( - 3./ 2.*mp2*mk2 + 10./9.*mp2*s3 - 7./18.*mp2*pow(s3,2)*propkp + 26./ 9.*pow(mp2,2)*s3*propkp - 28./9.*pow(mp2,2) - 8*pow(mp2,3)* propkp - 7./9.*mk2*s3 + 7./18.*pow(mk2,2) + 7./18.*pow(s3,2) ) ; K3Pp4Ix10331 += + C21b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G8 * ( - mp2 *mk2 + mp2*s3 - pow(mp2,2)*s3*propkp - 8./3.*pow(mp2,2) + 8./3. *pow(mp2,3)*propkp ); K3Pp4Ix10331 += + C21b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G27 * ( - 2./ 3.*mp2*mk2 + 2./3.*mp2*s3 - 2./3.*pow(mp2,2)*s3*propkp - 16./9. *pow(mp2,2) + 16./9.*pow(mp2,3)*propkp ); K3Pp4Ix10331 += + C22b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( 8./3.* mp2*mk2 + 4./3.*mp2*s3 + 4./3.*mp2*pow(s3,2)*propkp - 8./3.* pow(mp2,2) ); K3Pp4Ix10331 += + C22b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( 16./9. *mp2*mk2 + 8./9.*mp2*s3 + 8./9.*mp2*pow(s3,2)*propkp - 16./9.* pow(mp2,2) ); K3Pp4Ix10331 += + C22b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( 10./3. *mp2*mk2 - 1./3.*mp2*s3 - mp2*pow(s3,2)*propkp - 2./3.*pow( mp2,2) - 4./3.*mk2*s3 - 8./3.*pow(mk2,2) + 4./3.*pow(s3,2) ); K3Pp4Ix10331 += + C22b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( 50./9. *mp2*mk2 - 32./9.*mp2*s3 - 2./3.*mp2*pow(s3,2)*propkp - 34./9. *pow(mp2,2) + 22./9.*mk2*s3 - 16./9.*pow(mk2,2) + 8./9.*pow( s3,2) ); K3Pp4Ix10331 += + C22b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G8 * ( 2./3.* mp2*mk2 - 1./3.*mp2*pow(s3,2)*propkp + 1./3.*mk2*s3 - 2./3.* pow(mk2,2) + 1./3.*pow(s3,2) ); K3Pp4Ix10331 += + C22b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G27 * ( 4./9. *mp2*mk2 - 2./9.*mp2*pow(s3,2)*propkp + 2./9.*mk2*s3 - 4./9.* pow(mk2,2) + 2./9.*pow(s3,2) ); K3Pp4Ix10331 += + C23b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( 10./3. *mp2*mk2 - 4*mp2*s3 - 2./3.*mp2*pow(s3,2)*propkp - 8./3.*pow( mp2,2)*s3*propkp + 22./3.*pow(mp2,2) ); K3Pp4Ix10331 += + C23b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( 20./9. *mp2*mk2 - 8./3.*mp2*s3 - 4./9.*mp2*pow(s3,2)*propkp - 16./9.* pow(mp2,2)*s3*propkp + 44./9.*pow(mp2,2) ); K3Pp4Ix10331 += + C23b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( - 41./ 6.*mp2*mk2 - 11./3.*mp2*s3 + 1./6.*mp2*pow(s3,2)*propkp + 10./ 3.*pow(mp2,2)*s3*propkp + 13./6.*pow(mp2,2) + 4*mk2*s3 - 10./3. *pow(mk2,2) - 2./3.*pow(s3,2) ); K3Pp4Ix10331 += + C23b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( - 41. /9.*mp2*mk2 + 8./9.*mp2*s3 + 1./9.*mp2*pow(s3,2)*propkp + 20./ 9.*pow(mp2,2)*s3*propkp + 118./9.*pow(mp2,2) + 8./3.*mk2*s3 - 5./9.*pow(mk2,2) - 19./9.*pow(s3,2) ); K3Pp4Ix10331 += + C23b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G8 * ( - 13./ 6.*mp2*mk2 + 2./3.*mp2*s3 + 1./2.*mp2*pow(s3,2)*propkp - 2./3. *pow(mp2,2)*s3*propkp + mk2*s3 - 1./2.*pow(mk2,2) - 1./2.*pow( s3,2) ); K3Pp4Ix10331 += + C23b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G27 * ( - 13. /9.*mp2*mk2 + 4./9.*mp2*s3 + 1./3.*mp2*pow(s3,2)*propkp - 4./9. *pow(mp2,2)*s3*propkp + 2./3.*mk2*s3 - 1./3.*pow(mk2,2) - 1./3. *pow(s3,2) ); K3Pp4Ix10331 += + C24b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( 4*mp2* s3*propkp + 20./3.*mp2 - 32./3.*pow(mp2,2)*propkp ); K3Pp4Ix10331 += + C24b(mg2,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( 8./3. *mp2*s3*propkp + 40./9.*mp2 - 64./9.*pow(mp2,2)*propkp ); K3Pp4Ix10331 += + C24b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( - 3* mp2*s3*propkp - 29./3.*mp2 + 8*pow(mp2,2)*propkp - 20./3.*mk2 + 4*s3 ); K3Pp4Ix10331 += + C24b(mg2,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( 4./3. *mp2*s3*propkp - 28./9.*mp2 - 8*pow(mp2,2)*propkp - 10./9.*mk2 + 8./3.*s3 ); K3Pp4Ix10331 += + C24b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G8 * ( - mp2 *s3*propkp - 8./3.*mp2 + 8./3.*pow(mp2,2)*propkp - mk2 + s3 ); K3Pp4Ix10331 += + C24b(mg2,mk2,mk2,mp2,s3,mk2,mu2)*G27 * ( - 2./ 3.*mp2*s3*propkp - 16./9.*mp2 + 16./9.*pow(mp2,2)*propkp - 2./ 3.*mk2 + 2./3.*s3 ); K3Pp4Ix10331 += + dqsqBb(mg2,mp2,mp2,mu2)*G8 * ( - 8*mp2*s3 + 8 *pow(mp2,2) ); K3Pp4Ix10331 += + dqsqBb(mg2,mp2,mp2,mu2)*G27 * ( 20./3.*mp2*mk2 - 76./3.*mp2*s3 + 20*pow(mp2,2)*s3*propkp + 32*pow(mp2,2) - 80./3.*pow(mp2,3)*propkp ); K3Pp4Ix10331 += + dqsqBb(mg2,mk2,mk2,mu2)*G8 * ( 4*mp2*mk2 - 4* mk2*s3 ); K3Pp4Ix10331 += + dqsqBb(mg2,mk2,mk2,mu2)*G27 * ( 16*mp2*mk2 - 10*mp2*s3 + 10*pow(mp2,2)*s3*propkp + 40./3.*pow(mp2,2) - 40./ 3.*pow(mp2,3)*propkp - 38./3.*mk2*s3 + 10./3.*pow(mk2,2) ); K3Pp4Ix10331 += + dqsqB1b(mg2,mp2,mp2,mu2)*G8 * ( 8*mp2*s3 - 8* pow(mp2,2) ); K3Pp4Ix10331 += + dqsqB1b(mg2,mp2,mp2,mu2)*G27 * ( - 20./3.*mp2 *mk2 + 76./3.*mp2*s3 - 20*pow(mp2,2)*s3*propkp - 32*pow(mp2,2) + 80./3.*pow(mp2,3)*propkp ); K3Pp4Ix10331 += + dqsqB1b(mg2,mk2,mk2,mu2)*G8 * ( - 4*mp2*mk2 + 4*mk2*s3 ); K3Pp4Ix10331 += + dqsqB1b(mg2,mk2,mk2,mu2)*G27 * ( - 16*mp2*mk2 + 10*mp2*s3 - 10*pow(mp2,2)*s3*propkp - 40./3.*pow(mp2,2) + 40./3.*pow(mp2,3)*propkp + 38./3.*mk2*s3 - 10./3.*pow(mk2,2) ) ; return K3Pp4Ix10331*e2/(pow(fpip,3)*fkp); } dcomplex Ap4IKPppm(double s1, double s2,double s3){ using namespace inputpars; double mko2 = mko*mko, mkp2 = mkp*mkp,mpp2 = mpp*mpp,mpo2=mpo*mpo, me2 = me*me,mu2 = mu*mu; double mp2 = mpp2,mk2=mkp2; double propkp = prop(mk2,mp2); dcomplex K3Pp4Ix10112=dcomplex(0.); K3Pp4Ix10112 += + Ab(mp2,mu2)*G8 * ( 1./3. + 2*mp2*propkp ); K3Pp4Ix10112 += + Ab(mp2,mu2)*G27 * ( 92./9. + 4./3.*mp2*propkp ); K3Pp4Ix10112 += + Ab(mk2,mu2)*G8 * ( 7./3. - 2*mp2*propkp ); K3Pp4Ix10112 += + Ab(mk2,mu2)*G27 * ( 44./9. - 4./3.*mp2*propkp ); K3Pp4Ix10112 += + Bb(mp2,mp2,s1,mu2)*G8 * ( 5./2.*mp2*s1*propkp + 1./3.*mp2 - 1./2.*mk2 - 3./2.*s1 ); K3Pp4Ix10112 += + Bb(mp2,mp2,s1,mu2)*G27 * ( 5./3.*mp2*s1*propkp - 41./18.*mp2 - 7./6.*mk2 + 19./6.*s1 ); K3Pp4Ix10112 += + Bb(mp2,mp2,s2,mu2)*G8 * ( 5./2.*mp2*s2*propkp + 1./3.*mp2 - 1./2.*mk2 - 3./2.*s2 ); K3Pp4Ix10112 += + Bb(mp2,mp2,s2,mu2)*G27 * ( 5./3.*mp2*s2*propkp - 41./18.*mp2 - 7./6.*mk2 + 19./6.*s2 ); K3Pp4Ix10112 += + Bb(mp2,mp2,s3,mu2)*G8 * ( 20./3.*mp2*s3*propkp + 7./3.*mp2 - 8*pow(mp2,2)*propkp + mk2 - 5./3.*s3 ); K3Pp4Ix10112 += + Bb(mp2,mp2,s3,mu2)*G27 * ( 40./9.*mp2*s3* propkp - 31./9.*mp2 - 16./3.*pow(mp2,2)*propkp - mk2 + 65./9.* s3 ); K3Pp4Ix10112 += + Bb(mp2,mk2,s1,mu2)*G8 * ( - 5./3.*mp2*s1* propkp + 2./3.*mp2 - 2*mk2 + 2*s1 ); K3Pp4Ix10112 += + Bb(mp2,mk2,s1,mu2)*G27 * ( - 10./9.*mp2*s1* propkp - 7./18.*mp2 - 23./6.*mk2 + 11./2.*s1 ); K3Pp4Ix10112 += + Bb(mp2,mk2,s2,mu2)*G8 * ( - 5./3.*mp2*s2* propkp + 2./3.*mp2 - 2*mk2 + 2*s2 ); K3Pp4Ix10112 += + Bb(mp2,mk2,s2,mu2)*G27 * ( - 10./9.*mp2*s2* propkp - 7./18.*mp2 - 23./6.*mk2 + 11./2.*s2 ); K3Pp4Ix10112 += + Bb(mp2,mk2,s3,mu2)*G8 * ( - 20./3.*mp2*s3* propkp - 7*mp2 + 8*pow(mp2,2)*propkp - 3*mk2 + 17./3.*s3 ); K3Pp4Ix10112 += + Bb(mp2,mk2,s3,mu2)*G27 * ( - 40./9.*mp2*s3* propkp - 19./3.*mp2 + 16./3.*pow(mp2,2)*propkp - 7*mk2 + 109./ 9.*s3 ); K3Pp4Ix10112 += + Bb(mk2,mk2,s1,mu2)*G8 * ( - 5./6.*mp2*s1* propkp + 5./6.*s1 ); K3Pp4Ix10112 += + Bb(mk2,mk2,s1,mu2)*G27 * ( - 5./9.*mp2*s1* propkp + 5./9.*s1 ); K3Pp4Ix10112 += + Bb(mk2,mk2,s2,mu2)*G8 * ( - 5./6.*mp2*s2* propkp + 5./6.*s2 ); K3Pp4Ix10112 += + Bb(mk2,mk2,s2,mu2)*G27 * ( - 5./9.*mp2*s2* propkp + 5./9.*s2 ); K3Pp4Ix10112 += + Bb(0.,mp2,mp2,mu2)*G8 * ( - 6*mp2 + 8./3.* mk2 - 2*s3 ); K3Pp4Ix10112 += + Bb(0.,mp2,mp2,mu2)*G27 * ( 118./3.*mp2 - 14./ 9.*mk2 + 26./3.*s3 ); K3Pp4Ix10112 += + Bb(0.,mp2,mk2,mu2)*G8 * ( - 2*mp2*s3*propkp - 10./3.*mp2 + 4*pow(mp2,2)*propkp ); K3Pp4Ix10112 += + Bb(0.,mp2,mk2,mu2)*G27 * ( - 4./3.*mp2*s3* propkp - 20./9.*mp2 + 8./3.*pow(mp2,2)*propkp ); K3Pp4Ix10112 += + Bb(0.,mk2,mp2,mu2)*G8 * ( 2*mp2*s3*propkp + 4 *mp2 - 4*pow(mp2,2)*propkp - 4./3.*mk2 - 2*s3 ); K3Pp4Ix10112 += + Bb(0.,mk2,mp2,mu2)*G27 * ( 4./3.*mp2*s3* propkp + 8./3.*mp2 - 8./3.*pow(mp2,2)*propkp - 8./9.*mk2 - 4./ 3.*s3 ); K3Pp4Ix10112 += + Bb(0.,mk2,mk2,mu2)*G8 * ( - 4./3.*mp2 - 2./3. *mk2 + 2./3.*s3 ); K3Pp4Ix10112 += + Bb(0.,mk2,mk2,mu2)*G27 * ( - 38./9.*mp2 + 146./9.*mk2 + 34./9.*s3 ); K3Pp4Ix10112 += + B1b(mp2,mk2,s1,mu2)*G8 * ( - 4./3.*mp2 + 3* mk2 + 2./3.*s1 ); K3Pp4Ix10112 += + B1b(mp2,mk2,s1,mu2)*G27 * ( - 38./9.*mp2 + 16. /3.*mk2 + 4./9.*s1 ); K3Pp4Ix10112 += + B1b(mp2,mk2,s2,mu2)*G8 * ( - 4./3.*mp2 + 3* mk2 + 2./3.*s2 ); K3Pp4Ix10112 += + B1b(mp2,mk2,s2,mu2)*G27 * ( - 38./9.*mp2 + 16. /3.*mk2 + 4./9.*s2 ); K3Pp4Ix10112 += + B1b(mp2,mk2,s3,mu2)*G8 * ( 4./3.*mp2 + 4./3.* mk2 - 4./3.*s3 ); K3Pp4Ix10112 += + B1b(mp2,mk2,s3,mu2)*G27 * ( - 52./9.*mp2 + 68. /9.*mk2 - 8./9.*s3 ); K3Pp4Ix10112 += + B1b(0.,mp2,mp2,mu2)*G8 * ( 28./3.*mp2 + 4* pow(mp2,2)*s3*pow(propkp,2) - 8./3.*pow(mp2,2)*propkp - 8*pow( mp2,3)*pow(propkp,2) + 26./3.*mk2 - 28./3.*s3 ); K3Pp4Ix10112 += + B1b(0.,mp2,mp2,mu2)*G27 * ( - 769./9.*mp2 + 8./3.*pow(mp2,2)*s3*pow(propkp,2) - 16./9.*pow(mp2,2)*propkp - 16./3.*pow(mp2,3)*pow(propkp,2) - 23./9.*mk2 + 169./9.*s3 ) ; K3Pp4Ix10112 += + B1b(0.,mp2,mk2,mu2)*G8 * ( 4*mp2*s3*propkp + 20./3.*mp2 - 4*pow(mp2,2)*s3*pow(propkp,2) - 44./3.*pow(mp2,2) *propkp + 8*pow(mp2,3)*pow(propkp,2) ); K3Pp4Ix10112 += + B1b(0.,mp2,mk2,mu2)*G27 * ( 8./3.*mp2*s3* propkp + 40./9.*mp2 - 8./3.*pow(mp2,2)*s3*pow(propkp,2) - 88./ 9.*pow(mp2,2)*propkp + 16./3.*pow(mp2,3)*pow(propkp,2) ); K3Pp4Ix10112 += + B1b(0.,mk2,mp2,mu2)*G8 * ( 4*mp2*s3*propkp - 8./3.*mp2 - 4*pow(mp2,2)*s3*pow(propkp,2) - 16./3.*pow(mp2,2)* propkp + 8*pow(mp2,3)*pow(propkp,2) ); K3Pp4Ix10112 += + B1b(0.,mk2,mp2,mu2)*G27 * ( 8./3.*mp2*s3* propkp - 16./9.*mp2 - 8./3.*pow(mp2,2)*s3*pow(propkp,2) - 32./ 9.*pow(mp2,2)*propkp + 16./3.*pow(mp2,3)*pow(propkp,2) ); K3Pp4Ix10112 += + B1b(0.,mk2,mk2,mu2)*G8 * ( - 8*mp2*s3*propkp - 12*mp2 + 4*pow(mp2,2)*s3*pow(propkp,2) + 68./3.*pow(mp2,2)* propkp - 8*pow(mp2,3)*pow(propkp,2) - 6*mk2 + 4./3.*s3 ); K3Pp4Ix10112 += + B1b(0.,mk2,mk2,mu2)*G27 * ( - 16./3.*mp2*s3* propkp - 49./3.*mp2 + 8./3.*pow(mp2,2)*s3*pow(propkp,2) + 136./ 9.*pow(mp2,2)*propkp - 16./3.*pow(mp2,3)*pow(propkp,2) - 29* mk2 + 83./9.*s3 ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s1,mp2,mu2)*G8 * ( 2./3.* mp2*mk2 + mp2*s1*s2*propkp - mp2*s1*s3*propkp - 8./3.*mp2*s1 + 7./3.*mp2*s2 + mp2*pow(s2,2)*propkp - 5./3.*mp2*s3 - mp2* pow(s3,2)*propkp - 4*pow(mp2,2)*s2*propkp + 4*pow(mp2,2)*s3* propkp + pow(mp2,2) - 2*mk2*s3 + pow(mk2,2) + pow(s1,2) - pow( s2,2) + pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s1,mp2,mu2)*G27 * ( - 26./ 9.*mp2*mk2 + 2./3.*mp2*s1*s2*propkp - 2./3.*mp2*s1*s3*propkp + 74./9.*mp2*s1 - 76./9.*mp2*s2 + 2./3.*mp2*pow(s2,2)*propkp - 10./9.*mp2*s3 - 2./3.*mp2*pow(s3,2)*propkp - 8./3.*pow( mp2,2)*s2*propkp + 8./3.*pow(mp2,2)*s3*propkp + 17./3.*pow( mp2,2) + 10./3.*mk2*s1 - 10./3.*mk2*s2 + 16./3.*mk2*s3 - pow( mk2,2) - 13./3.*pow(s1,2) + 13./3.*pow(s2,2) - 13./3.*pow( s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s1,mk2,mu2)*G8 * ( - 1./3. *mp2*mk2 + 2*mp2*s1*s2*propkp - 2*mp2*s1*s3*propkp - mp2*s1 + 7./3.*mp2*s2 + 2*mp2*pow(s2,2)*propkp + 1./3.*mp2*s3 - 2*mp2* pow(s3,2)*propkp - 8*pow(mp2,2)*s2*propkp + 8*pow(mp2,2)*s3* propkp - 7./3.*pow(mp2,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s1,mk2,mu2)*G27 * ( - 2./9. *mp2*mk2 + 4./3.*mp2*s1*s2*propkp - 4./3.*mp2*s1*s3*propkp - 2. /3.*mp2*s1 + 14./9.*mp2*s2 + 4./3.*mp2*pow(s2,2)*propkp + 2./9. *mp2*s3 - 4./3.*mp2*pow(s3,2)*propkp - 16./3.*pow(mp2,2)*s2* propkp + 16./3.*pow(mp2,2)*s3*propkp - 14./9.*pow(mp2,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s2,mp2,mu2)*G8 * ( 2./3.* mp2*mk2 + mp2*s1*s2*propkp + 7./3.*mp2*s1 + mp2*pow(s1,2)* propkp - mp2*s2*s3*propkp - 8./3.*mp2*s2 - 5./3.*mp2*s3 - mp2* pow(s3,2)*propkp - 4*pow(mp2,2)*s1*propkp + 4*pow(mp2,2)*s3* propkp + pow(mp2,2) - 2*mk2*s3 + pow(mk2,2) - pow(s1,2) + pow( s2,2) + pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s2,mp2,mu2)*G27 * ( - 26./ 9.*mp2*mk2 + 2./3.*mp2*s1*s2*propkp - 76./9.*mp2*s1 + 2./3.* mp2*pow(s1,2)*propkp - 2./3.*mp2*s2*s3*propkp + 74./9.*mp2*s2 - 10./9.*mp2*s3 - 2./3.*mp2*pow(s3,2)*propkp - 8./3.*pow( mp2,2)*s1*propkp + 8./3.*pow(mp2,2)*s3*propkp + 17./3.*pow( mp2,2) - 10./3.*mk2*s1 + 10./3.*mk2*s2 + 16./3.*mk2*s3 - pow( mk2,2) + 13./3.*pow(s1,2) - 13./3.*pow(s2,2) - 13./3.*pow( s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s2,mk2,mu2)*G8 * ( - 1./3. *mp2*mk2 + 2*mp2*s1*s2*propkp + 7./3.*mp2*s1 + 2*mp2*pow(s1,2) *propkp - 2*mp2*s2*s3*propkp - mp2*s2 + 1./3.*mp2*s3 - 2*mp2* pow(s3,2)*propkp - 8*pow(mp2,2)*s1*propkp + 8*pow(mp2,2)*s3* propkp - 7./3.*pow(mp2,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s2,mk2,mu2)*G27 * ( - 2./9. *mp2*mk2 + 4./3.*mp2*s1*s2*propkp + 14./9.*mp2*s1 + 4./3.*mp2* pow(s1,2)*propkp - 4./3.*mp2*s2*s3*propkp - 2./3.*mp2*s2 + 2./ 9.*mp2*s3 - 4./3.*mp2*pow(s3,2)*propkp - 16./3.*pow(mp2,2)*s1* propkp + 16./3.*pow(mp2,2)*s3*propkp - 14./9.*pow(mp2,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( 4*mp2* mk2 - 6*mp2*s3 + 4*pow(mp2,2) - 2*mk2*s3 + 2*pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( - 4* mp2*mk2 + 26*mp2*s3 - 52./3.*pow(mp2,2) + 2*mk2*s3 - 26./3.* pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( 2*mp2* mk2 - 14./3.*mp2*s3 + 22./3.*pow(mp2,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( 4./3.* mp2*mk2 - 28./9.*mp2*s3 + 44./9.*pow(mp2,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s1,mp2,mu2)*G8 * ( - 5./6. *mp2*mk2 - mp2*s1*s2*propkp + mp2*s1*s3*propkp - 4*mp2*s2 - mp2*pow(s2,2)*propkp + 4*mp2*s3 + mp2*pow(s3,2)*propkp + 4* pow(mp2,2)*s2*propkp - 4*pow(mp2,2)*s3*propkp - 1./2.*mk2*s1 - 5./6.*mk2*s2 + 13./6.*mk2*s3 - 1./2.*pow(mk2,2) + s1*s2 - s1*s3 + pow(s2,2) - pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s1,mp2,mu2)*G27 * ( - 5./9. *mp2*mk2 - 2./3.*mp2*s1*s2*propkp + 2./3.*mp2*s1*s3*propkp - 8. /3.*mp2*s2 - 2./3.*mp2*pow(s2,2)*propkp + 8./3.*mp2*s3 + 2./3. *mp2*pow(s3,2)*propkp + 8./3.*pow(mp2,2)*s2*propkp - 8./3.* pow(mp2,2)*s3*propkp - 1./3.*mk2*s1 - 5./9.*mk2*s2 + 13./9.* mk2*s3 - 1./3.*pow(mk2,2) + 2./3.*s1*s2 - 2./3.*s1*s3 + 2./3.* pow(s2,2) - 2./3.*pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s1,mk2,mu2)*G8 * ( 11./3.* mp2*mk2 - mp2*s1*s2*propkp + mp2*s1*s3*propkp - 5./3.*mp2*s1 - 14./3.*mp2*s2 - mp2*pow(s2,2)*propkp + 16./3.*mp2*s3 + mp2* pow(s3,2)*propkp + 4*pow(mp2,2)*s2*propkp - 4*pow(mp2,2)*s3* propkp - 4./3.*mk2*s1 - 10./3.*mk2*s2 + 2./3.*mk2*s3 + 1./3.* pow(mk2,2) + 2*s1*s2 - 2*s1*s3 + pow(s1,2) + pow(s2,2) - pow( s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s1,mk2,mu2)*G27 * ( 52./9.* mp2*mk2 - 2./3.*mp2*s1*s2*propkp + 2./3.*mp2*s1*s3*propkp + 80. /9.*mp2*s1 - 118./9.*mp2*s2 - 2./3.*mp2*pow(s2,2)*propkp + 122. /9.*mp2*s3 + 2./3.*mp2*pow(s3,2)*propkp + 8./3.*pow(mp2,2)*s2* propkp - 8./3.*pow(mp2,2)*s3*propkp - 5*pow(mp2,2) + 22./9.* mk2*s1 - 50./9.*mk2*s2 - 26./9.*mk2*s3 + 17./9.*pow(mk2,2) + 4. /3.*s1*s2 - 4./3.*s1*s3 - 13./3.*pow(s1,2) + 17./3.*pow(s2,2) - 17./3.*pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s2,mp2,mu2)*G8 * ( - 5./6. *mp2*mk2 - mp2*s1*s2*propkp - 4*mp2*s1 - mp2*pow(s1,2)*propkp + mp2*s2*s3*propkp + 4*mp2*s3 + mp2*pow(s3,2)*propkp + 4*pow( mp2,2)*s1*propkp - 4*pow(mp2,2)*s3*propkp - 5./6.*mk2*s1 - 1./ 2.*mk2*s2 + 13./6.*mk2*s3 - 1./2.*pow(mk2,2) + s1*s2 + pow( s1,2) - s2*s3 - pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s2,mp2,mu2)*G27 * ( - 5./9. *mp2*mk2 - 2./3.*mp2*s1*s2*propkp - 8./3.*mp2*s1 - 2./3.*mp2* pow(s1,2)*propkp + 2./3.*mp2*s2*s3*propkp + 8./3.*mp2*s3 + 2./ 3.*mp2*pow(s3,2)*propkp + 8./3.*pow(mp2,2)*s1*propkp - 8./3.* pow(mp2,2)*s3*propkp - 5./9.*mk2*s1 - 1./3.*mk2*s2 + 13./9.* mk2*s3 - 1./3.*pow(mk2,2) + 2./3.*s1*s2 + 2./3.*pow(s1,2) - 2./ 3.*s2*s3 - 2./3.*pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s2,mk2,mu2)*G8 * ( 11./3.* mp2*mk2 - mp2*s1*s2*propkp - 14./3.*mp2*s1 - mp2*pow(s1,2)* propkp + mp2*s2*s3*propkp - 5./3.*mp2*s2 + 16./3.*mp2*s3 + mp2 *pow(s3,2)*propkp + 4*pow(mp2,2)*s1*propkp - 4*pow(mp2,2)*s3* propkp - 10./3.*mk2*s1 - 4./3.*mk2*s2 + 2./3.*mk2*s3 + 1./3.* pow(mk2,2) + 2*s1*s2 + pow(s1,2) - 2*s2*s3 + pow(s2,2) - pow( s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s2,mk2,mu2)*G27 * ( 52./9.* mp2*mk2 - 2./3.*mp2*s1*s2*propkp - 118./9.*mp2*s1 - 2./3.*mp2* pow(s1,2)*propkp + 2./3.*mp2*s2*s3*propkp + 80./9.*mp2*s2 + 122./9.*mp2*s3 + 2./3.*mp2*pow(s3,2)*propkp + 8./3.*pow(mp2,2) *s1*propkp - 8./3.*pow(mp2,2)*s3*propkp - 5*pow(mp2,2) - 50./9. *mk2*s1 + 22./9.*mk2*s2 - 26./9.*mk2*s3 + 17./9.*pow(mk2,2) + 4./3.*s1*s2 + 17./3.*pow(s1,2) - 4./3.*s2*s3 - 13./3.*pow( s2,2) - 17./3.*pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( 20./3.* mp2*mk2 - 14./3.*mk2*s3 + 8./3.*pow(mk2,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( 40./9.* mp2*mk2 - 28./9.*mk2*s3 + 16./9.*pow(mk2,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( 4*mp2* mk2 - 4*mp2*s3 + 2*pow(mp2,2) - 4*mk2*s3 + 2*pow(mk2,2) + 2* pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( - 32./ 3.*mp2*mk2 + 52./3.*mp2*s3 - 26./3.*pow(mp2,2) + 32./3.*mk2*s3 - 2*pow(mk2,2) - 26./3.*pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mk2,mk2,mp2,s1,mk2,mu2)*G8 * ( 3./2.* mp2*mk2 - mp2*s1*s2*propkp + mp2*s1*s3*propkp - 4*mp2*s2 - mp2 *pow(s2,2)*propkp + 4*mp2*s3 + mp2*pow(s3,2)*propkp + 4*pow( mp2,2)*s2*propkp - 4*pow(mp2,2)*s3*propkp - 3./2.*mk2*s1 - 11./ 6.*mk2*s2 + 7./6.*mk2*s3 + 7./6.*pow(mk2,2) + s1*s2 - s1*s3 + pow(s2,2) - pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mk2,mk2,mp2,s1,mk2,mu2)*G27 * ( mp2*mk2 - 2./3.*mp2*s1*s2*propkp + 2./3.*mp2*s1*s3*propkp - 8./3.*mp2 *s2 - 2./3.*mp2*pow(s2,2)*propkp + 8./3.*mp2*s3 + 2./3.*mp2* pow(s3,2)*propkp + 8./3.*pow(mp2,2)*s2*propkp - 8./3.*pow( mp2,2)*s3*propkp - mk2*s1 - 11./9.*mk2*s2 + 7./9.*mk2*s3 + 7./ 9.*pow(mk2,2) + 2./3.*s1*s2 - 2./3.*s1*s3 + 2./3.*pow(s2,2) - 2./3.*pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mk2,mk2,mp2,s2,mk2,mu2)*G8 * ( 3./2.* mp2*mk2 - mp2*s1*s2*propkp - 4*mp2*s1 - mp2*pow(s1,2)*propkp + mp2*s2*s3*propkp + 4*mp2*s3 + mp2*pow(s3,2)*propkp + 4*pow( mp2,2)*s1*propkp - 4*pow(mp2,2)*s3*propkp - 11./6.*mk2*s1 - 3./ 2.*mk2*s2 + 7./6.*mk2*s3 + 7./6.*pow(mk2,2) + s1*s2 + pow( s1,2) - s2*s3 - pow(s3,2) ); K3Pp4Ix10112 += + Cb(0.,mk2,mk2,mp2,s2,mk2,mu2)*G27 * ( mp2*mk2 - 2./3.*mp2*s1*s2*propkp - 8./3.*mp2*s1 - 2./3.*mp2*pow(s1,2) *propkp + 2./3.*mp2*s2*s3*propkp + 8./3.*mp2*s3 + 2./3.*mp2* pow(s3,2)*propkp + 8./3.*pow(mp2,2)*s1*propkp - 8./3.*pow( mp2,2)*s3*propkp - 11./9.*mk2*s1 - mk2*s2 + 7./9.*mk2*s3 + 7./ 9.*pow(mk2,2) + 2./3.*s1*s2 + 2./3.*pow(s1,2) - 2./3.*s2*s3 - 2./3.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G8 * ( 16./3. *mp2*mk2 - mp2*s1*s2*propkp + mp2*s1*s3*propkp + 47./6.*mp2*s1 + 7./6.*mp2*pow(s1,2)*propkp - 41./6.*mp2*s2 - 3./2.*mp2*pow( s2,2)*propkp - 17./6.*mp2*s3 + 3./2.*mp2*pow(s3,2)*propkp - 8./ 3.*pow(mp2,2)*s1*propkp + 6*pow(mp2,2)*s2*propkp - 4*pow( mp2,2)*s3*propkp - 1./6.*pow(mp2,2) - 4./3.*pow(mp2,3)*propkp - mk2*s1 + mk2*s2 + 4*mk2*s3 - 5./2.*pow(mk2,2) - 13./6.*pow( s1,2) + 3./2.*pow(s2,2) - 3./2.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G27 * ( 107./ 9.*mp2*mk2 - 2./3.*mp2*s1*s2*propkp + 2./3.*mp2*s1*s3*propkp - 43./9.*mp2*s1 + 7./9.*mp2*pow(s1,2)*propkp + 49./9.*mp2*s2 - mp2*pow(s2,2)*propkp + 28./9.*mp2*s3 + mp2*pow(s3,2)*propkp - 16./9.*pow(mp2,2)*s1*propkp + 4*pow(mp2,2)*s2*propkp - 8./3. *pow(mp2,2)*s3*propkp - 167./18.*pow(mp2,2) - 8./9.*pow(mp2,3) *propkp - 22./3.*mk2*s1 + 22./3.*mk2*s2 - 17./3.*mk2*s3 - 5./6. *pow(mk2,2) + 49./18.*pow(s1,2) - 13./2.*pow(s2,2) + 13./2.* pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G8 * ( - 2* mp2*s1*s2*propkp + 2*mp2*s1*s3*propkp + 7./3.*mp2*s1 + 7./3.* mp2*pow(s1,2)*propkp - 7./3.*mp2*s2 - 3*mp2*pow(s2,2)*propkp - 1./3.*mp2*s3 + 3*mp2*pow(s3,2)*propkp - 16./3.*pow(mp2,2)* s1*propkp + 12*pow(mp2,2)*s2*propkp - 8*pow(mp2,2)*s3*propkp - 8./3.*pow(mp2,3)*propkp ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G27 * ( - 4./ 3.*mp2*s1*s2*propkp + 4./3.*mp2*s1*s3*propkp + 14./9.*mp2*s1 + 14./9.*mp2*pow(s1,2)*propkp - 14./9.*mp2*s2 - 2*mp2*pow( s2,2)*propkp - 2./9.*mp2*s3 + 2*mp2*pow(s3,2)*propkp - 32./9.* pow(mp2,2)*s1*propkp + 8*pow(mp2,2)*s2*propkp - 16./3.*pow( mp2,2)*s3*propkp - 16./9.*pow(mp2,3)*propkp ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G8 * ( 16./3. *mp2*mk2 - mp2*s1*s2*propkp - 41./6.*mp2*s1 - 3./2.*mp2*pow( s1,2)*propkp + mp2*s2*s3*propkp + 47./6.*mp2*s2 + 7./6.*mp2* pow(s2,2)*propkp - 17./6.*mp2*s3 + 3./2.*mp2*pow(s3,2)*propkp + 6*pow(mp2,2)*s1*propkp - 8./3.*pow(mp2,2)*s2*propkp - 4* pow(mp2,2)*s3*propkp - 1./6.*pow(mp2,2) - 4./3.*pow(mp2,3)* propkp + mk2*s1 - mk2*s2 + 4*mk2*s3 - 5./2.*pow(mk2,2) + 3./2. *pow(s1,2) - 13./6.*pow(s2,2) - 3./2.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G27 * ( 107./ 9.*mp2*mk2 - 2./3.*mp2*s1*s2*propkp + 49./9.*mp2*s1 - mp2*pow( s1,2)*propkp + 2./3.*mp2*s2*s3*propkp - 43./9.*mp2*s2 + 7./9.* mp2*pow(s2,2)*propkp + 28./9.*mp2*s3 + mp2*pow(s3,2)*propkp + 4*pow(mp2,2)*s1*propkp - 16./9.*pow(mp2,2)*s2*propkp - 8./3.* pow(mp2,2)*s3*propkp - 167./18.*pow(mp2,2) - 8./9.*pow(mp2,3)* propkp + 22./3.*mk2*s1 - 22./3.*mk2*s2 - 17./3.*mk2*s3 - 5./6. *pow(mk2,2) - 13./2.*pow(s1,2) + 49./18.*pow(s2,2) + 13./2.* pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G8 * ( - 2* mp2*s1*s2*propkp - 7./3.*mp2*s1 - 3*mp2*pow(s1,2)*propkp + 2* mp2*s2*s3*propkp + 7./3.*mp2*s2 + 7./3.*mp2*pow(s2,2)*propkp - 1./3.*mp2*s3 + 3*mp2*pow(s3,2)*propkp + 12*pow(mp2,2)*s1* propkp - 16./3.*pow(mp2,2)*s2*propkp - 8*pow(mp2,2)*s3*propkp - 8./3.*pow(mp2,3)*propkp ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G27 * ( - 4./ 3.*mp2*s1*s2*propkp - 14./9.*mp2*s1 - 2*mp2*pow(s1,2)*propkp + 4./3.*mp2*s2*s3*propkp + 14./9.*mp2*s2 + 14./9.*mp2*pow( s2,2)*propkp - 2./9.*mp2*s3 + 2*mp2*pow(s3,2)*propkp + 8*pow( mp2,2)*s1*propkp - 32./9.*pow(mp2,2)*s2*propkp - 16./3.*pow( mp2,2)*s3*propkp - 16./9.*pow(mp2,3)*propkp ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( - 4* mp2*mk2 + 31./3.*mp2*s3 + 14./3.*mp2*pow(s3,2)*propkp - 56./3. *pow(mp2,2)*s3*propkp - 12*pow(mp2,2) + 56./3.*pow(mp2,3)* propkp + mk2*s3 - 5./3.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( 4*mp2 *mk2 - 343./9.*mp2*s3 + 28./9.*mp2*pow(s3,2)*propkp - 112./9.* pow(mp2,2)*s3*propkp + 116./3.*pow(mp2,2) + 112./9.*pow(mp2,3) *propkp - mk2*s3 + 65./9.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( - 8./ 3.*mp2*mk2 + 22./3.*mp2*s3 + 14./3.*mp2*pow(s3,2)*propkp - 56./ 3.*pow(mp2,2)*s3*propkp - 16*pow(mp2,2) + 56./3.*pow(mp2,3)* propkp ); K3Pp4Ix10112 += + C11b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( - 16. /9.*mp2*mk2 + 44./9.*mp2*s3 + 28./9.*mp2*pow(s3,2)*propkp - 112./9.*pow(mp2,2)*s3*propkp - 32./3.*pow(mp2,2) + 112./9.* pow(mp2,3)*propkp ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G8 * ( 3./2.* mp2*mk2 + mp2*s1*s2*propkp - mp2*s1*s3*propkp - 8./3.*mp2*s1 - 7./6.*mp2*pow(s1,2)*propkp + 6*mp2*s2 + 3./2.*mp2*pow(s2,2) *propkp - 4*mp2*s3 - 3./2.*mp2*pow(s3,2)*propkp + 8./3.*pow( mp2,2)*s1*propkp - 6*pow(mp2,2)*s2*propkp + 4*pow(mp2,2)*s3* propkp - 4./3.*pow(mp2,2) + 4./3.*pow(mp2,3)*propkp - 4./3.* mk2*s1 + 2./3.*mk2*s2 - 10./3.*mk2*s3 + 11./6.*pow(mk2,2) - s1 *s2 + s1*s3 + 7./6.*pow(s1,2) - 3./2.*pow(s2,2) + 3./2.*pow( s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G27 * ( mp2* mk2 + 2./3.*mp2*s1*s2*propkp - 2./3.*mp2*s1*s3*propkp - 16./9. *mp2*s1 - 7./9.*mp2*pow(s1,2)*propkp + 4*mp2*s2 + mp2*pow( s2,2)*propkp - 8./3.*mp2*s3 - mp2*pow(s3,2)*propkp + 16./9.* pow(mp2,2)*s1*propkp - 4*pow(mp2,2)*s2*propkp + 8./3.*pow( mp2,2)*s3*propkp - 8./9.*pow(mp2,2) + 8./9.*pow(mp2,3)*propkp - 8./9.*mk2*s1 + 4./9.*mk2*s2 - 20./9.*mk2*s3 + 11./9.*pow( mk2,2) - 2./3.*s1*s2 + 2./3.*s1*s3 + 7./9.*pow(s1,2) - pow( s2,2) + pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G8 * ( - 5./ 2.*mp2*mk2 + 1./2.*mp2*s1*s2*propkp - 1./2.*mp2*s1*s3*propkp - 9./2.*mp2*s1 - 7./3.*mp2*pow(s1,2)*propkp + 26./3.*mp2*s2 + 3./2.*mp2*pow(s2,2)*propkp - 13./3.*mp2*s3 - 3./2.*mp2*pow( s3,2)*propkp + 4*pow(mp2,2)*s1*propkp - 6*pow(mp2,2)*s2*propkp + 4*pow(mp2,2)*s3*propkp - 13./3.*pow(mp2,2) + 4./3.*pow( mp2,3)*propkp - 11./3.*mk2*s1 + 13./3.*mk2*s2 + 1./3.*mk2*s3 + 5./6.*pow(mk2,2) - 2*s1*s2 + 2*s1*s3 + 17./6.*pow(s1,2) - 3. /2.*pow(s2,2) + 3./2.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G27 * ( - 10. /3.*mp2*mk2 + 1./3.*mp2*s1*s2*propkp - 1./3.*mp2*s1*s3*propkp - 49./3.*mp2*s1 - 14./9.*mp2*pow(s1,2)*propkp + 172./9.*mp2* s2 + mp2*pow(s2,2)*propkp - 86./9.*mp2*s3 - mp2*pow(s3,2)* propkp + 8./3.*pow(mp2,2)*s1*propkp - 4*pow(mp2,2)*s2*propkp + 8./3.*pow(mp2,2)*s3*propkp - 7./18.*pow(mp2,2) + 8./9.*pow( mp2,3)*propkp - 52./9.*mk2*s1 + 56./9.*mk2*s2 + 32./9.*mk2*s3 - 5./18.*pow(mk2,2) - 4./3.*s1*s2 + 4./3.*s1*s3 + 109./18.* pow(s1,2) - 17./2.*pow(s2,2) + 17./2.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G8 * ( 3./2.* mp2*mk2 + mp2*s1*s2*propkp + 6*mp2*s1 + 3./2.*mp2*pow(s1,2)* propkp - mp2*s2*s3*propkp - 8./3.*mp2*s2 - 7./6.*mp2*pow(s2,2) *propkp - 4*mp2*s3 - 3./2.*mp2*pow(s3,2)*propkp - 6*pow(mp2,2) *s1*propkp + 8./3.*pow(mp2,2)*s2*propkp + 4*pow(mp2,2)*s3* propkp - 4./3.*pow(mp2,2) + 4./3.*pow(mp2,3)*propkp + 2./3.* mk2*s1 - 4./3.*mk2*s2 - 10./3.*mk2*s3 + 11./6.*pow(mk2,2) - s1 *s2 - 3./2.*pow(s1,2) + s2*s3 + 7./6.*pow(s2,2) + 3./2.*pow( s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G27 * ( mp2* mk2 + 2./3.*mp2*s1*s2*propkp + 4*mp2*s1 + mp2*pow(s1,2)*propkp - 2./3.*mp2*s2*s3*propkp - 16./9.*mp2*s2 - 7./9.*mp2*pow( s2,2)*propkp - 8./3.*mp2*s3 - mp2*pow(s3,2)*propkp - 4*pow( mp2,2)*s1*propkp + 16./9.*pow(mp2,2)*s2*propkp + 8./3.*pow( mp2,2)*s3*propkp - 8./9.*pow(mp2,2) + 8./9.*pow(mp2,3)*propkp + 4./9.*mk2*s1 - 8./9.*mk2*s2 - 20./9.*mk2*s3 + 11./9.*pow( mk2,2) - 2./3.*s1*s2 - pow(s1,2) + 2./3.*s2*s3 + 7./9.*pow( s2,2) + pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G8 * ( - 5./ 2.*mp2*mk2 + 1./2.*mp2*s1*s2*propkp + 26./3.*mp2*s1 + 3./2.* mp2*pow(s1,2)*propkp - 1./2.*mp2*s2*s3*propkp - 9./2.*mp2*s2 - 7./3.*mp2*pow(s2,2)*propkp - 13./3.*mp2*s3 - 3./2.*mp2*pow( s3,2)*propkp - 6*pow(mp2,2)*s1*propkp + 4*pow(mp2,2)*s2*propkp + 4*pow(mp2,2)*s3*propkp - 13./3.*pow(mp2,2) + 4./3.*pow( mp2,3)*propkp + 13./3.*mk2*s1 - 11./3.*mk2*s2 + 1./3.*mk2*s3 + 5./6.*pow(mk2,2) - 2*s1*s2 - 3./2.*pow(s1,2) + 2*s2*s3 + 17. /6.*pow(s2,2) + 3./2.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G27 * ( - 10. /3.*mp2*mk2 + 1./3.*mp2*s1*s2*propkp + 172./9.*mp2*s1 + mp2* pow(s1,2)*propkp - 1./3.*mp2*s2*s3*propkp - 49./3.*mp2*s2 - 14. /9.*mp2*pow(s2,2)*propkp - 86./9.*mp2*s3 - mp2*pow(s3,2)* propkp - 4*pow(mp2,2)*s1*propkp + 8./3.*pow(mp2,2)*s2*propkp + 8./3.*pow(mp2,2)*s3*propkp - 7./18.*pow(mp2,2) + 8./9.*pow( mp2,3)*propkp + 56./9.*mk2*s1 - 52./9.*mk2*s2 + 32./9.*mk2*s3 - 5./18.*pow(mk2,2) - 4./3.*s1*s2 - 17./2.*pow(s1,2) + 4./3.* s2*s3 + 109./18.*pow(s2,2) + 17./2.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( - 56./ 3.*mp2*s3 - 14./3.*mp2*pow(s3,2)*propkp + 56./3.*pow(mp2,2)*s3 *propkp + 56./3.*pow(mp2,2) - 56./3.*pow(mp2,3)*propkp - 2*mk2 *s3 + 14./3.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( - 112./9.*mp2*s3 - 28./9.*mp2*pow(s3,2)*propkp + 112./9.*pow( mp2,2)*s3*propkp + 112./9.*pow(mp2,2) - 112./9.*pow(mp2,3)* propkp - 4./3.*mk2*s3 + 28./9.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( 6*mp2* mk2 - 12*mp2*s3 - 14./3.*mp2*pow(s3,2)*propkp + 56./3.*pow( mp2,2)*s3*propkp + 15*pow(mp2,2) - 56./3.*pow(mp2,3)*propkp - 2*mk2*s3 + 1./3.*pow(mk2,2) + 5./3.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( 82./3. *mp2*mk2 - 124./3.*mp2*s3 - 28./9.*mp2*pow(s3,2)*propkp + 112./ 9.*pow(mp2,2)*s3*propkp + 35*pow(mp2,2) - 112./9.*pow(mp2,3)* propkp - 44./3.*mk2*s3 + 47./9.*pow(mk2,2) + 85./9.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G8 * ( - 5./ 3.*mp2*mk2 + 3./2.*mp2*s1*s2*propkp - 3./2.*mp2*s1*s3*propkp - 4./3.*mp2*s1 + 6*mp2*s2 + 3./2.*mp2*pow(s2,2)*propkp - 4* mp2*s3 - 3./2.*mp2*pow(s3,2)*propkp + 4./3.*pow(mp2,2)*s1* propkp - 6*pow(mp2,2)*s2*propkp + 4*pow(mp2,2)*s3*propkp - 4./ 3.*pow(mp2,2) + 4./3.*pow(mp2,3)*propkp + mk2*s1 + 5./2.*mk2* s2 - 3./2.*mk2*s3 - pow(mk2,2) - 3./2.*s1*s2 + 3./2.*s1*s3 - 3. /2.*pow(s2,2) + 3./2.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G27 * ( - 10. /9.*mp2*mk2 + mp2*s1*s2*propkp - mp2*s1*s3*propkp - 8./9.*mp2* s1 + 4*mp2*s2 + mp2*pow(s2,2)*propkp - 8./3.*mp2*s3 - mp2*pow( s3,2)*propkp + 8./9.*pow(mp2,2)*s1*propkp - 4*pow(mp2,2)*s2* propkp + 8./3.*pow(mp2,2)*s3*propkp - 8./9.*pow(mp2,2) + 8./9. *pow(mp2,3)*propkp + 2./3.*mk2*s1 + 5./3.*mk2*s2 - mk2*s3 - 2./ 3.*pow(mk2,2) - s1*s2 + s1*s3 - pow(s2,2) + pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G8 * ( - 5./ 3.*mp2*mk2 + 3./2.*mp2*s1*s2*propkp + 6*mp2*s1 + 3./2.*mp2* pow(s1,2)*propkp - 3./2.*mp2*s2*s3*propkp - 4./3.*mp2*s2 - 4* mp2*s3 - 3./2.*mp2*pow(s3,2)*propkp - 6*pow(mp2,2)*s1*propkp + 4./3.*pow(mp2,2)*s2*propkp + 4*pow(mp2,2)*s3*propkp - 4./3. *pow(mp2,2) + 4./3.*pow(mp2,3)*propkp + 5./2.*mk2*s1 + mk2*s2 - 3./2.*mk2*s3 - pow(mk2,2) - 3./2.*s1*s2 - 3./2.*pow(s1,2) + 3./2.*s2*s3 + 3./2.*pow(s3,2) ); K3Pp4Ix10112 += + C11b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G27 * ( - 10. /9.*mp2*mk2 + mp2*s1*s2*propkp + 4*mp2*s1 + mp2*pow(s1,2)* propkp - mp2*s2*s3*propkp - 8./9.*mp2*s2 - 8./3.*mp2*s3 - mp2* pow(s3,2)*propkp - 4*pow(mp2,2)*s1*propkp + 8./9.*pow(mp2,2)* s2*propkp + 8./3.*pow(mp2,2)*s3*propkp - 8./9.*pow(mp2,2) + 8./ 9.*pow(mp2,3)*propkp + 5./3.*mk2*s1 + 2./3.*mk2*s2 - mk2*s3 - 2./3.*pow(mk2,2) - s1*s2 - pow(s1,2) + s2*s3 + pow(s3,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G8 * ( - 19./ 3.*mp2*mk2 - 1./2.*mp2*s1*s2*propkp + 1./2.*mp2*s1*s3*propkp - 31./6.*mp2*s1 - 7./6.*mp2*pow(s1,2)*propkp + 10./3.*mp2*s2 + 10./3.*mp2*s3 + 4./3.*pow(mp2,2)*s1*propkp + 2*pow(mp2,2) + 3*mk2*s1 - 3*mk2*s2 - 3*mk2*s3 + 3*pow(mk2,2) + 4./3.*pow( s1,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G27 * ( - 128./9.*mp2*mk2 - 1./3.*mp2*s1*s2*propkp + 1./3.*mp2*s1*s3* propkp - 76./9.*mp2*s1 - 7./9.*mp2*pow(s1,2)*propkp + 65./9.* mp2*s2 + 65./9.*mp2*s3 + 8./9.*pow(mp2,2)*s1*propkp + 19./3.* pow(mp2,2) + 7*mk2*s1 - 7*mk2*s2 - 7*mk2*s3 + 7*pow(mk2,2) + 8. /9.*pow(s1,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G8 * ( 1./3.* mp2*mk2 - mp2*s1*s2*propkp + mp2*s1*s3*propkp - 2./3.*mp2*s1 - 7./3.*mp2*pow(s1,2)*propkp - 5./3.*mp2*s2 - 11./3.*mp2*s3 + 8./3.*pow(mp2,2)*s1*propkp + 23./3.*pow(mp2,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G27 * ( 2./9. *mp2*mk2 - 2./3.*mp2*s1*s2*propkp + 2./3.*mp2*s1*s3*propkp - 4. /9.*mp2*s1 - 14./9.*mp2*pow(s1,2)*propkp - 10./9.*mp2*s2 - 22./ 9.*mp2*s3 + 16./9.*pow(mp2,2)*s1*propkp + 46./9.*pow(mp2,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G8 * ( - 19./ 3.*mp2*mk2 - 1./2.*mp2*s1*s2*propkp + 10./3.*mp2*s1 + 1./2.* mp2*s2*s3*propkp - 31./6.*mp2*s2 - 7./6.*mp2*pow(s2,2)*propkp + 10./3.*mp2*s3 + 4./3.*pow(mp2,2)*s2*propkp + 2*pow(mp2,2) - 3*mk2*s1 + 3*mk2*s2 - 3*mk2*s3 + 3*pow(mk2,2) + 4./3.*pow( s2,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G27 * ( - 128./9.*mp2*mk2 - 1./3.*mp2*s1*s2*propkp + 65./9.*mp2*s1 + 1./ 3.*mp2*s2*s3*propkp - 76./9.*mp2*s2 - 7./9.*mp2*pow(s2,2)* propkp + 65./9.*mp2*s3 + 8./9.*pow(mp2,2)*s2*propkp + 19./3.* pow(mp2,2) - 7*mk2*s1 + 7*mk2*s2 - 7*mk2*s3 + 7*pow(mk2,2) + 8. /9.*pow(s2,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G8 * ( 1./3.* mp2*mk2 - mp2*s1*s2*propkp - 5./3.*mp2*s1 + mp2*s2*s3*propkp - 2./3.*mp2*s2 - 7./3.*mp2*pow(s2,2)*propkp - 11./3.*mp2*s3 + 8./3.*pow(mp2,2)*s2*propkp + 23./3.*pow(mp2,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G27 * ( 2./9. *mp2*mk2 - 2./3.*mp2*s1*s2*propkp - 10./9.*mp2*s1 + 2./3.*mp2* s2*s3*propkp - 4./9.*mp2*s2 - 14./9.*mp2*pow(s2,2)*propkp - 22. /9.*mp2*s3 + 16./9.*pow(mp2,2)*s2*propkp + 46./9.*pow(mp2,2) ) ; K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( - 4./ 3.*mp2*s3 - 14./3.*mp2*pow(s3,2)*propkp + 28./3.*pow(mp2,2)*s3 *propkp ); K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( - 8./ 9.*mp2*s3 - 28./9.*mp2*pow(s3,2)*propkp + 56./9.*pow(mp2,2)*s3 *propkp ); K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( - 14./ 3.*mp2*mk2 + 4*mp2*s3 - 14./3.*mp2*pow(s3,2)*propkp + 28./3.* pow(mp2,2)*s3*propkp - 14./3.*pow(mp2,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( - 28. /9.*mp2*mk2 + 8./3.*mp2*s3 - 28./9.*mp2*pow(s3,2)*propkp + 56./ 9.*pow(mp2,2)*s3*propkp - 28./9.*pow(mp2,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G8 * ( - 2* mp2*mk2 + 1./2.*mp2*s1*s2*propkp - 1./2.*mp2*s1*s3*propkp + 4./ 3.*mp2*s1 + 7./6.*mp2*pow(s1,2)*propkp - 4./3.*pow(mp2,2)*s1* propkp + 3./2.*mk2*s1 + mk2*s2 + mk2*s3 - 2*pow(mk2,2) - 1./2. *s1*s2 + 1./2.*s1*s3 - 7./6.*pow(s1,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G27 * ( - 4./ 3.*mp2*mk2 + 1./3.*mp2*s1*s2*propkp - 1./3.*mp2*s1*s3*propkp + 8./9.*mp2*s1 + 7./9.*mp2*pow(s1,2)*propkp - 8./9.*pow( mp2,2)*s1*propkp + mk2*s1 + 2./3.*mk2*s2 + 2./3.*mk2*s3 - 4./3. *pow(mk2,2) - 1./3.*s1*s2 + 1./3.*s1*s3 - 7./9.*pow(s1,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G8 * ( - 25./ 6.*mp2*mk2 + 3./2.*mp2*s1*s2*propkp - 3./2.*mp2*s1*s3*propkp + 4*mp2*s1 + 7./2.*mp2*pow(s1,2)*propkp - 11./6.*mp2*s2 - 5./ 6.*mp2*s3 - 4*pow(mp2,2)*s1*propkp + 7./2.*pow(mp2,2) + 7./3.* mk2*s1 + 4./3.*mk2*s2 + 4./3.*mk2*s3 + 2*pow(mk2,2) - s1*s2 + s1*s3 - 11./3.*pow(s1,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G27 * ( - 115./9.*mp2*mk2 + mp2*s1*s2*propkp - mp2*s1*s3*propkp + 13./3. *mp2*s1 + 7./3.*mp2*pow(s1,2)*propkp - 26./9.*mp2*s2 - 110./9. *mp2*s3 - 8./3.*pow(mp2,2)*s1*propkp + 14*pow(mp2,2) - 1./9.* mk2*s1 + 23./9.*mk2*s2 + 113./9.*mk2*s3 - 1./3.*pow(mk2,2) - 2. /3.*s1*s2 + 2./3.*s1*s3 - 22./9.*pow(s1,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G8 * ( - 2* mp2*mk2 + 1./2.*mp2*s1*s2*propkp - 1./2.*mp2*s2*s3*propkp + 4./ 3.*mp2*s2 + 7./6.*mp2*pow(s2,2)*propkp - 4./3.*pow(mp2,2)*s2* propkp + mk2*s1 + 3./2.*mk2*s2 + mk2*s3 - 2*pow(mk2,2) - 1./2. *s1*s2 + 1./2.*s2*s3 - 7./6.*pow(s2,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G27 * ( - 4./ 3.*mp2*mk2 + 1./3.*mp2*s1*s2*propkp - 1./3.*mp2*s2*s3*propkp + 8./9.*mp2*s2 + 7./9.*mp2*pow(s2,2)*propkp - 8./9.*pow( mp2,2)*s2*propkp + 2./3.*mk2*s1 + mk2*s2 + 2./3.*mk2*s3 - 4./3. *pow(mk2,2) - 1./3.*s1*s2 + 1./3.*s2*s3 - 7./9.*pow(s2,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G8 * ( - 25./ 6.*mp2*mk2 + 3./2.*mp2*s1*s2*propkp - 11./6.*mp2*s1 - 3./2.* mp2*s2*s3*propkp + 4*mp2*s2 + 7./2.*mp2*pow(s2,2)*propkp - 5./ 6.*mp2*s3 - 4*pow(mp2,2)*s2*propkp + 7./2.*pow(mp2,2) + 4./3.* mk2*s1 + 7./3.*mk2*s2 + 4./3.*mk2*s3 + 2*pow(mk2,2) - s1*s2 + s2*s3 - 11./3.*pow(s2,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G27 * ( - 115./9.*mp2*mk2 + mp2*s1*s2*propkp - 26./9.*mp2*s1 - mp2*s2*s3 *propkp + 13./3.*mp2*s2 + 7./3.*mp2*pow(s2,2)*propkp - 110./9. *mp2*s3 - 8./3.*pow(mp2,2)*s2*propkp + 14*pow(mp2,2) + 23./9.* mk2*s1 - 1./9.*mk2*s2 + 113./9.*mk2*s3 - 1./3.*pow(mk2,2) - 2./ 3.*s1*s2 + 2./3.*s2*s3 - 22./9.*pow(s2,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( 28./3. *mp2*s3 + 14./3.*mp2*pow(s3,2)*propkp - 28./3.*pow(mp2,2)*s3* propkp + 4*mk2*s3 - 14./3.*pow(s3,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( 56./9. *mp2*s3 + 28./9.*mp2*pow(s3,2)*propkp - 56./9.*pow(mp2,2)*s3* propkp + 8./3.*mk2*s3 - 28./9.*pow(s3,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( 6*mp2* mk2 + 10./3.*mp2*s3 + 14./3.*mp2*pow(s3,2)*propkp - 28./3.* pow(mp2,2)*s3*propkp + 10./3.*pow(mp2,2) - 2./3.*mk2*s3 - 2* pow(s3,2) ); K3Pp4Ix10112 += + C12b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( 32./3. *mp2*mk2 + 170./9.*mp2*s3 + 28./9.*mp2*pow(s3,2)*propkp - 56./ 9.*pow(mp2,2)*s3*propkp - 130./9.*pow(mp2,2) - 154./9.*mk2*s3 + 10*pow(mk2,2) - 4./3.*pow(s3,2) ); K3Pp4Ix10112 += + C12b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G8 * ( - 17./ 6.*mp2*mk2 - 1./2.*mp2*s1*s2*propkp + 1./2.*mp2*s1*s3*propkp - 4./3.*mp2*s1 - 7./6.*mp2*pow(s1,2)*propkp + 4./3.*pow( mp2,2)*s1*propkp - mk2*s1 + 5./6.*mk2*s2 + 11./6.*mk2*s3 - 7./ 6.*pow(mk2,2) + 1./2.*s1*s2 - 1./2.*s1*s3 + 7./6.*pow(s1,2) ); K3Pp4Ix10112 += + C12b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G27 * ( - 17. /9.*mp2*mk2 - 1./3.*mp2*s1*s2*propkp + 1./3.*mp2*s1*s3*propkp - 8./9.*mp2*s1 - 7./9.*mp2*pow(s1,2)*propkp + 8./9.*pow( mp2,2)*s1*propkp - 2./3.*mk2*s1 + 5./9.*mk2*s2 + 11./9.*mk2*s3 - 7./9.*pow(mk2,2) + 1./3.*s1*s2 - 1./3.*s1*s3 + 7./9.*pow( s1,2) ); K3Pp4Ix10112 += + C12b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G8 * ( - 17./ 6.*mp2*mk2 - 1./2.*mp2*s1*s2*propkp + 1./2.*mp2*s2*s3*propkp - 4./3.*mp2*s2 - 7./6.*mp2*pow(s2,2)*propkp + 4./3.*pow( mp2,2)*s2*propkp + 5./6.*mk2*s1 - mk2*s2 + 11./6.*mk2*s3 - 7./ 6.*pow(mk2,2) + 1./2.*s1*s2 - 1./2.*s2*s3 + 7./6.*pow(s2,2) ); K3Pp4Ix10112 += + C12b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G27 * ( - 17. /9.*mp2*mk2 - 1./3.*mp2*s1*s2*propkp + 1./3.*mp2*s2*s3*propkp - 8./9.*mp2*s2 - 7./9.*mp2*pow(s2,2)*propkp + 8./9.*pow( mp2,2)*s2*propkp + 5./9.*mk2*s1 - 2./3.*mk2*s2 + 11./9.*mk2*s3 - 7./9.*pow(mk2,2) + 1./3.*s1*s2 - 1./3.*s2*s3 + 7./9.*pow( s2,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G8 * ( - 37./ 6.*mp2*mk2 - 19./6.*mp2*s1 - 1./6.*mp2*pow(s1,2)*propkp + 19./ 6.*mp2*s2 + 1./2.*mp2*pow(s2,2)*propkp + 19./6.*mp2*s3 - 1./2. *mp2*pow(s3,2)*propkp + pow(mp2,2)*s1*propkp - 5./3.*pow( mp2,2)*s2*propkp + 1./3.*pow(mp2,2)*s3*propkp + 17./6.*pow( mp2,2) + 1./2.*mk2*s1 - 1./2.*mk2*s2 - 3./2.*mk2*s3 + pow( mk2,2) + 5./6.*pow(s1,2) - 1./2.*pow(s2,2) + 1./2.*pow(s3,2) ) ; K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G27 * ( - 97. /9.*mp2*mk2 - 23./18.*mp2*s1 - 1./9.*mp2*pow(s1,2)*propkp + 23. /18.*mp2*s2 + 1./3.*mp2*pow(s2,2)*propkp - 67./18.*mp2*s3 - 1./ 3.*mp2*pow(s3,2)*propkp + 2./3.*pow(mp2,2)*s1*propkp - 10./9.* pow(mp2,2)*s2*propkp + 2./9.*pow(mp2,2)*s3*propkp + 77./9.* pow(mp2,2) + 17./6.*mk2*s1 - 17./6.*mk2*s2 + 3./2.*mk2*s3 + 2./ 3.*pow(mk2,2) - 5./18.*pow(s1,2) + 13./6.*pow(s2,2) - 13./6.* pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G8 * ( - 1./ 3.*mp2*s1 - 1./3.*mp2*pow(s1,2)*propkp + 1./3.*mp2*s2 + mp2* pow(s2,2)*propkp + 1./3.*mp2*s3 - mp2*pow(s3,2)*propkp + 2* pow(mp2,2)*s1*propkp - 10./3.*pow(mp2,2)*s2*propkp + 2./3.* pow(mp2,2)*s3*propkp + 2./3.*pow(mp2,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G27 * ( - 2./ 9.*mp2*s1 - 2./9.*mp2*pow(s1,2)*propkp + 2./9.*mp2*s2 + 2./3.* mp2*pow(s2,2)*propkp + 2./9.*mp2*s3 - 2./3.*mp2*pow(s3,2)* propkp + 4./3.*pow(mp2,2)*s1*propkp - 20./9.*pow(mp2,2)*s2* propkp + 4./9.*pow(mp2,2)*s3*propkp + 4./9.*pow(mp2,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G8 * ( - 37./ 6.*mp2*mk2 + 19./6.*mp2*s1 + 1./2.*mp2*pow(s1,2)*propkp - 19./ 6.*mp2*s2 - 1./6.*mp2*pow(s2,2)*propkp + 19./6.*mp2*s3 - 1./2. *mp2*pow(s3,2)*propkp - 5./3.*pow(mp2,2)*s1*propkp + pow( mp2,2)*s2*propkp + 1./3.*pow(mp2,2)*s3*propkp + 17./6.*pow( mp2,2) - 1./2.*mk2*s1 + 1./2.*mk2*s2 - 3./2.*mk2*s3 + pow( mk2,2) - 1./2.*pow(s1,2) + 5./6.*pow(s2,2) + 1./2.*pow(s3,2) ) ; K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G27 * ( - 97. /9.*mp2*mk2 + 23./18.*mp2*s1 + 1./3.*mp2*pow(s1,2)*propkp - 23. /18.*mp2*s2 - 1./9.*mp2*pow(s2,2)*propkp - 67./18.*mp2*s3 - 1./ 3.*mp2*pow(s3,2)*propkp - 10./9.*pow(mp2,2)*s1*propkp + 2./3.* pow(mp2,2)*s2*propkp + 2./9.*pow(mp2,2)*s3*propkp + 77./9.* pow(mp2,2) - 17./6.*mk2*s1 + 17./6.*mk2*s2 + 3./2.*mk2*s3 + 2./ 3.*pow(mk2,2) + 13./6.*pow(s1,2) - 5./18.*pow(s2,2) - 13./6.* pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G8 * ( 1./3.* mp2*s1 + mp2*pow(s1,2)*propkp - 1./3.*mp2*s2 - 1./3.*mp2*pow( s2,2)*propkp + 1./3.*mp2*s3 - mp2*pow(s3,2)*propkp - 10./3.* pow(mp2,2)*s1*propkp + 2*pow(mp2,2)*s2*propkp + 2./3.*pow( mp2,2)*s3*propkp + 2./3.*pow(mp2,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G27 * ( 2./9. *mp2*s1 + 2./3.*mp2*pow(s1,2)*propkp - 2./9.*mp2*s2 - 2./9.* mp2*pow(s2,2)*propkp + 2./9.*mp2*s3 - 2./3.*mp2*pow(s3,2)* propkp - 20./9.*pow(mp2,2)*s1*propkp + 4./3.*pow(mp2,2)*s2* propkp + 4./9.*pow(mp2,2)*s3*propkp + 4./9.*pow(mp2,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( - 8./ 3.*mp2*s3 - 2./3.*mp2*pow(s3,2)*propkp + 20./3.*pow(mp2,2)*s3* propkp + 16./3.*pow(mp2,2) - 32./3.*pow(mp2,3)*propkp + 1./3.* pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( 104./ 9.*mp2*s3 - 4./9.*mp2*pow(s3,2)*propkp + 40./9.*pow(mp2,2)*s3* propkp - 208./9.*pow(mp2,2) - 64./9.*pow(mp2,3)*propkp - 13./9. *pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( 2./3.* mp2*mk2 - 4./3.*mp2*s3 - 2./3.*mp2*pow(s3,2)*propkp + 20./3.* pow(mp2,2)*s3*propkp + 6*pow(mp2,2) - 32./3.*pow(mp2,3)*propkp ); K3Pp4Ix10112 += + C21b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( 4./9. *mp2*mk2 - 8./9.*mp2*s3 - 4./9.*mp2*pow(s3,2)*propkp + 40./9.* pow(mp2,2)*s3*propkp + 4*pow(mp2,2) - 64./9.*pow(mp2,3)*propkp ); K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G8 * ( 1./6.* mp2*mk2 + mp2*s1 + 1./6.*mp2*pow(s1,2)*propkp - 5./3.*mp2*s2 - 1./2.*mp2*pow(s2,2)*propkp + 1./3.*mp2*s3 + 1./2.*mp2*pow( s3,2)*propkp - pow(mp2,2)*s1*propkp + 5./3.*pow(mp2,2)*s2* propkp - 1./3.*pow(mp2,2)*s3*propkp + mk2*s3 - 1./2.*pow( mk2,2) - 1./6.*pow(s1,2) + 1./2.*pow(s2,2) - 1./2.*pow(s3,2) ) ; K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G27 * ( 1./9. *mp2*mk2 + 2./3.*mp2*s1 + 1./9.*mp2*pow(s1,2)*propkp - 10./9.* mp2*s2 - 1./3.*mp2*pow(s2,2)*propkp + 2./9.*mp2*s3 + 1./3.*mp2 *pow(s3,2)*propkp - 2./3.*pow(mp2,2)*s1*propkp + 10./9.*pow( mp2,2)*s2*propkp - 2./9.*pow(mp2,2)*s3*propkp + 2./3.*mk2*s3 - 1./3.*pow(mk2,2) - 1./9.*pow(s1,2) + 1./3.*pow(s2,2) - 1./3. *pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G8 * ( - 1./ 6.*mp2*mk2 + 1./2.*mp2*s1*s2*propkp - 1./2.*mp2*s1*s3*propkp + 3*mp2*s1 + 1./3.*mp2*pow(s1,2)*propkp - 7./2.*mp2*s2 - 1./2. *mp2*pow(s2,2)*propkp - 1./2.*mp2*s3 + 1./2.*mp2*pow(s3,2)* propkp - 5./3.*pow(mp2,2)*s1*propkp + 5./3.*pow(mp2,2)*s2* propkp - 1./3.*pow(mp2,2)*s3*propkp + 11./6.*pow(mp2,2) + 5./6. *mk2*s1 - 5./6.*mk2*s2 - 5./6.*mk2*s3 - 5./6.*pow(s1,2) + 1./2. *pow(s2,2) - 1./2.*pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G27 * ( - 1./ 9.*mp2*mk2 + 1./3.*mp2*s1*s2*propkp - 1./3.*mp2*s1*s3*propkp + 9./2.*mp2*s1 + 2./9.*mp2*pow(s1,2)*propkp - 29./6.*mp2*s2 - 1./3.*mp2*pow(s2,2)*propkp - 17./6.*mp2*s3 + 1./3.*mp2*pow( s3,2)*propkp - 10./9.*pow(mp2,2)*s1*propkp + 10./9.*pow(mp2,2) *s2*propkp - 2./9.*pow(mp2,2)*s3*propkp + 11./9.*pow(mp2,2) + 25./18.*mk2*s1 - 25./18.*mk2*s2 - 25./18.*mk2*s3 - 25./18.* pow(s1,2) + 17./6.*pow(s2,2) - 17./6.*pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G8 * ( 1./6.* mp2*mk2 - 5./3.*mp2*s1 - 1./2.*mp2*pow(s1,2)*propkp + mp2*s2 + 1./6.*mp2*pow(s2,2)*propkp + 1./3.*mp2*s3 + 1./2.*mp2*pow( s3,2)*propkp + 5./3.*pow(mp2,2)*s1*propkp - pow(mp2,2)*s2* propkp - 1./3.*pow(mp2,2)*s3*propkp + mk2*s3 - 1./2.*pow( mk2,2) + 1./2.*pow(s1,2) - 1./6.*pow(s2,2) - 1./2.*pow(s3,2) ) ; K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G27 * ( 1./9. *mp2*mk2 - 10./9.*mp2*s1 - 1./3.*mp2*pow(s1,2)*propkp + 2./3.* mp2*s2 + 1./9.*mp2*pow(s2,2)*propkp + 2./9.*mp2*s3 + 1./3.*mp2 *pow(s3,2)*propkp + 10./9.*pow(mp2,2)*s1*propkp - 2./3.*pow( mp2,2)*s2*propkp - 2./9.*pow(mp2,2)*s3*propkp + 2./3.*mk2*s3 - 1./3.*pow(mk2,2) + 1./3.*pow(s1,2) - 1./9.*pow(s2,2) - 1./3. *pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G8 * ( - 1./ 6.*mp2*mk2 + 1./2.*mp2*s1*s2*propkp - 7./2.*mp2*s1 - 1./2.*mp2 *pow(s1,2)*propkp - 1./2.*mp2*s2*s3*propkp + 3*mp2*s2 + 1./3.* mp2*pow(s2,2)*propkp - 1./2.*mp2*s3 + 1./2.*mp2*pow(s3,2)* propkp + 5./3.*pow(mp2,2)*s1*propkp - 5./3.*pow(mp2,2)*s2* propkp - 1./3.*pow(mp2,2)*s3*propkp + 11./6.*pow(mp2,2) - 5./6. *mk2*s1 + 5./6.*mk2*s2 - 5./6.*mk2*s3 + 1./2.*pow(s1,2) - 5./6. *pow(s2,2) - 1./2.*pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G27 * ( - 1./ 9.*mp2*mk2 + 1./3.*mp2*s1*s2*propkp - 29./6.*mp2*s1 - 1./3.* mp2*pow(s1,2)*propkp - 1./3.*mp2*s2*s3*propkp + 9./2.*mp2*s2 + 2./9.*mp2*pow(s2,2)*propkp - 17./6.*mp2*s3 + 1./3.*mp2*pow( s3,2)*propkp + 10./9.*pow(mp2,2)*s1*propkp - 10./9.*pow(mp2,2) *s2*propkp - 2./9.*pow(mp2,2)*s3*propkp + 11./9.*pow(mp2,2) - 25./18.*mk2*s1 + 25./18.*mk2*s2 - 25./18.*mk2*s3 + 17./6.*pow( s1,2) - 25./18.*pow(s2,2) - 17./6.*pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( 20./3. *mp2*s3 + 2./3.*mp2*pow(s3,2)*propkp - 20./3.*pow(mp2,2)*s3* propkp - 32./3.*pow(mp2,2) + 32./3.*pow(mp2,3)*propkp - 2./3.* pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( 40./9. *mp2*s3 + 4./9.*mp2*pow(s3,2)*propkp - 40./9.*pow(mp2,2)*s3* propkp - 64./9.*pow(mp2,2) + 64./9.*pow(mp2,3)*propkp - 4./9.* pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( - 8./ 3.*mp2*mk2 + 10./3.*mp2*s3 + 2./3.*mp2*pow(s3,2)*propkp - 20./ 3.*pow(mp2,2)*s3*propkp - 29./3.*pow(mp2,2) + 32./3.*pow( mp2,3)*propkp - 2./3.*mk2*s3 + 1./3.*pow(mk2,2) + 1./3.*pow( s3,2) ); K3Pp4Ix10112 += + C21b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( - 106./9.*mp2*mk2 + 110./9.*mp2*s3 + 4./9.*mp2*pow(s3,2)*propkp - 40./9.*pow(mp2,2)*s3*propkp - 193./9.*pow(mp2,2) + 64./9.* pow(mp2,3)*propkp + 26./9.*mk2*s3 - 13./9.*pow(mk2,2) - 13./9. *pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G8 * ( 1./6.* mp2*mk2 - 1./2.*mp2*s1*s2*propkp + 1./2.*mp2*s1*s3*propkp + 1./ 3.*mp2*s1 - 5./3.*mp2*s2 - 1./2.*mp2*pow(s2,2)*propkp + 1./3.* mp2*s3 + 1./2.*mp2*pow(s3,2)*propkp - 1./3.*pow(mp2,2)*s1* propkp + 5./3.*pow(mp2,2)*s2*propkp - 1./3.*pow(mp2,2)*s3* propkp - 1./6.*mk2*s1 - 2./3.*mk2*s2 + 1./3.*mk2*s3 + 1./6.* pow(mk2,2) + 1./2.*s1*s2 - 1./2.*s1*s3 + 1./2.*pow(s2,2) - 1./ 2.*pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G27 * ( 1./9. *mp2*mk2 - 1./3.*mp2*s1*s2*propkp + 1./3.*mp2*s1*s3*propkp + 2. /9.*mp2*s1 - 10./9.*mp2*s2 - 1./3.*mp2*pow(s2,2)*propkp + 2./9. *mp2*s3 + 1./3.*mp2*pow(s3,2)*propkp - 2./9.*pow(mp2,2)*s1* propkp + 10./9.*pow(mp2,2)*s2*propkp - 2./9.*pow(mp2,2)*s3* propkp - 1./9.*mk2*s1 - 4./9.*mk2*s2 + 2./9.*mk2*s3 + 1./9.* pow(mk2,2) + 1./3.*s1*s2 - 1./3.*s1*s3 + 1./3.*pow(s2,2) - 1./ 3.*pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G8 * ( 1./6.* mp2*mk2 - 1./2.*mp2*s1*s2*propkp - 5./3.*mp2*s1 - 1./2.*mp2* pow(s1,2)*propkp + 1./2.*mp2*s2*s3*propkp + 1./3.*mp2*s2 + 1./ 3.*mp2*s3 + 1./2.*mp2*pow(s3,2)*propkp + 5./3.*pow(mp2,2)*s1* propkp - 1./3.*pow(mp2,2)*s2*propkp - 1./3.*pow(mp2,2)*s3* propkp - 2./3.*mk2*s1 - 1./6.*mk2*s2 + 1./3.*mk2*s3 + 1./6.* pow(mk2,2) + 1./2.*s1*s2 + 1./2.*pow(s1,2) - 1./2.*s2*s3 - 1./ 2.*pow(s3,2) ); K3Pp4Ix10112 += + C21b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G27 * ( 1./9. *mp2*mk2 - 1./3.*mp2*s1*s2*propkp - 10./9.*mp2*s1 - 1./3.*mp2* pow(s1,2)*propkp + 1./3.*mp2*s2*s3*propkp + 2./9.*mp2*s2 + 2./ 9.*mp2*s3 + 1./3.*mp2*pow(s3,2)*propkp + 10./9.*pow(mp2,2)*s1* propkp - 2./9.*pow(mp2,2)*s2*propkp - 2./9.*pow(mp2,2)*s3* propkp - 4./9.*mk2*s1 - 1./9.*mk2*s2 + 2./9.*mk2*s3 + 1./9.* pow(mk2,2) + 1./3.*s1*s2 + 1./3.*pow(s1,2) - 1./3.*s2*s3 - 1./ 3.*pow(s3,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G8 * ( mp2*s1 + 1./3.*mp2*pow(s1,2)*propkp ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G27 * ( 2./3. *mp2*s1 + 2./9.*mp2*pow(s1,2)*propkp ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G8 * ( 4./3.* mp2*mk2 + 2./3.*mp2*s1 + 2./3.*mp2*pow(s1,2)*propkp - 4./3.* pow(mp2,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G27 * ( 8./9. *mp2*mk2 + 4./9.*mp2*s1 + 4./9.*mp2*pow(s1,2)*propkp - 8./9.* pow(mp2,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G8 * ( mp2*s2 + 1./3.*mp2*pow(s2,2)*propkp ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G27 * ( 2./3. *mp2*s2 + 2./9.*mp2*pow(s2,2)*propkp ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G8 * ( 4./3.* mp2*mk2 + 2./3.*mp2*s2 + 2./3.*mp2*pow(s2,2)*propkp - 4./3.* pow(mp2,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G27 * ( 8./9. *mp2*mk2 + 4./9.*mp2*s2 + 4./9.*mp2*pow(s2,2)*propkp - 8./9.* pow(mp2,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( 4./3.* mp2*pow(s3,2)*propkp ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( 8./9. *mp2*pow(s3,2)*propkp ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( 8./3.* mp2*mk2 + 4./3.*mp2*s3 + 4./3.*mp2*pow(s3,2)*propkp - 8./3.* pow(mp2,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( 16./9. *mp2*mk2 + 8./9.*mp2*s3 + 8./9.*mp2*pow(s3,2)*propkp - 16./9.* pow(mp2,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G8 * ( - 1./ 3.*mp2*pow(s1,2)*propkp - mk2*s1 + 1./3.*pow(s1,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G27 * ( - 2./ 9.*mp2*pow(s1,2)*propkp - 2./3.*mk2*s1 + 2./9.*pow(s1,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G8 * ( 2*mp2* mk2 + 1./3.*mp2*s1 - 1./3.*mp2*pow(s1,2)*propkp - 2./3.*mp2*s2 - 2./3.*mp2*s3 + 2./3.*pow(mp2,2) - 4./3.*mk2*s1 + 2./3.*mk2* s2 + 2./3.*mk2*s3 - 8./3.*pow(mk2,2) + 2./3.*pow(s1,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G27 * ( 8*mp2 *mk2 + 2./9.*mp2*s1 - 2./9.*mp2*pow(s1,2)*propkp - 4./9.*mp2* s2 - 4./9.*mp2*s3 - 26./9.*pow(mp2,2) - 8./9.*mk2*s1 + 4./9.* mk2*s2 + 4./9.*mk2*s3 - 46./9.*pow(mk2,2) + 4./9.*pow(s1,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G8 * ( - 1./ 3.*mp2*pow(s2,2)*propkp - mk2*s2 + 1./3.*pow(s2,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G27 * ( - 2./ 9.*mp2*pow(s2,2)*propkp - 2./3.*mk2*s2 + 2./9.*pow(s2,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G8 * ( 2*mp2* mk2 - 2./3.*mp2*s1 + 1./3.*mp2*s2 - 1./3.*mp2*pow(s2,2)*propkp - 2./3.*mp2*s3 + 2./3.*pow(mp2,2) + 2./3.*mk2*s1 - 4./3.*mk2* s2 + 2./3.*mk2*s3 - 8./3.*pow(mk2,2) + 2./3.*pow(s2,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G27 * ( 8*mp2 *mk2 - 4./9.*mp2*s1 + 2./9.*mp2*s2 - 2./9.*mp2*pow(s2,2)* propkp - 4./9.*mp2*s3 - 26./9.*pow(mp2,2) + 4./9.*mk2*s1 - 8./ 9.*mk2*s2 + 4./9.*mk2*s3 - 46./9.*pow(mk2,2) + 4./9.*pow(s2,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( - 4./ 3.*mp2*pow(s3,2)*propkp + 4./3.*pow(s3,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( - 8./ 9.*mp2*pow(s3,2)*propkp + 8./9.*pow(s3,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( - 8./ 3.*mp2*s3 - 4./3.*mp2*pow(s3,2)*propkp + 4./3.*pow(mp2,2) + 4./ 3.*mk2*s3 - 4./3.*pow(mk2,2) + 4./3.*pow(s3,2) ); K3Pp4Ix10112 += + C22b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( 40./3. *mp2*mk2 - 16./9.*mp2*s3 - 8./9.*mp2*pow(s3,2)*propkp - 52./9. *pow(mp2,2) + 8./9.*mk2*s3 - 68./9.*pow(mk2,2) + 8./9.*pow( s3,2) ); K3Pp4Ix10112 += + C22b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G8 * ( 2./3.* mp2*mk2 - 1./3.*mp2*pow(s1,2)*propkp + 1./3.*mk2*s1 - 2./3.* pow(mk2,2) + 1./3.*pow(s1,2) ); K3Pp4Ix10112 += + C22b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G27 * ( 4./9. *mp2*mk2 - 2./9.*mp2*pow(s1,2)*propkp + 2./9.*mk2*s1 - 4./9.* pow(mk2,2) + 2./9.*pow(s1,2) ); K3Pp4Ix10112 += + C22b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G8 * ( 2./3.* mp2*mk2 - 1./3.*mp2*pow(s2,2)*propkp + 1./3.*mk2*s2 - 2./3.* pow(mk2,2) + 1./3.*pow(s2,2) ); K3Pp4Ix10112 += + C22b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G27 * ( 4./9. *mp2*mk2 - 2./9.*mp2*pow(s2,2)*propkp + 2./9.*mk2*s2 - 4./9.* pow(mk2,2) + 2./9.*pow(s2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G8 * ( 59./6. *mp2*mk2 + 1./2.*mp2*s1*s2*propkp - 1./2.*mp2*s1*s3*propkp + 11./6.*mp2*s1 - 1./6.*mp2*pow(s1,2)*propkp - 7./3.*mp2*s2 - 7./ 3.*mp2*s3 - 2./3.*pow(mp2,2)*s1*propkp - 3*pow(mp2,2) - 3./2.* mk2*s1 + 3./2.*mk2*s2 + 3./2.*mk2*s3 - 3./2.*pow(mk2,2) - 2./3. *pow(s1,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G27 * ( 194./ 9.*mp2*mk2 + 1./3.*mp2*s1*s2*propkp - 1./3.*mp2*s1*s3*propkp + 67./18.*mp2*s1 - 1./9.*mp2*pow(s1,2)*propkp - 73./18.*mp2* s2 - 73./18.*mp2*s3 - 4./9.*pow(mp2,2)*s1*propkp - 29./2.*pow( mp2,2) - 7./2.*mk2*s1 + 7./2.*mk2*s2 + 7./2.*mk2*s3 - 7./2.* pow(mk2,2) - 4./9.*pow(s1,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G8 * ( mp2*s1 *s2*propkp - mp2*s1*s3*propkp - 1./3.*mp2*s1 - 1./3.*mp2*pow( s1,2)*propkp + 2./3.*mp2*s2 + 8./3.*mp2*s3 - 4./3.*pow(mp2,2)* s1*propkp - 4./3.*pow(mp2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G27 * ( 2./3. *mp2*s1*s2*propkp - 2./3.*mp2*s1*s3*propkp - 2./9.*mp2*s1 - 2./ 9.*mp2*pow(s1,2)*propkp + 4./9.*mp2*s2 + 16./9.*mp2*s3 - 8./9. *pow(mp2,2)*s1*propkp - 8./9.*pow(mp2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G8 * ( 59./6. *mp2*mk2 + 1./2.*mp2*s1*s2*propkp - 7./3.*mp2*s1 - 1./2.*mp2* s2*s3*propkp + 11./6.*mp2*s2 - 1./6.*mp2*pow(s2,2)*propkp - 7./ 3.*mp2*s3 - 2./3.*pow(mp2,2)*s2*propkp - 3*pow(mp2,2) + 3./2.* mk2*s1 - 3./2.*mk2*s2 + 3./2.*mk2*s3 - 3./2.*pow(mk2,2) - 2./3. *pow(s2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G27 * ( 194./ 9.*mp2*mk2 + 1./3.*mp2*s1*s2*propkp - 73./18.*mp2*s1 - 1./3.* mp2*s2*s3*propkp + 67./18.*mp2*s2 - 1./9.*mp2*pow(s2,2)*propkp - 73./18.*mp2*s3 - 4./9.*pow(mp2,2)*s2*propkp - 29./2.*pow( mp2,2) + 7./2.*mk2*s1 - 7./2.*mk2*s2 + 7./2.*mk2*s3 - 7./2.* pow(mk2,2) - 4./9.*pow(s2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G8 * ( mp2*s1 *s2*propkp + 2./3.*mp2*s1 - mp2*s2*s3*propkp - 1./3.*mp2*s2 - 1./3.*mp2*pow(s2,2)*propkp + 8./3.*mp2*s3 - 4./3.*pow(mp2,2)* s2*propkp - 4./3.*pow(mp2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G27 * ( 2./3. *mp2*s1*s2*propkp + 4./9.*mp2*s1 - 2./3.*mp2*s2*s3*propkp - 2./ 9.*mp2*s2 - 2./9.*mp2*pow(s2,2)*propkp + 16./9.*mp2*s3 - 8./9. *pow(mp2,2)*s2*propkp - 8./9.*pow(mp2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( - 2./ 3.*mp2*pow(s3,2)*propkp - 8./3.*pow(mp2,2)*s3*propkp ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( - 4./ 9.*mp2*pow(s3,2)*propkp - 16./9.*pow(mp2,2)*s3*propkp ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( 10./3. *mp2*mk2 - 4*mp2*s3 - 2./3.*mp2*pow(s3,2)*propkp - 8./3.*pow( mp2,2)*s3*propkp + 22./3.*pow(mp2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( 20./9. *mp2*mk2 - 8./3.*mp2*s3 - 4./9.*mp2*pow(s3,2)*propkp - 16./9.* pow(mp2,2)*s3*propkp + 44./9.*pow(mp2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G8 * ( - 1./ 2.*mp2*mk2 - 1./2.*mp2*s1*s2*propkp + 1./2.*mp2*s1*s3*propkp - 2./3.*mp2*s1 + 1./6.*mp2*pow(s1,2)*propkp + 2./3.*pow( mp2,2)*s1*propkp + mk2*s1 - 1./2.*mk2*s2 - 1./2.*mk2*s3 + 1./2. *pow(mk2,2) + 1./2.*s1*s2 - 1./2.*s1*s3 - 1./6.*pow(s1,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G27 * ( - 1./ 3.*mp2*mk2 - 1./3.*mp2*s1*s2*propkp + 1./3.*mp2*s1*s3*propkp - 4./9.*mp2*s1 + 1./9.*mp2*pow(s1,2)*propkp + 4./9.*pow( mp2,2)*s1*propkp + 2./3.*mk2*s1 - 1./3.*mk2*s2 - 1./3.*mk2*s3 + 1./3.*pow(mk2,2) + 1./3.*s1*s2 - 1./3.*s1*s3 - 1./9.*pow( s1,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G8 * ( - 2./ 3.*mp2*mk2 - 3./2.*mp2*s1*s2*propkp + 3./2.*mp2*s1*s3*propkp - 11./3.*mp2*s1 - 1./6.*mp2*pow(s1,2)*propkp + 5./2.*mp2*s2 + 3./2.*mp2*s3 + 2*pow(mp2,2)*s1*propkp - 5./2.*pow(mp2,2) + 7./6.*mk2*s1 - 3./2.*mk2*s2 - 3./2.*mk2*s3 - 3./2.*pow(mk2,2) + s1*s2 - s1*s3 + 1./3.*pow(s1,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G27 * ( 11./9. *mp2*mk2 - mp2*s1*s2*propkp + mp2*s1*s3*propkp - 29./18.*mp2* s1 - 1./9.*mp2*pow(s1,2)*propkp + 5./6.*mp2*s2 + 61./6.*mp2*s3 + 4./3.*pow(mp2,2)*s1*propkp - 25./6.*pow(mp2,2) - 1./18.*mk2 *s1 - 1./6.*mk2*s2 - 61./6.*mk2*s3 - 1./6.*pow(mk2,2) + 2./3.* s1*s2 - 2./3.*s1*s3 + 2./9.*pow(s1,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G8 * ( - 1./ 2.*mp2*mk2 - 1./2.*mp2*s1*s2*propkp + 1./2.*mp2*s2*s3*propkp - 2./3.*mp2*s2 + 1./6.*mp2*pow(s2,2)*propkp + 2./3.*pow( mp2,2)*s2*propkp - 1./2.*mk2*s1 + mk2*s2 - 1./2.*mk2*s3 + 1./2. *pow(mk2,2) + 1./2.*s1*s2 - 1./2.*s2*s3 - 1./6.*pow(s2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G27 * ( - 1./ 3.*mp2*mk2 - 1./3.*mp2*s1*s2*propkp + 1./3.*mp2*s2*s3*propkp - 4./9.*mp2*s2 + 1./9.*mp2*pow(s2,2)*propkp + 4./9.*pow( mp2,2)*s2*propkp - 1./3.*mk2*s1 + 2./3.*mk2*s2 - 1./3.*mk2*s3 + 1./3.*pow(mk2,2) + 1./3.*s1*s2 - 1./3.*s2*s3 - 1./9.*pow( s2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G8 * ( - 2./ 3.*mp2*mk2 - 3./2.*mp2*s1*s2*propkp + 5./2.*mp2*s1 + 3./2.*mp2 *s2*s3*propkp - 11./3.*mp2*s2 - 1./6.*mp2*pow(s2,2)*propkp + 3. /2.*mp2*s3 + 2*pow(mp2,2)*s2*propkp - 5./2.*pow(mp2,2) - 3./2. *mk2*s1 + 7./6.*mk2*s2 - 3./2.*mk2*s3 - 3./2.*pow(mk2,2) + s1* s2 - s2*s3 + 1./3.*pow(s2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G27 * ( 11./9. *mp2*mk2 - mp2*s1*s2*propkp + 5./6.*mp2*s1 + mp2*s2*s3*propkp - 29./18.*mp2*s2 - 1./9.*mp2*pow(s2,2)*propkp + 61./6.*mp2*s3 + 4./3.*pow(mp2,2)*s2*propkp - 25./6.*pow(mp2,2) - 1./6.*mk2* s1 - 1./18.*mk2*s2 - 61./6.*mk2*s3 - 1./6.*pow(mk2,2) + 2./3.* s1*s2 - 2./3.*s2*s3 + 2./9.*pow(s2,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( - 8./ 3.*mp2*s3 + 2./3.*mp2*pow(s3,2)*propkp + 8./3.*pow(mp2,2)*s3* propkp - 2./3.*pow(s3,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( - 16. /9.*mp2*s3 + 4./9.*mp2*pow(s3,2)*propkp + 16./9.*pow(mp2,2)*s3 *propkp - 4./9.*pow(s3,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( - 22./ 3.*mp2*mk2 + 4*mp2*s3 + 2./3.*mp2*pow(s3,2)*propkp + 8./3.* pow(mp2,2)*s3*propkp - 8./3.*pow(mp2,2) + 8./3.*mk2*s3 - 2./3. *pow(mk2,2) - 2*pow(s3,2) ); K3Pp4Ix10112 += + C23b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( - 164./9.*mp2*mk2 - 4*mp2*s3 + 4./9.*mp2*pow(s3,2)*propkp + 16./ 9.*pow(mp2,2)*s3*propkp + 164./9.*pow(mp2,2) + 76./9.*mk2*s3 - 64./9.*pow(mk2,2) - 4./3.*pow(s3,2) ); K3Pp4Ix10112 += + C23b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G8 * ( 1./3.* mp2*mk2 + 1./2.*mp2*s1*s2*propkp - 1./2.*mp2*s1*s3*propkp + 2./ 3.*mp2*s1 + 1./2.*mp2*pow(s1,2)*propkp - 2./3.*pow(mp2,2)*s1* propkp + 1./6.*mk2*s1 - 1./3.*mk2*s2 - 4./3.*mk2*s3 + 1./3.* pow(mk2,2) - 1./2.*s1*s2 + 1./2.*s1*s3 - 1./2.*pow(s1,2) ); K3Pp4Ix10112 += + C23b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G27 * ( 2./9. *mp2*mk2 + 1./3.*mp2*s1*s2*propkp - 1./3.*mp2*s1*s3*propkp + 4. /9.*mp2*s1 + 1./3.*mp2*pow(s1,2)*propkp - 4./9.*pow(mp2,2)*s1* propkp + 1./9.*mk2*s1 - 2./9.*mk2*s2 - 8./9.*mk2*s3 + 2./9.* pow(mk2,2) - 1./3.*s1*s2 + 1./3.*s1*s3 - 1./3.*pow(s1,2) ); K3Pp4Ix10112 += + C23b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G8 * ( 1./3.* mp2*mk2 + 1./2.*mp2*s1*s2*propkp - 1./2.*mp2*s2*s3*propkp + 2./ 3.*mp2*s2 + 1./2.*mp2*pow(s2,2)*propkp - 2./3.*pow(mp2,2)*s2* propkp - 1./3.*mk2*s1 + 1./6.*mk2*s2 - 4./3.*mk2*s3 + 1./3.* pow(mk2,2) - 1./2.*s1*s2 + 1./2.*s2*s3 - 1./2.*pow(s2,2) ); K3Pp4Ix10112 += + C23b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G27 * ( 2./9. *mp2*mk2 + 1./3.*mp2*s1*s2*propkp - 1./3.*mp2*s2*s3*propkp + 4. /9.*mp2*s2 + 1./3.*mp2*pow(s2,2)*propkp - 4./9.*pow(mp2,2)*s2* propkp - 2./9.*mk2*s1 + 1./9.*mk2*s2 - 8./9.*mk2*s3 + 2./9.* pow(mk2,2) - 1./3.*s1*s2 + 1./3.*s2*s3 - 1./3.*pow(s2,2) ); K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G8 * ( - 1./ 3.*mp2*s1*propkp - 1./3.*mp2*s2*propkp - 7./3.*mp2*s3*propkp + 5./3.*mp2 + 8./3.*pow(mp2,2)*propkp - 2*mk2 + 4*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s1,mp2,mu2)*G27 * ( - 2./ 9.*mp2*s1*propkp - 2./9.*mp2*s2*propkp - 14./9.*mp2*s3*propkp + 70./9.*mp2 + 16./9.*pow(mp2,2)*propkp + 16./3.*mk2 - 10./3. *s1 + 10./3.*s2 - 14*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G8 * ( - 2./ 3.*mp2*s1*propkp - 2./3.*mp2*s2*propkp - 14./3.*mp2*s3*propkp + 2./3.*mp2 + 16./3.*pow(mp2,2)*propkp ); K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s1,mk2,mu2)*G27 * ( - 4./ 9.*mp2*s1*propkp - 4./9.*mp2*s2*propkp - 28./9.*mp2*s3*propkp + 4./9.*mp2 + 32./9.*pow(mp2,2)*propkp ); K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G8 * ( - 1./ 3.*mp2*s1*propkp - 1./3.*mp2*s2*propkp - 7./3.*mp2*s3*propkp + 5./3.*mp2 + 8./3.*pow(mp2,2)*propkp - 2*mk2 + 4*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s2,mp2,mu2)*G27 * ( - 2./ 9.*mp2*s1*propkp - 2./9.*mp2*s2*propkp - 14./9.*mp2*s3*propkp + 70./9.*mp2 + 16./9.*pow(mp2,2)*propkp + 16./3.*mk2 + 10./3. *s1 - 10./3.*s2 - 14*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G8 * ( - 2./ 3.*mp2*s1*propkp - 2./3.*mp2*s2*propkp - 14./3.*mp2*s3*propkp + 2./3.*mp2 + 16./3.*pow(mp2,2)*propkp ); K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s2,mk2,mu2)*G27 * ( - 4./ 9.*mp2*s1*propkp - 4./9.*mp2*s2*propkp - 28./9.*mp2*s3*propkp + 4./9.*mp2 + 32./9.*pow(mp2,2)*propkp ); K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G8 * ( 4*mp2* s3*propkp + 16./3.*mp2 - 32./3.*pow(mp2,2)*propkp - 4./3.*s3 ) ; K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s3,mp2,mu2)*G27 * ( 8./3. *mp2*s3*propkp - 208./9.*mp2 - 64./9.*pow(mp2,2)*propkp + 52./ 9.*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G8 * ( 4*mp2* s3*propkp + 20./3.*mp2 - 32./3.*pow(mp2,2)*propkp ); K3Pp4Ix10112 += + C24b(0.,mp2,mp2,mp2,s3,mk2,mu2)*G27 * ( 8./3. *mp2*s3*propkp + 40./9.*mp2 - 64./9.*pow(mp2,2)*propkp ); K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G8 * ( 1./3.* mp2*s1*propkp + 1./3.*mp2*s2*propkp + 7./3.*mp2*s3*propkp + 8./ 3.*mp2 - 8./3.*pow(mp2,2)*propkp + 1./3.*mk2 - 1./3.*s1 - 1./3. *s2 - 7./3.*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s1,mp2,mu2)*G27 * ( 2./9. *mp2*s1*propkp + 2./9.*mp2*s2*propkp + 14./9.*mp2*s3*propkp + 16./9.*mp2 - 16./9.*pow(mp2,2)*propkp + 2./9.*mk2 - 2./9.*s1 - 2./9.*s2 - 14./9.*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G8 * ( 1./3.* mp2*s1*propkp + 1./3.*mp2*s2*propkp + 7./3.*mp2*s3*propkp + 3* mp2 - 8./3.*pow(mp2,2)*propkp - 4*mk2 + 2./3.*s1 - 2*s2 - 2*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s1,mk2,mu2)*G27 * ( 2./9. *mp2*s1*propkp + 2./9.*mp2*s2*propkp + 14./9.*mp2*s3*propkp + 46./3.*mp2 - 16./9.*pow(mp2,2)*propkp - 8./3.*mk2 - 26./9.*s1 + 2*s2 - 18*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G8 * ( 1./3.* mp2*s1*propkp + 1./3.*mp2*s2*propkp + 7./3.*mp2*s3*propkp + 8./ 3.*mp2 - 8./3.*pow(mp2,2)*propkp + 1./3.*mk2 - 1./3.*s1 - 1./3. *s2 - 7./3.*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s2,mp2,mu2)*G27 * ( 2./9. *mp2*s1*propkp + 2./9.*mp2*s2*propkp + 14./9.*mp2*s3*propkp + 16./9.*mp2 - 16./9.*pow(mp2,2)*propkp + 2./9.*mk2 - 2./9.*s1 - 2./9.*s2 - 14./9.*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G8 * ( 1./3.* mp2*s1*propkp + 1./3.*mp2*s2*propkp + 7./3.*mp2*s3*propkp + 3* mp2 - 8./3.*pow(mp2,2)*propkp - 4*mk2 - 2*s1 + 2./3.*s2 - 2*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s2,mk2,mu2)*G27 * ( 2./9. *mp2*s1*propkp + 2./9.*mp2*s2*propkp + 14./9.*mp2*s3*propkp + 46./3.*mp2 - 16./9.*pow(mp2,2)*propkp - 8./3.*mk2 + 2*s1 - 26./ 9.*s2 - 18*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G8 * ( - 4* mp2*s3*propkp - 32./3.*mp2 + 32./3.*pow(mp2,2)*propkp + 4*s3 ) ; K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s3,mp2,mu2)*G27 * ( - 8./ 3.*mp2*s3*propkp - 64./9.*mp2 + 64./9.*pow(mp2,2)*propkp + 8./ 3.*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G8 * ( - 4* mp2*s3*propkp - 32./3.*mp2 + 32./3.*pow(mp2,2)*propkp - 4./3.* mk2 + 8./3.*s3 ); K3Pp4Ix10112 += + C24b(0.,mp2,mk2,mp2,s3,mk2,mu2)*G27 * ( - 8./ 3.*mp2*s3*propkp - 184./9.*mp2 + 64./9.*pow(mp2,2)*propkp - 128./9.*mk2 + 76./9.*s3 ); K3Pp4Ix10112 += + C24b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G8 * ( 1./3.* mp2*s1*propkp + 1./3.*mp2*s2*propkp + 7./3.*mp2*s3*propkp + 8./ 3.*mp2 - 8./3.*pow(mp2,2)*propkp + 1./3.*mk2 - 1./3.*s1 - 1./3. *s2 - 7./3.*s3 ); K3Pp4Ix10112 += + C24b(0.,mk2,mk2,mp2,s1,mk2,mu2)*G27 * ( 2./9. *mp2*s1*propkp + 2./9.*mp2*s2*propkp + 14./9.*mp2*s3*propkp + 16./9.*mp2 - 16./9.*pow(mp2,2)*propkp + 2./9.*mk2 - 2./9.*s1 - 2./9.*s2 - 14./9.*s3 ); K3Pp4Ix10112 += + C24b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G8 * ( 1./3.* mp2*s1*propkp + 1./3.*mp2*s2*propkp + 7./3.*mp2*s3*propkp + 8./ 3.*mp2 - 8./3.*pow(mp2,2)*propkp + 1./3.*mk2 - 1./3.*s1 - 1./3. *s2 - 7./3.*s3 ); K3Pp4Ix10112 += + C24b(0.,mk2,mk2,mp2,s2,mk2,mu2)*G27 * ( 2./9. *mp2*s1*propkp + 2./9.*mp2*s2*propkp + 14./9.*mp2*s3*propkp + 16./9.*mp2 - 16./9.*pow(mp2,2)*propkp + 2./9.*mk2 - 2./9.*s1 - 2./9.*s2 - 14./9.*s3 ); K3Pp4Ix10112 += + dqsqBb(0.,mp2,mp2,mu2)*G8 * ( - 12*mp2*mk2 + 12*mp2*s3 - 12*pow(mp2,2) ); K3Pp4Ix10112 += + dqsqBb(0.,mp2,mp2,mu2)*G27 * ( 12*mp2*mk2 - 52*mp2*s3 + 52*pow(mp2,2) ); K3Pp4Ix10112 += + dqsqBb(0.,mk2,mk2,mu2)*G8 * ( - 4*mp2*mk2 + 4*mk2*s3 - 4*pow(mk2,2) ); K3Pp4Ix10112 += + dqsqBb(0.,mk2,mk2,mu2)*G27 * ( 52./3.*mp2*mk2 - 52./3.*mk2*s3 + 4*pow(mk2,2) ); K3Pp4Ix10112 += + dqsqB1b(0.,mp2,mp2,mu2)*G8 * ( 12*mp2*mk2 - 12*mp2*s3 + 12*pow(mp2,2) ); K3Pp4Ix10112 += + dqsqB1b(0.,mp2,mp2,mu2)*G27 * ( - 12*mp2*mk2 + 52*mp2*s3 - 52*pow(mp2,2) ); K3Pp4Ix10112 += + dqsqB1b(0.,mk2,mk2,mu2)*G8 * ( 4*mp2*mk2 - 4* mk2*s3 + 4*pow(mk2,2) ); K3Pp4Ix10112 += + dqsqB1b(0.,mk2,mk2,mu2)*G27 * ( - 52./3.*mp2 *mk2 + 52./3.*mk2*s3 - 4*pow(mk2,2) ); return K3Pp4Ix10112*e2/(pow(fpip,3)*fkp); } dcomplex Ap4EKPppm(double s1, double s2,double s3){ using namespace inputpars; double mko2 = mko*mko, mkp2 = mkp*mkp,mpp2 = mpp*mpp,mpo2=mpo*mpo, me2 = me*me,mu2 = mu*mu; double mp2 = mpp2,mk2=mkp2; double props1 = 1./s1,props2 = 1./s2; dcomplex K3Pp4Ex10112 = dcomplex(0.); K3Pp4Ex10112 = + G8*props1 * ( - 2./3.*mp2*pi16*s2 + 2./3.*mp2*pi16*s3 - 2./3. *mk2*pi16*s2 + 2./3.*mk2*pi16*s3 + 4./3.*N14r*s1*s2 - 4./3.* N14r*s1*s3 - 4./3.*N15r*s1*s2 + 4./3.*N15r*s1*s3 + 4*L9r*s1*s2 - 4*L9r*s1*s3 + 2./9.*pi16*s1*s2 - 2./9.*pi16*s1*s3 ); K3Pp4Ex10112 += + G8*props2 * ( - 2./3.*mp2*pi16*s1 + 2./3.*mp2 *pi16*s3 - 2./3.*mk2*pi16*s1 + 2./3.*mk2*pi16*s3 + 4./3.*N14r* s1*s2 - 4./3.*N14r*s2*s3 - 4./3.*N15r*s1*s2 + 4./3.*N15r*s2*s3 + 4*L9r*s1*s2 - 4*L9r*s2*s3 + 2./9.*pi16*s1*s2 - 2./9.*pi16* s2*s3 ); K3Pp4Ex10112 += + G27*props1 * ( 26./9.*mp2*pi16*s2 - 26./9.*mp2 *pi16*s3 + 26./9.*mk2*pi16*s2 - 26./9.*mk2*pi16*s3 + 4./3.* D13r*s1*s2 - 4./3.*D13r*s1*s3 - 2./3.*D15r*s1*s2 + 2./3.*D15r* s1*s3 + 8./3.*L9r*s1*s2 - 8./3.*L9r*s1*s3 - 26./27.*pi16*s1*s2 + 26./27.*pi16*s1*s3 ); K3Pp4Ex10112 += + G27*props2 * ( 26./9.*mp2*pi16*s1 - 26./9.*mp2 *pi16*s3 + 26./9.*mk2*pi16*s1 - 26./9.*mk2*pi16*s3 + 4./3.* D13r*s1*s2 - 4./3.*D13r*s2*s3 - 2./3.*D15r*s1*s2 + 2./3.*D15r* s2*s3 + 8./3.*L9r*s1*s2 - 8./3.*L9r*s2*s3 - 26./27.*pi16*s1*s2 + 26./27.*pi16*s2*s3 ); K3Pp4Ex10112 += + Ab(mp2,mu2)*G8*props1 * ( 2./3.*s2 - 2./3.*s3 ); K3Pp4Ex10112 += + Ab(mp2,mu2)*G8*props2 * ( 2./3.*s1 - 2./3.*s3 ); K3Pp4Ex10112 += + Ab(mp2,mu2)*G27*props1 * ( - 26./9.*s2 + 26./ 9.*s3 ); K3Pp4Ex10112 += + Ab(mp2,mu2)*G27*props2 * ( - 26./9.*s1 + 26./ 9.*s3 ); K3Pp4Ex10112 += + Ab(mk2,mu2)*G8*props1 * ( 2./3.*s2 - 2./3.*s3 ); K3Pp4Ex10112 += + Ab(mk2,mu2)*G8*props2 * ( 2./3.*s1 - 2./3.*s3 ); K3Pp4Ex10112 += + Ab(mk2,mu2)*G27*props1 * ( - 26./9.*s2 + 26./ 9.*s3 ); K3Pp4Ex10112 += + Ab(mk2,mu2)*G27*props2 * ( - 26./9.*s1 + 26./ 9.*s3 ); K3Pp4Ex10112 += + Bb(mp2,mp2,s1,mu2)*G8*props1 * ( - 2./3.*mp2* s2 + 2./3.*mp2*s3 + 1./6.*s1*s2 - 1./6.*s1*s3 ); K3Pp4Ex10112 += + Bb(mp2,mp2,s1,mu2)*G27*props1 * ( 26./9.*mp2* s2 - 26./9.*mp2*s3 - 13./18.*s1*s2 + 13./18.*s1*s3 ); K3Pp4Ex10112 += + Bb(mp2,mp2,s2,mu2)*G8*props2 * ( - 2./3.*mp2* s1 + 2./3.*mp2*s3 + 1./6.*s1*s2 - 1./6.*s2*s3 ); K3Pp4Ex10112 += + Bb(mp2,mp2,s2,mu2)*G27*props2 * ( 26./9.*mp2* s1 - 26./9.*mp2*s3 - 13./18.*s1*s2 + 13./18.*s2*s3 ); K3Pp4Ex10112 += + Bb(mk2,mk2,s1,mu2)*G8*props1 * ( - 2./3.*mk2* s2 + 2./3.*mk2*s3 + 1./6.*s1*s2 - 1./6.*s1*s3 ); K3Pp4Ex10112 += + Bb(mk2,mk2,s1,mu2)*G27*props1 * ( 26./9.*mk2* s2 - 26./9.*mk2*s3 - 13./18.*s1*s2 + 13./18.*s1*s3 ); K3Pp4Ex10112 += + Bb(mk2,mk2,s2,mu2)*G8*props2 * ( - 2./3.*mk2* s1 + 2./3.*mk2*s3 + 1./6.*s1*s2 - 1./6.*s2*s3 ); K3Pp4Ex10112 += + Bb(mk2,mk2,s2,mu2)*G27*props2 * ( 26./9.*mk2* s1 - 26./9.*mk2*s3 - 13./18.*s1*s2 + 13./18.*s2*s3 ); return K3Pp4Ex10112*e2/(pow(fpip,3)*fkp); } dcomplex AfullemKLooo(double s1, double s2, double s3){ using namespace photoncutoffs; return Ap2KLooo(s1,s2,s3)+Ap4KLooo(s1,s2,s3)+Ap4IKLooo(s1,s2,s3)+ Ap4KLoooIR(s1,s2,s3,eg3pi,egfpi); } dcomplex AfullemKLpmo(double s1, double s2, double s3){ using namespace photoncutoffs; return Ap2KLpmo(s1,s2,s3)+Ap4KLpmo(s1,s2,s3)+Ap4IKLpmo(s1,s2,s3)+ Ap4KLpmoIR(s1,s2,s3,eg3pi,egfpi); } dcomplex AfullemKSpmo(double s1, double s2, double s3){ using namespace photoncutoffs; return Ap2KSpmo(s1,s2,s3)+Ap4KSpmo(s1,s2,s3)+Ap4IKSpmo(s1,s2,s3)+ Ap4EKSpmo(s1,s2,s3)+Ap4KSpmoIR(s1,s2,s3,eg3pi,egfpi); } dcomplex AfullemKPoop(double s1, double s2, double s3){ using namespace photoncutoffs; return Ap2KPoop(s1,s2,s3)+Ap4KPoop(s1,s2,s3)+Ap4IKPoop(s1,s2,s3)+ Ap4KPoopIR(s1,s2,s3,eg3pi,egfpi); } dcomplex AfullemKPppm(double s1, double s2, double s3){ using namespace photoncutoffs; using namespace inputpars; double thresh = 2.*mpp+egcoulomb; if( (sqrt(s1) < thresh) ||(sqrt(s2) < thresh)||(sqrt(s3) < thresh) ) return dcomplex(0.); return Ap2KPppm(s1,s2,s3)+Ap4KPppm(s1,s2,s3)+Ap4IKPppm(s1,s2,s3)+ Ap4EKPppm(s1,s2,s3)+Ap4KPppmIR(s1,s2,s3,eg3pi,egfpi); } dcomplex AfullcutKPppm(double s1, double s2, double s3){ using namespace photoncutoffs; using namespace inputpars; double thresh = 2.*mpp+egcoulomb; if( (sqrt(s1) < thresh) ||(sqrt(s2) < thresh)||(sqrt(s3) < thresh) ) return dcomplex(0.); return Ap2KPppm(s1,s2,s3)+Ap4KPppm(s1,s2,s3); }