*FORM by J.Vermaseren,version 3.1(Oct 17 2002) Run at: Mon Sep 19 15:32:14 2005 * * This FORM code processes the result produced by the (PQint) code: * * Simplifications and sums over sea and valence quark indices are * * implemented. Includes the possiblilty of producing output intended for * * fortran functions and LaTeX. * * ************************************************************************** * * comp 1 : Simplification of the pseudoscalar meson mass to NNLO * * comp 2 : Simplification of the pseudoscalar meson decay constant to NNLO * * #define comp "1" * #define nfl "3" * * * The parameters (mode) and (type) select the values of (dval) and (dsea), * * according to the tables below: * * * nfl = 2: * * dval + dsea Mode Type * ***************************** * * 1 + 1 1 1 * * 1 + 2 2 1 * * 2 + 1 4 1 * * 2 + 2 5 1 * * * nfl = 3: * * dval + dsea Mode Type * ***************************** * * 1 + 1 2 1 * * 1 + 2 4 1 * * 1 + 3 5 1 * * 2 + 1 2 2 * * 2 + 2 4 2 * * 2 + 3 5 2 * * #define mode "4" * #define type "2" * * #define pout "0" * #define nsim "0" * * #- * Sigma2 loaded * Sigma4 loaded * Sigma6 loaded * *Time = 3.27 sec Generated terms = 1 * M0p2 Terms in output = 1 * Bytes used = 24 * *Time = 3.27 sec Generated terms = 6 * M0p4 Terms in output = 6 * Bytes used = 250 * *Time = 3.27 sec Generated terms = 349 * M0p6 Terms in output = 349 * Bytes used = 14950 L M0p2timo22 = + X13 ; L M0p4timo22 = + F0^-2*L4r * ( - 24*Xb(1)*X13 ) + F0^-2*L5r * ( - 8*X13^2 ) + F0^-2*L6r * ( + 48*Xb(1)*X13 ) + F0^-2*L8r * ( + 16*X13^2 ) + Ab(Xpp)*F0^-2 * ( - 1/3*R(Xpp,Xqq,Xe0)*X13 ) + Ab(Xe0)*F0^-2 * ( - 1/3*R(Xe0,X11,X33)*X13 ); L M0p6timo22 = + F0^-4*pi16*L0r * ( - 8/9*Xe0*X13*X44 + X11*X13*X33 - 2*X13^3 + 26/3*Xb(1)*X13^2 + 35/9*Xb(2)*X13 ) + F0^-4*pi16*L1r * ( + 4*X13^3 ) + F0^-4*pi16*L2r * ( + 22/3*Xe0*X13*X44 + 2*X13^3 + 26/3*Xb(2)*X13 ) + F0^-4*pi16*L3r * ( - 8/9*Xe0*X13*X44 + 7/6*X11*X13*X33 - 11/3*X13^3 + 17/3*Xb(1)*X13^2 + 43/18*Xb(2)*X13 ) + F0^-4*pi16^2 * ( + 15/32*Xe0*X13*X44 + 59/192*X11*X13*X33 - 77/192*X13^3 + 15/32*Xb(1)*X13^2 + 43/64*Xb(2)*X13 ) + F0^-4*L4r*L5r * ( + 384*Xb(1)*X13^2 ) + F0^-4*L4r*L6r * ( - 1152*Xb(1)^2*X13 ) + F0^-4*L4r*L8r * ( - 384*Xb(1)*X13^2 ) + F0^-4*L4r^2 * ( + 576*Xb(1)^2*X13 ) + F0^-4*L5r*L6r * ( - 384*Xb(1)*X13^2 ) + F0^-4*L5r*L8r * ( - 128*X13^3 ) + F0^-4*L5r^2 * ( + 64*X13^3 ) + F0^-4*K17r * ( - 32*X13^3 ) + F0^-4*K18r * ( - 96*Xb(1)*X13^2 ) + F0^-4*K19r * ( - 8*Xpp^2*X13 ) + F0^-4*K20r * ( - 48*Xb(1)*X13^2 ) + F0^-4*K21r * ( - 48*Xb(2)*X13 ) + F0^-4*K22r * ( - 144*Xb(1)^2*X13 ) + F0^-4*K23r * ( - 16*X11*X13*X33 ) + F0^-4*K25r * ( + 24*Xpp^2*X13 ) + F0^-4*K26r * ( + 96*Xb(1)*X13^2 + 48*Xb(2)*X13 ) + F0^-4*K27r * ( + 432*Xb(1)^2*X13 ) + F0^-4*K39r * ( + 32*X13^3 ) + F0^-4*K40r * ( + 96*Xb(1)*X13^2 ) + Ab(Xpp)*F0^-4*pi16 * ( + 1/4*Xb(1)*R(Xpp,Xqq,Xe0)*X13 - 1/8*Xb(1)*RC(Xpp)*X13 + 1/12*R(Xpp,Xqq,Xe0)*X13^2 - 1/24*RC(Xpp)*Xqq*X13 ) + Ab(Xpp)*F0^-4*L0r * ( - 16/3*R(Xpp,Xqq,Xe0)*Xpp*X13 - 4/3*RC(Xpp)*Xpp*X13 - 4/3*RD(Xpp)*X13 ) + Ab(Xpp)*F0^-4*L3r * ( - 4/3*R(Xpp,Xqq,Xe0)*Xpp*X13 - 10/3*RC(Xpp)*Xpp*X13 - 10/3*RD(Xpp)*X13 ) + Ab(Xpp)*F0^-4*L4r * ( + 16*Xb(1)*R(Xpp,Xqq,Xe0)*X13 ) + Ab(Xpp)*F0^-4*L5r * ( + 16/3*R(Xpp,Xqq,Xe0)*Xpp*X13 + 8/3*RC(Xpp)*Xpp*X13 ) + Ab(Xpp)*F0^-4*L6r * ( - 16*Xb(1)*R(Xpp,Xqq,Xe0)*X13 ) + Ab(Xpp)*F0^-4*L7r * ( + 32/3*RD(Xpp)*Xpp + 16/3*RD(Xpp)*Xqq - 8/3*RD(Xpp)^2 ) + Ab(Xpp)*F0^-4*L8r * ( - 32/9*R(Xpp,Xqq,Xe0)*X11*X33 - 64/9*R(Xpp,Xqq,Xe0)*X13^2 - 16/3*RC(Xpp)*Xpp*X13 - 8/9*RD(Xpp)^2 ) + Ab(Xpp)^2*F0^-4 * ( + 1/8*X13 + 1/18*R(Xpp,Xqq,Xe0)^2*X13 + 1/54*R(Xpp,Xqq,Xe0)*RC(Xpp)*Xpp - 1/54*RC(Xpp)^2*Xpp + 1/72*RC(Xpp)^2*X13 ) + Ab(Xpp)*Ab(Xps)*F0^-4 * ( + 1/72*R(Xpp,Xqq,Xe0)*Xpp + 19/216*R(Xpp,Xqq,Xe0)*Xqq - 5/216*R(Xpp,Xss,Xe0)*Xpp - 1/24*R(Xpp,Xss,Xe0)*Xqq + 1/27*RC(Xpp)*Xpp ) + Ab(Xpp)*Ab(Xqs)*F0^-4 * ( + 1/18*R(Xpp,Xqq,Xe0)*Xqq - 1/72*R(Xpp,Xss,Xe0)*RS(Xqq,Xss,Xpp)*Xpp + 1/216*R(Xpp,Xss,Xe0)*RS(Xqq,Xss,Xpp)*Xqq + 1/12*RC(Xpp)*X13 ) + Ab(Xpp)*Ab(Xe0)*F0^-4 * ( + 1/54*R(Xpp,Xqq,Xe0)*R(Xe0,Xpp,Xpp)*Xpp - 1/90*R(Xpp,Xqq,Xe0)*R(Xe0,Xqq,Xqq)*Xqq + 13/270*R(Xpp,Xqq,Xe0)*R(Xe0,X11,X33)*Xpp + 1/18*R(Xpp,Xqq,Xe0)*R(Xe0,X11,X33)*Xqq - 1/135*R(Xpp,Xe0,Xe0)*R(Xe0,Xqq,Xqq)*Xe0 + 1/162*R(Xpp,Xss,Xe0)^2*R(Xe0,Xss,Xss)*Xss - 1/54*R(Xe0,Xqq,Xqq)*RC(Xpp)*Xpp + 1/216*RC(Xpp)*RV(Xe0,X11,X33)*Xpp + 1/72*RC(Xpp)*RV(Xe0,X11,X33)*Xqq - 1/36*RC(Xpp)*RV(Xe0,X11,X33)*Xe0 ) + Ab(Xpp)*Ab(X13)*F0^-4 * ( + 1/18*Xpp - 1/36*X13 ) + Ab(Xpp)*Ab(X44)*F0^-4 * ( - 1/18*R(Xpp,X44,Xe0)^2*X44 ) + Ab(Xpp)*Ab(X46)*F0^-4 * ( - 2/27*R(Xpp,Xe0,Xe0)*X46 + 1/9*R(Xe0,Xpp,Xpp)*Xpp ) + Ab(Xpp)*Bb(Xpp,Xpp,0)*F0^-4 * ( + 11/54*R(Xpp,Xqq,Xe0)*RC(Xpp)*Xpp*X13 - 1/27*R(Xpp,Xqq,Xe0)*RC(Xpp)*Xpp^2 - 1/18*R(Xpp,Xqq,Xe0)*RC(Xpp)*X13^2 + 1/54*R(Xpp,Xqq,Xe0)*RD(Xpp)*Xpp ) + Ab(Xpp)*Bb(Xpp,Xe0,0)*F0^-4 * ( + 2/27*R(Xe0,X11,X33)*RC(Xpp)*Xpp^2 + 1/27*R(Xe0,X11,X33)*RC(Xpp)*X11*X33 ) + Ab(Xpp)*Bb(Xqq,Xqq,0)*F0^-4 * ( - 1/54*R(Xpp,Xqq,Xe0)*RD(Xqq)*Xqq - 1/36*RC(Xpp)*RD(Xqq)*X13 ) + Ab(Xpp)*Bb(X11,X33,0)*F0^-4 * ( + 2/27*R(Xqq,Xpp,Xe0)*RC(Xpp)*Xpp*X13 ) + Ab(Xpp)*Bb(X11,X33,0,kk)*F0^-4 * ( + 1/27*R(Xqq,Xpp,Xe0)*RC(Xpp)*Xpp ) + Ab(Xpp,eps)*F0^-4*pi16 * ( - 1/4*Xb(1)*R(Xpp,Xqq,Xe0)*X13 + 1/8*Xb(1)*RC(Xpp)*X13 + 5/18*R(Xpp,Xqq,Xe0)*Xpp*X13 + 7/36*RC(Xpp)*Xpp*X13 + 5/72*RC(Xpp)*Xqq*X13 - 1/72*RD(Xpp)*X13 ) + Ab(Xps)*F0^-4*pi16 * ( - 1/4*Xpp*X13 + 1/8*X13*Xss + 1/4*X13^2 + 3/8*Xb(1)*X13 ) + Ab(Xps)*F0^-4*L0r * ( + 4*Xps*X13 ) + Ab(Xps)*F0^-4*L3r * ( + 10*Xps*X13 ) + Ab(Xps)*F0^-4*L5r * ( - 8*Xps*X13 ) + Ab(Xps)*F0^-4*L8r * ( + 16*Xps*X13 ) + Ab(Xps)*Ab(Xe0)*F0^-4 * ( + 1/27*R(Xe0,Xpp,Xpp)*Xpp - 5/216*R(Xe0,Xpp,Xss)*Xpp - 1/24*R(Xe0,Xpp,Xss)*Xqq + 1/12*R(Xe0,Xqq,Xqq)*X13 - 5/72*R(Xe0,Xqq,Xss)*Xpp - 19/216*R(Xe0,Xqq,Xss)*RS(Xpp,Xss,Xe0)*Xqq + 1/72*R(Xe0,X11,X33)*Xpp ) + Ab(Xps)*Bb(Xpp,Xe0,0)*F0^-4 * ( - 2/9*R(Xe0,Xqq,Xss)*Xpp*X13 - 2/27*R(Xe0,Xqq,Xss)*Xpp*Xss - 1/27*R(Xe0,Xqq,Xss)*Xqq*Xss - 1/27*R(Xe0,X11,X33)*Xpp^2 + 1/27*R(Xe0,X11,X33)*X11*X33 ) + Ab(Xps)*Bb(Xqq,Xqq,0)*F0^-4 * ( + 1/12*RD(Xqq)*X13 ) + Ab(Xps)*Bb(X11,X33,0)*F0^-4 * ( - 1/27*R(Xqq,Xpp,Xe0)*Xpp^2 - 4/27*R(Xqq,Xss,Xe0)*Xps*X13 - 1/27*R(Xqq,Xss,Xe0)*X11*X33 ) + Ab(Xps)*Bb(X11,X33,0,kk)*F0^-4 * ( + 1/27*R(Xqq,Xpp,Xe0)*Xqq - 2/27*R(Xqq,Xss,Xe0)*Xps - 2/27*R(Xqq,Xss,Xe0)*X13 ) + Ab(Xps,eps)*F0^-4*pi16 * ( - 3/8*X13*Xss - 5/12*X13^2 - 3/8*Xb(1)*X13 ) + Ab(Xp4)*Ab(Xq6)*F0^-4 * ( - 1/2*X13 ) + Ab(Xe0)*F0^-4*pi16 * ( - 1/8*Xb(1)*RV(Xe0,X11,X33)*X13 + 1/24*RV(Xe0,X11,X33)*Xe0*X13 - 1/12*RV(Xe0,X11,X33)*X13^2 ) + Ab(Xe0)*F0^-4*L0r * ( - 8*R(Xe0,X11,X33)*Xe0*X13 - 4/3*RV(Xe0,X11,X33)*Xe0*X13 ) + Ab(Xe0)*F0^-4*L1r * ( + 16*Xe0*X13 ) + Ab(Xe0)*F0^-4*L2r * ( + 4*Xe0*X13 ) + Ab(Xe0)*F0^-4*L3r * ( - 8*R(Xe0,X11,X33)*Xe0*X13 - 10/3*RV(Xe0,X11,X33)*Xe0*X13 ) + Ab(Xe0)*F0^-4*L4r * ( - 16*Xe0*X13 + 16*Xb(1)*R(Xe0,X11,X33)*X13 ) + Ab(Xe0)*F0^-4*L5r * ( + 8/3*R(Xe0,Xpp,Xpp)*Xe0*X13 + 16/3*R(Xe0,X11,X33)*X13^2 ) + Ab(Xe0)*F0^-4*L6r * ( + 16*Xe0*X13 - 16*Xb(1)*R(Xe0,X11,X33)*X13 ) + Ab(Xe0)*F0^-4*L7r * ( + 64/3*R(Xe0,X11,X33)*Xe0*X13 - 64/3*R(Xe0,X11,X33)*X13^2 - 4/3*RD(Xe0,X44)^2*RV(Xe0,X11,X33) + 16/3*RD(Xe0,X44,X66,Xpp)*Xqq ) + Ab(Xe0)*F0^-4*L8r * ( - 64/9*R(Xe0,X11,X33)*X11*X33 - 128/9*R(Xe0,X11,X33)*X13^2 - 16/9*RV(Xe0,X11,X33)*Xe0^2 - 8/3*RV(Xe0,X11,X33)*X11*X33 - 8/27*RV(Xe0,X11,X33)*X44^2 - 16/27*RV(Xe0,X11,X33)*X66^2 ) + Ab(Xe0)^2*F0^-4 * ( + 1/54*Xb(1)*RV(Xe0,X11,X33) + 1/18*R(Xe0,X11,X33)^2*X13 - 1/54*R(Xe0,X11,X33)*RV(Xe0,X11,X33)*Xe0 + 1/54*R(Xe0,X11,X33)*RV(Xe0,X11,X33)*X13 - 5/108*RV(Xe0,X11,X33)*Xe0 - 1/54*RV(Xe0,X11,X33)^2*Xe0 + 1/72*RV(Xe0,X11,X33)^2*X13 ) + Ab(Xe0)*Ab(X44)*F0^-4 * ( - 1/36*RV(Xe0,X11,X33)*X44 ) + Ab(Xe0)*Ab(X46)*F0^-4 * ( + 1/9*R(Xe0,X11,X33)*X13 + 2/27*RV(Xe0,X11,X33)*X46 ) + Ab(Xe0)*Bb(Xpp,Xpp,0)*F0^-4 * ( + 1/9*R(Xpp,Xqq,Xe0)*R(Xe0,Xpp,Xpp)*Xpp^2 - 1/27*R(Xe0,Xpp,Xpp)*RD(Xpp)*Xpp - 1/54*R(Xe0,Xqq,Xqq)*RD(Xpp)*Xpp + 1/216*RD(Xpp)*RV(Xe0,X11,X33)*Xpp + 1/72*RD(Xpp)*RV(Xe0,X11,X33)*Xqq - 1/36*RD(Xpp)*RV(Xe0,X11,X33)*Xe0 ) + Ab(Xe0)*Bb(Xpp,Xe0,0)*F0^-4 * ( - 2/27*R(Xpp,Xe0,Xe0)*R(Xe0,Xqq,X44)*Xpp*Xe0 - 1/27*R(Xpp,Xe0,Xe0)*R(Xe0,Xqq,X44)*Xqq*Xe0 + 2/27*R(Xpp,Xe0,Xe0)*R(Xe0,X11,X33)*Xpp^2 + 1/27*R(Xpp,Xe0,Xe0)*R(Xe0,X11,X33)*X11*X33 ) + Ab(Xe0)*Bb(Xe0,Xe0,0)*F0^-4 * ( + 1/9*Xb(1)*R(Xe0,X11,X33)*X13 - 5/18*R(Xe0,X11,X33)*Xe0*X13 ) + Ab(Xe0)*Bb(X11,X33,0)*F0^-4 * ( + 2/27*R(X11,Xe0,Xe0)*R(X33,Xe0,Xe0)*Xe0*X13 ) + Ab(Xe0)*Bb(X11,X33,0,kk)*F0^-4 * ( + 1/27*R(X11,Xe0,Xe0)*R(X33,Xe0,Xe0)*Xe0 ) + Ab(Xe0,eps)*F0^-4*pi16 * ( - 1/4*Xe0*X13 + 1/8*Xb(1)*RV(Xe0,X11,X33)*X13 + 1/4*R(Xe0,X11,X33)*Xe0*X13 + 5/9*R(Xe0,X11,X33)*X13^2 + 1/8*RV(Xe0,X11,X33)*Xe0*X13 + 5/36*RV(Xe0,X11,X33)*X13^2 ) + Ab(X11)*Ab(X33)*F0^-4 * ( - 1/54*R(Xpp,Xqq,Xe0)*RC(Xqq)*Xqq + 1/9*R(X11,X33,Xe0)*R(X33,X11,Xe0)*X13 - 1/36*RC(X11)*RC(X33)*X13 ) + Ab(X13)*F0^-4*L1r * ( + 8*X13^2 ) + Ab(X13)*F0^-4*L2r * ( + 20*X13^2 ) + Ab(X13)*F0^-4*L6r * ( + 32*X13^2 ) + Ab(X13)^2*F0^-4 * ( - 1/4*X13 ) + Ab(X13)*Bb(X11,X33,0)*F0^-4 * ( + 1/9*X11*X33 + 4/9*X13^2 ) + Ab(X13)*Bb(X11,X33,0,kk)*F0^-4 * ( + 4/9*X13 ) + Ab(X13,eps)*F0^-4*pi16 * ( - 31/18*X13^2 ) + Ab(X14)*Ab(X34)*F0^-4 * ( - X13 ) + Ab(X16)*Ab(X36)*F0^-4 * ( - 1/4*X13 ) + Ab(X44)*F0^-4*L1r * ( + 48*X13*X44 ) + Ab(X44)*F0^-4*L2r * ( + 12*X13*X44 ) + Ab(X44)*F0^-4*L4r * ( - 48*X13*X44 ) + Ab(X44)*F0^-4*L6r * ( + 48*X13*X44 ) + Ab(X44)*Bb(Xpp,Xe0,0)*F0^-4 * ( + 2/9*R(Xpp,X44,Xe0)*R(Xe0,Xqq,X44)*Xpp*X44 + 1/9*R(Xpp,X44,Xe0)*R(Xe0,Xqq,X44)*Xqq*X44 ) + Ab(X44)*Bb(Xe0,Xe0,0)*F0^-4 * ( - 1/6*R(Xe0,X11,X33)*X13*X44 ) + Ab(X44)*Bb(X11,X33,0)*F0^-4 * ( + 2/9*R(X11,X44,Xe0)*R(X33,X44,Xe0)*X13*X44 ) + Ab(X44)*Bb(X11,X33,0,kk)*F0^-4 * ( + 1/9*R(X11,X44,Xe0)*R(X33,X44,Xe0)*X44 ) + Ab(X44,eps)*F0^-4*pi16 * ( - 3/4*X13*X44 ) + Ab(X46)*F0^-4*L1r * ( + 64*X13*X46 ) + Ab(X46)*F0^-4*L2r * ( + 16*X13*X46 ) + Ab(X46)*F0^-4*L4r * ( - 64*X13*X46 ) + Ab(X46)*F0^-4*L6r * ( + 64*X13*X46 ) + Ab(X46)*Bb(Xpp,Xe0,0)*F0^-4 * ( + 2/81*RU(Xpp,Xe0,Xpp,X44,X66)*Xpp*Xe0 - 2/81*RU(Xpp,Xe0,Xqq,X44,X66)*Xpp*Xe0 - 2/81*RU(Xpp,Xe0,Xqq,X44,X66)*X11*X33 + 4/81*RW(Xpp,Xe0,Xqq,X44,X66)*Xpp + 2/81*RW(Xpp,Xe0,Xqq,X44,X66)*Xqq ) + Ab(X46)*Bb(Xpp,Xe0,0,kk)*F0^-4 * ( - 2/81*RU(Xpp,Xe0,Xpp,X44,X66)*Xpp - 2/81*RU(Xpp,Xe0,Xqq,X44,X66)*Xpp ) + Ab(X46)*Bb(Xe0,Xe0,0)*F0^-4 * ( + 4/9*R(Xe0,X11,X33)*X13*X46 ) + Ab(X46)*Bb(X11,X33,0)*F0^-4 * ( - 2/81*RU(X11,X33,X44,X66)*X11*X33 + 4/81*RW(X11,X33,X44,X66)*X13 ) + Ab(X46)*Bb(X11,X33,0,kk)*F0^-4 * ( - 4/81*RU(X11,X33,X44,X66)*X13 + 2/81*RW(X11,X33,X44,X66) ) + Ab(X46,eps)*F0^-4*pi16 * ( - X13*X46 ) + Bb(Xpp,Xpp,0)*F0^-4*pi16 * ( - 1/8*Xb(1)*RD(Xpp)*X13 - 1/24*RD(Xpp)*Xqq*X13 ) + Bb(Xpp,Xpp,0)*F0^-4*L0r * ( - 4/3*RD(Xpp)*Xpp*X13 ) + Bb(Xpp,Xpp,0)*F0^-4*L3r * ( - 10/3*RD(Xpp)*Xpp*X13 ) + Bb(Xpp,Xpp,0)*F0^-4*L4r * ( + 8*Xb(1)*R(Xpp,Xqq,Xe0)*Xpp*X13 ) + Bb(Xpp,Xpp,0)*F0^-4*L5r * ( + 16/3*R(Xpp,Xqq,Xe0)*Xpp^2*X13 - 8/3*R(Xpp,Xqq,Xe0)*X11*X13*X33 ) + Bb(Xpp,Xpp,0)*F0^-4*L6r * ( - 16*Xb(1)*R(Xpp,Xqq,Xe0)*Xpp*X13 ) + Bb(Xpp,Xpp,0)*F0^-4*L8r * ( - 32/3*R(Xpp,Xqq,Xe0)*Xpp^2*X13 + 16/3*R(Xpp,Xqq,Xe0)*X11*X13*X33 ) + Bb(Xpp,Xpp,0)^2*F0^-4 * ( + 1/8*R(Xpp,Xqq,Xe0)*RD(Xpp)*Xpp*X13 - 1/72*R(Xpp,Xqq,Xe0)*RD(Xpp)*Xqq*X13 ) + Bb(Xpp,Xpp,0)*Bb(Xpp,Xe0,0)*F0^-4 * ( + 2/27*R(Xe0,X11,X33)*RD(Xpp)*Xpp^2 + 1/27*R(Xe0,X11,X33)*RD(Xpp)*X11*X33 ) + Bb(Xpp,Xpp,0)*Bb(X11,X33,0)*F0^-4 * ( + 2/27*R(Xqq,Xpp,Xe0)*RD(Xpp)*Xpp*X13 ) + Bb(Xpp,Xpp,0)*Bb(X11,X33,0,kk)*F0^-4 * ( + 1/27*R(Xqq,Xpp,Xe0)*RD(Xpp)*Xpp ) + Bb(Xpp,Xpp,0,eps)*F0^-4*pi16 * ( + 1/8*Xb(1)*RD(Xpp)*X13 + 7/36*RD(Xpp)*Xpp*X13 + 5/72*RD(Xpp)*Xqq*X13 ) + Bb(Xpp,Xe0,0)*F0^-4*L7r * ( + 32/3*RD(Xpp)*RD(Xe0,X44,X66,Xqq)*Xpp + 16/3*RD(Xpp)*RD(Xe0,X44,X66,Xqq)*Xqq ) + Bb(Xpp,Xe0,0)*F0^-4*L8r * ( + 32/9*RD(Xpp)*RD(Xe0,X44,X66,Xqq)*Xpp + 16/9*RD(Xpp)*RD(Xe0,X44,X66,Xqq)*Xqq ) + Bb(Xe0,Xe0,0)*F0^-4*L4r * ( + 8*Xb(1)*R(Xe0,X11,X33)*Xe0*X13 ) + Bb(Xe0,Xe0,0)*F0^-4*L5r * ( + 8/3*R(Xe0,X11,X33)*Xe0^2*X13 ) + Bb(Xe0,Xe0,0)*F0^-4*L6r * ( - 16*Xb(1)*R(Xe0,X11,X33)*Xe0*X13 ) + Bb(Xe0,Xe0,0)*F0^-4*L7r * ( - 8*R(Xe0,X11,X33)*RD(Xe0,X44)^2*X13 ) + Bb(Xe0,Xe0,0)*F0^-4*L8r * ( - 16/9*R(Xe0,X11,X33)*X13*X44^2 - 32/9*R(Xe0,X11,X33)*X13*X66^2 ) + Bb(X11,X11,0)*Bb(X33,X33,0)*F0^-4 * ( - 1/36*RD(X11)*RD(X33)*X13 ) + Bb(X11,X33,0)*F0^-4*L7r * ( + 32/3*RD(X11)*RD(X33)*X13 ) + Bb(X11,X33,0)*F0^-4*L8r * ( + 32/9*RD(X11)*RD(X33)*X13 ) + Bb(X11,X33,0,kk)*F0^-4*L7r * ( + 16/3*RD(X11)*RD(X33) ) + Bb(X11,X33,0,kk)*F0^-4*L8r * ( + 16/9*RD(X11)*RD(X33) ) + Hbb(1,Xpp,Xpp,X13,X13)*F0^-4 * ( - 1/4*Xpp*X13 - 3/8*X13^2 - 1/18*R(Xpp,Xqq,Xe0)^2*X13^2 + 1/6*R(Xpp,Xqq,Xe0)*RC(Xpp)*X13^2 - 1/24*RC(Xpp)^2*X13^2 ) + Hbb(1,Xpp,X1s,X3s,X13)*F0^-4 * ( - 5/12*R(Xpp,Xqq,Xe0)*Xqq*X13 + 1/6*R(Xpp,Xqq,Xe0)*Xps*X13 + 1/8*R(Xpp,Xss,Xe0)*Xqq*X13 + 7/24*R(Xpp,Xss,Xe0)*X13*Xss - 1/12*RC(Xpp)*Xps*X13 ) + Hbb(1,Xp4,Xq6,X46,X13)*F0^-4 * ( + 1/2*X13*X46 ) + Hbb(1,Xe0,Xpp,X13,X13)*F0^-4 * ( + 2/9*R(Xpp,Xqq,Xe0)*R(Xe0,X11,X33)*X13^2 - 1/18*R(Xpp,Xqq,Xe0)*RV(Xe0,X11,X33)*X13^2 + 1/9*R(Xe0,Xpp,Xpp)*RC(Xpp)*X13^2 + 1/36*RC(Xpp)*RV(Xe0,X11,X33)*X13^2 ) + Hbb(1,Xe0,Xe0,X13,X13)*F0^-4 * ( + 2/9*R(Xe0,X11,X33)^2*X13^2 + 2/9*R(Xe0,X11,X33)*RV(Xe0,X11,X33)*X13^2 + 5/72*RV(Xe0,X11,X33)^2*X13^2 ) + Hbb(1,Xe0,X1s,X3s,X13)*F0^-4 * ( - 1/8*R(Xe0,Xpp,Xss)*Xpp*X13 - 5/24*R(Xe0,Xpp,Xss)*Xqq*X13 + 1/3*R(Xe0,X11,X33)*X13^2 - 1/24*RV(Xe0,Xpp,Xss)*Xe0*X13 - 1/24*RV(Xe0,X11,X33)*X13*Xss ) + Hbb(1,X11,X13,X33,X13)*F0^-4 * ( - 1/18*R(Xpp,Xqq,Xe0)*RC(Xqq)*X13^2 + 1/3*R(X11,X33,Xe0)*R(X33,X11,Xe0)*X13^2 + 1/36*RC(X11)*RC(X33)*X13^2 ) + Hbb(1,X13,X13,X13,X13)*F0^-4 * ( + 1/4*X13^2 ) + Hbb(1,X14,X34,X44,X13)*F0^-4 * ( + 3/4*X13*X44 ) + Hbb(2,Xpp,Xpp,X13,X13)*F0^-4 * ( - 1/18*R(Xpp,Xqq,Xe0)*RD(Xpp)*X13^2 + 5/36*RC(Xpp)*RD(Xpp)*X13^2 ) + Hbb(2,Xpp,Xe0,X13,X13)*F0^-4 * ( + 1/9*R(Xe0,Xpp,Xpp)*RD(Xpp)*X13^2 + 1/36*RD(Xpp)*RV(Xe0,X11,X33)*X13^2 ) + Hbb(2,Xpp,X13,Xqq,X13)*F0^-4 * ( - 1/18*R(Xqq,Xpp,Xe0)*RD(Xpp)*X13^2 + 1/36*RC(Xqq)*RD(Xpp)*X13^2 ) + Hbb(2,Xpp,X1s,X3s,X13)*F0^-4 * ( - 1/12*RD(Xpp)*Xps*X13 ) + Hbb(5,Xpp,Xpp,X13,X13)*F0^-4 * ( + 5/72*RD(Xpp)^2*X13^2 ) + Hbb(5,X11,X33,X13,X13)*F0^-4 * ( + 1/36*RD(X11)*RD(X33)*X13^2 ) + H1bb(1,Xpp,Xpp,X13,X13)*F0^-4 * ( + 2*X13^2 + 4/9*R(Xpp,Xqq,Xe0)^2*X13^2 - 4/9*R(Xpp,Xqq,Xe0)*RC(Xpp)*X13^2 + 2/9*RC(Xpp)^2*X13^2 ) + H1bb(1,Xpp,X1s,X3s,X13)*F0^-4 * ( + 2/3*R(Xpp,Xqq,Xe0)*RS(Xss,Xqq,Xpp)*X13^2 ) + H1bb(1,Xps,Xqs,Xpp,X13)*F0^-4 * ( - 2/3*R(Xpp,Xss,Xe0)*X13^2 ) + H1bb(1,Xps,Xqs,Xe0,X13)*F0^-4 * ( - 2/3*R(Xe0,X11,X33)*RS(Xss,Xpp,Xe0)*X13^2 ) + H1bb(1,X13,Xpp,Xe0,X13)*F0^-4 * ( + 2/9*R(Xpp,Xqq,Xe0)*RV(Xe0,X11,X33)*X13^2 - 2/9*R(Xe0,Xpp,Xpp)*RS(Xqq,Xpp,Xe0)*RC(Xpp)*X13^2 ) + H1bb(1,X13,Xe0,Xe0,X13)*F0^-4 * ( - 2/9*R(Xe0,X11,X33)*RV(Xe0,X11,X33)*X13^2 - 1/9*RV(Xe0,X11,X33)^2*X13^2 ) + H1bb(1,X13,X11,X33,X13)*F0^-4 * ( + 2/9*R(Xpp,Xqq,Xe0)*RC(Xqq)*X13^2 - 4/9*R(X11,X33,Xe0)*R(X33,X11,Xe0)*X13^2 ) + H1bb(3,X13,Xpp,Xpp,X13)*F0^-4 * ( + 2/9*R(Xpp,Xqq,Xe0)*RD(Xpp)*X13^2 - 2/9*RC(Xpp)*RD(Xpp)*X13^2 ) + H1bb(3,X13,Xpp,Xqq,X13)*F0^-4 * ( + 2/9*R(Xqq,Xpp,Xe0)*RD(Xpp)*X13^2 ) + H1bb(3,X13,Xpp,Xe0,X13)*F0^-4 * ( + 2/9*R(Xe0,X11,X33)*RS(Xpp,Xqq,Xe0)*RD(Xpp)*X13^2 ) + H1bb(7,X13,Xpp,Xpp,X13)*F0^-4 * ( - 1/9*RD(Xpp)^2*X13^2 ) + H21bb(1,Xpp,Xpp,X13,X13)*F0^-4 * ( - 3/4*X13^2 ) + H21bb(1,Xpp,X1s,X3s,X13)*F0^-4 * ( - 1/4*R(Xpp,Xqq,Xe0)*RS(Xss,Xqq,Xpp)*X13^2 ) + H21bb(1,Xps,Xqs,Xpp,X13)*F0^-4 * ( + 1/4*R(Xpp,Xqq,Xe0)*X13^2 + 1/4*R(Xpp,Xss,Xe0)*X13^2 - 1/4*RC(Xpp)*X13^2 ) + H21bb(1,Xps,Xqs,Xqq,X13)*F0^-4 * ( + 1/4*R(Xqq,Xpp,Xe0)*RS(Xss,Xpp,Xqq)*X13^2 ) + H21bb(1,Xps,Xqs,Xe0,X13)*F0^-4 * ( + 1/4*R(Xe0,X11,X33)*RS(Xpp,Xqq,Xe0)*RS(Xss,Xpp,Xe0)*X13^2 ) + H21bb(1,Xe0,X1s,X3s,X13)*F0^-4 * ( - 1/4*R(Xe0,X11,X33)*RS(Xss,X11,Xe0)*RS(Xss,X33,Xe0)*X13^2 ) + H21bb(1,X13,Xpp,Xpp,X13)*F0^-4 * ( + 3/8*X13^2 + 1/6*R(Xpp,Xqq,Xe0)^2*X13^2 - 1/6*R(Xpp,Xqq,Xe0)*RC(Xpp)*X13^2 + 1/24*RC(Xpp)^2*X13^2 ) + H21bb(1,X13,Xpp,Xe0,X13)*F0^-4 * ( - 1/6*R(Xpp,Xqq,Xe0)*RV(Xe0,X11,X33)*X13^2 + 1/12*RC(Xpp)*RV(Xe0,X11,X33)*X13^2 ) + H21bb(1,X13,Xe0,Xe0,X13)*F0^-4 * ( + 1/24*RV(Xe0,X11,X33)^2*X13^2 ) + H21bb(1,X13,X11,X33,X13)*F0^-4 * ( - 1/6*R(Xpp,Xqq,Xe0)*RC(Xqq)*X13^2 + 1/3*R(X11,X33,Xe0)*R(X33,X11,Xe0)*X13^2 + 1/12*RC(X11)*RC(X33)*X13^2 ) + H21bb(1,X13,X13,X13,X13)*F0^-4 * ( + 3/4*X13^2 ) + H21bb(1,X44,X14,X34,X13)*F0^-4 * ( + 9/4*X13^2 ) + H21bb(1,X46,Xp4,Xq6,X13)*F0^-4 * ( + 3/2*X13^2 ) + H21bb(3,Xps,Xpp,Xqs,X13)*F0^-4 * ( - 1/4*RD(Xpp)*X13^2 ) + H21bb(3,X13,Xpp,Xpp,X13)*F0^-4 * ( - 1/6*R(Xpp,Xqq,Xe0)*RD(Xpp)*X13^2 + 1/12*RC(Xpp)*RD(Xpp)*X13^2 ) + H21bb(3,X13,Xpp,Xqq,X13)*F0^-4 * ( - 1/6*R(Xqq,Xpp,Xe0)*RD(Xpp)*X13^2 + 1/12*RC(Xqq)*RD(Xpp)*X13^2 ) + H21bb(3,X13,Xpp,Xe0,X13)*F0^-4 * ( + 1/12*RD(Xpp)*RV(Xe0,X11,X33)*X13^2 ) + H21bb(7,X13,Xpp,Xpp,X13)*F0^-4 * ( + 1/24*RD(Xpp)^2*X13^2 ) + H21bb(7,X13,X11,X33,X13)*F0^-4 * ( + 1/12*RD(X11)*RD(X33)*X13^2 );