*FORM by J.Vermaseren,version 3.1(Oct 17 2002) Run at: Thu Nov 25 06:39:59 2004
*    * This program loads the result for the finite part of the meson mass
*    * produced by PQint.frm and simplifies it as much as possible. Should
*    * also include a facility to print out the result it fortran format.
*    
*    * NB! Definition of mode must be as in the save file:
*    
*    * Mode 1 : full 3+3 theory, all quarks nondegenerate
*    * Mode 2 : 1+1 theory, all valence and 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
*    
*    * Choice of quark indices to be assigned to the external legs:
*    * The choices available are (in the notation ijkl):
*    
*    * Type 1 : pi^+/- meson - 1221
*    * Type 2 : K^+/- meson - 1331
*    
*    #define comp "2"
*    #define mode "2"
*    #define type "2"
*    
*    #-
* Sigma2 loaded
* Sigma4 loaded
* Sigma6 loaded
*
*Time =       0.24 sec    Generated terms =          1
*             F0p2        Terms in output =          1
*                         Bytes used      =         24
*
*Time =       0.24 sec    Generated terms =          6
*             F0p4        Terms in output =          6
*                         Bytes used      =        196
*
*Time =       0.24 sec    Generated terms =        182
*             F0p6        Terms in output =        182
*                         Bytes used      =       7022

L   F0p2timo21 =

       + F0 * (
          + 1
          );

L   F0p4timo21 =

       + F0^-1*L4r * (
          + 12*X44
          )

       + F0^-1*L5r * (
          + 4*X13
          )

       + Ab(Xpp)*F0^-1 * (
          - 1/12
          + 1/6*R(Xpp,Xqq)
          )

       + Ab(Xp4)*F0^-1 * (
          + 3/4
          )

       + Bb(Xpp,Xpp,0)*F0^-1 * (
          - 1/12*RD(Xpp)
          );

