double complex function Fsc5(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 L0,L1,Lsm1_2me,Lnrat,Ls1
      double precision t  
      Fsc5= 
     .(Lsm1_2me(t(j1,j2,j3),t(j2,j3,j4),s(j2,j3),s(j5,j6))*za(j1,j3)*za(
     .j1,j5)**2*
     .za(j3,j4))/(za(j1,j2)*za(j1,j4)**3*za(j2,j3)*za(j5,j6))-
     .(za(j1,j3)*za(j3,j5)**2)/
     .(2d0*za(j1,j2)*za(j1,j4)*za(j2,j3)*za(j3,j4)*za(j5,j6))+
     .(Lnrat(-t(j2,j3,j4),-s(j5,j6))*za(j1,j3)*za(j3,j5)**2)/
     .(2d0*za(j1,j2)*za(j1,j4)*za(j2,j3)*za(j3,j4)*za(j5,j6))-
     .(L0(-t(j1,j2,j3),-s(j2,j3))*za(j1,j3)*za(j1,j5)**2*zb(j1,j2))/
     .(s(j2,j3)*za(j1,j2)*za(j1,j4)**2*za(j5,j6))+
     .(L0(-t(j2,j3,j4),-s(j5,j6))*za(j1,j3)*za(j1,j5)**2*
     .(za(j2,j3)*zb(j1,j2)-za(j3,j4)*zb(j1,j4)))/
     .(s(j5,j6)*za(j1,j2)*za(j1,j4)**2*za(j2,j3)*za(j5,j6))+
     .(L0(-t(j2,j3,j4),-s(j5,j6))*za(j1,j3)**2*za(j3,j5)*zb(j1,j6))/
     .(s(j5,j6)*za(j1,j2)*za(j1,j4)*za(j2,j3)*za(j3,j4))+
     .(L1(-t(j2,j3,j4),-s(j5,j6))*za(j1,j3)**3*za(j5,j6)*zb(j1,j6)**2)/
     .(2d0*s(j5,j6)**2*za(j1,j2)*za(j1,j4)*za(j2,j3)*za(j3,j4))+
     .(L0(-t(j2,j3,j4),-s(j2,j3))*za(j1,j5)**2*za(j3,j4)*zb(j2,j4))/
     .(s(j2,j3)*za(j1,j2)*za(j1,j4)**2*za(j5,j6))-
     .(L0(-t(j1,j2,j3),-s(j5,j6))*za(j1,j3)*za(j1,j5)*za(j4,j5)*
     .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4)))/
     .(s(j5,j6)*za(j1,j2)*za(j1,j4)**2*za(j2,j3)*za(j5,j6))-
     .(L0(-t(j1,j2,j3),-s(j5,j6))*za(j1,j3)*za(j3,j5)*zb(j4,j6))/
     .(s(j5,j6)*za(j1,j2)*za(j1,j4)*za(j2,j3))+
     .(L1(-t(j1,j2,j3),-s(j5,j6))*za(j1,j3)*za(j3,j4)*za(j5,j6)*zb(j4,j6
     .)**2)/
     .(s(j5,j6)**2*za(j1,j2)*za(j1,j4)*za(j2,j3))+
     .(L1(-s(j3,j4),-t(j2,j3,j4))*za(j2,j5)**2*za(j3,j4)*zb(j2,j4)**2)/
     .(2d0*za(j1,j2)*za(j2,j4)*za(j5,j6)*t(j2,j3,j4)**2)+
     .(Ls1(-s(j2,j3),-t(j2,j3,j4),-s(j3,j4),-t(j2,j3,j4))*za(j2,j5)**2*z
     .a(j3,j4)*
     .zb(j2,j4)**2)/(za(j1,j2)*za(j2,j4)*za(j5,j6)*t(j2,j3,j4)**2)-
     .(L1(-s(j2,j3),-t(j2,j3,j4))*za(j2,j3)*za(j4,j5)**2*zb(j2,j4)**2)/
     .(2d0*za(j1,j4)*za(j2,j4)*za(j5,j6)*t(j2,j3,j4)**2)-
     .(za(j1,j5)*za(j2,j4)*za(j3,j5)*zb(j2,j4)**2)/
     .(2d0*za(j1,j2)*za(j1,j4)*za(j5,j6)*t(j1,j2,j4)*t(j2,j3,j4))-
     .(zb(j2,j4)*(za(j1,j2)*zb(j2,j6)+za(j1,j4)*zb(j4,j6))*
     .(-(za(j2,j3)*zb(j2,j6))+za(j3,j4)*zb(j4,j6)))/
     .(2d0*za(j1,j2)*za(j1,j4)*zb(j5,j6)*t(j1,j2,j4)*t(j2,j3,j4))
      return
      end