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

L   M0p2timo12 =
       + X11
         ;

L   M0p4timo12 =

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

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

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

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

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

       + Ab(X11)*F0^-2 * (
          - 1/2*RC(X11)*X11
          )

       + Bb(X11,X11,0)*F0^-2 * (
          - 1/2*RD(X11)*X11
          );

L   M0p6timo12 =

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

       + F0^-4*K19r * (
          - 16*X11^3
          )

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

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

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

       + F0^-4*K23r * (
          - 16*X11^3
          )

       + F0^-4*K25r * (
          + 48*X11^3
          )

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

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

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

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

       + Ab(Xp0)*F0^-4*L0r * (
          - 12*R(Xp0,X11,X11)*Xp0*X11
          )

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

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

       + Ab(Xp0)*F0^-4*L3r * (
          - 12*R(Xp0,X11,X11)*Xp0*X11
          )

       + Ab(Xp0)*F0^-4*L4r * (
          - 16*RC(X11)*Xp0*X11
          )

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

       + Ab(Xp0)*F0^-4*L6r * (
          + 16*RC(X11)*Xp0*X11
          )

       + Ab(Xp0)*F0^-4*L7r * (
          + 32*RD(Xp0,X33,X44,X11)*X11
          )

       + Ab(Xp0)*F0^-4*L8r * (
          - 32*R(Xp0,X11,X11)*X11^2
          )

       + Ab(Xp0)^2*F0^-4 * (
          + 1/8*R(Xp0,X11,X11)^2*X11
          )

       + Ab(Xp0)*Ab(X11)*F0^-4 * (
          + 1/4*R(Xp0,X11,X11)*RC(X11)*X11
          )

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

       + Ab(Xp0)*Ab(X1s)*F0^-4 * (
          - 2/3*R(Xp0,X11,Xss)*RS(X11,Xss,Xp0)*X11
          )

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

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

       + Ab(Xp0)*Bb(X11,X11,0)*F0^-4 * (
          + 3/4*RC(X11)*X11^2
          - 1/2*RC(X11)^2*X11^2
          - 1/4*RC(X11)*RD(X11)*X11
          )

       + Ab(Xp0)*Cb(X11,X11,X11,0)*F0^-4 * (
          + 1/2*R(Xp0,X11,X11)*RD(X11)*X11^2
          )

       + Ab(Xp0,eps)*F0^-4*pi16 * (
          - 1/4*RC(X11)*Xp0*X11
          )

       + Ab(X11)*F0^-4*L0r * (
          - 12*RC(X11)*X11^2
          - 12*RD(X11)*X11
          )

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

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

       + Ab(X11)*F0^-4*L3r * (
          - 12*RC(X11)*X11^2
          - 12*RD(X11)*X11
          )

       + Ab(X11)*F0^-4*L4r * (
          + 16*RC(X11)*X11*X34
          )

       + Ab(X11)*F0^-4*L5r * (
          + 16*RC(X11)*X11^2
          + 8*RD(X11)*X11
          )

       + Ab(X11)*F0^-4*L6r * (
          + 32*X11^2
          - 16*RC(X11)*X11*X34
          )

       + Ab(X11)*F0^-4*L7r * (
          + 32*RD(X11)*X11
          )

       + Ab(X11)*F0^-4*L8r * (
          - 32*RC(X11)*X11^2
          )

       + Ab(X11)^2*F0^-4 * (
          + 1/2*X11
          + 1/8*RC(X11)^2*X11
          )

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

       + Ab(X11)*Ab(X1s)*F0^-4 * (
          + 2/3*R(Xp0,X11,Xss)*RS(X11,Xss,Xp0)*X11
          )

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

       + Ab(X11)*Bb(X11,X11,0)*F0^-4 * (
          + X11^2
          + 1/2*RC(X11)^2*X11^2
          + 1/4*RC(X11)*RD(X11)*X11
          )

       + Ab(X11)*Cb(X11,X11,X11,0)*F0^-4 * (
          + 1/2*RC(X11)*RD(X11)*X11^2
          )

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

       + Ab(X13)*Ab(X14)*F0^-4 * (
          - 1/2*X11
          )

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

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

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

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

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

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

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

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

       + Ab(X1s)*F0^-4*pi16 * (
          + 1/2*X11*X34
          + 1/6*X11*Xss
          )

       + Ab(X1s)*F0^-4*L0r * (
          + 8*X11*X1s
          )

       + Ab(X1s)*F0^-4*L3r * (
          + 20*X11*X1s
          )

       + Ab(X1s)*F0^-4*L5r * (
          - 16*X11*X1s
          )

       + Ab(X1s)*F0^-4*L8r * (
          + 32*X11*X1s
          )

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

       + Ab(X1s)*Bb(X11,Xp0,0)*F0^-4 * (
          - 2/3*R(Xp0,X11,Xss)*X11*X1s
          - 1/3*R(Xp0,X11,Xss)*X11^2
          )

       + Ab(X1s)*Bb(X11,X11,0)*F0^-4 * (
          - R(X11,Xss,Xp0)*X11*Xss
          )

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

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

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

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

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

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

       + Bb(X11,Xp0,0)*F0^-4*L7r * (
          + 32*RD(Xp0,X33,X44,X11)*RD(X11)*X11
          )

       + Bb(X11,Xp0,0)*F0^-4*L8r * (
          + 16*RD(Xp0,X33,X44,X11)*RD(X11)*X11
          )

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

       + Bb(X11,X11,0)*F0^-4*L0r * (
          - 12*RD(X11)*X11^2
          )

       + Bb(X11,X11,0)*F0^-4*L3r * (
          - 12*RD(X11)*X11^2
          )

       + Bb(X11,X11,0)*F0^-4*L4r * (
          + 8*RC(X11)*X11^2*X34
          + 24*RD(X11)*X11*X34
          )

       + Bb(X11,X11,0)*F0^-4*L5r * (
          + 4*RC(X11)*X11^3
          + 24*RD(X11)*X11^2
          )

       + Bb(X11,X11,0)*F0^-4*L6r * (
          - 16*RC(X11)*X11^2*X34
          - 32*RD(X11)*X11*X34
          )

       + Bb(X11,X11,0)*F0^-4*L7r * (
          + 16*RD(X11)^2*X11
          )

       + Bb(X11,X11,0)*F0^-4*L8r * (
          - 8*RC(X11)*X11^3
          - 48*RD(X11)*X11^2
          + 8*RD(X11)^2*X11
          )

       + Bb(X11,X11,0)^2*F0^-4 * (
          + 1/2*RC(X11)*RD(X11)*X11^2
          + 1/8*RD(X11)^2*X11
          )

       + Bb(X11,X11,0)*Cb(X11,X11,X11,0)*F0^-4 * (
          + 1/2*RD(X11)^2*X11^2
          )

       + Bb(X11,X11,0,eps)*F0^-4*pi16 * (
          + 1/4*RD(X11)*X11^2
          )

       + Cb(X11,X11,X11,0)*F0^-4*L4r * (
          + 16*RD(X11)*X11^2*X34
          )

       + Cb(X11,X11,X11,0)*F0^-4*L5r * (
          + 8*RD(X11)*X11^3
          )

       + Cb(X11,X11,X11,0)*F0^-4*L6r * (
          - 32*RD(X11)*X11^2*X34
          )

       + Cb(X11,X11,X11,0)*F0^-4*L8r * (
          - 16*RD(X11)*X11^3
          )

       + Hbb(1,Xp0,Xp0,X11,X11)*F0^-4 * (
          + 1/2*R(Xp0,X11,X11)^2*X11^2
          )

       + Hbb(1,Xp0,X11,X11,X11)*F0^-4 * (
          + R(Xp0,X11,X11)*RC(X11)*X11^2
          )

       + Hbb(1,Xp0,X1s,X1s,X11)*F0^-4 * (
          - 1/2*R(Xp0,X11,X11)*X11^2
          - 1/8*RV(Xp0,X11,Xss)*Xp0*X11
          )

       + Hbb(1,X11,X11,X11,X11)*F0^-4 * (
          + 1/3*X11^2
          + 1/2*RC(X11)^2*X11^2
          )

       + Hbb(1,X11,X1s,X1s,X11)*F0^-4 * (
          - 3/4*R(X11,Xss,Xp0)*X11^2
          + 3/8*RC(X11)*X11^2
          - 1/8*RD(X11)*X11
          )

       + Hbb(1,X13,X14,X34,X11)*F0^-4 * (
          + 1/2*X11*X34
          )

       + Hbb(2,X11,Xp0,X11,X11)*F0^-4 * (
          + R(Xp0,X11,X11)*RD(X11)*X11^2
          )

       + Hbb(2,X11,X11,X11,X11)*F0^-4 * (
          + RC(X11)*RD(X11)*X11^2
          )

       + Hbb(2,X11,X1s,X1s,X11)*F0^-4 * (
          + 3/8*RD(X11)*X11^2
          )

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

       + H1bb(1,Xp0,X1s,X1s,X11)*F0^-4 * (
          - R(Xp0,X11,Xss)*RS(X11,Xss,Xp0)*X11^2
          )

       + H1bb(1,X1s,X1s,X11,X11)*F0^-4 * (
          - 2*R(Xp0,X11,Xss)*RS(X11,Xss,Xp0)*X11^2
          )

       + H1bb(3,X1s,X11,X1s,X11)*F0^-4 * (
          - 2*RD(X11)*X11^2
          )

       + H21bb(1,Xp0,X1s,X1s,X11)*F0^-4 * (
          - 3/8*RV(Xp0,X11,Xss)*X11^2
          )

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

       + H21bb(1,X34,X13,X14,X11)*F0^-4 * (
          + 3/2*X11^2
          )

       + H21bb(2,X11,X1s,X1s,X11)*F0^-4 * (
          - 3/8*RD(X11)*X11^2
          );