*FORM by J.Vermaseren,version 3.1(Oct 17 2002) Run at: Sun Mar 13 15:18:49 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 "2"
*    
*    * 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 "5"
*    #define type "1"
*    
*    #define pout "0"
*    #define nsim "0"
*    
*    #-
* Sigma2 loaded
* Sigma4 loaded
* Sigma6 loaded
*
*Time =       2.75 sec    Generated terms =          1
*             M0p2        Terms in output =          1
*                         Bytes used      =         24
*
*Time =       2.75 sec    Generated terms =          6
*             M0p4        Terms in output =          6
*                         Bytes used      =        238
*
*Time =       2.75 sec    Generated terms =        323
*             M0p6        Terms in output =        323
*                         Bytes used      =      13392

L   M0p2timo22 =
       + X12
         ;

L   M0p4timo22 =

       + F0^-2*L4r * (
          - 16*X12*X34
          )

       + F0^-2*L5r * (
          - 8*X12^2
          )

       + F0^-2*L6r * (
          + 32*X12*X34
          )

       + F0^-2*L8r * (
          + 16*X12^2
          )

       + Ab(Xpp)*F0^-2 * (
          - 1/2*R(Xpp,Xqq,Xp0)*X12
          )

       + Ab(Xp0)*F0^-2 * (
          - 1/2*R(Xp0,X11,X22)*X12
          );

