double complex function Fcc_qpqmgpgm(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_2mh,Lnrat,I3m
      double precision t  

      Fcc_qpqmgpgm=
     .Lsm1(-s(j1,j2),-t(j1,j2,j4),-s(j1,j4),-t(j1,j2,j4))*
     .(-((zb(j1,j2)**2*(-(za(j1,j5)*zb(j1,j4))-za(j2,j5)*zb(j2,j4))**2)/
     .(za(j5,j6)*zb(j1,j4)*zb(j2,j4)**2*
     .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))))+
     .(zb(j1,j4)*(-(za(j1,j5)*zb(j1,j2))+za(j4,j5)*zb(j2,j4))**2)/
     .(za(j5,j6)*zb(j2,j4)**2*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4
     .))*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))))-
     .(2d0*L0(-t(j2,j3,j4),-s(j5,j6))*za(j1,j5)*zb(j1,j3)*
     .(za(j1,j2)*zb(j2,j3)-za(j1,j4)*zb(j3,j4))*
     .(-(za(j2,j5)*zb(j2,j3))+za(j4,j5)*zb(j3,j4)))/
     .(s(j5,j6)*za(j5,j6)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*
     .zb(j3,j4)*(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4)))+
     .Lsm1(-s(j1,j2),-t(j1,j2,j3),-s(j2,j3),-t(j1,j2,j3))*
     .((za(j1,j2)**2*(-(za(j1,j3)*zb(j1,j6))-za(j2,j3)*zb(j2,j6))**2)/
     .(za(j1,j3)**2*za(j2,j3)*(-(za(j1,j3)*zb(j1,j4))-
     .za(j2,j3)*zb(j2,j4))*(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))*
     .zb(j5,j6))-(za(j2,j3)*
     .(za(j1,j2)*zb(j2,j6)+za(j1,j3)*zb(j3,j6))**2)/
     .(za(j1,j3)**2*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))*zb(j5,j6)))-
     .(2d0*L0(-t(j1,j3,j4),-s(j5,j6))*za(j2,j4)*
     .(-(za(j1,j4)*zb(j1,j2))+za(j3,j4)*zb(j2,j3))*zb(j2,j6)*
     .(-(za(j1,j4)*zb(j1,j6))-za(j3,j4)*zb(j3,j6)))/
     .(s(j5,j6)*za(j3,j4)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*zb(j5,j6))
     
      Fcc_qpqmgpgm=Fcc_qpqmgpgm+
     .I3m(s(j1,j4),s(j2,j3),s(j5,j6))*
     .((za(j1,j5)*za(j2,j4)*zb(j1,j3)*zb(j2,j6)+
     .((-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))*
     .(za(j1,j3)*za(j2,j5)*zb(j1,j6)*zb(j2,j4)+
     .za(j1,j3)*zb(j1,j4)*
     .(-(za(j2,j5)*zb(j2,j6))+za(j3,j5)*zb(j3,j6))))/
     .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))+
     .((-(za(j1,j4)*zb(j1,j3))+za(j2,j4)*zb(j2,j3))*
     .(-(za(j1,j5)*zb(j1,j6))-za(j4,j5)*zb(j4,j6)))/2)/
     .((-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4)))-
     .(za(j4,j5)*zb(j1,j3)*(-(za(j1,j2)*zb(j1,j6))+za(j2,j3)*zb(j3,j6)))
     ./
     .((-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))*t(j1,j2,j3)))+
     .Lsm1_2mh(s(j1,j4),t(j1,j2,j3),s(j2,j3),s(j5,j6))*
     .(-(((-(za(j1,j2)*zb(j1,j4))+za(j2,j3)*zb(j3,j4))**2*
     .(za(j1,j2)*zb(j2,j6)+za(j1,j3)*zb(j3,j6))**2)/
     .(za(j2,j3)*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))**3*zb(j5,j6)))+
     .(za(j1,j2)**2*zb(j4,j6)**2*t(j1,j2,j3)**2)/
     .(za(j2,j3)*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))**3*zb(j5,j6)))+
     .Lsm1_2mh(s(j3,j4),t(j1,j2,j3),s(j1,j2),s(j5,j6))*
     .(((-(za(j1,j3)*zb(j1,j6))-za(j2,j3)*zb(j2,j6))**2*
     .(-(za(j1,j2)*zb(j1,j4))+za(j2,j3)*zb(j3,j4))**2)/
     .(za(j2,j3)*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))**3*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))*zb(j5,j6))-
     .(za(j2,j3)*zb(j4,j6)**2*t(j1,j2,j3)**2)/
     .((-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))**3*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))*zb(j5,j6)))+
     .(2d0*Lnrat(-s(j2,j3),-s(j5,j6))*(-2d0*(s(j1,j2)-s(j3,j4))*za(j2,j4
     .)*zb(j2,j3)*
     .(-(za(j2,j5)*zb(j2,j6))-za(j3,j5)*zb(j3,j6))+
     .(za(j2,j5)*zb(j1,j3)*((s(j1,j2)-s(j3,j4))*
     .((-s(j1,j4)+s(j2,j3)-s(j5,j6))*za(j3,j5)*zb(j2,j3)+
     .(-s(j1,j4)-s(j2,j3)+s(j5,j6))*za(j5,j6)*zb(j2,j6))+
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*
     .(-((-s(j1,j4)+s(j2,j3)-s(j5,j6))*za(j2,j5)*zb(j2,j3))+
     .(-s(j1,j4)-s(j2,j3)+s(j5,j6))*za(j5,j6)*zb(j3,j6))))/
     .(za(j5,j6)*zb(j1,j4))-
     .(za(j2,j4)*(zb(j1,j3)*
     .(za(j1,j2)*zb(j2,j6)+za(j1,j3)*zb(j3,j6))-
     .zb(j3,j4)*(-(za(j2,j4)*zb(j2,j6))-za(j3,j4)*zb(j3,j6))-
     .zb(j2,j3)*(-(za(j1,j2)*zb(j1,j6))+za(j2,j4)*zb(j4,j6)))*
     .(za(j3,j5)*zb(j2,j3)*zb(j5,j6)+zb(j2,j6)*t(j1,j2,j3)))/zb(j5,j6)
     .))/((s(j1,j4)**2-2d0*s(j1,j4)*s(j2,j3)+s(j2,j3)**2-2d0*s(j1,j4)*s(
     .j5,j6)-
     .2d0*s(j2,j3)*s(j5,j6)+s(j5,j6)**2)*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4)))
      
      Fcc_qpqmgpgm=Fcc_qpqmgpgm+
     .I3m(s(j2,j3),s(j1,j4),s(j5,j6))*
     .((za(j1,j5)*za(j2,j4)*zb(j1,j3)*zb(j2,j6)+
     .((za(j1,j4)*zb(j1,j3)-za(j2,j4)*zb(j2,j3))*
     .(-(za(j2,j5)*zb(j2,j6))-za(j3,j5)*zb(j3,j6)))/2+
     .((-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))*
     .(za(j1,j3)*za(j2,j5)*zb(j1,j6)*zb(j2,j4)+
     .za(j2,j3)*zb(j2,j4)*
     .(-(za(j1,j5)*zb(j1,j6))+za(j4,j5)*zb(j4,j6))))/
     .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4)))/
     .((-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4)))-
     .(za(j2,j4)*(za(j2,j5)*zb(j1,j2)+za(j4,j5)*zb(j1,j4))*zb(j3,j6))/
     .((-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))*t(j1,j2,j4)))+
     .Lsm1_2mh(s(j2,j3),t(j1,j2,j4),s(j1,j4),s(j5,j6))*
     .(((za(j2,j3)*zb(j1,j2)-za(j3,j4)*zb(j1,j4))**2*
     .(-(za(j1,j5)*zb(j1,j2))+za(j4,j5)*zb(j2,j4))**2)/
     .(za(j5,j6)*zb(j1,j4)*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))**3)-
     .(za(j3,j5)**2*zb(j1,j2)**2*t(j1,j2,j4)**2)/
     .(za(j5,j6)*zb(j1,j4)*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))**3))+
     .Lsm1_2mh(s(j3,j4),t(j1,j2,j4),s(j1,j2),s(j5,j6))*
     .(-(((za(j2,j3)*zb(j1,j2)-za(j3,j4)*zb(j1,j4))**2*
     .(-(za(j1,j5)*zb(j1,j4))-za(j2,j5)*zb(j2,j4))**2)/
     .(za(j5,j6)*zb(j1,j4)*(-(za(j1,j3)*zb(j1,j4))-
     .za(j2,j3)*zb(j2,j4))**3*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))))+
     .(za(j3,j5)**2*zb(j1,j4)*t(j1,j2,j4)**2)/
     .(za(j5,j6)*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))**3*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))))+
     .(2d0*Lnrat(-s(j1,j4),-s(j5,j6))*(-2d0*(s(j1,j2)-s(j3,j4))*za(j1,j4
     .)*zb(j1,j3)*
     .(-(za(j1,j5)*zb(j1,j6))-za(j4,j5)*zb(j4,j6))-
     .(za(j2,j4)*zb(j1,j6)*((s(j1,j2)-s(j3,j4))*
     .((s(j1,j4)-s(j2,j3)-s(j5,j6))*za(j1,j4)*zb(j4,j6)-
     .(-s(j1,j4)-s(j2,j3)+s(j5,j6))*za(j1,j5)*zb(j5,j6))+
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))*
     .(-((s(j1,j4)-s(j2,j3)-s(j5,j6))*za(j1,j4)*zb(j1,j6))-
     .(-s(j1,j4)-s(j2,j3)+s(j5,j6))*za(j4,j5)*zb(j5,j6))))/
     .(za(j2,j3)*zb(j5,j6))+
     .(zb(j1,j3)*(-(za(j1,j4)*
     .(za(j2,j5)*zb(j1,j2)+za(j3,j5)*zb(j1,j3)))+
     .za(j2,j4)*(-(za(j1,j5)*zb(j1,j2))+za(j4,j5)*zb(j2,j4))+
     .za(j3,j4)*(-(za(j1,j5)*zb(j1,j3))+za(j4,j5)*zb(j3,j4)))*
     .(-(za(j1,j4)*za(j5,j6)*zb(j4,j6))+za(j1,j5)*t(j1,j2,j4)))/
     .za(j5,j6)))/
     .((s(j1,j4)**2-2d0*s(j1,j4)*s(j2,j3)+s(j2,j3)**2-2d0*s(j1,j4)*s(j5,
     .j6)-
     .2d0*s(j2,j3)*s(j5,j6)+s(j5,j6)**2)*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4)))
     
      Fcc_qpqmgpgm=Fcc_qpqmgpgm+
     .Lsm1(-s(j1,j4),-t(j1,j3,j4),-s(j3,j4),-t(j1,j3,j4))*
     .((za(j3,j4)**2*(za(j1,j3)*zb(j3,j6)+za(j1,j4)*zb(j4,j6))**2)/
     .(za(j1,j3)**3*(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*
     .zb(j5,j6)*t(j1,j3,j4))-
     .(za(j1,j4)**2*(-(za(j1,j3)*zb(j1,j6))+za(j3,j4)*zb(j4,j6))**2)/
     .(za(j1,j3)**3*(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*zb(j5,j
     .6)*
     .t(j1,j3,j4)))+(2d0*L0(-t(j1,j3,j4),-s(j3,j4))*za(j1,j4)*zb(j1,j3)*
     .(-(za(j1,j4)*zb(j1,j6))-za(j3,j4)*zb(j3,j6))*
     .(za(j1,j3)*zb(j3,j6)+za(j1,j4)*zb(j4,j6)))/
     .(s(j3,j4)*za(j1,j3)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*
     .zb(j5,j6)*t(j1,j3,j4))+(2d0*L0(-t(j1,j3,j4),-s(j1,j4))*za(j1,j4)*z
     .b(j1,j3)*
     .(-(za(j1,j4)*zb(j1,j6))-za(j3,j4)*zb(j3,j6))*
     .(-(za(j1,j3)*zb(j1,j6))+za(j3,j4)*zb(j4,j6)))/
     .(s(j1,j4)*za(j1,j3)*(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*
     .zb(j5,j6)*t(j1,j3,j4))+Lsm1_2mh(s(j1,j2),t(j1,j3,j4),s(j3,j4),
     .s(j5,j6))*(-((za(j2,j5)**2*zb(j1,j3)**3)/
     .(za(j5,j6)*zb(j1,j4)*(-(za(j2,j3)*zb(j1,j3))-
     .za(j2,j4)*zb(j1,j4))*zb(j3,j4)*t(j1,j3,j4)))-
     .((-(za(j1,j4)*zb(j1,j2))+za(j3,j4)*zb(j2,j3))**3*
     .(za(j1,j3)*zb(j3,j6)+za(j1,j4)*zb(j4,j6))**2)/
     .(za(j3,j4)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))**3*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*zb(j5,j6)*t(j1,j3,j4)
     .)+(za(j1,j4)**2*(-(za(j1,j4)*zb(j1,j2))+za(j3,j4)*zb(j2,j3))*zb(j2
     .,j6)**2*
     .t(j1,j3,j4))/
     .(za(j3,j4)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))**3*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*zb(j5,j6)))+
     .(I3m(s(j1,j4),s(j2,j3),s(j5,j6))*zb(j1,j3)*
     .(-(s(j5,j6)*za(j2,j4)*((-s(j1,j4)-s(j2,j3)+s(j5,j6))*za(j1,j5)*
     .zb(j2,j6)-2d0*za(j1,j4)*za(j3,j5)*zb(j2,j3)*zb(j4,j6)))+
     .2d0*za(j1,j4)*za(j1,j5)*za(j2,j5)*
     .((-s(j1,j4)-s(j2,j3)+s(j5,j6))*zb(j1,j2)-
     .2d0*za(j3,j4)*zb(j1,j4)*zb(j2,j3))*zb(j5,j6)-
     .2d0*za(j1,j4)*za(j4,j5)*zb(j4,j6)*
     .(s(j5,j6)*(-s(j1,j4)-s(j2,j3)+s(j5,j6))+
     .(s(j1,j4)-s(j2,j3)-s(j5,j6))*t(j1,j3,j4))))/
     .((s(j1,j4)**2-2d0*s(j1,j4)*s(j2,j3)+s(j2,j3)**2-2d0*s(j1,j4)*s(j5,
     .j6)-
     .2d0*s(j2,j3)*s(j5,j6)+s(j5,j6)**2)*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4)))
     
      Fcc_qpqmgpgm=Fcc_qpqmgpgm+
     .Lsm1(-s(j2,j3),-t(j2,j3,j4),-s(j3,j4),-t(j2,j3,j4))*
     .(-(((za(j3,j5)*zb(j2,j3)+za(j4,j5)*zb(j2,j4))**2*zb(j3,j4)**2)/
     .(za(j5,j6)*zb(j2,j4)**3*(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))*
     .t(j2,j3,j4)))+(zb(j2,j3)**2*
     .(-(za(j2,j5)*zb(j2,j4))-za(j3,j5)*zb(j3,j4))**2)/
     .(za(j5,j6)*zb(j2,j4)**3*(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))*
     .t(j2,j3,j4)))-(2d0*L0(-t(j2,j3,j4),-s(j3,j4))*za(j2,j4)*zb(j2,j3)*
     .(za(j3,j5)*zb(j2,j3)+za(j4,j5)*zb(j2,j4))*
     .(-(za(j2,j5)*zb(j2,j3))+za(j4,j5)*zb(j3,j4)))/
     .(s(j3,j4)*za(j5,j6)*zb(j2,j4)*
     .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*t(j2,j3,j4))-
     .(2d0*L0(-t(j2,j3,j4),-s(j2,j3))*za(j2,j4)*zb(j2,j3)*
     .(-(za(j2,j5)*zb(j2,j4))-za(j3,j5)*zb(j3,j4))*
     .(-(za(j2,j5)*zb(j2,j3))+za(j4,j5)*zb(j3,j4)))/
     .(s(j2,j3)*za(j5,j6)*zb(j2,j4)*(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3
     .,j4))*
     .t(j2,j3,j4))+Lsm1_2mh(s(j1,j2),t(j2,j3,j4),s(j3,j4),s(j5,j6))*
     .(-(((za(j3,j5)*zb(j2,j3)+za(j4,j5)*zb(j2,j4))**2*
     .(za(j1,j2)*zb(j2,j3)-za(j1,j4)*zb(j3,j4))**3)/
     .(za(j5,j6)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))**3*
     .zb(j3,j4)*(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))*t(j2,j3,j4))
     .)-(za(j2,j4)**3*zb(j1,j6)**2)/
     .(za(j2,j3)*za(j3,j4)*(-(za(j2,j3)*zb(j1,j3))-za(j2,j4)*zb(j1,j4))*
     .zb(j5,j6)*t(j2,j3,j4))+
     .(za(j1,j5)**2*zb(j2,j3)**2*(za(j1,j2)*zb(j2,j3)-za(j1,j4)*zb(j3,j4
     .))*
     .t(j2,j3,j4))/
     .(za(j5,j6)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))**3*zb(j3,j
     .4)*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))))+
     .(I3m(s(j2,j3),s(j1,j4),s(j5,j6))*za(j2,j4)*
     .(-2d0*za(j5,j6)*zb(j1,j6)*zb(j2,j3)*zb(j2,j6)*
     .(-((-s(j1,j4)-s(j2,j3)+s(j5,j6))*za(j1,j2))+
     .2d0*za(j1,j4)*za(j2,j3)*zb(j3,j4))-
     .s(j5,j6)*zb(j1,j3)*((-s(j1,j4)-s(j2,j3)+s(j5,j6))*za(j1,j5)*
     .zb(j2,j6)-2d0*za(j1,j4)*za(j3,j5)*zb(j2,j3)*zb(j4,j6))-
     .2d0*za(j3,j5)*zb(j2,j3)*zb(j3,j6)*
     .(s(j5,j6)*(-s(j1,j4)-s(j2,j3)+s(j5,j6))+
     .(-s(j1,j4)+s(j2,j3)-s(j5,j6))*t(j2,j3,j4))))/
     .((s(j1,j4)**2-2d0*s(j1,j4)*s(j2,j3)+s(j2,j3)**2-2d0*s(j1,j4)*s(j5,
     .j6)-
     .2d0*s(j2,j3)*s(j5,j6)+s(j5,j6)**2)*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4)))
     
      Fcc_qpqmgpgm=Fcc_qpqmgpgm+
     .I3m(s(j1,j2),s(j3,j4),s(j5,j6))*
     .((za(j3,j5)*zb(j2,j3)*(s(j1,j4)*za(j1,j4)*zb(j1,j6)-
     .s(j2,j4)*za(j2,j4)*zb(j2,j6))-
     .5d0*za(j1,j3)*za(j2,j4)*za(j3,j5)*zb(j1,j2)*zb(j2,j3)*zb(j3,j6)-
     .(-s(j1,j3)+s(j2,j4))*za(j3,j4)*
     .(-(za(j1,j5)*zb(j1,j2))+za(j4,j5)*zb(j2,j4))*zb(j3,j6)-
     .za(j1,j4)*zb(j1,j2)*(-(s(j3,j4)*za(j1,j5)*zb(j1,j6))-
     .2d0*(-((s(j1,j3)+s(j2,j3))*za(j2,j5)*zb(j2,j6))-
     .za(j2,j3)*za(j4,j5)*zb(j2,j6)*zb(j3,j4))-
     .s(j2,j4)*za(j3,j5)*zb(j3,j6)+
     .(2d0*za(j2,j3)*zb(j1,j2)*zb(j3,j4)*
     .(-(za(j2,j5)*zb(j2,j6))-za(j4,j5)*zb(j4,j6)))/zb(j1,j4)+
     .3d0*(za(j2,j5)*za(j3,j4)*zb(j2,j4)*zb(j3,j6)+
     .s(j1,j3)*za(j4,j5)*zb(j4,j6))))/
     .(2d0*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*
     .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4)))-
     .(za(j2,j4)*zb(j1,j2)*(za(j2,j5)*zb(j1,j2)+za(j4,j5)*zb(j1,j4))*
     .zb(j3,j6))/
     .(zb(j1,j4)*(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))*t(j1,j2,j4
     .))
     .-(za(j2,j5)*zb(j1,j3)**2*(-(za(j1,j4)*zb(j1,j6))-za(j3,j4)*zb(j3,j
     .6)))/
     .(zb(j1,j4)*t(j1,j3,j4)**2)+
     .(za(j2,j5)*zb(j1,j3)**2*(-(za(j2,j5)*zb(j2,j3))+
     .za(j4,j5)*zb(j3,j4)))/
     .(2d0*za(j5,j6)*zb(j1,j4)*zb(j3,j4)*t(j1,j3,j4))+
     .(za(j2,j4)*zb(j1,j3)*zb(j1,j6)*
     .(-(za(j1,j4)*zb(j1,j6))-za(j3,j4)*zb(j3,j6)))/
     .(2d0*za(j3,j4)*zb(j1,j4)*zb(j5,j6)*t(j1,j3,j4))+
     .((za(j1,j5)**2*zb(j1,j3)**3*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4)))/
     .(za(j5,j6)*zb(j1,j4)*zb(j3,j4))+
     .(zb(j1,j2)*((-s(j1,j2)+2d0*s(j1,j3)+s(j3,j4)-s(j5,j6))*
     .za(j4,j5)-2d0*za(j2,j4)*za(j5,j6)*zb(j2,j6))*zb(j3,j6))/
     .zb(j1,j4)+(zb(j1,j3)*
     .(-(za(j1,j3)*za(j1,j5)*zb(j1,j2)*zb(j1,j3))+
     .za(j1,j3)*za(j2,j5)*zb(j1,j2)*zb(j2,j3)-
     .za(j1,j5)*za(j3,j4)*zb(j1,j3)*zb(j2,j4)+
     .za(j2,j5)*za(j3,j4)*zb(j2,j3)*zb(j2,j4))*zb(j3,j6))/
     .(zb(j1,j4)*zb(j3,j4))+
     .(za(j2,j4)**2*(-((s(j1,j2)-s(j3,j4)-s(j5,j6))*zb(j1,j2))+
     .za(j1,j3)*zb(j1,j2)*zb(j1,j3)+
     .za(j3,j4)*zb(j1,j3)*zb(j2,j4))*zb(j2,j6)**2)/
     .(za(j3,j4)*zb(j1,j4)*zb(j5,j6))+
     .(4d0*za(j1,j4)*za(j2,j5)*zb(j1,j3)*
     .(-(za(j1,j4)*zb(j1,j2))+za(j3,j4)*zb(j2,j3))*zb(j2,j6))/
     .t(j1,j3,j4)+((-(za(j1,j4)*zb(j1,j2))+za(j3,j4)*zb(j2,j3))*
     .(-(za(j1,j4)*zb(j1,j6))-za(j3,j4)*zb(j3,j6))**2*
     .(2d0*s(j3,j4)*s(j5,j6)+(s(j1,j2)-s(j3,j4)-s(j5,j6))*t(j1,j3,j4))
     .)/(za(j3,j4)*zb(j5,j6)*t(j1,j3,j4)**2)+
     .(za(j4,j5)*zb(j1,j2)*(za(j1,j4)*zb(j1,j6)-za(j2,j4)*zb(j2,j6))*
     .(s(j1,j4)+t(j2,j3,j4)))/(za(j3,j4)*zb(j1,j4)))/
     .(2d0*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*
     .(-(za(j1,j3)*zb(j1,j2))-za(j3,j4)*zb(j2,j4))))
     
      Fcc_qpqmgpgm=Fcc_qpqmgpgm+
     .I3m(s(j1,j2),s(j3,j4),s(j5,j6))*
     .(((s(j1,j3)-s(j2,j4))*za(j4,j5)*zb(j3,j4)*
     .(za(j1,j2)*zb(j2,j6)+za(j1,j3)*zb(j3,j6))+
     .za(j1,j4)*(-(s(j1,j3)*za(j1,j5)*zb(j1,j3))+
     .s(j2,j3)*za(j2,j5)*zb(j2,j3))*zb(j4,j6)+
     .5d0*za(j1,j2)*za(j1,j4)*za(j4,j5)*zb(j1,j3)*zb(j2,j4)*zb(j4,j6)+
     .za(j1,j2)*zb(j2,j3)*(-(s(j3,j4)*za(j2,j5)*zb(j2,j6))+
     .(2d0*za(j1,j2)*za(j3,j4)*zb(j1,j4)*
     .(-(za(j1,j5)*zb(j1,j6))-za(j3,j5)*zb(j3,j6)))/za(j2,j3)+
     .3d0*(-(za(j1,j3)*za(j4,j5)*zb(j1,j6)*zb(j3,j4))+
     .s(j2,j4)*za(j3,j5)*zb(j3,j6))-
     .2d0*(-((s(j1,j4)+s(j2,j4))*za(j1,j5)*zb(j1,j6))+
     .za(j1,j5)*za(j3,j4)*zb(j1,j4)*zb(j3,j6))-
     .s(j1,j3)*za(j4,j5)*zb(j4,j6)))/
     .(2d0*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*
     .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4)))+
     .(za(j1,j2)*za(j4,j5)*zb(j1,j3)*
     .(-(za(j1,j2)*zb(j1,j6))+za(j2,j3)*zb(j3,j6)))/
     .(za(j2,j3)*(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))*t(j1,j2,j3))-
     .(za(j2,j4)**2*zb(j1,j6)*(-(za(j2,j5)*zb(j2,j3))+za(j4,j5)*zb(j3,j4
     .)))/
     .(za(j2,j3)*t(j2,j3,j4)**2)+
     .(za(j2,j4)*za(j2,j5)*zb(j1,j3)*
     .(-(za(j2,j5)*zb(j2,j3))+za(j4,j5)*zb(j3,j4)))/
     .(2d0*za(j2,j3)*za(j5,j6)*zb(j3,j4)*t(j2,j3,j4))+
     .(za(j2,j4)**2*zb(j1,j6)*(-(za(j1,j4)*zb(j1,j6))-za(j3,j4)*zb(j3,j6
     .)))/
     .(2d0*za(j2,j3)*za(j3,j4)*zb(j5,j6)*t(j2,j3,j4))+
     .((za(j1,j5)**2*zb(j1,j3)**2*
     .((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j1,j2)-
     .za(j1,j2)*za(j2,j4)*zb(j2,j4)-za(j1,j3)*za(j2,j4)*zb(j3,j4)
     .))/(za(j2,j3)*za(j5,j6)*zb(j3,j4))-
     .(za(j2,j4)*za(j4,j5)*(-(za(j1,j2)*za(j1,j4)*zb(j1,j6)*
     .zb(j2,j4))+za(j1,j2)*za(j2,j4)*zb(j2,j4)*zb(j2,j6)-
     .za(j1,j3)*za(j1,j4)*zb(j1,j6)*zb(j3,j4)+
     .za(j1,j3)*za(j2,j4)*zb(j2,j6)*zb(j3,j4)))/
     .(za(j2,j3)*za(j3,j4))+
     .(za(j2,j4)**3*zb(j2,j6)**2*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4)))/
     .(za(j2,j3)*za(j3,j4)*zb(j5,j6))-
     .(za(j1,j2)*za(j4,j5)*((-s(j1,j2)+2d0*s(j2,j4)+s(j3,j4)-s(j5,j6))*
     .zb(j3,j6)+2d0*za(j1,j5)*zb(j1,j3)*zb(j5,j6)))/za(j2,j3)+
     .(za(j1,j2)*(-(za(j1,j5)*zb(j1,j3))+za(j2,j5)*zb(j2,j3))*
     .zb(j3,j6)*(s(j2,j3)+t(j1,j3,j4)))/(za(j2,j3)*zb(j3,j4))+
     .(4d0*za(j1,j5)*za(j2,j4)*zb(j1,j6)*zb(j2,j3)*
     .(za(j1,j2)*zb(j2,j3)-za(j1,j4)*zb(j3,j4)))/t(j2,j3,j4)+
     .((za(j1,j2)*zb(j2,j3)-za(j1,j4)*zb(j3,j4))*
     .(-(za(j2,j5)*zb(j2,j3))+za(j4,j5)*zb(j3,j4))**2*
     .(2d0*s(j3,j4)*s(j5,j6)+(s(j1,j2)-s(j3,j4)-s(j5,j6))*t(j2,j3,j4)))/
     .(za(j5,j6)*zb(j3,j4)*t(j2,j3,j4)**2))/
     .(2d0*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))))

      return
      end