*FORM by J.Vermaseren,version 3.1(Oct 17 2002) Run at: Tue Dec 7 19:24:49 2004 * * This FORM code processes the result produced by the (PQint) code and * * implements the simplifications defined in the (decsimpl) and (massimpl) * * procedures, as well as sums over valence and sea quark indices. Includes * * the option of producing raw output intended for fortran functions and * * LaTeX. These have to be further processed with perl. * * ************************************************************************** * * The parameters (comp), (mode) and (type) define the following things: * * Type of calculation, degree of quark mass degeneracy, and flavor indices * * of the external legs (ijkl). These parameters are used to find the * * relevant files to be loaded. The parameter (pout) selects whether output * * for fortran and LaTeX is to be produced. * * * Comp 1 : Pseudoscalar meson mass (M_ps^2) to NNLO * * Comp 2 : Pseudoscalar meson decay constant F_ps to NNLO * * * Mode 1 : full 3+3 theory, all quarks nondegenerate * * Mode 2 : SU(2) + 1 theory, all sea quarks degenerate * * Mode 3 : unquenched SU(2) symmetric limit * * Mode 4 : SU(2) + 2 theory, sea quarks 4 and 6 nondegenerate * * Mode 5 : SU(2) + 3 theory, all sea quarks nondegenerate * * * Type 1 : pi^+/- meson - 1221 * * Type 2 : K^+/- meson - 1331 * * #define comp "2" * #define mode "4" * #define type "2" * * #define pout "0" * * #- * Sigma2 loaded * Sigma4 loaded * Sigma6 loaded * *Time = 1.57 sec Generated terms = 1 * F0p2 Terms in output = 1 * Bytes used = 24 * *Time = 1.57 sec Generated terms = 7 * F0p4 Terms in output = 7 * Bytes used = 256 * *Time = 1.58 sec Generated terms = 338 * F0p6 Terms in output = 338 * Bytes used = 14574 L F0p2timo22 = + F0 * ( + 1 ); L F0p4timo22 = + F0^-1*L4r * ( + 12*Xb(1) ) + F0^-1*L5r * ( + 4*X13 ) + Ab(Xpp)*F0^-1 * ( + 1/6*R(Xpp,Xqq,Xe0) - 1/12*RC(Xpp) ) + Ab(Xps)*F0^-1 * ( + 1/4 ) + Ab(Xe0)*F0^-1 * ( - 1/12*RV(Xe0,X11,X33) ) + Bb(Xpp,Xpp,0)*F0^-1 * ( - 1/12*RD(Xpp) ); L F0p6timo22 = + F0^-3*pi16*L0r * ( + 4/9*Xe0*X44 - 1/2*X11*X33 + X13^2 - 13/3*Xb(1)*X13 - 35/18*Xb(2) ) + F0^-3*pi16*L1r * ( - 2*X13^2 ) + F0^-3*pi16*L2r * ( - 11/3*Xe0*X44 - X13^2 - 13/3*Xb(2) ) + F0^-3*pi16*L3r * ( + 4/9*Xe0*X44 - 7/12*X11*X33 + 11/6*X13^2 - 17/6*Xb(1)*X13 - 43/36*Xb(2) ) + F0^-3*pi16^2 * ( - 15/64*Xe0*X44 - 59/384*X11*X33 + 65/384*X13^2 - 1/2*Xb(1)*X13 - 43/128*Xb(2) ) + F0^-3*L4r*L5r * ( - 48*Xb(1)*X13 ) + F0^-3*L4r^2 * ( - 72*Xb(1)^2 ) + F0^-3*L5r^2 * ( - 8*X13^2 ) + F0^-3*K19r * ( + 4*Xpp^2 ) + F0^-3*K20r * ( + 24*Xb(1)*X13 ) + F0^-3*K21r * ( + 24*Xb(2) ) + F0^-3*K22r * ( + 72*Xb(1)^2 ) + F0^-3*K23r * ( + 8*X11*X33 ) + Ab(Xpp)*F0^-3*pi16 * ( - 1/24*Xpp + 1/48*Xb(1) - 1/8*Xb(1)*R(Xpp,Xqq,Xe0) + 1/16*Xb(1)*RC(Xpp) - 1/48*R(Xpp,Xqq,Xe0)*Xpp - 1/16*R(Xpp,Xqq,Xe0)*Xqq + 1/48*R(Xe0,Xpp,Xpp)*Xe0 + 1/16*RC(Xpp)*X13 ) + Ab(Xpp)*F0^-3*L0r * ( + 8/3*R(Xpp,Xqq,Xe0)*Xpp + 2/3*RC(Xpp)*Xpp + 2/3*RD(Xpp) ) + Ab(Xpp)*F0^-3*L3r * ( + 2/3*R(Xpp,Xqq,Xe0)*Xpp + 5/3*RC(Xpp)*Xpp + 5/3*RD(Xpp) ) + Ab(Xpp)*F0^-3*L4r * ( - 2*Xb(1)*Xb(Xpp,Xe0,Xpp,Xe0,0) - 2*Xb(1)*R(Xpp,Xqq,Xe0) + 3*Xb(1)*RC(Xpp) ) + Ab(Xpp)*F0^-3*L5r * ( - 2/3*Xb(Xpp,Xe0,Xpp,Xe0,1) - R(Xpp,Xqq,Xe0)*Xpp + 1/3*R(Xpp,Xqq,Xe0)*Xqq + 1/2*RC(Xpp)*Xpp - 1/6*RC(Xpp)*Xqq ) + Ab(Xpp)^2*F0^-3 * ( + 1/16 + 1/72*R(Xpp,Xqq,Xe0)^2 - 1/72*R(Xpp,Xqq,Xe0)*RC(Xpp) + 1/288*RC(Xpp)^2 ) + Ab(Xpp)*Ab(Xps)*F0^-3 * ( - 1/36*R(Xpp,Xqq,Xe0) - 5/72*R(Xpp,Xss,Xe0) + 7/144*RC(Xpp) ) + Ab(Xpp)*Ab(Xqs)*F0^-3 * ( - 1/36*R(Xpp,Xqq,Xe0) - 1/24*R(Xpp,Xss,Xe0) - 1/48*RC(Xpp) ) + Ab(Xpp)*Ab(Xe0)*F0^-3 * ( - 1/72*R(Xpp,Xqq,Xe0)*RV(Xe0,X11,X33) + 1/144*RC(Xpp)*RV(Xe0,X11,X33) ) + Ab(Xpp)*Ab(X13)*F0^-3 * ( + 1/8 ) + Ab(Xpp)*Ab(X46)*F0^-3 * ( + 1/12*R(Xe0,Xpp,Xpp) ) + Ab(Xpp)*Bb(Xpp,Xpp,0)*F0^-3 * ( + 1/4*Xpp - 1/18*R(Xpp,Xqq,Xe0)*RC(Xpp)*Xpp - 1/72*R(Xpp,Xqq,Xe0)*RD(Xpp) + 1/18*RC(Xpp)^2*Xpp + 1/144*RC(Xpp)*RD(Xpp) ) + Ab(Xpp)*Bb(Xpp,Xe0,0)*F0^-3 * ( + 1/18*R(Xe0,Xpp,Xpp)*RC(Xpp)*Xpp - 1/18*R(Xe0,X11,X33)*RC(Xpp)*Xpp ) + Ab(Xpp)*Bb(Xqq,Xqq,0)*F0^-3 * ( - 1/72*R(Xpp,Xqq,Xe0)*RD(Xqq) + 1/144*RC(Xpp)*RD(Xqq) ) + Ab(Xpp)*Bb(Xps,Xps,0)*F0^-3 * ( - 1/12*R(Xpp,Xss,Xe0)*Xps ) + Ab(Xpp)*Bb(X11,X33,0)*F0^-3 * ( - 1/18*R(Xqq,Xpp,Xe0)*RC(Xpp)*Xpp ) + Ab(Xpp)*Cb(Xpp,Xpp,Xpp,0)*F0^-3 * ( + 1/18*RC(Xpp)*RD(Xpp)*Xpp ) + Ab(Xpp,eps)*F0^-3*pi16 * ( + 1/8*Xb(1)*R(Xpp,Xqq,Xe0) - 1/16*Xb(1)*RC(Xpp) - 1/16*RC(Xpp)*Xpp - 1/16*RD(Xpp) ) + Ab(Xps)*F0^-3*pi16 * ( + 1/16*Xps - 3/16*Xqs - 3/16*Xb(1) ) + Ab(Xps)*F0^-3*L0r * ( - 2*Xps ) + Ab(Xps)*F0^-3*L3r * ( - 5*Xps ) + Ab(Xps)*F0^-3*L4r * ( - 3*Xb(1) ) + Ab(Xps)*F0^-3*L5r * ( + X13 ) + Ab(Xps)*Ab(Xe0)*F0^-3 * ( + 7/144*R(Xe0,Xpp,Xpp) - 5/72*R(Xe0,Xpp,Xss) - 1/48*R(Xe0,Xqq,Xqq) + 5/72*R(Xe0,Xqq,Xss) - 1/36*R(Xe0,X11,X33) ) + Ab(Xps)*Bb(Xpp,Xpp,0)*F0^-3 * ( + 1/24*R(Xpp,Xss,Xe0)*Xpp - 5/24*R(Xpp,Xss,Xe0)*Xps ) + Ab(Xps)*Bb(Xpp,Xe0,0)*F0^-3 * ( - 1/18*R(Xe0,Xpp,Xss)*RS(Xqq,Xpp,Xe0)*Xpp - 1/9*R(Xe0,Xpp,Xss)*RS(Xqq,Xpp,Xe0)*Xps ) + Ab(Xps)*Bb(Xqq,Xqq,0)*F0^-3 * ( - 1/48*RD(Xqq) ) + Ab(Xps)*Bb(X11,X33,0)*F0^-3 * ( + 1/18*R(Xqq,Xss,Xe0)*Xss ) + Ab(Xps)*Bb(X11,X33,0,kk)*F0^-3 * ( + 1/9*R(Xqq,Xss,Xe0) ) + Ab(Xps,eps)*F0^-3*pi16 * ( + 3/16*Xss + 3/16*Xb(1) ) + Ab(Xp4)^2*F0^-3 * ( - 1/8 ) + Ab(Xp4)*Ab(Xp6)*F0^-3 * ( - 1/8 ) + Ab(Xp4)*Ab(Xq6)*F0^-3 * ( + 1/8 ) + Ab(Xp6)^2*F0^-3 * ( - 1/32 ) + Ab(Xe0)*F0^-3*pi16 * ( + 1/16*Xb(1)*RV(Xe0,X11,X33) - 1/48*RV(Xe0,X11,X33)*Xe0 + 1/16*RV(Xe0,X11,X33)*X13 ) + Ab(Xe0)*F0^-3*L0r * ( + 4*R(Xe0,X11,X33)*Xe0 + 2/3*RV(Xe0,X11,X33)*Xe0 ) + Ab(Xe0)*F0^-3*L1r * ( - 8*Xe0 ) + Ab(Xe0)*F0^-3*L2r * ( - 2*Xe0 ) + Ab(Xe0)*F0^-3*L3r * ( + 4*R(Xe0,X11,X33)*Xe0 + 5/3*RV(Xe0,X11,X33)*Xe0 ) + Ab(Xe0)*F0^-3*L4r * ( + 4*Xe0 + Xb(1)*RV(Xe0,X11,X33) ) + Ab(Xe0)*F0^-3*L5r * ( - 1/6*R(Xe0,Xpp,Xpp)*Xqq - R(Xe0,X11,X33)*X13 - 1/6*RV(Xe0,X11,X33)*Xe0 ) + Ab(Xe0)^2*F0^-3 * ( + 1/288*RV(Xe0,X11,X33)^2 ) + Ab(Xe0)*Ab(X46)*F0^-3 * ( + 1/12*RV(Xe0,X11,X33) ) + Ab(Xe0)*Bb(Xpp,Xpp,0)*F0^-3 * ( - 1/36*Xb(Xpp,Xe0,Xpp,Xe0,Xe0,1) - 1/18*R(Xpp,Xqq,Xe0)*R(Xe0,Xpp,Xpp)*Xpp + 1/18*R(Xe0,Xpp,Xpp)*RC(Xpp)*Xpp + 1/144*RD(Xpp)*RV(Xe0,X11,X33) ) + Ab(Xe0)*Bb(Xpp,Xe0,0)*F0^-3 * ( - 1/18*Xb(Xe0,Xpp,Xpp,Xe0,Xe0,1) + 1/18*Xb(Xe0,Xqq,Xpp,Xe0,Xe0,1) + 1/18*R(Xe0,Xpp,Xpp)^2*RS(Xqq,Xpp,Xe0)*Xpp ) + Ab(Xe0)*Bb(Xps,Xps,0)*F0^-3 * ( - 1/12*R(Xe0,Xpp,Xss)*Xps ) + Ab(Xe0)*Bb(Xe0,Xe0,0)*F0^-3 * ( - 1/216*RV(Xe0,X11,X33)*X44 - 1/27*RV(Xe0,X11,X33)*X66 ) + Ab(Xe0)*Bb(X11,X33,0)*F0^-3 * ( - 1/18*R(X11,Xe0,Xe0)*R(X33,Xe0,Xe0)*Xe0 ) + Ab(Xe0)*Cb(Xpp,Xpp,Xpp,0)*F0^-3 * ( + 1/18*R(Xe0,Xpp,Xpp)*RD(Xpp)*Xpp ) + Ab(Xe0,eps)*F0^-3*pi16 * ( + 1/8*Xe0 - 1/16*Xb(1)*RV(Xe0,X11,X33) - 1/8*R(Xe0,X11,X33)*Xe0 - 1/16*RV(Xe0,X11,X33)*Xe0 ) + Ab(X11)*Ab(X33)*F0^-3 * ( - 1/72*R(Xpp,Xqq,Xe0)*RC(Xqq) + 1/36*R(X11,X33,Xe0)*R(X33,X11,Xe0) + 1/144*RC(X11)*RC(X33) ) + Ab(X13)*F0^-3*L1r * ( - 4*X13 ) + Ab(X13)*F0^-3*L2r * ( - 10*X13 ) + Ab(X13)^2*F0^-3 * ( + 1/8 ) + Ab(X13)*Bb(X11,X33,0,kk)*F0^-3 * ( - 1/2 ) + Ab(X13,eps)*F0^-3*pi16 * ( + 1/4*X13 ) + Ab(X14)*Ab(X34)*F0^-3 * ( + 1/4 ) + Ab(X16)*Ab(X36)*F0^-3 * ( + 1/16 ) + Ab(X44)*F0^-3*L1r * ( - 24*X44 ) + Ab(X44)*F0^-3*L2r * ( - 6*X44 ) + Ab(X44)*F0^-3*L4r * ( + 12*X44 ) + Ab(X44)*Bb(Xpp,Xpp,0)*F0^-3 * ( + 1/12*R(Xpp,X44,Xe0)^2*X44 ) + Ab(X44)*Bb(Xpp,Xe0,0)*F0^-3 * ( + 1/6*R(Xpp,X44,Xe0)*R(Xe0,Xpp,X44)*X44 - 1/6*R(Xpp,X44,Xe0)*R(Xe0,Xqq,X44)*X44 ) + Ab(X44)*Bb(Xe0,Xe0,0)*F0^-3 * ( - 1/24*RV(Xe0,X11,X33)*X44 ) + Ab(X44)*Bb(X11,X33,0)*F0^-3 * ( - 1/6*R(X11,X44,Xe0)*R(X33,X44,Xe0)*X44 ) + Ab(X44,eps)*F0^-3*pi16 * ( + 3/8*X44 ) + Ab(X46)*F0^-3*L1r * ( - 32*X46 ) + Ab(X46)*F0^-3*L2r * ( - 8*X46 ) + Ab(X46)*F0^-3*L4r * ( + 16*X46 ) + Ab(X46)*Bb(Xpp,Xpp,0)*F0^-3 * ( + 1/9*X46 + 1/12*R(Xe0,Xpp,Xpp)*Xpp + 1/36*R(Xe0,Xpp,Xpp)*X44 + 1/9*R(Xe0,Xpp,X44)*X66 ) + Ab(X46)*Bb(Xpp,Xe0,0)*F0^-3 * ( - 1/18*R(Xe0,Xpp,Xpp)*X44 - 1/9*R(Xe0,Xpp,X44)*X66 + 1/9*R(Xe0,Xqq,X44)*X66 + 1/18*R(Xe0,X11,X33)*X44 ) + Ab(X46)*Bb(Xpp,Xe0,0,kk)*F0^-3 * ( - 1/6*R(Xe0,Xpp,Xpp) + 1/6*R(Xe0,X11,X33) ) + Ab(X46)*Bb(Xe0,Xe0,0)*F0^-3 * ( + 1/9*RV(Xe0,X11,X33)*X46 ) + Ab(X46)*Bb(X11,X33,0)*F0^-3 * ( - 2/9*X46 - 1/9*R(Xe0,Xpp,X44)*X66 - 1/18*R(Xe0,X11,X33)*X44 ) + Ab(X46)*Bb(X11,X33,0,kk)*F0^-3 * ( - 1/6*R(Xe0,X11,X33) ) + Ab(X46,eps)*F0^-3*pi16 * ( + 1/2*X46 ) + Bb(Xpp,Xpp,0)*F0^-3*pi16 * ( + 1/16*Xb(1)*RD(Xpp) + 1/96*RD(Xpp)*Xpp + 1/32*RD(Xpp)*Xqq ) + Bb(Xpp,Xpp,0)*F0^-3*L0r * ( + 2/3*RD(Xpp)*Xpp ) + Bb(Xpp,Xpp,0)*F0^-3*L3r * ( + 5/3*RD(Xpp)*Xpp ) + Bb(Xpp,Xpp,0)*F0^-3*L4r * ( - 2*Xb(1)*Xb(Xpp,Xe0,Xpp,Xe0,0)*Xpp - 4*Xb(1)*R(Xpp,Xqq,Xe0)*Xpp + 4*Xb(1)*RC(Xpp)*Xpp + 3*Xb(1)*RD(Xpp) ) + Bb(Xpp,Xpp,0)*F0^-3*L5r * ( - 2/3*Xb(Xpp,Xe0,Xpp,Xe0,1)*Xpp - 4/3*R(Xpp,Xqq,Xe0)*Xpp^2 + 4/3*RC(Xpp)*Xpp^2 + 1/2*RD(Xpp)*Xpp - 1/6*RD(Xpp)*Xqq ) + Bb(Xpp,Xpp,0)*F0^-3*L6r * ( + 4*Xb(1)*Xb(Xpp,Xe0,Xpp,Xe0,1) + 8*Xb(1)*R(Xpp,Xqq,Xe0)*Xpp - 8*Xb(1)*RC(Xpp)*Xpp ) + Bb(Xpp,Xpp,0)*F0^-3*L7r * ( + 4*RD(Xpp)^2 ) + Bb(Xpp,Xpp,0)*F0^-3*L8r * ( + 4/3*Xb(Xpp,Xe0,Xpp,Xe0,2) + 8/3*R(Xpp,Xqq,Xe0)*Xpp^2 - 8/3*RC(Xpp)*Xpp^2 ) + Bb(Xpp,Xpp,0)^2*F0^-3 * ( - 1/18*R(Xpp,Xqq,Xe0)*RD(Xpp)*Xpp + 1/18*RC(Xpp)*RD(Xpp)*Xpp + 1/288*RD(Xpp)^2 ) + Bb(Xpp,Xpp,0)*Bb(Xpp,Xe0,0)*F0^-3 * ( + 1/18*R(Xe0,Xpp,Xpp)*RD(Xpp)*Xpp - 1/18*R(Xe0,X11,X33)*RD(Xpp)*Xpp ) + Bb(Xpp,Xpp,0)*Bb(X11,X33,0)*F0^-3 * ( - 1/18*R(Xqq,Xpp,Xe0)*RD(Xpp)*Xpp ) + Bb(Xpp,Xpp,0)*Cb(Xpp,Xpp,Xpp,0)*F0^-3 * ( + 1/18*RD(Xpp)^2*Xpp ) + Bb(Xpp,Xpp,0,eps)*F0^-3*pi16 * ( - 1/16*Xb(1)*RD(Xpp) - 1/16*RD(Xpp)*Xpp ) + Bb(Xpp,Xe0,0)*F0^-3*L7r * ( + 8*RS(Xqq,Xpp,Xe0)*RD(Xpp)*RD(Xe0,X44,X66,Xpp) ) + Bb(Xpp,Xe0,0)*F0^-3*L8r * ( + 8/3*RS(Xqq,Xpp,Xe0)*RD(Xpp)*RD(Xe0,X44,X66,Xpp) ) + Bb(Xps,Xps,0)*F0^-3*L4r * ( - 6*Xb(1)*Xps ) + Bb(Xps,Xps,0)*F0^-3*L5r * ( - 2*Xps^2 ) + Bb(Xps,Xps,0)*F0^-3*L6r * ( + 12*Xb(1)*Xps ) + Bb(Xps,Xps,0)*F0^-3*L8r * ( + 4*Xps^2 ) + Bb(Xe0,Xe0,0)*F0^-3*L4r * ( + 2*Xb(1)*RV(Xe0,X11,X33)*Xe0 ) + Bb(Xe0,Xe0,0)*F0^-3*L5r * ( + 2/3*RV(Xe0,X11,X33)*Xe0^2 ) + Bb(Xe0,Xe0,0)*F0^-3*L6r * ( - 4*Xb(1)*RV(Xe0,X11,X33)*Xe0 ) + Bb(Xe0,Xe0,0)*F0^-3*L7r * ( + 4*RS(X33,X11,X11,Xe0,Xe0)*RD(Xe0,X44,X66,X11)^2 ) + Bb(Xe0,Xe0,0)*F0^-3*L8r * ( - 4/9*RV(Xe0,X11,X33)*X44^2 - 8/9*RV(Xe0,X11,X33)*X66^2 ) + Bb(X11,X11,0)*Bb(X33,X33,0)*F0^-3 * ( + 1/144*RD(X11)*RD(X33) ) + Bb(X11,X33,0)*F0^-3*L7r * ( - 8*RD(X11)*RD(X33) ) + Bb(X11,X33,0)*F0^-3*L8r * ( - 8/3*RD(X11)*RD(X33) ) + Cb(Xpp,Xpp,Xpp,0)*F0^-3*L4r * ( + 4*Xb(1)*RD(Xpp)*Xpp ) + Cb(Xpp,Xpp,Xpp,0)*F0^-3*L5r * ( + 4/3*RD(Xpp)*Xpp^2 ) + Cb(Xpp,Xpp,Xpp,0)*F0^-3*L6r * ( - 8*Xb(1)*RD(Xpp)*Xpp ) + Cb(Xpp,Xpp,Xpp,0)*F0^-3*L8r * ( - 8/3*RD(Xpp)*Xpp^2 ) + Hbb(1,Xpp,Xpp,X13,X13)*F0^-3 * ( + 1/8*Xpp - 1/16*X13 - 1/36*R(Xpp,Xqq,Xe0)^2*X13 + 1/36*R(Xpp,Xqq,Xe0)*RC(Xpp)*X13 - 1/144*RC(Xpp)^2*X13 ) + Hbb(1,Xpp,X1s,X3s,X13)*F0^-3 * ( - 1/12*R(Xpp,Xqq,Xe0)*Xqs + 1/24*R(Xpp,Xqq,Xe0)*X13 - 1/16*R(Xpp,Xss,Xe0)*Xpp + 1/48*R(Xpp,Xss,Xe0)*Xqq + 1/24*RC(Xpp)*Xps ) + Hbb(1,Xp4,Xq6,X46,X13)*F0^-3 * ( - 1/4*X46 ) + Hbb(1,Xe0,Xpp,X13,X13)*F0^-3 * ( + 1/36*R(Xpp,Xqq,Xe0)*RV(Xe0,X11,X33)*X13 - 1/72*RC(Xpp)*RV(Xe0,X11,X33)*X13 ) + Hbb(1,Xe0,Xe0,X13,X13)*F0^-3 * ( - 1/144*RV(Xe0,X11,X33)^2*X13 ) + Hbb(1,Xe0,X1s,X3s,X13)*F0^-3 * ( - 1/48*R(Xe0,Xpp,Xss)*Xpp + 1/48*R(Xe0,Xpp,Xss)*Xqq + 1/48*RV(Xe0,Xpp,Xss)*Xe0 + 1/48*RV(Xe0,X11,X33)*Xss ) + Hbb(1,X11,X13,X33,X13)*F0^-3 * ( + 1/36*R(Xpp,Xqq,Xe0)*RC(Xqq)*X13 - 1/18*R(X11,X33,Xe0)*R(X33,X11,Xe0)*X13 - 1/72*RC(X11)*RC(X33)*X13 ) + Hbb(1,X13,X13,X13,X13)*F0^-3 * ( - 1/8*X13 ) + Hbb(1,X14,X34,X44,X13)*F0^-3 * ( - 3/8*X44 ) + Hbb(2,Xpp,Xpp,X13,X13)*F0^-3 * ( + 1/36*R(Xpp,Xqq,Xe0)*RD(Xpp)*X13 - 1/72*RC(Xpp)*RD(Xpp)*X13 ) + Hbb(2,Xpp,Xe0,X13,X13)*F0^-3 * ( - 1/72*RD(Xpp)*RV(Xe0,X11,X33)*X13 ) + Hbb(2,Xpp,X13,Xqq,X13)*F0^-3 * ( + 1/36*R(Xqq,Xpp,Xe0)*RD(Xpp)*X13 - 1/72*RC(Xqq)*RD(Xpp)*X13 ) + Hbb(2,Xpp,X1s,X3s,X13)*F0^-3 * ( + 1/24*RD(Xpp)*Xps ) + Hbb(5,Xpp,Xpp,X13,X13)*F0^-3 * ( - 1/144*RD(Xpp)^2*X13 ) + Hbb(5,X11,X33,X13,X13)*F0^-3 * ( - 1/72*RD(X11)*RD(X33)*X13 ) + dHbb(1,Xpp,X1s,X3s,X13)*F0^-3 * ( - 5/24*R(Xpp,Xqq,Xe0)*Xqq*X13 + 1/12*R(Xpp,Xqq,Xe0)*Xps*X13 + 1/16*R(Xpp,Xss,Xe0)*Xqq*X13 + 7/48*R(Xpp,Xss,Xe0)*X13*Xss - 1/24*RC(Xpp)*Xps*X13 ) + dHbb(1,Xp4,Xq6,X46,X13)*F0^-3 * ( + 1/4*X13*X46 ) + dHbb(1,Xe0,Xpp,X13,X13)*F0^-3 * ( + 1/9*R(Xpp,Xqq,Xe0)*R(Xe0,X11,X33)*X13^2 - 1/36*R(Xpp,Xqq,Xe0)*RV(Xe0,X11,X33)*X13^2 + 1/18*R(Xe0,Xpp,Xpp)*RC(Xpp)*X13^2 + 1/72*RC(Xpp)*RV(Xe0,X11,X33)*X13^2 ) + dHbb(1,Xe0,Xe0,X13,X13)*F0^-3 * ( + 1/9*R(Xe0,X11,X33)^2*X13^2 + 1/9*R(Xe0,X11,X33)*RV(Xe0,X11,X33)*X13^2 + 5/144*RV(Xe0,X11,X33)^2*X13^2 ) + dHbb(1,Xe0,X1s,X3s,X13)*F0^-3 * ( - 1/16*R(Xe0,Xpp,Xss)*Xpp*X13 - 5/48*R(Xe0,Xpp,Xss)*Xqq*X13 + 1/6*R(Xe0,X11,Xss)*RS(Xe0,Xss,X33)*X13^2 - 1/48*RV(Xe0,Xpp,Xss)*Xe0*X13 - 1/48*RV(Xe0,X11,X33)*X13*Xss ) + dHbb(1,X11,X13,X33,X13)*F0^-3 * ( - 1/36*R(Xpp,Xqq,Xe0)*RC(Xqq)*X13^2 + 1/6*R(X11,X33,Xe0)*R(X33,X11,Xe0)*X13^2 + 1/72*RC(X11)*RC(X33)*X13^2 ) + dHbb(1,X13,Xpp,Xpp,X13)*F0^-3 * ( - 1/8*Xpp*X13 - 3/16*X13^2 - 1/36*R(Xpp,Xqq,Xe0)^2*X13^2 + 1/12*R(Xpp,Xqq,Xe0)*RC(Xpp)*X13^2 - 1/48*RC(Xpp)^2*X13^2 ) + dHbb(1,X13,X13,X13,X13)*F0^-3 * ( + 1/8*X13^2 ) + dHbb(1,X14,X34,X44,X13)*F0^-3 * ( + 3/8*X13*X44 ) + dHbb(2,Xpp,Xpp,X13,X13)*F0^-3 * ( - 1/36*R(Xpp,Xqq,Xe0)*RD(Xpp)*X13^2 + 5/72*RC(Xpp)*RD(Xpp)*X13^2 ) + dHbb(2,Xpp,Xe0,X13,X13)*F0^-3 * ( + 1/18*R(Xe0,Xpp,Xpp)*RD(Xpp)*X13^2 + 1/72*RD(Xpp)*RV(Xe0,X11,X33)*X13^2 ) + dHbb(2,Xpp,X13,Xqq,X13)*F0^-3 * ( - 1/36*R(Xqq,Xpp,Xe0)*RD(Xpp)*X13^2 + 1/72*RC(Xqq)*RD(Xpp)*X13^2 ) + dHbb(2,Xpp,X1s,X3s,X13)*F0^-3 * ( - 1/24*RD(Xpp)*Xps*X13 ) + dHbb(5,Xpp,Xpp,X13,X13)*F0^-3 * ( + 5/144*RD(Xpp)^2*X13^2 ) + dHbb(5,X11,X33,X13,X13)*F0^-3 * ( + 1/72*RD(X11)*RD(X33)*X13^2 ) + dH1bb(1,Xpp,Xpp,X13,X13)*F0^-3 * ( + X13^2 + 2/9*R(Xpp,Xqq,Xe0)^2*X13^2 - 2/9*R(Xpp,Xqq,Xe0)*RC(Xpp)*X13^2 + 1/9*RC(Xpp)^2*X13^2 ) + dH1bb(1,Xpp,X1s,X3s,X13)*F0^-3 * ( + 1/3*R(Xpp,Xqq,Xe0)*RS(Xss,Xqq,Xpp)*X13^2 ) + dH1bb(1,Xps,Xqs,Xpp,X13)*F0^-3 * ( - 1/3*R(Xpp,Xss,Xe0)*X13^2 ) + dH1bb(1,Xps,Xqs,Xe0,X13)*F0^-3 * ( - 1/3*R(Xe0,X11,X33)*RS(Xss,Xpp,Xe0)*X13^2 ) + dH1bb(1,X13,Xpp,Xe0,X13)*F0^-3 * ( + 1/9*R(Xpp,Xqq,Xe0)*RV(Xe0,X11,X33)*X13^2 - 1/9*R(Xe0,Xpp,Xpp)*RS(Xqq,Xpp,Xe0)*RC(Xpp)*X13^2 ) + dH1bb(1,X13,Xe0,Xe0,X13)*F0^-3 * ( - 1/9*R(Xe0,X11,X33)*RV(Xe0,X11,X33)*X13^2 - 1/18*RV(Xe0,X11,X33)^2*X13^2 ) + dH1bb(1,X13,X11,X33,X13)*F0^-3 * ( + 1/9*R(Xpp,Xqq,Xe0)*RC(Xqq)*X13^2 - 2/9*R(X11,X33,Xe0)*R(X33,X11,Xe0)*X13^2 ) + dH1bb(3,X13,Xpp,Xpp,X13)*F0^-3 * ( + 1/9*R(Xpp,Xqq,Xe0)*RD(Xpp)*X13^2 - 1/9*RC(Xpp)*RD(Xpp)*X13^2 ) + dH1bb(3,X13,Xpp,Xqq,X13)*F0^-3 * ( + 1/9*R(Xqq,Xpp,Xe0)*RD(Xpp)*X13^2 ) + dH1bb(3,X13,Xpp,Xe0,X13)*F0^-3 * ( + 1/9*R(Xe0,X11,X33)*RS(Xpp,Xqq,Xe0)*RD(Xpp)*X13^2 ) + dH1bb(7,X13,Xpp,Xpp,X13)*F0^-3 * ( - 1/18*RD(Xpp)^2*X13^2 ) + dH21bb(1,Xpp,Xpp,X13,X13)*F0^-3 * ( - 3/8*X13^2 ) + dH21bb(1,Xpp,X1s,X3s,X13)*F0^-3 * ( - 1/8*R(Xpp,Xqq,Xe0)*RS(Xss,Xqq,Xpp)*X13^2 ) + dH21bb(1,Xps,Xqs,Xpp,X13)*F0^-3 * ( + 1/8*R(Xpp,Xqq,Xe0)*X13^2 + 1/8*R(Xpp,Xss,Xe0)*X13^2 - 1/8*RC(Xpp)*X13^2 ) + dH21bb(1,Xps,Xqs,Xqq,X13)*F0^-3 * ( + 1/8*R(Xqq,Xpp,Xe0)*RS(Xss,Xpp,Xqq)*X13^2 ) + dH21bb(1,Xps,Xqs,Xe0,X13)*F0^-3 * ( + 1/8*R(Xe0,X11,X33)*RS(Xpp,Xqq,Xe0)*RS(Xss,Xpp,Xe0)*X13^2 ) + dH21bb(1,Xe0,X1s,X3s,X13)*F0^-3 * ( - 1/8*R(Xe0,X11,X33)*RS(Xss,X11,Xe0)*RS(Xss,X33,Xe0)*X13^2 ) + dH21bb(1,X13,Xpp,Xpp,X13)*F0^-3 * ( + 3/16*X13^2 + 1/12*R(Xpp,Xqq,Xe0)^2*X13^2 - 1/12*R(Xpp,Xqq,Xe0)*RC(Xpp)*X13^2 + 1/48*RC(Xpp)^2*X13^2 ) + dH21bb(1,X13,Xpp,Xe0,X13)*F0^-3 * ( - 1/12*R(Xpp,Xqq,Xe0)*RV(Xe0,X11,X33)*X13^2 + 1/24*RC(Xpp)*RV(Xe0,X11,X33)*X13^2 ) + dH21bb(1,X13,Xe0,Xe0,X13)*F0^-3 * ( + 1/48*RV(Xe0,X11,X33)^2*X13^2 ) + dH21bb(1,X13,X11,X33,X13)*F0^-3 * ( - 1/12*R(Xpp,Xqq,Xe0)*RC(Xqq)*X13^2 + 1/6*R(X11,X33,Xe0)*R(X33,X11,Xe0)*X13^2 + 1/24*RC(X11)*RC(X33)*X13^2 ) + dH21bb(1,X13,X13,X13,X13)*F0^-3 * ( + 3/8*X13^2 ) + dH21bb(1,X44,X14,X34,X13)*F0^-3 * ( + 9/8*X13^2 ) + dH21bb(1,X46,Xp4,Xq6,X13)*F0^-3 * ( + 3/4*X13^2 ) + dH21bb(3,Xps,Xpp,Xqs,X13)*F0^-3 * ( - 1/8*RD(Xpp)*X13^2 ) + dH21bb(3,X13,Xpp,Xpp,X13)*F0^-3 * ( - 1/12*R(Xpp,Xqq,Xe0)*RD(Xpp)*X13^2 + 1/24*RC(Xpp)*RD(Xpp)*X13^2 ) + dH21bb(3,X13,Xpp,Xqq,X13)*F0^-3 * ( - 1/12*R(Xqq,Xpp,Xe0)*RD(Xpp)*X13^2 + 1/24*RC(Xqq)*RD(Xpp)*X13^2 ) + dH21bb(3,X13,Xpp,Xe0,X13)*F0^-3 * ( + 1/24*RD(Xpp)*RV(Xe0,X11,X33)*X13^2 ) + dH21bb(7,X13,Xpp,Xpp,X13)*F0^-3 * ( + 1/48*RD(Xpp)^2*X13^2 ) + dH21bb(7,X13,X11,X33,X13)*F0^-3 * ( + 1/24*RD(X11)*RD(X33)*X13^2 );