L   M0p6timo22 =

       + F0^-4*pi16*L0r * (
          + 8*Xp0*X12^2
          + 6*Xp0^2*X12
          + 4*X11*X12*X22
          - 4*X12*X33*X44
          - 8*X12^3
          )

       + F0^-4*pi16*L1r * (
          + 4*X12^3
          )

       + F0^-4*pi16*L2r * (
          + 6*Xp0^2*X12
          + 2*X12^3
          )

       + F0^-4*pi16*L3r * (
          + 6*Xp0*X12^2
          + 4*Xp0^2*X12
          + 3*X11*X12*X22
          - 3*X12*X33*X44
          - 8*X12^3
          )

       + F0^-4*pi16^2 * (
          + 35/144*Xp0*X12^2
          + 2/3*Xp0^2*X12
          + 17/96*X11*X12*X22
          - 17/96*X12*X33*X44
          - 17/72*X12^3
          )

       + F0^-4*L4r*L5r * (
          + 256*X12^2*X34
          )

       + F0^-4*L4r*L6r * (
          - 512*X12*X34^2
          )

       + F0^-4*L4r*L8r * (
          - 256*X12^2*X34
          )

       + F0^-4*L4r*L11r * (
          + 128*X12^2*X34
          )

       + F0^-4*L4r^2 * (
          + 256*X12*X34^2
          )

       + F0^-4*L5r*L6r * (
          - 256*X12^2*X34
          )

       + F0^-4*L5r*L8r * (
          - 128*X12^3
          )

       + F0^-4*L5r*L11r * (
          + 64*X12^3
          )

       + F0^-4*L5r^2 * (
          + 64*X12^3
          )

       + F0^-4*L6r*L11r * (
          - 256*X12^2*X34
          )

       + F0^-4*L8r*L11r * (
          - 128*X12^3
          )

       + F0^-4*K17r * (
          - 32*X12^3
          )

       + F0^-4*K18r * (
          - 64*X12^2*X34
          )

       + F0^-4*K19r * (
          - 8*Xpp^2*X12
          )

       + F0^-4*K20r * (
          - 32*X12^2*X34
          )

       + F0^-4*K21r * (
          - 16*X12*Xss^2
          )

       + F0^-4*K22r * (
          - 64*X12*X34^2
          )

       + F0^-4*K23r * (
          - 16*X11*X12*X22
          )

       + F0^-4*K25r * (
          + 24*Xpp^2*X12
          )

       + F0^-4*K26r * (
          + 16*X12*Xss^2
          + 64*X12^2*X34
          )

       + F0^-4*K27r * (
          + 192*X12*X34^2
          )

       + F0^-4*K39r * (
          + 32*X12^3
          )

       + F0^-4*K40r * (
          + 64*X12^2*X34
          )

       + Ab(Xpp)*F0^-4*pi16 * (
          + 1/4*R(Xpp,Xqq,Xp0)*X12*X34
          + 1/12*R(Xpp,Xqq,Xp0)*X12^2
          - 1/24*RC(Xpp)*Xqq*X12
          - 1/8*RC(Xpp)*X12*X34
          )

       + Ab(Xpp)*F0^-4*L0r * (
          - 8*R(Xpp,Xqq,Xp0)*Xpp*X12
          - 2*RC(Xpp)*Xpp*X12
          - 2*RD(Xpp)*X12
          )

       + Ab(Xpp)*F0^-4*L3r * (
          - 2*R(Xpp,Xqq,Xp0)*Xpp*X12
          - 5*RC(Xpp)*Xpp*X12
          - 5*RD(Xpp)*X12
          )

       + Ab(Xpp)*F0^-4*L4r * (
          + 16*R(Xpp,Xqq,Xp0)*X12*X34
          )

       + Ab(Xpp)*F0^-4*L5r * (
          + 8*R(Xpp,Xqq,Xp0)*Xpp*X12
          + 4*RC(Xpp)*Xpp*X12
          )

       + Ab(Xpp)*F0^-4*L6r * (
          - 16*R(Xpp,Xqq,Xp0)*X12*X34
          )

       + 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 * (
          - 16/3*R(Xpp,Xqq,Xp0)*X11*X22
          - 32/3*R(Xpp,Xqq,Xp0)*X12^2
          - 8*RC(Xpp)*Xpp*X12
          - 4/3*RD(Xpp)^2
          )

       + Ab(Xpp)^2*F0^-4 * (
          + 1/8*X12
          + 1/8*R(Xpp,Xqq,Xp0)^2*X12
          + 1/24*R(Xpp,Xqq,Xp0)*RC(Xpp)*Xpp
          - 1/24*RC(Xpp)^2*Xpp
          + 1/32*RC(Xpp)^2*X12
          )

       + Ab(Xpp)*Ab(Xps)*F0^-4 * (
          + 1/48*R(Xpp,Xqq,Xp0)*Xpp
          + 19/144*R(Xpp,Xqq,Xp0)*Xqq
          - 5/144*R(Xpp,Xss,Xp0)*Xpp
          - 1/16*R(Xpp,Xss,Xp0)*Xqq
          + 1/18*RC(Xpp)*Xpp
          )

       + Ab(Xpp)*Ab(Xqs)*F0^-4 * (
          + 1/12*R(Xpp,Xqq,Xp0)*Xqq
          - 1/48*R(Xpp,Xss,Xp0)*RS(Xqq,Xss,Xpp)*Xpp
          + 1/144*R(Xpp,Xss,Xp0)*RS(Xqq,Xss,Xpp)*Xqq
          + 1/8*RC(Xpp)*X12
          )

       + Ab(Xpp)*Ab(Xp0)*F0^-4 * (
          + 1/24*R(Xpp,Xqq,Xp0)*R(Xp0,Xpp,Xpp)*Xpp
          - 1/40*R(Xpp,Xqq,Xp0)*R(Xp0,Xqq,Xqq)*Xqq
          + 13/120*R(Xpp,Xqq,Xp0)*R(Xp0,X11,X22)*Xpp
          + 1/8*R(Xpp,Xqq,Xp0)*R(Xp0,X11,X22)*Xqq
          - 1/60*R(Xpp,Xp0,Xp0)*R(Xp0,Xqq,Xqq)*Xp0
          - 1/24*R(Xp0,Xqq,Xqq)*RC(Xpp)*Xpp
          - 1/24*RC(Xpp)*Xp0
          + 1/96*RC(Xpp)*RV(Xp0,X11,X22)*Xpp
          + 1/32*RC(Xpp)*RV(Xp0,X11,X22)*Xqq
          - 1/16*RC(Xpp)*RV(Xp0,X11,X22)*Xp0
          - 1/12*RD(Xp0,X33,X44,Xpp)
          )

       + Ab(Xpp)*Ab(X12)*F0^-4 * (
          + 1/18*Xpp
          - 1/36*X12
          )

       + Ab(Xpp)*Ab(X34)*F0^-4 * (
          - 1/12*R(Xpp,Xp0,Xp0)*X34
          + 1/9*R(Xp0,Xpp,Xpp)*Xpp
          )

       + Ab(Xpp)*Bb(Xpp,Xpp,0)*F0^-4 * (
          + 11/24*R(Xpp,Xqq,Xp0)*RC(Xpp)*Xpp*X12
          - 1/12*R(Xpp,Xqq,Xp0)*RC(Xpp)*Xpp^2
          - 1/8*R(Xpp,Xqq,Xp0)*RC(Xpp)*X12^2
          + 1/24*R(Xpp,Xqq,Xp0)*RD(Xpp)*Xpp
          )

       + Ab(Xpp)*Bb(Xpp,Xp0,0)*F0^-4 * (
          + 1/6*R(Xp0,X11,X22)*RC(Xpp)*Xpp^2
          + 1/12*R(Xp0,X11,X22)*RC(Xpp)*X11*X22
          )

       + Ab(Xpp)*Bb(Xqq,Xqq,0)*F0^-4 * (
          - 1/24*R(Xpp,Xqq,Xp0)*RD(Xqq)*Xqq
          - 1/16*RC(Xpp)*RD(Xqq)*X12
          )

       + Ab(Xpp)*Bb(X11,X22,0)*F0^-4 * (
          + 1/6*R(Xqq,Xpp,Xp0)*RC(Xpp)*Xpp*X12
          )

       + Ab(Xpp)*Bb(X11,X22,0,kk)*F0^-4 * (
          + 1/12*R(Xqq,Xpp,Xp0)*RC(Xpp)*Xpp
          )

       + Ab(Xpp,eps)*F0^-4*pi16 * (
          + 1/4*Xpp*X12
          - 1/4*R(Xpp,Xqq,Xp0)*X12*X34
          - 1/4*R(Xp0,Xpp,Xpp)*X12*X34
          )

       + Ab(Xps)*F0^-4*pi16 * (
          + 1/6*Xqq*X12
          + 1/4*X12*X34
          + 1/12*X12*Xss
          - 1/6*X12^2
          )

       + Ab(Xps)*F0^-4*L0r * (
          + 4*Xps*X12
          )

       + Ab(Xps)*F0^-4*L3r * (
          + 10*Xps*X12
          )

       + Ab(Xps)*F0^-4*L5r * (
          - 8*Xps*X12
          )

       + Ab(Xps)*F0^-4*L8r * (
          + 16*Xps*X12
          )

       + Ab(Xps)*Ab(Xp0)*F0^-4 * (
          + 1/18*R(Xp0,Xpp,Xpp)*Xpp
          - 5/144*R(Xp0,Xpp,Xss)*Xpp
          - 1/16*R(Xp0,Xpp,Xss)*Xqq
          + 1/8*R(Xp0,Xqq,Xqq)*X12
          - 5/48*R(Xp0,Xqq,Xss)*Xpp
          - 19/144*R(Xp0,Xqq,Xss)*RS(Xpp,Xss,Xp0)*Xqq
          + 1/48*R(Xp0,X11,X22)*Xpp
          )

       + Ab(Xps)*Bb(Xpp,Xp0,0)*F0^-4 * (
          - 2/9*R(Xp0,Xqq,Xss)*Xpp*Xps
          - 1/9*R(Xp0,Xqq,Xss)*Xpp*X12
          - 1/9*R(Xp0,Xqq,Xss)*Xqq*Xps
          - 1/18*R(Xp0,Xqq,Xss)*RS(Xpp,Xss,Xp0)*X11*X22
          - 1/18*R(Xp0,X11,X22)*Xpp^2
          )

       + Ab(Xps)*Bb(Xqq,Xqq,0)*F0^-4 * (
          + 1/8*RD(Xqq)*X12
          )

       + Ab(Xps)*Bb(X11,X22,0)*F0^-4 * (
          - 1/18*R(Xqq,Xpp,Xp0)*Xpp^2
          - 2/9*R(Xqq,Xss,Xp0)*Xps*X12
          - 1/18*R(Xqq,Xss,Xp0)*X11*X22
          )

       + Ab(Xps)*Bb(X11,X22,0,kk)*F0^-4 * (
          + 1/18*R(Xqq,Xpp,Xp0)*Xqq
          - 1/9*R(Xqq,Xss,Xp0)*Xps
          - 1/9*R(Xqq,Xss,Xp0)*X12
          )

       + Ab(Xps,eps)*F0^-4*pi16 * (
          - 1/4*X12*X34
          - 1/4*X12*Xss
          )

       + Ab(Xp3)*Ab(Xq4)*F0^-4 * (
          - 1/4*X12
          )

       + Ab(Xp0)*F0^-4*pi16 * (
          - 1/12*RV(Xp0,X11,X22)*Xp0*X12
          - 1/12*RV(Xp0,X11,X22)*X12^2
          )

       + Ab(Xp0)*F0^-4*L0r * (
          - 12*R(Xp0,X11,X22)*Xp0*X12
          - 2*RV(Xp0,X11,X22)*Xp0*X12
          )

       + Ab(Xp0)*F0^-4*L1r * (
          + 16*Xp0*X12
          )

       + Ab(Xp0)*F0^-4*L2r * (
          + 4*Xp0*X12
          )

       + Ab(Xp0)*F0^-4*L3r * (
          - 12*R(Xp0,X11,X22)*Xp0*X12
          - 5*RV(Xp0,X11,X22)*Xp0*X12
          )

       + Ab(Xp0)*F0^-4*L4r * (
          - 16*Xp0*X12
          + 16*R(Xp0,X11,X22)*Xp0*X12
          )

       + Ab(Xp0)*F0^-4*L5r * (
          + 4*R(Xp0,Xpp,Xpp)*Xp0*X12
          + 8*R(Xp0,X11,X22)*X12^2
          )

       + Ab(Xp0)*F0^-4*L6r * (
          + 16*Xp0*X12
          - 16*R(Xp0,X11,X22)*Xp0*X12
          )

       + Ab(Xp0)*F0^-4*L7r * (
          + 64/3*R(Xp0,X11,X22)*Xp0*X12
          - 64/3*R(Xp0,X11,X22)*X12^2
          + 16/3*RD(Xp0,X33,X44,Xpp)*Xqq
          - 8/3*RD(Xp0,X44)^2*RV(Xp0,X11,X22)
          )

       + Ab(Xp0)*F0^-4*L8r * (
          - 32/3*R(Xp0,X11,X22)*X11*X22
          - 64/3*R(Xp0,X11,X22)*X12^2
          - 16/3*RV(Xp0,X11,X22)*Xp0^2
          - 4*RV(Xp0,X11,X22)*X11*X22
          + 4/3*RV(Xp0,X11,X22)*X33*X44
          )

       + Ab(Xp0)^2*F0^-4 * (
          + 1/8*R(Xp0,X11,X22)^2*X12
          - 1/24*R(Xp0,X11,X22)*RV(Xp0,X11,X22)*Xp0
          + 1/24*R(Xp0,X11,X22)*RV(Xp0,X11,X22)*X12
          - 1/24*RV(Xp0,X11,X22)*Xp0
          - 1/24*RV(Xp0,X11,X22)^2*Xp0
          + 1/32*RV(Xp0,X11,X22)^2*X12
          )

       + Ab(Xp0)*Ab(X34)*F0^-4 * (
          - 1/3*R(Xp0,X11,X22)*X12
          - 1/36*RV(Xp0,X11,X22)*Xp0
          )

       + Ab(Xp0)*Bb(Xpp,Xpp,0)*F0^-4 * (
          + 1/4*R(Xpp,Xqq,Xp0)*R(Xp0,Xpp,Xpp)*Xpp^2
          - 1/12*R(Xp0,Xpp,Xpp)*RD(Xpp)*Xpp
          - 1/24*R(Xp0,Xqq,Xqq)*RD(Xpp)*Xpp
          + 1/96*RD(Xpp)*RV(Xp0,X11,X22)*Xpp
          + 1/32*RD(Xpp)*RV(Xp0,X11,X22)*Xqq
          - 1/16*RD(Xpp)*RV(Xp0,X11,X22)*Xp0
          )

       + Ab(Xp0)*Bb(Xpp,Xp0,0)*F0^-4 * (
          + 1/6*R(Xpp,Xp0,Xp0)*R(Xp0,X11,X22)*Xpp*X12
          + 1/12*R(Xpp,Xp0,Xp0)*R(Xp0,X11,X22)*Xpp^2
          )

       + Ab(Xp0)*Bb(Xp0,Xp0,0)*F0^-4 * (
          - 1/4*R(Xp0,X11,X22)*X12*X34
          )

       + Ab(Xp0)*Bb(X11,X22,0)*F0^-4 * (
          + 1/6*R(X11,Xp0,Xp0)*R(X22,Xp0,Xp0)*X12*X34
          )

       + Ab(Xp0)*Bb(X11,X22,0,kk)*F0^-4 * (
          + 1/12*R(X11,Xp0,Xp0)*R(X22,Xp0,Xp0)*X34
          )

       + Ab(Xp0,eps)*F0^-4*pi16 * (
          - 1/4*Xp0*X12
          + 1/4*R(Xp0,X11,X22)*Xp0*X12
          + 1/4*RV(Xp0,X11,X22)*Xp0*X12
          )

       + Ab(X11)*Ab(X22)*F0^-4 * (
          - 1/24*R(Xpp,Xqq,Xp0)*RC(Xqq)*Xqq
          + 1/4*R(X11,X22,Xp0)*R(X22,X11,Xp0)*X12
          - 1/16*RC(X11)*RC(X22)*X12
          )

       + Ab(X12)*F0^-4*L1r * (
          + 8*X12^2
          )

       + Ab(X12)*F0^-4*L2r * (
          + 20*X12^2
          )

       + Ab(X12)*F0^-4*L6r * (
          + 32*X12^2
          )

       + Ab(X12)^2*F0^-4 * (
          - 1/4*X12
          )

       + Ab(X12)*Bb(X11,X22,0)*F0^-4 * (
          + 1/9*X11*X22
          + 4/9*X12^2
          )

       + Ab(X12)*Bb(X11,X22,0,kk)*F0^-4 * (
          + 4/9*X12
          )

       + Ab(X12,eps)*F0^-4*pi16 * (
          - 2*X12^2
          )

       + Ab(X34)*F0^-4*L1r * (
          + 32*X12*X34
          )

       + Ab(X34)*F0^-4*L2r * (
          + 8*X12*X34
          )

       + Ab(X34)*F0^-4*L4r * (
          - 32*X12*X34
          )

       + Ab(X34)*F0^-4*L6r * (
          + 32*X12*X34
          )

       + Ab(X34)*Bb(Xpp,Xp0,0)*F0^-4 * (
          - 2/9*R(Xp0,X11,X22)*Xpp*Xp0
          - 1/9*R(Xp0,X11,X22)*Xqq*Xp0
          - 1/6*R(Xp0,X11,X22)*X11*X22
          )

       + Ab(X34)*Bb(Xp0,Xp0,0)*F0^-4 * (
          + 1/2*R(Xp0,X11,X22)*Xp0*X12
          )

       + Ab(X34)*Bb(X11,X22,0)*F0^-4 * (
          + 1/3*X12*X34
          + 1/9*R(Xp0,X11,X22)*X11*X22
          + 1/3*R(Xp0,X11,X22)*X12*X34
          - 2/9*R(Xp0,X11,X22)*X12^2
          )

       + Ab(X34)*Bb(X11,X22,0,kk)*F0^-4 * (
          + 1/6*X34
          + 1/9*R(Xp0,X11,X22)*X12
          + 1/6*R(Xp0,X11,X22)*X34
          )

       + Ab(X34,eps)*F0^-4*pi16 * (
          - 1/2*X12*X34
          )

       + Ab(X1s)*Ab(X2s)*F0^-4 * (
          - 1/4*X12
          )

       + Bb(Xpp,Xpp,0)*F0^-4*pi16 * (
          - 1/24*RD(Xpp)*Xqq*X12
          - 1/8*RD(Xpp)*X12*X34
          )

       + Bb(Xpp,Xpp,0)*F0^-4*L0r * (
          - 2*RD(Xpp)*Xpp*X12
          )

       + Bb(Xpp,Xpp,0)*F0^-4*L3r * (
          - 5*RD(Xpp)*Xpp*X12
          )

       + Bb(Xpp,Xpp,0)*F0^-4*L4r * (
          + 8*R(Xpp,Xqq,Xp0)*Xpp*X12*X34
          )

       + Bb(Xpp,Xpp,0)*F0^-4*L5r * (
          + 4*R(Xpp,Xqq,Xp0)*Xpp^3
          + 2*RD(Xpp)*X11*X22
          )

       + Bb(Xpp,Xpp,0)*F0^-4*L6r * (
          - 16*R(Xpp,Xqq,Xp0)*Xpp*X12*X34
          )

       + Bb(Xpp,Xpp,0)*F0^-4*L8r * (
          - 8*R(Xpp,Xqq,Xp0)*Xpp^3
          - 4*RD(Xpp)*X11*X22
          )

       + Bb(Xpp,Xpp,0)^2*F0^-4 * (
          + 9/32*R(Xpp,Xqq,Xp0)*RD(Xpp)*Xpp*X12
          - 1/32*R(Xpp,Xqq,Xp0)*RD(Xpp)*Xqq*X12
          )

       + Bb(Xpp,Xpp,0)*Bb(Xpp,Xp0,0)*F0^-4 * (
          + 1/6*R(Xp0,X11,X22)*RD(Xpp)*Xpp^2
          + 1/12*R(Xp0,X11,X22)*RD(Xpp)*X11*X22
          )

       + Bb(Xpp,Xpp,0)*Bb(X11,X22,0)*F0^-4 * (
          + 1/6*R(Xqq,Xpp,Xp0)*RD(Xpp)*Xpp*X12
          )

       + Bb(Xpp,Xpp,0)*Bb(X11,X22,0,kk)*F0^-4 * (
          + 1/12*R(Xqq,Xpp,Xp0)*RD(Xpp)*Xpp
          )

       + Bb(Xpp,Xpp,0,eps)*F0^-4*pi16 * (
          + 1/8*RD(Xpp)*Xpp*X12
          + 1/8*RD(Xpp)*X12*X34
          )

       + Bb(Xpp,Xp0,0)*F0^-4*L7r * (
          + 16/3*RD(Xpp)*RD(Xp0,X33,X44,Xqq)*Xpp
          + 32/3*RD(Xpp)*RD(Xp0,X33,X44,Xqq)*X12
          )

       + Bb(Xpp,Xp0,0)*F0^-4*L8r * (
          + 8/3*RD(Xpp)*RD(Xp0,X33,X44,Xqq)*Xpp
          + 16/3*RD(Xpp)*RD(Xp0,X33,X44,Xqq)*X12
          )

       + Bb(Xp0,Xp0,0)*F0^-4*L4r * (
          + 8*R(Xp0,X11,X22)*Xp0^2*X12
          )

       + Bb(Xp0,Xp0,0)*F0^-4*L5r * (
          + 4*R(Xp0,X11,X22)*Xp0^2*X12
          )

       + Bb(Xp0,Xp0,0)*F0^-4*L6r * (
          - 16*R(Xp0,X11,X22)*Xp0^2*X12
          )

       + Bb(Xp0,Xp0,0)*F0^-4*L7r * (
          - 16*R(Xp0,X11,X22)*RD(Xp0,X44)^2*X12
          )

       + Bb(Xp0,Xp0,0)*F0^-4*L8r * (
          - 4*R(Xp0,X11,X22)*X12*X33^2
          - 4*R(Xp0,X11,X22)*X12*X44^2
          )

       + Bb(X11,X11,0)*Bb(X22,X22,0)*F0^-4 * (
          - 1/16*RD(X11)*RD(X22)*X12
          )

       + Bb(X11,X22,0)*F0^-4*L7r * (
          + 32/3*RD(X11)*RD(X22)*X12
          )

       + Bb(X11,X22,0)*F0^-4*L8r * (
          + 16/3*RD(X11)*RD(X22)*X12
          )

       + Bb(X11,X22,0,kk)*F0^-4*L7r * (
          + 16/3*RD(X11)*RD(X22)
          )

       + Bb(X11,X22,0,kk)*F0^-4*L8r * (
          + 8/3*RD(X11)*RD(X22)
          )

       + Hbb(1,Xpp,Xpp,X12,X12)*F0^-4 * (
          - 1/4*Xpp*X12
          - 3/8*X12^2
          - 1/8*R(Xpp,Xqq,Xp0)^2*X12^2
          + 3/8*R(Xpp,Xqq,Xp0)*RC(Xpp)*X12^2
          - 3/32*RC(Xpp)^2*X12^2
          )

       + Hbb(1,Xpp,X1s,X2s,X12)*F0^-4 * (
          - 5/8*R(Xpp,Xqq,Xp0)*Xqq*X12
          + 1/4*R(Xpp,Xqq,Xp0)*Xps*X12
          + 3/16*R(Xpp,Xss,Xp0)*Xqq*X12
          + 7/16*R(Xpp,Xss,Xp0)*X12*Xss
          - 1/8*RC(Xpp)*Xps*X12
          )

       + Hbb(1,Xp3,Xq4,X34,X12)*F0^-4 * (
          + 1/4*X12*X34
          )

       + Hbb(1,Xp0,Xpp,X12,X12)*F0^-4 * (
          + 1/2*R(Xpp,Xqq,Xp0)*R(Xp0,X11,X22)*X12^2
          - 1/8*R(Xpp,Xqq,Xp0)*RV(Xp0,X11,X22)*X12^2
          + 1/4*R(Xp0,Xpp,Xpp)*RC(Xpp)*X12^2
          + 1/16*RC(Xpp)*RV(Xp0,X11,X22)*X12^2
          )

       + Hbb(1,Xp0,Xp0,X12,X12)*F0^-4 * (
          + 1/2*R(Xp0,X11,X22)^2*X12^2
          + 1/2*R(Xp0,X11,X22)*RV(Xp0,X11,X22)*X12^2
          + 5/32*RV(Xp0,X11,X22)^2*X12^2
          )

       + Hbb(1,Xp0,X1s,X2s,X12)*F0^-4 * (
          - 3/16*R(Xp0,Xpp,Xss)*Xpp*X12
          - 5/16*R(Xp0,Xpp,Xss)*Xqq*X12
          + 1/2*R(Xp0,X11,X22)*X12^2
          - 1/16*RV(Xp0,Xpp,Xss)*Xp0*X12
          - 1/16*RV(Xp0,X11,X22)*X12*Xss
          )

       + Hbb(1,X11,X12,X22,X12)*F0^-4 * (
          - 1/8*R(Xpp,Xqq,Xp0)*RC(Xqq)*X12^2
          + 3/4*R(X11,X22,Xp0)*R(X22,X11,Xp0)*X12^2
          + 1/16*RC(X11)*RC(X22)*X12^2
          )

       + Hbb(1,X12,X12,X12,X12)*F0^-4 * (
          + 1/4*X12^2
          )

       + Hbb(2,Xpp,Xpp,X12,X12)*F0^-4 * (
          - 1/8*R(Xpp,Xqq,Xp0)*RD(Xpp)*X12^2
          + 5/16*RC(Xpp)*RD(Xpp)*X12^2
          )

       + Hbb(2,Xpp,Xp0,X12,X12)*F0^-4 * (
          + 1/4*R(Xp0,Xpp,Xpp)*RD(Xpp)*X12^2
          + 1/16*RD(Xpp)*RV(Xp0,X11,X22)*X12^2
          )

       + Hbb(2,Xpp,X12,Xqq,X12)*F0^-4 * (
          - 1/8*R(Xqq,Xpp,Xp0)*RD(Xpp)*X12^2
          + 1/16*RC(Xqq)*RD(Xpp)*X12^2
          )

       + Hbb(2,Xpp,X1s,X2s,X12)*F0^-4 * (
          - 1/8*RD(Xpp)*Xps*X12
          )

       + Hbb(5,Xpp,Xpp,X12,X12)*F0^-4 * (
          + 5/32*RD(Xpp)^2*X12^2
          )

       + Hbb(5,X11,X22,X12,X12)*F0^-4 * (
          + 1/16*RD(X11)*RD(X22)*X12^2
          )

       + H1bb(1,Xpp,Xpp,X12,X12)*F0^-4 * (
          + 2*X12^2
          + R(Xpp,Xqq,Xp0)^2*X12^2
          - R(Xpp,Xqq,Xp0)*RC(Xpp)*X12^2
          + 1/2*RC(Xpp)^2*X12^2
          )

       + H1bb(1,Xpp,X1s,X2s,X12)*F0^-4 * (
          + R(Xpp,Xqq,Xp0)*RS(Xss,Xqq,Xpp)*X12^2
          )

       + H1bb(1,Xps,Xqs,Xpp,X12)*F0^-4 * (
          - R(Xpp,Xss,Xp0)*X12^2
          )

       + H1bb(1,Xps,Xqs,Xp0,X12)*F0^-4 * (
          - R(Xp0,X11,X22)*RS(Xss,Xpp,Xp0)*X12^2
          )

       + H1bb(1,X12,Xpp,Xp0,X12)*F0^-4 * (
          + 1/2*R(Xpp,Xqq,Xp0)*RV(Xp0,X11,X22)*X12^2
          - 1/2*R(Xp0,Xpp,Xpp)*RS(Xqq,Xpp,Xp0)*RC(Xpp)*X12^2
          )

       + H1bb(1,X12,Xp0,Xp0,X12)*F0^-4 * (
          - 1/2*R(Xp0,X11,X22)*RV(Xp0,X11,X22)*X12^2
          - 1/4*RV(Xp0,X11,X22)^2*X12^2
          )

       + H1bb(1,X12,X11,X22,X12)*F0^-4 * (
          + 1/2*R(Xpp,Xqq,Xp0)*RC(Xqq)*X12^2
          - R(X11,X22,Xp0)*R(X22,X11,Xp0)*X12^2
          )

       + H1bb(3,X12,Xpp,Xpp,X12)*F0^-4 * (
          + 1/2*R(Xpp,Xqq,Xp0)*RD(Xpp)*X12^2
          - 1/2*RC(Xpp)*RD(Xpp)*X12^2
          )

       + H1bb(3,X12,Xpp,Xqq,X12)*F0^-4 * (
          + 1/2*R(Xqq,Xpp,Xp0)*RD(Xpp)*X12^2
          )

       + H1bb(3,X12,Xpp,Xp0,X12)*F0^-4 * (
          + 1/2*R(Xp0,X11,X22)*RS(Xpp,Xqq,Xp0)*RD(Xpp)*X12^2
          )

       + H1bb(7,X12,Xpp,Xpp,X12)*F0^-4 * (
          - 1/4*RD(Xpp)^2*X12^2
          )

       + H21bb(1,Xpp,Xpp,X12,X12)*F0^-4 * (
          - 3/4*X12^2
          )

       + H21bb(1,Xpp,X1s,X2s,X12)*F0^-4 * (
          - 3/8*R(Xpp,Xqq,Xp0)*RS(Xss,Xqq,Xpp)*X12^2
          )

       + H21bb(1,Xps,Xqs,Xpp,X12)*F0^-4 * (
          + 3/8*R(Xpp,Xqq,Xp0)*X12^2
          + 3/8*R(Xpp,Xss,Xp0)*X12^2
          - 3/8*RC(Xpp)*X12^2
          )

       + H21bb(1,Xps,Xqs,Xqq,X12)*F0^-4 * (
          + 3/8*R(Xqq,Xpp,Xp0)*RS(Xss,Xpp,Xqq)*X12^2
          )

       + H21bb(1,Xps,Xqs,Xp0,X12)*F0^-4 * (
          + 3/8*R(Xp0,X11,X22)*RS(Xpp,Xqq,Xp0)*RS(Xss,Xpp,Xp0)*X12^2
          )

       + H21bb(1,Xp0,X1s,X2s,X12)*F0^-4 * (
          - 3/8*R(Xp0,X11,X22)*RS(Xss,X11,Xp0)*RS(Xss,X22,Xp0)*X12^2
          )

       + H21bb(1,X12,Xpp,Xpp,X12)*F0^-4 * (
          + 3/8*X12^2
          + 3/8*R(Xpp,Xqq,Xp0)^2*X12^2
          - 3/8*R(Xpp,Xqq,Xp0)*RC(Xpp)*X12^2
          + 3/32*RC(Xpp)^2*X12^2
          )

       + H21bb(1,X12,Xpp,Xp0,X12)*F0^-4 * (
          - 3/8*R(Xpp,Xqq,Xp0)*RV(Xp0,X11,X22)*X12^2
          + 3/16*RC(Xpp)*RV(Xp0,X11,X22)*X12^2
          )

       + H21bb(1,X12,Xp0,Xp0,X12)*F0^-4 * (
          + 3/32*RV(Xp0,X11,X22)^2*X12^2
          )

       + H21bb(1,X12,X11,X22,X12)*F0^-4 * (
          - 3/8*R(Xpp,Xqq,Xp0)*RC(Xqq)*X12^2
          + 3/4*R(X11,X22,Xp0)*R(X22,X11,Xp0)*X12^2
          + 3/16*RC(X11)*RC(X22)*X12^2
          )

       + H21bb(1,X12,X12,X12,X12)*F0^-4 * (
          + 3/4*X12^2
          )

       + H21bb(1,X34,Xp3,Xq4,X12)*F0^-4 * (
          + 3/4*X12^2
          )

       + H21bb(3,Xps,Xpp,Xqs,X12)*F0^-4 * (
          - 3/8*RD(Xpp)*X12^2
          )

       + H21bb(3,X12,Xpp,Xpp,X12)*F0^-4 * (
          - 3/8*R(Xpp,Xqq,Xp0)*RD(Xpp)*X12^2
          + 3/16*RC(Xpp)*RD(Xpp)*X12^2
          )

       + H21bb(3,X12,Xpp,Xqq,X12)*F0^-4 * (
          - 3/8*R(Xqq,Xpp,Xp0)*RD(Xpp)*X12^2
          + 3/16*RC(Xqq)*RD(Xpp)*X12^2
          )

       + H21bb(3,X12,Xpp,Xp0,X12)*F0^-4 * (
          + 3/16*RD(Xpp)*RV(Xp0,X11,X22)*X12^2
          )

       + H21bb(7,X12,Xpp,Xpp,X12)*F0^-4 * (
          + 3/32*RD(Xpp)^2*X12^2
          )

       + H21bb(7,X12,X11,X22,X12)*F0^-4 * (
          + 3/16*RD(X11)*RD(X22)*X12^2
          );


