double complex function Faxsl(st,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' include 'masses.f' character*9 st double complex L1 double precision t,mtsq mtsq=mt**2 if(st.eq.'q+qb-g+g-') then Faxsl= .(2d0*(1/(24d0*mtsq)-L1(-t(j1,j2,j4),-s(j5,j6))/(2d0*s(j5,j6)))*zb( .j1,j3)* .(za(j2,j5)*zb(j1,j2)+za(j4,j5)*zb(j1,j4))*zb(j3,j6))/ .(s(j5,j6)*zb(j1,j4)*zb(j2,j4))+ .(2d0*(1/(24d0*mtsq)-L1(-t(j1,j2,j3),-s(j5,j6))/(2d0*s(j5,j6)))*za( .j2,j4)*za(j4,j5)* .(-(za(j1,j2)*zb(j1,j6))+za(j2,j3)*zb(j3,j6)))/ .(s(j5,j6)*za(j1,j3)*za(j2,j3)) elseif(st.eq.'q+qb-g+g+') then Faxsl= .(2d0*(1/(24d0*mtsq)-L1(-t(j1,j2,j4),-s(j5,j6))/(2d0*s(j5,j6)))*za( .j2,j5)* .(-(za(j1,j2)*zb(j1,j3))-za(j2,j4)*zb(j3,j4))*zb(j3,j6))/ .(s(j5,j6)*za(j1,j4)*za(j2,j4))+ .(2d0*(1/(24d0*mtsq)-L1(-t(j1,j2,j3),-s(j5,j6))/(2d0*s(j5,j6)))*za( .j2,j5)* .(-(za(j1,j2)*zb(j1,j4))+za(j2,j3)*zb(j3,j4))*zb(j4,j6))/ .(s(j5,j6)*za(j1,j3)*za(j2,j3)) endif return end