double complex function bub6(j1,j2,j3,j4,j5,j6,za,zb)
!----- C. Williams June 2011
!----- This function returns the coefficient of the log(s_gg) with j1=p and j2=m 
      implicit none 
      include 'constants.f'
      include 'sprods_com.f'
      include 'zprods_decl.f'
      integer j1,j2,j3,j4,j5,j6
      integer del1,del2
      double complex sym,antisym
      parameter(del1=7,del2=8)
      

      bub6= 
     &+sym(j1,j2,j3,j4,j5,j6,del1,del2,za,zb)
     &+sym(j1,j2,j3,j4,j5,j6,del2,del1,za,zb)
     &+antisym(j1,j2,j3,j4,j5,j6,del1,del2,za,zb)
     &-antisym(j1,j2,j3,j4,j5,j6,del2,del1,za,zb)
      
!---- W propagators 
    
      bub6=bub6/(s(j3,j4)*s(j5,j6))
      return 
      end


      double complex function sym(j1,j2,j3,j4,j5,j6,j7,j8,za,zb) 
!---- symmetric piece of coefficient 
      implicit none 
      include 'constants.f' 
      include 'zprods_decl.f' 
      include 'masses.f'
    
      integer j1,j2,j3,j4,j5,j6,j7,j8 
      double complex mtsq 
      double complex X 

      mtsq=mt**2

      X=-mtsq+za(j5,j6)*zb(j6,j5) 

      sym=  -((X*za(j2,j3)*za(j2,j5)*za(j2,j7)*za(j7,j8)**2*zb(j2,j1)*
     &     zb(j7,j4)*zb(j7,j6)*
     &      (-(za(j1,j2)*zb(j7,j2)) + (za(j1,j2)*zb(j2,j1)*(-(za(j1,j5)*
     &     zb(j7,j5)) - za(j1,j6)*zb(j7,j6)))/X)*zb(j8,j1)*zb(j8,j7))/
     &     (za(j1,j2)*zb(j7,j2)*(za(j1,j7)*zb(j7,j1) + za(j2,j7)*
     &     zb(j7,j2))*(za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5) 
     &     - za(j7,j8)**2*zb(j8,j7)**2)**2)) + 
     &(X*za(j2,j3)*za(j2,j5)*za(j2,j7)*za(j7,j8)**2*zb(j2,j1)*zb(j7,j6)*
     &     (za(j1,j2)*zb(j7,j1) + (za(j1,j2)*zb(j2,j1)*(-(za(j2,j5)*
     &    zb(j7,j5)) - za(j2,j6)*zb(j7,j6)))/X)*zb(j8,j4)*zb(j8,j7))/
     &   (za(j1,j2)*(za(j1,j7)*zb(j7,j1) + za(j2,j7)*zb(j7,j2))*
     & (za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5) 
     &    - za(j7,j8)**2*zb(j8,j7)**2)**2) + 
     &(X*za(j1,j3)*za(j2,j5)*za(j2,j7)*za(j7,j8)**2*zb(j2,j1)*zb(j7,j1)*
     &     (za(j1,j2)*zb(j7,j1) + (za(j1,j2)*zb(j2,j1)*(-(za(j2,j5)*
     &   zb(j7,j5)) - za(j2,j6)*zb(j7,j6)))/X)*(zb(j7,j6)*zb(j8,j4) 
     &   + zb(j7,j4)*zb(j8,j6))*
     &     zb(j8,j7))/(za(j1,j2)*zb(j7,j2)*(za(j1,j7)*zb(j7,j1) 
     &   + za(j2,j7)*zb(j7,j2))*
     &     (za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5)
     &     - za(j7,j8)**2*zb(j8,j7)**2)**2) + 
     &     (X*za(j1,j5)*za(j2,j7)*za(j7,j8)**2*zb(j2,j1)*zb(j7,j1)*
     &     (za(j1,j3)*zb(j7,j1) + za(j2,j3)*zb(j7,j2))*
     &     (za(j1,j2)*zb(j7,j1) + (za(j1,j2)*zb(j2,j1)*
     & (-(za(j2,j5)*zb(j7,j5)) - za(j2,j6)*zb(j7,j6)))/X)*
     &(zb(j7,j6)*zb(j8,j4) + zb(j7,j4)*zb(j8,j6))*
     &     zb(j8,j7))/(za(j1,j2)*zb(j7,j2)**2*(za(j1,j7)*zb(j7,j1) 
     &     + za(j2,j7)*zb(j7,j2))*
     &     (za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5) 
     &- za(j7,j8)**2*zb(j8,j7)**2)**2) - 
     &(X*za(j2,j5)*za(j2,j7)*za(j7,j8)**2*zb(j2,j1)*(za(j1,j3)*zb(j2,j1)
     &+ (za(j1,j2)*zb(j2,j1)*(za(j3,j5)*zb(j5,j2)
     & + za(j3,j6)*zb(j6,j2)))/X)*zb(j7,j1)*
     &     zb(j7,j4)*zb(j7,j6)*zb(j8,j7)**2)/
     &   (zb(j7,j2)**2*(za(j1,j7)*zb(j7,j1) + za(j2,j7)*zb(j7,j2))*
     &(za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5) - za(j7,j8)**2*
     & zb(j8,j7)**2)**2)  
     
      sym=sym
     & -(X*za(j1,j5)*za(j2,j7)*za(j7,j8)**2*zb(j2,j1)*(za(j1,j2)*
     & zb(j2,j1) + (za(j1,j2)*zb(j2,j1)*(za(j2,j5)*zb(j5,j2)
     & + za(j2,j6)*zb(j6,j2)))/X)*zb(j7,j1)*
     &     (za(j1,j3)*zb(j7,j1) + za(j2,j3)*zb(j7,j2))*
     &  zb(j7,j4)*zb(j7,j6)*zb(j8,j7)**2)/
     &   (za(j1,j2)*zb(j7,j2)**3*(za(j1,j7)*zb(j7,j1) + za(j2,j7)*
     & zb(j7,j2))*(za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5) 
     & - za(j7,j8)**2*zb(j8,j7)**2)**2) + 
     &  (X*za(j2,j3)*za(j2,j5)*za(j2,j7)*za(j7,j8)**2*
     & zb(j2,j1)*zb(j6,j2)*zb(j7,j4)*
     &     (za(j1,j2)*zb(j7,j1) + (za(j1,j2)*zb(j2,j1)*
     &(-(za(j2,j5)*zb(j7,j5)) - za(j2,j6)*zb(j7,j6)))/X)*zb(j8,j7)**2)/
     &   (za(j1,j2)*zb(j7,j2)*(za(j1,j7)*zb(j7,j1) + za(j2,j7)*
     &     zb(j7,j2))*(za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5) 
     & - za(j7,j8)**2*zb(j8,j7)**2)**2) - 
     &  (X*za(j1,j5)*za(j2,j7)*za(j7,j8)**2*zb(j2,j1)**2*(za(j1,j3)*
     & zb(j7,j1) + za(j2,j3)*zb(j7,j2))*zb(j7,j4)*zb(j7,j6)*
     &     (za(j1,j2)*zb(j7,j1) + (za(j1,j2)*zb(j2,j1)*(-(za(j2,j5)*
     & zb(j7,j5)) - za(j2,j6)*zb(j7,j6)))/X)*zb(j8,j7)**2)/
     &   (za(j1,j2)*zb(j7,j2)**3*(za(j1,j7)*zb(j7,j1) + za(j2,j7)*
     & zb(j7,j2))*(za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5)
     & - za(j7,j8)**2*zb(j8,j7)**2)**2) - 
     &  (X*za(j2,j7)*za(j3,j7)*za(j7,j8)**2*zb(j2,j1)**2*(za(j1,j5)*
     & zb(j7,j1) + za(j2,j5)*zb(j7,j2))*zb(j7,j4)*zb(j7,j6)*
     &     (za(j1,j2)*zb(j7,j1) + (za(j1,j2)*zb(j2,j1)*(-(za(j2,j5)*
     & zb(j7,j5)) - za(j2,j6)*zb(j7,j6)))/X)*zb(j8,j7)**2)/
     &   (zb(j7,j2)**2*(za(j1,j7)*zb(j7,j1) + za(j2,j7)*zb(j7,j2))**2*
     & (-(za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5)) + za(j7,j8)**2*
     &zb(j8,j7)**2)**2)

      return 
      end
      

      double complex function antisym(j1,j2,j3,j4,j5,j6,j7,j8,za,zb) 
