double complex function ttbqqbrqpp(k1,k2,k3,k4,k5,k6,k7)
      implicit none
      include 'constants.f'
      include 'sprods_com.f'
      include 'zprods_com.f'
      include 'masses.f'
      integer k1,k2,k3,k4,k5,k6,k7
      double precision s3459,s6789,mtsq
      s3459=s(k4,k3)+s(k5,k3)
      s6789=s(k3,k6)+s(k3,k7)
      mtsq=mt**2
      ttbqqbrqpp =  + xn**(-1)*s6789**(-1) * ( 1/(za(k1,k2))/(za(k5,k3)
     &    )/(zb(k1,k2))*za(k1,k5)*za(k5,k6)*za(k6,k7)*zb(k2,k6)*zb(k4,
     &    k5)*zb(k6,k3) + 1/(za(k1,k2))/(za(k5,k3))/(zb(k1,k2))*za(k1,
     &    k5)*za(k5,k7)*za(k6,k7)*zb(k2,k7)*zb(k4,k5)*zb(k6,k3) + 1/(
     &    za(k1,k2))/(zb(k1,k2))*za(k1,k5)*za(k6,k7)*zb(k2,k3)*zb(k4,k5
     &    )*zb(k6,k3) )
      ttbqqbrqpp = ttbqqbrqpp + xn**(-1)*mtsq*s6789**(-1) * ( 1/(za(k1,
     &    k2))/(za(k5,k3))/(zb(k1,k2))*za(k1,k5)*za(k5,k7)*zb(k2,k3)*
     &    zb(k4,k5) + 1/(za(k1,k2))/(za(k5,k3))/(zb(k1,k2))*za(k1,k5)*
     &    za(k6,k7)*zb(k2,k4)*zb(k6,k3) - 1/(za(k1,k2))/(za(k5,k3))/(
     &    zb(k1,k2))*za(k1,k6)*za(k5,k7)*zb(k2,k4)*zb(k6,k3) - 1/(za(k1
     &    ,k2))/(za(k5,k3))/(zb(k1,k2))*za(k1,k7)*za(k5,k7)*zb(k2,k4)*
     &    zb(k7,k3) )
      ttbqqbrqpp = ttbqqbrqpp + xn**(-1)*mtsq*s3459**(-1) * (  - 1/(za(
     &    k1,k2))/(za(k5,k3))/(zb(k1,k2))*za(k1,k5)*za(k6,k7)*zb(k2,k6)
     &    *zb(k4,k3) - 1/(za(k1,k2))/(za(k5,k3))/(zb(k1,k2))*za(k1,k7)*
     &    za(k4,k5)*zb(k2,k4)*zb(k4,k3) + 1/(za(k1,k2))/(zb(k1,k2))*za(
     &    k1,k7)*zb(k2,k3)*zb(k4,k3) )

      return
      end