*FORM by J.Vermaseren,version 3.1(Oct 17 2002) Run at: Mon Oct  4 00:49:13 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 "1"
*    
*    #-
* Sigma2 loaded
* Sigma4 loaded
* Sigma6 loaded
*
*Time =       0.01 sec    Generated terms =          1
*             F0p2        Terms in output =          1
*                         Bytes used      =         24
*
*Time =       0.01 sec    Generated terms =          3
*             F0p4        Terms in output =          3
*                         Bytes used      =         96
*
*Time =       0.01 sec    Generated terms =         63
*             F0p6        Terms in output =         63
*                         Bytes used      =       2290

L   F0p2timo11 =

       + F0 * (
          + 1
          );

L   F0p4timo11 =

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

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

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

L   F0p6timo11 =

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

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

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

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

       + F0^-3*pi16^2 * (
          - 1/2*X11*X44
          + 1/64*X11^2
          - 73/128*X44^2
          )

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

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

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

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

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

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

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

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

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

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

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

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

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

       + Ab(X11)*Bb(X14,X14,0)*F0^-3 * (
          - 1/2*X14
          )

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

       + Ab(X14)*F0^-3*pi16 * (
          - 3/4*X14
          - 9/8*X44
          )

       + Ab(X14)*F0^-3*L0r * (
          - 12*X14
          )

       + Ab(X14)*F0^-3*L3r * (
          - 30*X14
          )

       + Ab(X14)*F0^-3*L4r * (
          - 18*X44
          )

       + Ab(X14)*F0^-3*L5r * (
          + 6*X11
          )

       + Ab(X14,eps)*F0^-3*pi16 * (
          + 9/4*X44
          )

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

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

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

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

       + 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(X14,X14,0)*F0^-3*L4r * (
          - 36*X14*X44
          )

       + Bb(X14,X14,0)*F0^-3*L5r * (
          - 12*X14^2
          )

       + Bb(X14,X14,0)*F0^-3*L6r * (
          + 72*X14*X44
          )

       + Bb(X14,X14,0)*F0^-3*L8r * (
          + 24*X14^2
          )

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

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

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

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

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

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

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

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

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

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

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

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

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