subroutine dksrwzc(p1,p2,p3,p4,p5,p6,p7,cc11,cc12) implicit none include 'constants.f' include 'sprods_com.f' include 'zprods_com.f' integer p1,p2,p3,p4,p5,p6,p7 double precision s12,s34,s126,s345 double complex iza,izb,cc11,cc12 iza(p1,p2)=cone/za(p1,p2) izb(p1,p2)=cone/zb(p1,p2) s126=s(p1,p2)+s(p1,p6)+s(p2,p6) s345=s(p3,p4)+s(p3,p5)+s(p4,p5) s34=s(p3,p4) s12=s(p1,p2) cc11= + s12**(-1)*s34**(-1) * ( - za(p1,p2)*za(p3,p5)*za(p3,p7)* & zb(p1,p6)**2*zb(p3,p4)*izb(p6,p7)*s126**(-1)*s345**(-1) + za( & p1,p2)*za(p3,p5)*za(p5,p7)*zb(p1,p6)**2*zb(p4,p5)*izb(p6,p7)* & s126**(-1)*s345**(-1) + za(p1,p2)*za(p3,p5)*zb(p1,p4)*zb(p1, & p6)*zb(p5,p6)*izb(p5,p7)*izb(p6,p7)*s126**(-1) - za(p1,p2)* & za(p3,p7)*zb(p1,p4)*zb(p1,p6)*izb(p5,p7)*s126**(-1) + za(p2, & p6)*za(p3,p5)*zb(p1,p6)*zb(p4,p6)*zb(p5,p6)*izb(p5,p7)*izb(p6 & ,p7)*s126**(-1) - za(p2,p6)*za(p3,p7)*zb(p1,p6)*zb(p4,p6)* & izb(p5,p7)*s126**(-1) ) cc12= + s12**(-1)*s34**(-1) * ( za(p1,p2)*za(p3,p5)**2*zb(p1,p6)* & zb(p1,p7)*zb(p3,p4)*iza(p5,p7)*s126**(-1)*s345**(-1) - za(p2, & p3)*za(p3,p5)*za(p5,p6)*zb(p1,p6)*zb(p3,p4)*iza(p5,p7)*iza(p6 & ,p7)*s345**(-1) - za(p2,p3)*za(p3,p5)*zb(p1,p7)*zb(p3,p4)* & iza(p6,p7)*s345**(-1) + za(p2,p5)*za(p3,p5)*za(p5,p6)*zb(p1, & p6)*zb(p4,p5)*iza(p5,p7)*iza(p6,p7)*s345**(-1) + za(p2,p5)* & za(p3,p5)*zb(p1,p7)*zb(p4,p5)*iza(p6,p7)*s345**(-1) - za(p2, & p6)*za(p3,p5)**2*zb(p1,p6)*zb(p3,p4)*zb(p6,p7)*iza(p5,p7)* & s126**(-1)*s345**(-1) ) return end