*FORM by J.Vermaseren,version 3.1(Oct 17 2002) Run at: Mon Oct  4 00:50:30 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 "4"
*    #define type "1"
*    
*    #-
* Sigma2 loaded
* Sigma4 loaded
* Sigma6 loaded
*
*Time =       0.05 sec    Generated terms =          1
*             F0p2        Terms in output =          1
*                         Bytes used      =         24
*
*Time =       0.05 sec    Generated terms =          3
*             F0p4        Terms in output =          3
*                         Bytes used      =         98
*
*Time =       0.05 sec    Generated terms =         98
*             F0p6        Terms in output =         98
*                         Bytes used      =       4052

L   F0p2timo12 =

       + F0 * (
          + 1
          );

L   F0p4timo12 =

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

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

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

L   F0p6timo12 =

       + F0^-3*pi16*L0r * (
          + 4/9*Xe0*X44
          + 1/2*X11^2
          - 13/3*Xb(1)*X11
          - 35/18*Xb(2)
          )

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

       + F0^-3*pi16*L2r * (
          - 11/3*Xe0*X44
          - X11^2
          - 13/3*Xb(2)
          )

       + F0^-3*pi16*L3r * (
          + 4/9*Xe0*X44
          + 5/4*X11^2
          - 17/6*Xb(1)*X11
          - 43/36*Xb(2)
          )

       + F0^-3*pi16^2 * (
          - 15/64*Xe0*X44
          + 1/64*X11^2
          - 1/2*Xb(1)*X11
          - 43/128*Xb(2)
          )

       + F0^-3*L4r*L5r * (
          - 48*Xb(1)*X11
          )

       + F0^-3*L4r^2 * (
          - 72*Xb(1)^2
          )

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

       + F0^-3*K19r * (
          + 8*X11^2
          )

       + F0^-3*K20r * (
          + 24*Xb(1)*X11
          )

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

       + F0^-3*K22r * (
          + 72*Xb(1)^2
          )

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

       + Ab(Xe0)*F0^-3*L0r * (
          + 4*R(Xe0,X11,X11)*Xe0
          )

       + Ab(Xe0)*F0^-3*L1r * (
          - 8*Xe0
          )

       + Ab(Xe0)*F0^-3*L2r * (
          - 2*Xe0
          )

       + Ab(Xe0)*F0^-3*L3r * (
          + 4*R(Xe0,X11,X11)*Xe0
          )

       + Ab(Xe0)*F0^-3*L4r * (
          + 4*Xe0
          )

       + Ab(Xe0)*F0^-3*L5r * (
          - 4/3*R(Xe0,X11,X11)*X11
          )

       + Ab(Xe0)*Bb(X1s,X1s,0)*F0^-3 * (
          - 1/6*R(Xe0,X11,Xss)*X1s
          )

       + Ab(Xe0,eps)*F0^-3*pi16 * (
          + 1/8*RC(X11)*Xe0
          )

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

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

       + Ab(X11)*F0^-3*L2r * (
          - 10*X11
          )

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

       + Ab(X11)*F0^-3*L5r * (
          - 4/3*RC(X11)*X11
          )

       + Ab(X11)*Bb(X1s,X1s,0)*F0^-3 * (
          - 1/6*R(X11,Xss,Xe0)*X1s
          )

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

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

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

       + Ab(X44)*F0^-3*L4r * (
          + 12*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,eps)*F0^-3*pi16 * (
          + 1/2*X46
          )

       + Ab(X1s)*F0^-3*pi16 * (
          - 1/4*X1s
          - 3/8*Xb(1)
          )

       + Ab(X1s)*F0^-3*L0r * (
          - 4*X1s
          )

       + Ab(X1s)*F0^-3*L3r * (
          - 10*X1s
          )

       + Ab(X1s)*F0^-3*L4r * (
          - 6*Xb(1)
          )

       + Ab(X1s)*F0^-3*L5r * (
          + 2*X11
          )

       + Ab(X1s,eps)*F0^-3*pi16 * (
          + 3/8*Xss
          + 3/8*Xb(1)
          )

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

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

       + Bb(X11,X11,0)*F0^-3*L5r * (
          - 4/3*RD(X11)*X11
          )

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

       + Bb(X1s,X1s,0)*F0^-3*L4r * (
          - 12*Xb(1)*X1s
          )

       + Bb(X1s,X1s,0)*F0^-3*L5r * (
          - 4*X1s^2
          )

       + Bb(X1s,X1s,0)*F0^-3*L6r * (
          + 24*Xb(1)*X1s
          )

       + Bb(X1s,X1s,0)*F0^-3*L8r * (
          + 8*X1s^2
          )

       + Hbb(1,Xe0,X1s,X1s,X11)*F0^-3 * (
          + 1/24*RV(Xe0,X11,Xss)*Xe0
          )

       + Hbb(1,X11,X1s,X1s,X11)*F0^-3 * (
          - 1/12*R(X11,Xss,Xe0)*X11
          + 1/24*RC(X11)*X11
          + 1/24*RD(X11)
          )

       + Hbb(1,X14,X14,X44,X11)*F0^-3 * (
          - 3/8*X44
          )

       + Hbb(1,X14,X16,X46,X11)*F0^-3 * (
          - 1/2*X46
          )

       + Hbb(2,X11,X1s,X1s,X11)*F0^-3 * (
          + 1/24*RD(X11)*X11
          )

       + dHbb(1,Xe0,Xe0,X11,X11)*F0^-3 * (
          + 1/9*R(Xe0,X11,X11)^2*X11^2
          )

       + dHbb(1,Xe0,X11,X11,X11)*F0^-3 * (
          + 2/9*R(Xe0,X11,X11)*RC(X11)*X11^2
          )

       + dHbb(1,Xe0,X1s,X1s,X11)*F0^-3 * (
          - 1/6*R(Xe0,X11,X11)*X11^2
          - 1/24*RV(Xe0,X11,Xss)*Xe0*X11
          )

       + dHbb(1,X11,X11,X11,X11)*F0^-3 * (
          + 1/6*X11^2
          + 1/9*RC(X11)^2*X11^2
          )

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

       + dHbb(1,X14,X14,X44,X11)*F0^-3 * (
          + 3/8*X11*X44
          )

       + dHbb(1,X14,X16,X46,X11)*F0^-3 * (
          + 1/2*X11*X46
          )

       + dHbb(2,X11,Xe0,X11,X11)*F0^-3 * (
          + 2/9*R(Xe0,X11,X11)*RD(X11)*X11^2
          )

       + dHbb(2,X11,X11,X11,X11)*F0^-3 * (
          + 2/9*RC(X11)*RD(X11)*X11^2
          )

       + dHbb(2,X11,X1s,X1s,X11)*F0^-3 * (
          + 1/8*RD(X11)*X11^2
          )

       + dHbb(5,X11,X11,X11,X11)*F0^-3 * (
          + 1/9*RD(X11)^2*X11^2
          )

       + dH1bb(1,Xe0,X1s,X1s,X11)*F0^-3 * (
          - 1/3*R(Xe0,X11,Xss)*RS(X11,Xss,Xe0)*X11^2
          )

       + dH1bb(1,X1s,X1s,X11,X11)*F0^-3 * (
          - 2/3*R(Xe0,X11,Xss)*RS(X11,Xss,Xe0)*X11^2
          )

       + dH1bb(3,X1s,X11,X1s,X11)*F0^-3 * (
          - 2/3*RD(X11)*X11^2
          )

       + dH21bb(1,Xe0,X1s,X1s,X11)*F0^-3 * (
          - 1/8*RV(Xe0,X11,Xss)*X11^2
          )

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

       + dH21bb(1,X44,X14,X14,X11)*F0^-3 * (
          + 9/8*X11^2
          )

       + dH21bb(1,X46,X14,X16,X11)*F0^-3 * (
          + 3/2*X11^2
          )

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