double complex function Fcc_qpqmgpgp(j1,j2,j3,j4,j5,j6,za,zb) 
      implicit none
      integer j1,j2,j3,j4,j5,j6
      include 'constants.f'
      include 'zprods_decl.f'
      include 'sprods_com.f'
      double complex L0,Lsm1,Lsm1_2me
      double precision t  

      Fcc_qpqmgpgp=
     .((Lsm1(-s(j1,j4),-t(j1,j2,j4),-s(j1,j2),-t(j1,j2,j4))+
     .Lsm1_2me(t(j1,j3,j4),t(j1,j2,j4),s(j1,j4),s(j5,j6)))*za(j2,j5)**2)
     ./
     .(za(j1,j4)*za(j2,j3)*za(j3,j4)*za(j5,j6))-
     .(Lsm1(-s(j1,j2),-t(j1,j2,j3),-s(j2,j3),-t(j1,j2,j3))*za(j2,j5)*
     .(za(j1,j5)*za(j2,j3)-za(j1,j2)*za(j3,j5)))/
     .(za(j1,j3)*za(j1,j4)*za(j2,j3)*za(j3,j4)*za(j5,j6))+
     .(Lsm1_2me(t(j1,j2,j3),t(j2,j3,j4),s(j2,j3),s(j5,j6))*za(j2,j5)*
     .(-(za(j1,j5)*za(j2,j4))+za(j1,j2)*za(j4,j5)))/
     .(za(j1,j4)**2*za(j2,j3)*za(j3,j4)*za(j5,j6))+
     .(Lsm1_2me(t(j1,j2,j4),t(j1,j2,j3),s(j1,j2),s(j5,j6))*za(j2,j5)*
     .(za(j2,j4)*za(j3,j5)+za(j2,j3)*za(j4,j5)))/
     .(za(j1,j4)*za(j2,j3)*za(j3,j4)**2*za(j5,j6))-
     .(2d0*L0(-t(j2,j3,j4),-s(j5,j6))*za(j1,j5)*za(j2,j5)*
     .(-(za(j2,j3)*zb(j1,j3))-za(j2,j4)*zb(j1,j4)))/
     .(s(j5,j6)*za(j1,j4)*za(j2,j3)*za(j3,j4)*za(j5,j6))-
     .(2d0*L0(-t(j2,j3,j4),-s(j2,j3))*za(j2,j5)*za(j4,j5)*zb(j3,j4))/
     .(s(j2,j3)*za(j1,j4)*za(j3,j4)*za(j5,j6))

      return
      end