* contains final results allowing for later redefinitions * here all masses and Fpi at p^4 are the physical masses and decay constants * see also Mpi.frm and Fpi.frm where the p^4 quantities are not renormalized * pion mass : mpi^2 = P2x11 + P4realx11 + P6finx11 P2x11 = + mhat; P4realx11 = + L4r*Fpi^-2 * ( - 8*mp2^2 - 16*mp2*mk2 ) + L5r*Fpi^-2 * ( - 8*mp2^2 ) + L6r*Fpi^-2 * ( 16*mp2^2 + 32*mp2*mk2 ) + L8r*Fpi^-2 * ( 16*mp2^2 ) + Ab(mp2)*Fpi^-2 * ( - 1/2*mp2 ) + Ab(me2)*Fpi^-2 * ( 1/6*mp2 ); P6finx11 = + F^-2 * ( - 32*mp2^3*CC12 - 32*mp2^3*CC13 - 16*mp2^3*CC14 - 16*mp2^3* CC15 - 48*mp2^3*CC16 - 16*mp2^3*CC17 + 48*mp2^3*CC19 + 80*mp2^3*CC20 + 48*mp2^3*CC21 + 32*mp2^3*CC31 + 32*mp2^3*CC32 - 64*mp2^2*mk2*CC13 - 32*mp2^2*mk2*CC15 + 64*mp2^2*mk2*CC16 + 192*mp2^2*mk2*CC21 + 64* mp2^2*mk2*CC32 - 64*mp2*mk2^2*CC16 + 64*mp2*mk2^2*CC20 + 192*mp2* mk2^2*CC21 ) + F^-4*pi16^2 * ( - 527/1296*mp2^3*pi^2 - 3217/1728*mp2^3 - 37/324* mp2^2*mk2*pi^2 - 139/216*mp2^2*mk2 - 11/36*mp2*mk2^2*pi^2 - 15/16*mp2 *mk2^2 ) + F^-4*pi16*L1r * ( 4*mp2^3 ) + F^-4*pi16*L2r * ( 74/9*mp2^3 - 16/9*mp2^2*mk2 + 104/9*mp2*mk2^2 ) + F^-4*pi16*L3r * ( 56/27*mp2^3 - 16/27*mp2^2*mk2 + 86/27*mp2*mk2^2 ) + F^-4*L4r^2 * ( - 128*mp2^3 - 512*mp2^2*mk2 - 512*mp2*mk2^2 ) + F^-4*L4r*L5r * ( - 256*mp2^3 - 384*mp2^2*mk2 - 128*mp2*mk2^2 ) + F^-4*L4r*L6r * ( 512*mp2^3 + 2048*mp2^2*mk2 + 2048*mp2*mk2^2 ) + F^-4*L4r*L8r * ( 512*mp2^3 + 768*mp2^2*mk2 + 256*mp2*mk2^2 ) + F^-4*L5r^2 * ( - 128*mp2^3 ) + F^-4*L5r*L6r * ( 512*mp2^3 + 768*mp2^2*mk2 + 256*mp2*mk2^2 ) + F^-4*L5r*L8r * ( 512*mp2^3 ) + F^-4*L6r^2 * ( - 512*mp2^3 - 2048*mp2^2*mk2 - 2048*mp2*mk2^2 ) + F^-4*L6r*L8r * ( - 1024*mp2^3 - 1536*mp2^2*mk2 - 512*mp2*mk2^2 ) + F^-4*L8r^2 * ( - 512*mp2^3 ) + F^-4 * ( 5/6*Hbb(mp2,mp2,mp2,mp2,mu2)*mp2^2 - 5/8*Hbb(mp2,mk2,mk2,mp2 ,mu2)*mp2^2 + 1/18*Hbb(mp2,me2,me2,mp2,mu2)*mp2^2 + Hbb(mk2,mp2,mk2, mp2,mu2)*mp2*mk2 - 5/6*Hbb(mk2,mk2,me2,mp2,mu2)*mp2^2 - 1/8*Hbb(me2, mk2,mk2,mp2,mu2)*mp2^2 + 1/2*Hbb(me2,mk2,mk2,mp2,mu2)*mp2*mk2 + H1bb( mp2,mk2,mk2,mp2,mu2)*mp2^2 + 2*H1bb(mk2,mk2,me2,mp2,mu2)*mp2^2 + 3* H21bb(mp2,mp2,mp2,mp2,mu2)*mp2^2 - 3/8*H21bb(mp2,mk2,mk2,mp2,mu2)* mp2^2 + 3*H21bb(mk2,mp2,mk2,mp2,mu2)*mp2^2 + 9/8*H21bb(me2,mk2,mk2, mp2,mu2)*mp2^2 ) + Ab(mp2)*F^-4*pi16 * ( mp2^2 + mp2*mk2 ) + Ab(mp2)*F^-4*L1r * ( 56*mp2^2 ) + Ab(mp2)*F^-4*L2r * ( 32*mp2^2 ) + Ab(mp2)*F^-4*L3r * ( 28*mp2^2 ) + Ab(mp2)*F^-4*L4r * ( - 72*mp2^2 - 40*mp2*mk2 ) + Ab(mp2)*F^-4*L5r * ( - 48*mp2^2 ) + Ab(mp2)*F^-4*L6r * ( 128*mp2^2 + 80*mp2*mk2 ) + Ab(mp2)*F^-4*L8r * ( 88*mp2^2 ) + Ab(mp2)*Ab(mp2)*F^-4 * ( - 589/144*mp2 - 1/2*mk2 ) + Ab(mp2)*Ab(mk2)*F^-4 * ( - mp2 ) + Ab(mp2)*Ab(me2)*F^-4 * ( 5/12*mp2 ) + Ab(mk2)*F^-4*pi16 * ( mp2*mk2 ) + Ab(mk2)*F^-4*L1r * ( 64*mp2*mk2 ) + Ab(mk2)*F^-4*L2r * ( 16*mp2*mk2 ) + Ab(mk2)*F^-4*L3r * ( 20*mp2*mk2 ) + Ab(mk2)*F^-4*L4r * ( - 8*mp2^2 - 80*mp2*mk2 ) + Ab(mk2)*F^-4*L5r * ( - 8*mp2^2 - 16*mp2*mk2 ) + Ab(mk2)*F^-4*L6r * ( 16*mp2^2 + 96*mp2*mk2 ) + Ab(mk2)*F^-4*L8r * ( 16*mp2^2 + 32*mp2*mk2 ) + Ab(mk2)*Ab(mk2)*F^-4 * ( - 1/6*mp2^2*mk2^-1 - 7/4*mp2 ) + Ab(mk2)*Ab(me2)*F^-4 * ( - 1/3*mp2 ) + Ab(me2)*F^-4*L1r * ( - 16/3*mp2^2 + 64/3*mp2*mk2 ) + Ab(me2)*F^-4*L2r * ( - 4/3*mp2^2 + 16/3*mp2*mk2 ) + Ab(me2)*F^-4*L3r * ( - 4/3*mp2^2 + 16/3*mp2*mk2 ) + Ab(me2)*F^-4*L4r * ( 8*mp2^2 - 24*mp2*mk2 ) + Ab(me2)*F^-4*L5r * ( 32/9*mp2^2 - 32/9*mp2*mk2 ) + Ab(me2)*F^-4*L6r * ( - 32/3*mp2^2 + 80/3*mp2*mk2 ) + Ab(me2)*F^-4*L7r * ( 64/3*mp2^2 - 64/3*mp2*mk2 ) + Ab(me2)*F^-4*L8r * ( 8/3*mp2^2 ) + Ab(me2)*Ab(me2)*F^-4 * ( - 5/36*mp2^2*me2^-1 - 29/144*mp2 ); * Kaon mass mK^2 = P2x44 + P4realx44 + P6finx44 P2x44 = + 1/2*mhat + 1/2*mstrange; P4realx44 = + L4r*Fpi^-2 * ( - 8*mp2*mk2 - 16*mk2^2 ) + L5r*Fpi^-2 * ( - 8*mk2^2 ) + L6r*Fpi^-2 * ( 16*mp2*mk2 + 32*mk2^2 ) + L8r*Fpi^-2 * ( 16*mk2^2 ) + Ab(me2)*Fpi^-2 * ( - 1/12*mp2 - 1/4*me2 ); P6finx44 = + F^-2 * ( - 16*mp2^2*mk2*CC14 - 48*mp2^2*mk2*CC16 + 16*mp2^2*mk2*CC17 + 48*mp2^2*mk2*CC19 + 48*mp2^2*mk2*CC20 + 48*mp2^2*mk2*CC21 - 32*mp2 *mk2^2*CC13 + 32*mp2*mk2^2*CC14 - 16*mp2*mk2^2*CC15 + 64*mp2*mk2^2* CC16 - 32*mp2*mk2^2*CC17 - 96*mp2*mk2^2*CC19 - 32*mp2*mk2^2*CC20 + 192*mp2*mk2^2*CC21 + 32*mp2*mk2^2*CC32 - 32*mk2^3*CC12 - 64*mk2^3* CC13 - 32*mk2^3*CC14 - 32*mk2^3*CC15 - 64*mk2^3*CC16 + 96*mk2^3*CC19 + 128*mk2^3*CC20 + 192*mk2^3*CC21 + 32*mk2^3*CC31 + 64*mk2^3*CC32 ) + F^-4*pi16^2 * ( - 1/8*mp2^2*mk2*pi^2 - 13/24*mp2^2*mk2 - 73/648*mp2* mk2^2*pi^2 - 19/108*mp2*mk2^2 - 763/1296*mk2^3*pi^2 - 4709/1728*mk2^3 ) + F^-4*pi16*L1r * ( 4*mk2^3 ) + F^-4*pi16*L2r * ( 56/9*mp2^2*mk2 - 16/9*mp2*mk2^2 + 122/9*mk2^3 ) + F^-4*pi16*L3r * ( 41/27*mp2^2*mk2 - 4/27*mp2*mk2^2 + 89/27*mk2^3 ) + F^-4*L4r^2 * ( - 128*mp2^2*mk2 - 512*mp2*mk2^2 - 512*mk2^3 ) + F^-4*L4r*L5r * ( - 128*mp2^2*mk2 - 256*mp2*mk2^2 - 384*mk2^3 ) + F^-4*L4r*L6r * ( 512*mp2^2*mk2 + 2048*mp2*mk2^2 + 2048*mk2^3 ) + F^-4*L4r*L8r * ( 128*mp2^2*mk2 + 384*mp2*mk2^2 + 1024*mk2^3 ) + F^-4*L5r^2 * ( - 64*mp2*mk2^2 - 64*mk2^3 ) + F^-4*L5r*L6r * ( 256*mp2^2*mk2 + 512*mp2*mk2^2 + 768*mk2^3 ) + F^-4*L5r*L8r * ( 128*mp2*mk2^2 + 384*mk2^3 ) + F^-4*L6r^2 * ( - 512*mp2^2*mk2 - 2048*mp2*mk2^2 - 2048*mk2^3 ) + F^-4*L6r*L8r * ( - 256*mp2^2*mk2 - 768*mp2*mk2^2 - 2048*mk2^3 ) + F^-4*L8r^2 * ( - 512*mk2^3 ) + F^-4 * ( 3/4*Hbb(mp2,mp2,mk2,mk2,mu2)*mp2*mk2 + 3/8*Hbb(mp2,mp2,mk2, mk2,mu2)*mk2^2 + 1/4*Hbb(mp2,mk2,me2,mk2,mu2)*mk2^2 - 3/32*Hbb(mk2, mp2,mp2,mk2,mu2)*mk2^2 + 9/16*Hbb(mk2,mp2,me2,mk2,mu2)*mk2^2 + 3/4* Hbb(mk2,mk2,mk2,mk2,mu2)*mk2^2 + 181/288*Hbb(mk2,me2,me2,mk2,mu2)* mk2^2 - 3/2*H1bb(mp2,mp2,mk2,mk2,mu2)*mk2^2 - 3/2*H1bb(mk2,mp2,me2, mk2,mu2)*mk2^2 - 5/4*H1bb(mk2,me2,me2,mk2,mu2)*mk2^2 + 9/4*H21bb(mp2, mp2,mk2,mk2,mu2)*mk2^2 - 9/32*H21bb(mk2,mp2,mp2,mk2,mu2)*mk2^2 + 27/ 16*H21bb(mk2,mp2,me2,mk2,mu2)*mk2^2 + 9/4*H21bb(mk2,mk2,mk2,mk2,mu2)* mk2^2 + 27/32*H21bb(mk2,me2,me2,mk2,mu2)*mk2^2 ) + Ab(mp2)*F^-4*pi16 * ( 3/4*mk2^2 ) + Ab(mp2)*F^-4*L1r * ( 48*mp2*mk2 ) + Ab(mp2)*F^-4*L2r * ( 12*mp2*mk2 ) + Ab(mp2)*F^-4*L3r * ( 15*mp2*mk2 ) + Ab(mp2)*F^-4*L4r * ( - 68*mp2*mk2 - 32*mk2^2 ) + Ab(mp2)*F^-4*L5r * ( - 12*mp2*mk2 - 16*mk2^2 ) + Ab(mp2)*F^-4*L6r * ( 88*mp2*mk2 + 64*mk2^2 ) + Ab(mp2)*F^-4*L8r * ( 24*mp2*mk2 + 32*mk2^2 ) + Ab(mp2)*Ab(mp2)*F^-4 * ( - 1/2*mp2^-1*mk2^2 - 27/32*mk2 ) + Ab(mp2)*Ab(mk2)*F^-4 * ( - 3/4*mk2 ) + Ab(mp2)*Ab(me2)*F^-4 * ( 1/12*mp2 - 41/48*mk2 ) + Ab(mk2)*F^-4*pi16 * ( 3/4*mp2*mk2 + 3/4*mk2^2 ) + Ab(mk2)*F^-4*L1r * ( 72*mk2^2 ) + Ab(mk2)*F^-4*L2r * ( 36*mk2^2 ) + Ab(mk2)*F^-4*L3r * ( 30*mk2^2 ) + Ab(mk2)*F^-4*L4r * ( - 8*mp2*mk2 - 80*mk2^2 ) + Ab(mk2)*F^-4*L5r * ( - 32*mk2^2 ) + Ab(mk2)*F^-4*L6r * ( 16*mp2*mk2 + 128*mk2^2 ) + Ab(mk2)*F^-4*L8r * ( 64*mk2^2 ) + Ab(mk2)*Ab(mk2)*F^-4 * ( - 3/8*mp2 - 251/72*mk2 ) + Ab(mk2)*Ab(me2)*F^-4 * ( - 2/3*mk2 ) + Ab(me2)*F^-4*pi16 * ( 1/4*mp2*mk2 + 1/2*mk2^2 ) + Ab(me2)*F^-4*L1r * ( - 16/3*mp2*mk2 + 64/3*mk2^2 ) + Ab(me2)*F^-4*L2r * ( - 4/3*mp2*mk2 + 16/3*mk2^2 ) + Ab(me2)*F^-4*L3r * ( - 7/3*mp2*mk2 + 28/3*mk2^2 ) + Ab(me2)*F^-4*L4r * ( 4*mp2*mk2 - 32*mk2^2 ) + Ab(me2)*F^-4*L5r * ( - 8/9*mp2^2 + 4/3*mp2*mk2 - 112/9*mk2^2 ) + Ab(me2)*F^-4*L6r * ( - 8/3*mp2*mk2 + 128/3*mk2^2 ) + Ab(me2)*F^-4*L7r * ( 32/3*mp2^2 - 32*mp2*mk2 + 64/3*mk2^2 ) + Ab(me2)*F^-4*L8r * ( 16/3*mp2^2 - 56/3*mp2*mk2 + 32*mk2^2 ) + Ab(me2)*Ab(me2)*F^-4 * ( - 43/1152*mp2^2*me2^-1 - 25/128*mp2 - 5/36* mk2 ); * eta mass meta^2 = P2x88 + P4realx88 + P6finx88 P2x88 = + 1/3*mhat + 2/3*mstrange - qext.qext; P4realx88 = + L4r*Fpi^-2 * ( - 8*mp2*me2 - 16*mk2*me2 ) + L5r*Fpi^-2 * ( 8/3*mp2*me2 - 32/3*mk2*me2 ) + L6r*Fpi^-2 * ( - 16/3*mp2^2 + 32/3*mp2*mk2 + 128/3*mk2^2 ) + L7r*Fpi^-2 * ( 128/3*mp2^2 - 256/3*mp2*mk2 + 128/3*mk2^2 ) + L8r*Fpi^-2 * ( 16*mp2^2 - 128/3*mp2*mk2 + 128/3*mk2^2 ) + Ab(mp2)*Fpi^-2 * ( 1/2*mp2 ) + Ab(mk2)*Fpi^-2 * ( - 1/3*mp2 - me2 ) + Ab(me2)*Fpi^-2 * ( - 7/18*mp2 + 8/9*mk2 ); P6finx88 = + F^-2 * ( 32/27*mp2^3*CC12 - 32/9*mp2^3*CC13 + 16/3*mp2^3*CC14 - 16/9* mp2^3*CC15 + 16*mp2^3*CC16 + 16/3*mp2^3*CC17 + 128/9*mp2^3*CC18 - 16* mp2^3*CC19 + 16*mp2^3*CC20 - 16*mp2^3*CC21 - 32/3*mp2^3*CC31 + 32* mp2^3*CC32 - 128/9*mp2^2*mk2*CC12 + 64/3*mp2^2*mk2*CC13 - 320/9*mp2^2 *mk2*CC14 + 32/3*mp2^2*mk2*CC15 - 256/3*mp2^2*mk2*CC16 - 320/9*mp2^2* mk2*CC17 - 256/3*mp2^2*mk2*CC18 + 192*mp2^2*mk2*CC19 + 64*mp2^2*mk2* CC20 + 128*mp2^2*mk2*CC31 - 64/3*mp2^2*mk2*CC32 + 512/3*mp2^2*mk2* CC33 + 512/9*mp2*mk2^2*CC12 + 640/9*mp2*mk2^2*CC14 + 320/3*mp2*mk2^2* CC16 + 640/9*mp2*mk2^2*CC17 + 128*mp2*mk2^2*CC18 - 384*mp2*mk2^2*CC19 - 192*mp2*mk2^2*CC20 + 192*mp2*mk2^2*CC21 - 256*mp2*mk2^2*CC31 - 256/ 3*mp2*mk2^2*CC32 - 1024/3*mp2*mk2^2*CC33 - 2048/27*mk2^3*CC12 - 1024/ 9*mk2^3*CC13 - 512/9*mk2^3*CC14 - 512/9*mk2^3*CC15 - 256/3*mk2^3*CC16 - 512/9*mk2^3*CC17 - 512/9*mk2^3*CC18 + 256*mk2^3*CC19 + 256*mk2^3* CC20 + 256*mk2^3*CC21 + 512/3*mk2^3*CC31 + 512/3*mk2^3*CC32 + 512/3* mk2^3*CC33 ) + F^-4*pi16^2 * ( - 91/11664*mp2^3*pi^2 - 1781/15552*mp2^3 - 269/486* mp2^2*mk2*pi^2 - 4133/1296*mp2^2*mk2 + 133/108*mp2*mk2^2*pi^2 + 367/ 48*mp2*mk2^2 - 1091/729*mk2^3*pi^2 - 7567/972*mk2^3 ) + F^-4*pi16*L1r * ( - 4/27*mp2^3 + 16/9*mp2^2*mk2 - 64/9*mp2*mk2^2 + 256/27*mk2^3 ) + F^-4*pi16*L2r * ( - 58/27*mp2^3 + 88/9*mp2^2*mk2 - 88/9*mp2*mk2^2 + 544/27*mk2^3 ) + F^-4*pi16*L3r * ( - 20/27*mp2^3 + 32/9*mp2^2*mk2 - 34/9*mp2*mk2^2 + 152/27*mk2^3 ) + F^-4*L4r^2 * ( 128/3*mp2^3 - 512*mp2*mk2^2 - 2048/3*mk2^3 ) + F^-4*L4r*L5r * ( 256/9*mp2^3 - 128/3*mp2^2*mk2 - 128*mp2*mk2^2 - 5632/ 9*mk2^3 ) + F^-4*L4r*L6r * ( - 512/3*mp2^3 + 2048*mp2*mk2^2 + 8192/3*mk2^3 ) + F^-4*L4r*L7r * ( 1024*mp2^3 - 3072*mp2*mk2^2 + 2048*mk2^3 ) + F^-4*L4r*L8r * ( 1024/3*mp2^3 - 256/3*mp2^2*mk2 - 3328/3*mp2*mk2^2 + 7168/3*mk2^3 ) + F^-4*L5r^2 * ( - 128/9*mp2^3 + 256/3*mp2^2*mk2 - 256/3*mp2*mk2^2 - 1024/9*mk2^3 ) + F^-4*L5r*L6r * ( - 1024/9*mp2^3 + 256/3*mp2^2*mk2 + 1280/3*mp2*mk2^2 + 10240/9*mk2^3 ) + F^-4*L5r*L7r * ( 1024*mp2^3 - 4096/3*mp2^2*mk2 - 1024/3*mp2*mk2^2 + 2048/3*mk2^3 ) + F^-4*L5r*L8r * ( 3584/9*mp2^3 - 6656/9*mp2^2*mk2 - 512/9*mp2*mk2^2 + 8192/9*mk2^3 ) + F^-4*L6r^2 * ( 512/3*mp2^3 - 2048*mp2*mk2^2 - 8192/3*mk2^3 ) + F^-4*L6r*L7r * ( - 4096/3*mp2^3 + 4096*mp2*mk2^2 - 8192/3*mk2^3 ) + F^-4*L6r*L8r * ( - 1024/3*mp2^3 + 512/3*mp2^2*mk2 + 3584/3*mp2*mk2^2 - 4096*mk2^3 ) + F^-4*L7r*L8r * ( - 4096/3*mp2^3 + 4096/3*mp2^2*mk2 + 4096/3*mp2* mk2^2 - 4096/3*mk2^3 ) + F^-4*L8r^2 * ( - 512*mp2^3 + 2048/3*mp2^2*mk2 + 2048/3*mp2*mk2^2 - 4096/3*mk2^3 ) + F^-4 * ( 1/6*Hbb(mp2,mp2,me2,me2,mu2)*mp2^2 + 1/8*Hbb(mp2,mk2,mk2,me2 ,mu2)*mp2^2 + 3/2*Hbb(mp2,mk2,mk2,me2,mu2)*mp2*mk2 - 1/2*Hbb(mk2,mk2, me2,me2,mu2)*mp2^2 + 4*Hbb(mk2,mk2,me2,me2,mu2)*mp2*mk2 - 64/9*Hbb( mk2,mk2,me2,me2,mu2)*mk2^2 + 1/8*Hbb(me2,mk2,mk2,me2,mu2)*mp2^2 - Hbb(me2,mk2,mk2,me2,mu2)*mp2*mk2 + 2*Hbb(me2,mk2,mk2,me2,mu2)*mk2^2 + 49/486*Hbb(me2,me2,me2,me2,mu2)*mp2^2 - 112/243*Hbb(me2,me2,me2, me2,mu2)*mp2*mk2 + 128/243*Hbb(me2,me2,me2,me2,mu2)*mk2^2 + H1bb(mp2, mk2,mk2,me2,mu2)*mp2^2 - 4*H1bb(mp2,mk2,mk2,me2,mu2)*mp2*mk2 + 2* H1bb(mk2,mk2,me2,me2,mu2)*mp2^2 - 40/3*H1bb(mk2,mk2,me2,me2,mu2)*mp2* mk2 + 64/3*H1bb(mk2,mk2,me2,me2,mu2)*mk2^2 + 3/8*H21bb(mp2,mk2,mk2, me2,mu2)*mp2^2 - 3*H21bb(mp2,mk2,mk2,me2,mu2)*mp2*mk2 + 6*H21bb(mp2, mk2,mk2,me2,mu2)*mk2^2 + 3/8*H21bb(me2,mk2,mk2,me2,mu2)*mp2^2 - 3* H21bb(me2,mk2,mk2,me2,mu2)*mp2*mk2 + 6*H21bb(me2,mk2,mk2,me2,mu2)* mk2^2 ) + Ab(mp2)*F^-4*L1r * ( - 16*mp2^2 + 64*mp2*mk2 ) + Ab(mp2)*F^-4*L2r * ( - 4*mp2^2 + 16*mp2*mk2 ) + Ab(mp2)*F^-4*L3r * ( - 4*mp2^2 + 16*mp2*mk2 ) + Ab(mp2)*F^-4*L4r * ( 80/3*mp2^2 - 72*mp2*mk2 - 128/3*mk2^2 ) + Ab(mp2)*F^-4*L5r * ( - 8/9*mp2^2 + 32/3*mp2*mk2 - 256/9*mk2^2 ) + Ab(mp2)*F^-4*L6r * ( - 32*mp2^2 + 272/3*mp2*mk2 + 256/3*mk2^2 ) + Ab(mp2)*F^-4*L7r * ( 192*mp2^2 - 832/3*mp2*mk2 + 256/3*mk2^2 ) + Ab(mp2)*F^-4*L8r * ( 72*mp2^2 - 320/3*mp2*mk2 + 256/3*mk2^2 ) + Ab(mp2)*Ab(mp2)*F^-4 * ( 43/48*mp2 - 3/4*mk2 ) + Ab(mp2)*Ab(mk2)*F^-4 * ( - 2/3*mp2 - 8/3*mk2 ) + Ab(mp2)*Ab(me2)*F^-4 * ( - 35/36*mp2 + 16/9*mk2 ) + Ab(mk2)*F^-4*pi16 * ( - 1/3*mp2^2 + 2/3*mp2*mk2 + 8/3*mk2^2 ) + Ab(mk2)*F^-4*L1r * ( - 64/3*mp2*mk2 + 256/3*mk2^2 ) + Ab(mk2)*F^-4*L2r * ( - 16/3*mp2*mk2 + 64/3*mk2^2 ) + Ab(mk2)*F^-4*L3r * ( - 28/3*mp2*mk2 + 112/3*mk2^2 ) + Ab(mk2)*F^-4*L4r * ( 8/3*mp2^2 + 16/3*mp2*mk2 - 128*mk2^2 ) + Ab(mk2)*F^-4*L5r * ( - 40/9*mp2^2 + 16*mp2*mk2 - 512/9*mk2^2 ) + Ab(mk2)*F^-4*L6r * ( - 16/3*mp2^2 - 32/3*mp2*mk2 + 128*mk2^2 ) + Ab(mk2)*F^-4*L7r * ( 256/3*mp2^2 - 640/3*mp2*mk2 + 128*mk2^2 ) + Ab(mk2)*F^-4*L8r * ( 112/3*mp2^2 - 352/3*mp2*mk2 + 128*mk2^2 ) + Ab(mk2)*Ab(mk2)*F^-4 * ( - mp2^2*mk2^-1 + 43/12*mp2 - 53/9*mk2 ) + Ab(mk2)*Ab(me2)*F^-4 * ( 7/9*mp2 - 16/9*mk2 ) + Ab(me2)*F^-4*L1r * ( 8/3*mp2^2 - 64/3*mp2*mk2 + 128/3*mk2^2 ) + Ab(me2)*F^-4*L2r * ( 8/3*mp2^2 - 64/3*mp2*mk2 + 128/3*mk2^2 ) + Ab(me2)*F^-4*L3r * ( 4/3*mp2^2 - 32/3*mp2*mk2 + 64/3*mk2^2 ) + Ab(me2)*F^-4*L4r * ( - 16/3*mp2^2 + 56/3*mp2*mk2 - 64/3*mk2^2 ) + Ab(me2)*F^-4*L5r * ( - 248/27*mp2^2 + 736/27*mp2*mk2 - 704/27*mk2^2 ) + Ab(me2)*F^-4*L6r * ( 64/9*mp2^2 - 368/9*mp2*mk2 + 1024/9*mk2^2 ) + Ab(me2)*F^-4*L7r * ( 512/9*mp2^2 - 1792/9*mp2*mk2 + 1280/9*mk2^2 ) + Ab(me2)*F^-4*L8r * ( 296/9*mp2^2 - 1088/9*mp2*mk2 + 128*mk2^2 ) + Ab(me2)*Ab(me2)*F^-4 * ( - 5/36*mp2^2*me2^-1 + 55/48*mp2 - 25/12*mk2 ); * Fpi Fpi = F*(F2x11+F4x11+F6x11) F2x11 = + 1; F4x11 = + L4r*Fpi^-2 * ( 4*mp2 + 8*mk2 ) + L5r*Fpi^-2 * ( 4*mp2 ) + Ab(mp2)*Fpi^-2 * ( 1 ) + Ab(mk2)*Fpi^-2 * ( 1/2 ); F6x11 = + F^-2 * ( 8*mp2^2*CC14 + 8*mp2^2*CC15 + 24*mp2^2*CC16 + 8*mp2^2 *CC17 + 16*mp2*mk2*CC15 - 32*mp2*mk2*CC16 + 32*mk2^2*CC16 ) + F^-4*pi16^2 * ( 35/288*mp2^2*pi^2 + 41/128*mp2^2 + 1/144*mp2* mk2*pi^2 - 5/32*mp2*mk2 + 11/72*mk2^2*pi^2 + 15/32*mk2^2 ) + F^-4*pi16*L1r * ( - 2*mp2^2 ) + F^-4*pi16*L2r * ( - 37/9*mp2^2 + 8/9*mp2*mk2 - 52/9*mk2^2 ) + F^-4*pi16*L3r * ( - 28/27*mp2^2 + 8/27*mp2*mk2 - 43/27*mk2^2 ) + F^-4*L4r^2 * ( 56*mp2^2 + 224*mp2*mk2 + 224*mk2^2 ) + F^-4*L4r*L5r * ( 112*mp2^2 + 160*mp2*mk2 + 64*mk2^2 ) + F^-4*L4r*L6r * ( - 64*mp2^2 - 256*mp2*mk2 - 256*mk2^2 ) + F^-4*L4r*L8r * ( - 64*mp2^2 - 128*mk2^2 ) + F^-4*L5r^2 * ( 56*mp2^2 ) + F^-4*L5r*L6r * ( - 64*mp2^2 - 128*mp2*mk2 ) + F^-4*L5r*L8r * ( - 64*mp2^2 ) + F^-4 * ( - 1/2*Hbb(mp2,mp2,mp2,mp2)*mp2 + 1/16*Hbb(mp2,mk2, mk2,mp2)*mp2 - 1/2*Hbb(mk2,mp2,mk2,mp2)*mk2 + 1/16*Hbb(me2,mk2 ,mk2,mp2)*mp2 - 1/4*Hbb(me2,mk2,mk2,mp2)*mk2 + 5/12*Hd(mp2,mp2 ,mp2,mp2)*mp2^2 - 5/16*Hd(mp2,mk2,mk2,mp2)*mp2^2 + 1/36*Hd(mp2 ,me2,me2,mp2)*mp2^2 + 1/2*Hd(mk2,mp2,mk2,mp2)*mp2*mk2 - 5/12* Hd(mk2,mk2,me2,mp2)*mp2^2 - 1/16*Hd(me2,mk2,mk2,mp2)*mp2^2 + 1/ 4*Hd(me2,mk2,mk2,mp2)*mp2*mk2 + 1/2*H1d(mp2,mk2,mk2,mp2)*mp2^2 + H1d(mk2,mk2,me2,mp2)*mp2^2 + 3/2*H21d(mp2,mp2,mp2,mp2)* mp2^2 - 3/16*H21d(mp2,mk2,mk2,mp2)*mp2^2 + 3/2*H21d(mk2,mp2, mk2,mp2)*mp2^2 + 9/16*H21d(me2,mk2,mk2,mp2)*mp2^2 ) + Ab(mp2)*F^-4*pi16 * ( - 3/4*mp2 - 1/2*mk2 ) + Ab(mp2)*F^-4*L1r * ( - 28*mp2 ) + Ab(mp2)*F^-4*L2r * ( - 16*mp2 ) + Ab(mp2)*F^-4*L3r * ( - 14*mp2 ) + Ab(mp2)*F^-4*L4r * ( 26*mp2 + 24*mk2 ) + Ab(mp2)*F^-4*L5r * ( 20*mp2 ) + Ab(mp2)*Ab(mp2)*F^-4 * ( 87/32 + 1/4*mp2^-1*mk2 ) + Ab(mp2)*Ab(mk2)*F^-4 * ( 2 ) + Ab(mk2)*F^-4*pi16 * ( - 1/8*mp2 - 1/2*mk2 ) + Ab(mk2)*F^-4*L1r * ( - 32*mk2 ) + Ab(mk2)*F^-4*L2r * ( - 8*mk2 ) + Ab(mk2)*F^-4*L3r * ( - 10*mk2 ) + Ab(mk2)*F^-4*L4r * ( 6*mp2 + 28*mk2 ) + Ab(mk2)*F^-4*L5r * ( 10*mp2 ) + Ab(mk2)*Ab(mk2)*F^-4 * ( 5/4 - 1/8*mp2*mk2^-1 ) + Ab(me2)*F^-4*L1r * ( 8/3*mp2 - 32/3*mk2 ) + Ab(me2)*F^-4*L2r * ( 2/3*mp2 - 8/3*mk2 ) + Ab(me2)*F^-4*L3r * ( 2/3*mp2 - 8/3*mk2 ) + Ab(me2)*F^-4*L4r * ( - 2*mp2 + 8*mk2 ) + Ab(me2)*Ab(me2)*F^-4 * ( 3/32 ); * FK FK = F*(F2x44+F4x44+F6x44) F2x44 = + 1; F4x44 = + L4r*Fpi^-2 * ( 4*mp2 + 8*mk2 ) + L5r*Fpi^-2 * ( 4*mk2 ) + Ab(mp2)*Fpi^-2 * ( 3/8 ) + Ab(mk2)*Fpi^-2 * ( 3/4 ) + Ab(me2)*Fpi^-2 * ( 3/8 ); F6x44 = + F^-2 * ( 8*mp2^2*CC14 + 24*mp2^2*CC16 - 8*mp2^2*CC17 - 16*mp2* mk2*CC14 + 8*mp2*mk2*CC15 - 32*mp2*mk2*CC16 + 16*mp2*mk2*CC17 + 16*mk2^2*CC14 + 16*mk2^2*CC15 + 32*mk2^2*CC16 ) + F^-4*pi16^2 * ( 1/16*mp2^2*pi^2 + 13/48*mp2^2 + 1/32*mp2*mk2* pi^2 - 29/192*mp2*mk2 + 3/16*mk2^2*pi^2 + 197/384*mk2^2 ) + F^-4*pi16*L1r * ( - 2*mk2^2 ) + F^-4*pi16*L2r * ( - 28/9*mp2^2 + 8/9*mp2*mk2 - 61/9*mk2^2 ) + F^-4*pi16*L3r * ( - 41/54*mp2^2 + 2/27*mp2*mk2 - 89/54*mk2^2 ) + F^-4*L4r^2 * ( 56*mp2^2 + 224*mp2*mk2 + 224*mk2^2 ) + F^-4*L4r*L5r * ( 64*mp2^2 + 112*mp2*mk2 + 160*mk2^2 ) + F^-4*L4r*L6r * ( - 64*mp2^2 - 256*mp2*mk2 - 256*mk2^2 ) + F^-4*L4r*L8r * ( - 64*mp2^2 - 128*mk2^2 ) + F^-4*L5r^2 * ( 32*mp2*mk2 + 24*mk2^2 ) + F^-4*L5r*L6r * ( - 64*mp2*mk2 - 128*mk2^2 ) + F^-4*L5r*L8r * ( - 64*mk2^2 ) + F^-4 * ( - 3/8*Hbb(mp2,mp2,mk2,mk2)*mp2 + 3/64*Hbb(mk2,mp2, mp2,mk2)*mk2 - 9/32*Hbb(mk2,mp2,me2,mk2)*mk2 - 3/8*Hbb(mk2,mk2 ,mk2,mk2)*mk2 - 9/64*Hbb(mk2,me2,me2,mk2)*mk2 + 3/8*Hd(mp2,mp2 ,mk2,mk2)*mp2*mk2 + 3/16*Hd(mp2,mp2,mk2,mk2)*mk2^2 + 1/8*Hd( mp2,mk2,me2,mk2)*mk2^2 - 3/64*Hd(mk2,mp2,mp2,mk2)*mk2^2 + 9/32 *Hd(mk2,mp2,me2,mk2)*mk2^2 + 3/8*Hd(mk2,mk2,mk2,mk2)*mk2^2 + 181/576*Hd(mk2,me2,me2,mk2)*mk2^2 - 3/4*H1d(mp2,mp2,mk2,mk2)* mk2^2 - 3/4*H1d(mk2,mp2,me2,mk2)*mk2^2 - 5/8*H1d(mk2,me2,me2, mk2)*mk2^2 + 9/8*H21d(mp2,mp2,mk2,mk2)*mk2^2 - 9/64*H21d(mk2, mp2,mp2,mk2)*mk2^2 + 27/32*H21d(mk2,mp2,me2,mk2)*mk2^2 + 9/8* H21d(mk2,mk2,mk2,mk2)*mk2^2 + 27/64*H21d(mk2,me2,me2,mk2)* mk2^2 ) + Ab(mp2)*F^-4*pi16 * ( - 15/32*mk2 ) + Ab(mp2)*F^-4*L1r * ( - 24*mp2 ) + Ab(mp2)*F^-4*L2r * ( - 6*mp2 ) + Ab(mp2)*F^-4*L3r * ( - 15/2*mp2 ) + Ab(mp2)*F^-4*L4r * ( 47/2*mp2 + 19*mk2 ) + Ab(mp2)*F^-4*L5r * ( 3*mp2 + 19/2*mk2 ) + Ab(mp2)*Ab(mp2)*F^-4 * ( 129/128 + 3/32*mp2^-1*mk2 ) + Ab(mp2)*Ab(mk2)*F^-4 * ( 63/32 ) + Ab(mp2)*Ab(me2)*F^-4 * ( 57/64 ) + Ab(mk2)*F^-4*pi16 * ( - 3/8*mp2 - 9/16*mk2 ) + Ab(mk2)*F^-4*L1r * ( - 36*mk2 ) + Ab(mk2)*F^-4*L2r * ( - 18*mk2 ) + Ab(mk2)*F^-4*L3r * ( - 15*mk2 ) + Ab(mk2)*F^-4*L4r * ( 7*mp2 + 30*mk2 ) + Ab(mk2)*F^-4*L5r * ( 6*mp2 + 7*mk2 ) + Ab(mk2)*Ab(mk2)*F^-4 * ( 51/32 + 3/16*mp2*mk2^-1 ) + Ab(mk2)*Ab(me2)*F^-4 * ( 3/32 ) + Ab(me2)*F^-4*pi16 * ( - 1/8*mp2 - 11/32*mk2 ) + Ab(me2)*F^-4*L1r * ( 8/3*mp2 - 32/3*mk2 ) + Ab(me2)*F^-4*L2r * ( 2/3*mp2 - 8/3*mk2 ) + Ab(me2)*F^-4*L3r * ( 7/6*mp2 - 14/3*mk2 ) + Ab(me2)*F^-4*L4r * ( - 1/2*mp2 + 11*mk2 ) + Ab(me2)*F^-4*L5r * ( 3*mp2 + 3/2*mk2 ) + Ab(me2)*Ab(me2)*F^-4 * ( 9/32 + 9/128*mp2*me2^-1 ); * Feta (octet) Feta = F (F2x88+ F4x88+ F6x88) F2x88 = + 1; F4x88 = + L4r*Fpi^-2 * ( 4*mp2 + 8*mk2 ) + L5r*Fpi^-2 * ( - 4/3*mp2 + 16/3*mk2 ) + Ab(mk2)*Fpi^-2 * ( 3/2 ); F6x88 = + F^-2 * ( 8*mp2^2*CC14 - 8/3*mp2^2*CC15 + 24*mp2^2*CC16 + 8* mp2^2*CC17 + 64/3*mp2^2*CC18 - 64/3*mp2*mk2*CC14 + 16/3*mp2* mk2*CC15 - 32*mp2*mk2*CC16 - 64/3*mp2*mk2*CC17 - 128/3*mp2*mk2 *CC18 + 64/3*mk2^2*CC14 + 64/3*mk2^2*CC15 + 32*mk2^2*CC16 + 64/ 3*mk2^2*CC17 + 64/3*mk2^2*CC18 ) + F^-4*pi16^2 * ( 5/96*mp2^2*pi^2 + 91/384*mp2^2 + 1/48*mp2*mk2* pi^2 - 11/96*mp2*mk2 + 5/24*mk2^2*pi^2 + 49/96*mk2^2 ) + F^-4*pi16*L1r * ( - 2/9*mp2^2 + 16/9*mp2*mk2 - 32/9*mk2^2 ) + F^-4*pi16*L2r * ( - 29/9*mp2^2 + 16/9*mp2*mk2 - 68/9*mk2^2 ) + F^-4*pi16*L3r * ( - 10/9*mp2^2 + 8/9*mp2*mk2 - 19/9*mk2^2 ) + F^-4*L4r^2 * ( 56*mp2^2 + 224*mp2*mk2 + 224*mk2^2 ) + F^-4*L4r*L5r * ( 48*mp2^2 + 96*mp2*mk2 + 192*mk2^2 ) + F^-4*L4r*L6r * ( - 64*mp2^2 - 256*mp2*mk2 - 256*mk2^2 ) + F^-4*L4r*L8r * ( - 64*mp2^2 - 128*mk2^2 ) + F^-4*L5r^2 * ( - 200/9*mp2^2 + 448/9*mp2*mk2 + 256/9*mk2^2 ) + F^-4*L5r*L6r * ( 64/3*mp2^2 - 128/3*mp2*mk2 - 512/3*mk2^2 ) + F^-4*L5r*L8r * ( 64/3*mp2^2 - 256/3*mk2^2 ) + F^-4 * ( - 9/16*Hbb(mp2,mk2,mk2,me2)*mp2 + 3/16*Hbb(me2,mk2, mk2,me2)*mp2 - 3/4*Hbb(me2,mk2,mk2,me2)*mk2 + 1/12*Hd(mp2,mp2, me2,me2)*mp2^2 + 1/16*Hd(mp2,mk2,mk2,me2)*mp2^2 + 3/4*Hd(mp2, mk2,mk2,me2)*mp2*mk2 - 1/4*Hd(mk2,mk2,me2,me2)*mp2^2 + 2*Hd( mk2,mk2,me2,me2)*mp2*mk2 - 32/9*Hd(mk2,mk2,me2,me2)*mk2^2 + 1/ 16*Hd(me2,mk2,mk2,me2)*mp2^2 - 1/2*Hd(me2,mk2,mk2,me2)*mp2*mk2 + Hd(me2,mk2,mk2,me2)*mk2^2 + 49/972*Hd(me2,me2,me2,me2)* mp2^2 - 56/243*Hd(me2,me2,me2,me2)*mp2*mk2 + 64/243*Hd(me2,me2 ,me2,me2)*mk2^2 + 1/2*H1d(mp2,mk2,mk2,me2)*mp2^2 - 2*H1d(mp2, mk2,mk2,me2)*mp2*mk2 + H1d(mk2,mk2,me2,me2)*mp2^2 - 20/3*H1d( mk2,mk2,me2,me2)*mp2*mk2 + 32/3*H1d(mk2,mk2,me2,me2)*mk2^2 + 3/ 16*H21d(mp2,mk2,mk2,me2)*mp2^2 - 3/2*H21d(mp2,mk2,mk2,me2)*mp2 *mk2 + 3*H21d(mp2,mk2,mk2,me2)*mk2^2 + 3/16*H21d(me2,mk2,mk2, me2)*mp2^2 - 3/2*H21d(me2,mk2,mk2,me2)*mp2*mk2 + 3*H21d(me2, mk2,mk2,me2)*mk2^2 ) + Ab(mp2)*F^-4*L1r * ( - 24*mp2 ) + Ab(mp2)*F^-4*L2r * ( - 6*mp2 ) + Ab(mp2)*F^-4*L3r * ( - 6*mp2 ) + Ab(mp2)*F^-4*L4r * ( 22*mp2 + 16*mk2 ) + Ab(mp2)*F^-4*L5r * ( - 4/3*mp2 + 32/3*mk2 ) + Ab(mp2)*Ab(mp2)*F^-4 * ( 9/32 ) + Ab(mp2)*Ab(mk2)*F^-4 * ( 3 ) + Ab(mk2)*F^-4*pi16 * ( - 3/8*mp2 - 3/2*mk2 ) + Ab(mk2)*F^-4*L1r * ( - 32*mk2 ) + Ab(mk2)*F^-4*L2r * ( - 8*mk2 ) + Ab(mk2)*F^-4*L3r * ( - 14*mk2 ) + Ab(mk2)*F^-4*L4r * ( 10*mp2 + 36*mk2 ) + Ab(mk2)*F^-4*L5r * ( 26/3*mp2 + 8*mk2 ) + Ab(mk2)*Ab(mk2)*F^-4 * ( 9/4 + 3/8*mp2*mk2^-1 ) + Ab(me2)*F^-4*L1r * ( 4*mp2 - 16*mk2 ) + Ab(me2)*F^-4*L2r * ( 4*mp2 - 16*mk2 ) + Ab(me2)*F^-4*L3r * ( 2*mp2 - 8*mk2 ) + Ab(me2)*F^-4*L4r * ( - 2*mp2 + 8*mk2 ) + Ab(me2)*F^-4*L5r * ( - 4/3*mp2 + 16/3*mk2 ) + Ab(me2)*Ab(me2)*F^-4 * ( 9/32 );