!---- antisymmetric piece of coefficient 
      implicit none 
      include 'constants.f' 
      include 'zprods_decl.f' 
    
      integer j1,j2,j3,j4,j5,j6,j7,j8 

      antisym=
     &-((za(j2,j7)*za(j3,j5)*za(j7,j8)*zb(j2,j1)*zb(j6,j5)*
     & zb(j7,j1)*(za(j1,j5)*zb(j7,j1) + za(j2,j5)*zb(j7,j2))*
     & zb(j7,j4)*zb(j8,j7))/
     &     (zb(j7,j2)**2*(za(j1,j7)*zb(j7,j1) + za(j2,j7)*zb(j7,j2))*
     &     (za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5) - 
     & za(j7,j8)**2*zb(j8,j7)**2))) - 
     &  (za(j2,j7)*za(j7,j8)*zb(j2,j1)*zb(j4,j1)*(za(j1,j3)*zb(j7,j1) + 
     & za(j2,j3)*zb(j7,j2))*(za(j1,j5)*zb(j7,j1) + za(j2,j5)*zb(j7,j2))*
     & zb(j7,j6)*zb(j8,j7))/
     &   (zb(j7,j2)**2*(za(j1,j7)*zb(j7,j1) + za(j2,j7)*zb(j7,j2))*
     & (za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5) - 
     & za(j7,j8)**2*zb(j8,j7)**2)) + 
     &  (za(j2,j7)**2*za(j7,j8)*zb(j2,j1)**2*(za(j1,j3)*zb(j7,j1) 
     &+ za(j2,j3)*zb(j7,j2))*(za(j1,j5)*zb(j7,j1) + za(j2,j5)*zb(j7,j2))
     & *zb(j7,j4)*zb(j7,j6)*
     &     zb(j8,j7))/(2d0*zb(j7,j2)**2*(za(j1,j7)*zb(j7,j1)
     &+ za(j2,j7)*zb(j7,j2))**2*(za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5) 
     &     - za(j7,j8)**2*zb(j8,j7)**2))
     &   + (za(j1,j3)*za(j1,j5)*za(j2,j7)*za(j7,j8)*zb(j2,j1)**2*
     & zb(j7,j1)**2*zb(j7,j4)*zb(j7,j6)*zb(j8,j7))/
     &   (zb(j7,j2)**3*(za(j1,j7)*zb(j7,j1) + za(j2,j7)*zb(j7,j2))*
     &(za(j1,j2)*za(j5,j6)*zb(j2,j1)*zb(j6,j5) 
     &- za(j7,j8)**2*zb(j8,j7)**2))

      return 
      end