double complex function Fsc2(j1,j2,j3,j4,j5,j6,za,zb) 
      implicit none
      include 'constants.f'
      include 'zprods_decl.f'
      include 'sprods_com.f'
      integer j1,j2,j3,j4,j5,j6
      double complex L1
      double precision t  
      Fsc2= 
     .za(j3,j5)**2/(2d0*za(j1,j2)*za(j2,j3)*za(j5,j6)*zb(j3,j4))+
     .zb(j2,j6)**2/(2d0*za(j1,j2)*zb(j2,j3)*zb(j3,j4)*zb(j5,j6))-
     .(za(j1,j3)*((L1(-t(j1,j2,j3),-s(j2,j3))*za(j2,j3)**2*zb(j1,j2)**2*
     .(za(j1,j2)*zb(j2,j6)+za(j1,j3)*zb(j3,j6))**2)/s(j2,j3)**2+
     .L1(-s(j5,j6),-t(j1,j2,j3))*za(j3,j4)**2*zb(j4,j6)**2))/
     .(2d0*za(j1,j2)*za(j2,j3)*(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))
     .*
     .zb(j5,j6)*t(j1,j2,j3))-(zb(j2,j4)*
     .(L1(-s(j5,j6),-t(j2,j3,j4))*za(j1,j5)**2*zb(j1,j2)**2+
     .(L1(-t(j2,j3,j4),-s(j2,j3))*za(j3,j4)**2*zb(j2,j3)**2*
     .(-(za(j2,j5)*zb(j2,j4))-za(j3,j5)*zb(j3,j4))**2)/s(j2,j3)**2))/
     .(2d0*za(j5,j6)*zb(j2,j3)*zb(j3,j4)*
     .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4))*t(j2,j3,j4))
      return
      end