L   F0p6timo21 =

       + F0^-3*pi16*L0r * (
          - 1/2*X11*X33
          - 13/3*X13*X44
          + X13^2
          - 3/2*X44^2
          )

       + F0^-3*pi16*L1r * (
          - 2*X13^2
          )

       + F0^-3*pi16*L2r * (
          - X13^2
          - 8*X44^2
          )

       + F0^-3*pi16*L3r * (
          - 7/12*X11*X33
          - 17/6*X13*X44
          + 11/6*X13^2
          - 3/4*X44^2
          )

       + F0^-3*pi16^2 * (
          - 59/384*X11*X33
          - 1/2*X13*X44
          + 65/384*X13^2
          - 73/128*X44^2
          )

       + F0^-3*L4r*L5r * (
          - 48*X13*X44
          )

       + F0^-3*L4r^2 * (
          - 72*X44^2
          )

       + F0^-3*L5r^2 * (
          - 8*X13^2
          )

       + F0^-3*K19r * (
          + 4*Xpp^2
          )

       + F0^-3*K20r * (
          + 24*X13*X44
          )

       + F0^-3*K21r * (
          + 24*X44^2
          )

       + F0^-3*K22r * (
          + 72*X44^2
          )

       + F0^-3*K23r * (
          + 8*X11*X33
          )

       + Ab(Xpp)*F0^-3*pi16 * (
          + 5/96*Xpp
          + 1/32*Xqq
          + 1/48*X44
          - 1/12*R(Xpp,Xqq)*Xpp
          - 1/8*R(Xpp,Xqq)*X44
          )

       + Ab(Xpp)*F0^-3*L0r * (
          + 2/3*Xpp
          + 8/3*R(Xpp,Xqq)*Xpp
          + 2/3*RD(Xpp)
          )

       + Ab(Xpp)*F0^-3*L3r * (
          + 5/3*Xpp
          + 2/3*R(Xpp,Xqq)*Xpp
          + 5/3*RD(Xpp)
          )

       + Ab(Xpp)*F0^-3*L4r * (
          + X44
          - 2*R(Xpp,Xqq)*X44
          )

       + Ab(Xpp)*F0^-3*L5r * (
          + 1/6*Xpp
          - 1/6*Xqq
          - 1/3*X44
          - 2/3*R(Xpp,Xqq)*Xpp
          )

       + Ab(Xpp)^2*F0^-3 * (
          + 19/288
          - 1/72*R(Xpp,Xqq)*R(Xqq,Xpp)
          )

       + Ab(Xpp)*Ab(Xp4)*F0^-3 * (
          - 1/16
          - 1/12*R(Xpp,Xqq)
          )

       + Ab(Xpp)*Ab(Xq4)*F0^-3 * (
          - 3/16
          - 1/12*R(Xpp,Xqq)
          )

       + Ab(Xpp)*Ab(X13)*F0^-3 * (
          + 1/8
          )

       + Ab(Xpp)*Bb(Xpp,Xpp,0)*F0^-3 * (
          + 11/36*Xpp
          - 1/18*R(Xpp,Xqq)*Xpp
          - 1/72*R(Xpp,Xqq)*RD(Xpp)
          + 1/144*RD(Xpp)
          )

       + Ab(Xpp)*Bb(Xqq,Xqq,0)*F0^-3 * (
          - 1/72*R(Xpp,Xqq)*RD(Xqq)
          + 1/144*RD(Xqq)
          )

       + Ab(Xpp)*Bb(Xp4,Xp4,0)*F0^-3 * (
          - 1/4*Xp4
          )

       + Ab(Xpp)*Bb(X11,X33,0)*F0^-3 * (
          - 1/18*R(Xqq,Xpp)*Xpp
          )

       + Ab(Xpp)*Cb(Xpp,Xpp,Xpp,0)*F0^-3 * (
          + 1/18*RD(Xpp)*Xpp
          )

       + Ab(Xpp,eps)*F0^-3*pi16 * (
          - 1/8*Xpp
          + 1/8*R(Xpp,Xqq)*X44
          )

       + Ab(Xp4)*F0^-3*pi16 * (
          + 3/16*Xp4
          - 9/16*Xq4
          - 9/16*X44
          )

       + Ab(Xp4)*F0^-3*L0r * (
          - 6*Xp4
          )

       + Ab(Xp4)*F0^-3*L3r * (
          - 15*Xp4
          )

       + Ab(Xp4)*F0^-3*L4r * (
          - 9*X44
          )

       + Ab(Xp4)*F0^-3*L5r * (
          + 3*X13
          )

       + Ab(Xp4)^2*F0^-3 * (
          - 9/32
          )

       + Ab(Xp4)*Bb(Xpp,Xpp,0)*F0^-3 * (
          + 1/8*Xpp
          - 5/8*Xp4
          )

       + Ab(Xp4)*Bb(Xqq,Xqq,0)*F0^-3 * (
          - 1/16*RD(Xqq)
          )

       + Ab(Xp4)*Bb(X11,X33,0)*F0^-3 * (
          + 1/6*X44
          )

       + Ab(Xp4)*Bb(X11,X33,0,kk)*F0^-3 * (
          + 1/3
          )

       + Ab(Xp4,eps)*F0^-3*pi16 * (
          + 9/8*X44
          )

       + Ab(X11)*Ab(X33)*F0^-3 * (
          - 1/144
          + 1/36*R(X11,X33)*R(X33,X11)
          )

       + 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 * (
          + 9/16
          )

       + Ab(X44)*F0^-3*L1r * (
          - 64*X44
          )

       + Ab(X44)*F0^-3*L2r * (
          - 16*X44
          )

       + Ab(X44)*F0^-3*L4r * (
          + 32*X44
          )

       + Ab(X44)*Bb(Xpp,Xpp,0)*F0^-3 * (
          + 2/9*X44
          )

       + Ab(X44)*Bb(X11,X33,0)*F0^-3 * (
          - 4/9*X44
          )

       + Ab(X44,eps)*F0^-3*pi16 * (
          + X44
          )

       + Bb(Xpp,Xpp,0)*F0^-3*pi16 * (
          + 1/96*RD(Xpp)*Xpp
          + 1/32*RD(Xpp)*Xqq
          + 1/16*RD(Xpp)*X44
          )

       + 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*Xpp*X44
          - 4*R(Xpp,Xqq)*Xpp*X44
          + 3*RD(Xpp)*X44
          )

       + Bb(Xpp,Xpp,0)*F0^-3*L5r * (
          - 2/3*Xpp*X44
          + 4/3*Xpp^2
          - 4/3*R(Xpp,Xqq)*Xpp*X13
          - 1/3*RD(Xpp)*X13
          )

       + Bb(Xpp,Xpp,0)*F0^-3*L6r * (
          + 4*X44^2
          - 8*R(Xqq,Xpp)*Xpp*X44
          )

       + Bb(Xpp,Xpp,0)*F0^-3*L7r * (
          + 4*RD(Xpp)^2
          )

       + Bb(Xpp,Xpp,0)*F0^-3*L8r * (
          + 4/3*X44^2
          - 8/3*R(Xqq,Xpp)*Xpp^2
          )

       + Bb(Xpp,Xpp,0)^2*F0^-3 * (
          + 1/18*R(Xqq,Xpp)*RD(Xpp)*Xpp
          + 1/288*RD(Xpp)^2
          )

       + Bb(Xpp,Xpp,0)*Bb(X11,X33,0)*F0^-3 * (
          - 1/18*R(Xqq,Xpp)*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/8*RD(Xpp)*Xp4
          )

       + Bb(Xp4,Xp4,0)*F0^-3*L4r * (
          - 18*Xp4*X44
          )

       + Bb(Xp4,Xp4,0)*F0^-3*L5r * (
          - 6*Xp4^2
          )

       + Bb(Xp4,Xp4,0)*F0^-3*L6r * (
          + 36*Xp4*X44
          )

       + Bb(Xp4,Xp4,0)*F0^-3*L8r * (
          + 12*Xp4^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*RD(Xpp)*Xpp*X44
          )

       + Cb(Xpp,Xpp,Xpp,0)*F0^-3*L5r * (
          + 4/3*RD(Xpp)*Xpp^2
          )

       + Cb(Xpp,Xpp,Xpp,0)*F0^-3*L6r * (
          - 8*RD(Xpp)*Xpp*X44
          )

       + 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
          - 5/72*X13
          + 1/36*R(Xpp,Xqq)*R(Xqq,Xpp)*X13
          )

       + Hbb(1,Xpp,X14,X34,X13)*F0^-3 * (
          - 1/16*Xpp
          + 1/16*Xqq
          - 1/8*R(Xpp,Xqq)*X44
          )

       + Hbb(1,X11,X13,X33,X13)*F0^-3 * (
          + 1/72*X13
          - 1/18*R(X11,X33)*R(X33,X11)*X13
          )

       + Hbb(1,X13,X13,X13,X13)*F0^-3 * (
          - 1/8*X13
          )

       + Hbb(1,X14,X34,X44,X13)*F0^-3 * (
          - X44
          )

       + Hbb(2,Xpp,Xpp,X13,X13)*F0^-3 * (
          + 1/36*R(Xpp,Xqq)*RD(Xpp)*X13
          - 1/72*RD(Xpp)*X13
          )

       + Hbb(2,Xpp,X13,Xqq,X13)*F0^-3 * (
          + 1/36*R(Xqq,Xpp)*RD(Xpp)*X13
          - 1/72*RD(Xpp)*X13
          )

       + Hbb(2,Xpp,X14,X34,X13)*F0^-3 * (
          + 1/8*RD(Xpp)*Xp4
          )

       + 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,Xpp,X13,X13)*F0^-3 * (
          - 1/8*Xpp*X13
          - 5/24*X13^2
          + 1/12*R(Xpp,Xqq)*X13^2
          - 1/36*R(Xpp,Xqq)^2*X13^2
          )

       + dHbb(1,Xpp,X14,X34,X13)*F0^-3 * (
          + 3/16*Xpp*X13
          - 1/16*X13*X44
          - 1/8*R(Xpp,Xqq)*RD(Xpp)*X13
          + 3/8*R(Xqq,Xpp)*X13^2
          )

       + dHbb(1,X11,X13,X33,X13)*F0^-3 * (
          - 1/72*X13^2
          + 1/6*R(X11,X33)*R(X33,X11)*X13^2
          )

       + dHbb(1,X13,X13,X13,X13)*F0^-3 * (
          + 1/8*X13^2
          )

       + dHbb(1,X14,X34,X44,X13)*F0^-3 * (
          + X13*X44
          )

       + dHbb(2,Xpp,Xpp,X13,X13)*F0^-3 * (
          - 1/36*R(Xpp,Xqq)*RD(Xpp)*X13^2
          + 5/72*RD(Xpp)*X13^2
          )

       + dHbb(2,Xpp,X13,Xqq,X13)*F0^-3 * (
          - 1/36*R(Xqq,Xpp)*RD(Xpp)*X13^2
          + 1/72*RD(Xpp)*X13^2
          )

       + dHbb(2,Xpp,X14,X34,X13)*F0^-3 * (
          - 1/8*RD(Xpp)*Xp4*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 * (
          + 10/9*X13^2
          - 2/9*R(Xpp,Xqq)*R(Xqq,Xpp)*X13^2
          )

       + dH1bb(1,Xpp,X14,X34,X13)*F0^-3 * (
          - R(Xqq,Xpp)*X13^2
          )

       + dH1bb(1,Xp4,Xq4,Xpp,X13)*F0^-3 * (
          - X13^2
          )

       + dH1bb(1,X13,X11,X33,X13)*F0^-3 * (
          + 1/9*X13^2
          - 2/9*R(X11,X33)*R(X33,X11)*X13^2
          )

       + dH1bb(3,X13,Xpp,Xpp,X13)*F0^-3 * (
          + 1/9*R(Xpp,Xqq)*RD(Xpp)*X13^2
          - 1/9*RD(Xpp)*X13^2
          )

       + dH1bb(3,X13,Xpp,Xqq,X13)*F0^-3 * (
          + 1/9*R(Xqq,Xpp)*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,X14,X34,X13)*F0^-3 * (
          + 3/8*R(Xqq,Xpp)*X13^2
          )

       + dH21bb(1,Xp4,Xpp,Xq4,X13)*F0^-3 * (
          + 3/8*R(Xpp,Xqq)*X13^2
          )

       + dH21bb(1,Xp4,Xqq,Xq4,X13)*F0^-3 * (
          - 3/8*R(Xpp,Xqq)*X13^2
          )

       + dH21bb(1,X13,Xpp,Xpp,X13)*F0^-3 * (
          + 5/24*X13^2
          - 1/12*R(Xpp,Xqq)*R(Xqq,Xpp)*X13^2
          )

       + dH21bb(1,X13,X11,X33,X13)*F0^-3 * (
          - 1/24*X13^2
          + 1/6*R(X11,X33)*R(X33,X11)*X13^2
          )

       + dH21bb(1,X13,X13,X13,X13)*F0^-3 * (
          + 3/8*X13^2
          )

       + dH21bb(1,X44,X14,X34,X13)*F0^-3 * (
          + 3*X13^2
          )

       + dH21bb(3,Xp4,Xpp,Xq4,X13)*F0^-3 * (
          - 3/8*RD(Xpp)*X13^2
          )

       + dH21bb(3,X13,Xpp,Xpp,X13)*F0^-3 * (
          - 1/12*R(Xpp,Xqq)*RD(Xpp)*X13^2
          + 1/24*RD(Xpp)*X13^2
          )

       + dH21bb(3,X13,Xpp,Xqq,X13)*F0^-3 * (
          - 1/12*R(Xqq,Xpp)*RD(Xpp)*X13^2
          + 1/24*RD(Xpp)*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
          );