double complex function Fsc6(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_2mh,I3m,Lnrat double precision t Fsc6= .-(za(j3,j5)*(-(za(j2,j5)*za(j3,j4))+za(j2,j3)*za(j4,j5)))/ .(2d0*za(j1,j4)*za(j3,j4)*za(j5,j6)* .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3)))+ .(Lsm1_2mh(s(j3,j4),t(j1,j2,j3),s(j1,j2),s(j5,j6))* .(za(j2,j4)*zb(j1,j2)+za(j3,j4)*zb(j1,j3))* .(-(za(j1,j5)*zb(j1,j3))-za(j2,j5)*zb(j2,j3))**2)/ .(za(j5,j6)*zb(j2,j3)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))* .*3)- .(L0(-t(j1,j2,j3),-s(j1,j2))*za(j2,j3)*zb(j1,j2)* .(-(za(j1,j5)*zb(j1,j3))-za(j2,j5)*zb(j2,j3))**2)/ .(s(j1,j2)*za(j5,j6)*zb(j2,j3)* .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))**2)+ .(Lnrat(-s(j1,j2),-s(j5,j6))*za(j4,j5)* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j2,j5)*zb(j2,j4))-za(j3,j5)*zb(j3,j4)))/ .(2d0*s(j3,j4)*za(j5,j6)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3 .))* .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4)))- .(za(j3,j5)*((-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j2,j4)+ .2d0*za(j1,j2)*za(j3,j4)*zb(j1,j3))*zb(j4,j6))/ .((s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5, .j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*za(j1,j4)* .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3)))- .(za(j3,j5)*((-s(j1,j2)-s(j3,j4)+s(j5,j6))*zb(j1,j3)+ .2d0*za(j2,j4)*zb(j1,j2)*zb(j3,j4))*zb(j4,j6))/ .((s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5, .j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*zb(j2,j3)* .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))) Fsc6=Fsc6+ .(za(j1,j3)*(-((za(j3,j5)*(-((za(j2,j5)* .(-((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j1,j3)* .zb(j1,j2))+ .(-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j3,j4)*zb(j2,j4) .))/za(j5,j6))+zb(j1,j6)*((-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j1,j3)+ .2d0*za(j1,j2)*za(j3,j4)*zb(j2,j4))))/za(j3,j4))+ .(zb(j4,j6)*(-(za(j2,j5)* .((-s(j1,j2)-s(j3,j4)+s(j5,j6))*zb(j2,j4)+ .2d0*za(j1,j3)*zb(j1,j2)*zb(j3,j4)))- .(zb(j1,j6)*((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j1,j2)* .zb(j2,j4)- .(-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j1,j3)*zb(j3,j4)))/ .zb(j5,j6)))/zb(j3,j4)))/ .(2d0*(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s .(j5,j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*za(j1,j4)* .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4)))+ .(zb(j2,j4)*(-((za(j3,j5)*((za(j2,j5)* .(-((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j1,j3)* .zb(j1,j2))+ .(-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j3,j4)*zb(j2,j4)))/ .za(j5,j6)-zb(j1,j6)* .((-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j1,j3)+ .2d0*za(j1,j2)*za(j3,j4)*zb(j2,j4))))/za(j3,j4))+ .(zb(j4,j6)*(za(j2,j5)* .((-s(j1,j2)-s(j3,j4)+s(j5,j6))*zb(j2,j4)+ .2d0*za(j1,j3)*zb(j1,j2)*zb(j3,j4))+ .(zb(j1,j6)*((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j1,j2)* .zb(j2,j4)- .(-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j1,j3)*zb(j3,j4)))/ .zb(j5,j6)))/zb(j3,j4)))/ .(2d0*(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s .(j5,j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*zb(j2,j3)* .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))) Fsc6=Fsc6- .(((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j2,j4)*zb(j1,j2)- .(-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j3,j4)*zb(j1,j3))* .((za(j3,j5)**2*zb(j3,j4))/za(j5,j6)+ .(za(j3,j4)*zb(j4,j6)**2)/zb(j5,j6)))/ .(2d0*(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s .(j5,j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*za(j3,j4)*zb(j2,j3)* .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3)))+ .((-((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j1,j2)*zb(j1,j3))+ .(-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j2,j4)*zb(j3,j4))* .((za(j3,j5)**2*zb(j3,j4))/za(j5,j6)+ .(za(j3,j4)*zb(j4,j6)**2)/zb(j5,j6)))/ .(2d0*(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s .(j5,j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*za(j1,j4)* .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))*zb(j3,j4))- .(L0(-t(j1,j2,j4),-s(j1,j2))*za(j1,j2)*zb(j1,j4)* .(-(za(j1,j4)*zb(j1,j6))-za(j2,j4)*zb(j2,j6))**2)/ .(s(j1,j2)*za(j1,j4)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))** .2* .zb(j5,j6))-(Lsm1_2mh(s(j3,j4),t(j1,j2,j4),s(j1,j2),s(j5,j6))* .(-(za(j1,j4)*zb(j1,j6))-za(j2,j4)*zb(j2,j6))**2* .(-(za(j1,j2)*zb(j1,j3))-za(j2,j4)*zb(j3,j4)))/ .(za(j1,j4)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))**3*zb(j5,j .6))- .((zb(j1,j6)*zb(j3,j4)+zb(j1,j4)*zb(j3,j6))*zb(j4,j6))/ .(2d0*zb(j2,j3)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))*zb(j3, .j4)* .zb(j5,j6))-(Lnrat(-s(j1,j2),-s(j5,j6))* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))*zb(j3,j6)* .(-(za(j1,j3)*zb(j1,j6))+za(j3,j4)*zb(j4,j6)))/ .(2d0*s(j3,j4)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))* .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*zb(j5,j6)) Fsc6=Fsc6- .(Lnrat(-s(j1,j2),-s(j5,j6))*za(j1,j2)* .((-6d0*zb(j1,j2)*((-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j2,j4)+ .2d0*za(j1,j2)*za(j3,j4)*zb(j1,j3))* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6)))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)**2+ .((-(za(j1,j4)*zb(j1,j6))-za(j2,j4)*zb(j2,j6))* .(-(zb(j1,j6)*zb(j3,j4))-zb(j1,j4)*zb(j3,j6)))/ .((-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))*zb(j3,j4)* .zb(j5,j6))+(-((za(j3,j5)* .(2d0*((-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j4,j5)* .zb(j1,j4)+ .(-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j5,j6)*zb(j1,j6)) .-za(j3,j5)*(za(j2,j4)*zb(j1,j2)+za(j3,j4)*zb(j1,j3))*zb(j3,j4)))/ .za(j5,j6))-((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j2,j4)* .zb(j1,j2)*(za(j3,j5)*zb(j3,j6)-za(j4,j5)*zb(j4,j6)))/ .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))- .(za(j2,j3)*zb(j1,j2)* .(-(za(j1,j4)*zb(j1,j6))-za(j2,j4)*zb(j2,j6))*zb(j4,j6))/ .zb(j5,j6)+(zb(j1,j2)*zb(j3,j6)* .(((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j2,j4)* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .zb(j3,j6))/ .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))- .(2d0*(s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j2,j3)- .za(j1,j2)*za(j3,j4)*zb(j1,j4))*zb(j4,j6)))/ .(zb(j3,j4)*zb(j5,j6))- .((s(j1,j2)-s(j3,j4)-s(j5,j6))*zb(j1,j3)*zb(j4,j6)* .(-(za(j1,j3)*zb(j1,j6))-za(j2,j3)*zb(j2,j6)- .za(j3,j5)*zb(j5,j6)))/(zb(j3,j4)*zb(j5,j6))+ .za(j3,j5)*(-(za(j2,j4)*zb(j1,j2)*zb(j4,j6))+ .za(j3,j5)*zb(j1,j3)*zb(j5,j6))+ .4d0*za(j4,j5)*zb(j1,j2)* .(za(j2,j3)*zb(j4,j6)+ .(za(j2,j4)*za(j3,j5)*zb(j3,j4)*zb(j5,j6))/ .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)))/ .(2d0*za(j1,j4)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))) Fsc6=Fsc6+ .(Lnrat(-s(j1,j2),-s(j5,j6))*zb(j1,j2)* .(((za(j2,j5)*za(j3,j4)-za(j2,j3)*za(j4,j5))* .(-(za(j1,j5)*zb(j1,j3))-za(j2,j5)*zb(j2,j3)))/ .(za(j3,j4)*za(j5,j6)*(-(za(j1,j4)*zb(j1,j3))- .za(j2,j4)*zb(j2,j3)))+ .(6d0*za(j1,j2)*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))* .((-s(j1,j2)-s(j3,j4)+s(j5,j6))*zb(j1,j3)+ .2d0*za(j2,j4)*zb(j1,j2)*zb(j3,j4)))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)**2+ .(-((za(j1,j2)*za(j3,j5)*zb(j1,j4)* .(-(za(j1,j5)*zb(j1,j3))-za(j2,j5)*zb(j2,j3)))/za(j5,j6)) .-(za(j1,j2)*za(j4,j5)*(((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j4,j5)* .zb(j1,j3)* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4)))/ .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))- .za(j3,j5)*(2d0*(s(j1,j2)-s(j3,j4)-s(j5,j6))*zb(j1,j4)- .za(j2,j3)*zb(j1,j2)*zb(j3,j4))))/(za(j3,j4)*za(j5,j6)) .+((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j1,j2)*zb(j1,j3)* .(-(za(j3,j5)*zb(j3,j6))+za(j4,j5)*zb(j4,j6)))/ .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))- .((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j2,j4)*za(j3,j5)* .(-(za(j1,j5)*zb(j1,j4))-za(j2,j5)*zb(j2,j4)+ .za(j5,j6)*zb(j4,j6)))/(za(j3,j4)*za(j5,j6))+ .zb(j4,j6)*(za(j1,j2)*za(j3,j5)*zb(j1,j3)- .za(j2,j4)*za(j5,j6)*zb(j4,j6))- .4d0*za(j1,j2)*zb(j3,j6)* .(za(j3,j5)*zb(j1,j4)+ .(za(j3,j4)*za(j5,j6)*zb(j1,j3)*zb(j4,j6))/ .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3)))+ .(zb(j4,j6)*(za(j3,j4)* .(-(za(j1,j2)*zb(j1,j3))-za(j2,j4)*zb(j3,j4))*zb(j4,j6) .+2d0*((-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j2,j3)*zb(j3,j6)- .(-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j2,j5)*zb(j5,j6))))/ .zb(j5,j6))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)))/ .(2d0*zb(j2,j3)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3)))+ .(L1(-s(j5,j6),-t(j1,j2,j3))*za(j4,j5)*zb(j1,j4)*zb(j4,j6))/ .(zb(j2,j3)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))*t(j1,j2,j3 .)) Fsc6=Fsc6- .(L0(-t(j1,j2,j3),-s(j5,j6))*za(j4,j5)*zb(j1,j4)* .(-(za(j1,j5)*zb(j1,j3))-za(j2,j5)*zb(j2,j3))*t(j1,j2,j3))/ .(s(j5,j6)*za(j5,j6)*zb(j2,j3)* .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))**2)- .(Lnrat(-s(j3,j4),-s(j5,j6))*(-((za(j3,j5)* .(2d0*((-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j4,j5)*zb(j1,j4)+ .(-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j5,j6)*zb(j1,j6))+ .(za(j2,j4)*zb(j1,j2)+za(j3,j4)*zb(j1,j3))* .(-(za(j1,j5)*zb(j1,j4))-za(j2,j5)*zb(j2,j4)))*zb(j3,j4))/ .za(j5,j6))+(6d0*zb(j1,j2)* .((-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j2,j4)+ .2d0*za(j1,j2)*za(j3,j4)*zb(j1,j3))* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))*zb(j3,j4))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)- .(zb(j4,j6)*(-(za(j2,j4)*zb(j1,j2)* .(-(za(j1,j3)*zb(j1,j6))-za(j2,j3)*zb(j2,j6))*zb(j3,j4)) .-za(j3,j4)*zb(j1,j3)*(-(za(j1,j3)*zb(j1,j6))-za(j2,j3)*zb(j2,j6))* .zb(j3,j4)+3d0*(-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j3,j4)* .zb(j1,j3)*zb(j4,j6)+ .2d0*za(j2,j3)*zb(j1,j2)* .((s(j1,j2)-s(j3,j4)-s(j5,j6))*zb(j3,j6)- .2d0*za(j4,j5)*zb(j3,j4)*zb(j5,j6))))/zb(j5,j6)+ .(zb(j1,j3)*(-((za(j3,j5)*za(j4,j5)*zb(j3,j4))/za(j5,j6))- .(za(j3,j4)*zb(j3,j6)*zb(j4,j6))/zb(j5,j6))* .(-s(j1,j2)**2+2d0*s(j1,j2)*s(j3,j4)-s(j3,j4)**2+ .2d0*s(j1,j2)*s(j5,j6)+2d0*s(j3,j4)*s(j5,j6)-s(j5,j6)**2+ .(-s(j1,j2)+s(j3,j4)-s(j5,j6))*(t(j1,j2,j3)-t(j1,j2,j4))))/ .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))))/ .(2d0*(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s .(j5,j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*zb(j2,j3)* .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))) Fsc6=Fsc6- .(I3m(s(j1,j2),s(j3,j4),s(j5,j6))*za(j1,j2)* .((3d0*za(j3,j4)*zb(j1,j2)* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))* .(-((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j1,j2)*zb(j1,j3))+ .(-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j2,j4)*zb(j3,j4)))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)**2- .((-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(za(j3,j4)*zb(j1,j3)* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))- .3d0*za(j2,j4)*za(j3,j5)*zb(j1,j2)*zb(j3,j6))- .za(j2,j3)*zb(j1,j2)* .((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j4,j5)- .za(j3,j5)*(-(za(j1,j4)*zb(j1,j3))- .za(j2,j4)*zb(j2,j3))-2d0*za(j3,j4)*za(j5,j6)*zb(j3,j6))* .zb(j4,j6)+za(j3,j5)*zb(j1,j4)* .(-((-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j3,j4)*zb(j3,j6))- .(-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j4,j5)*zb(j5,j6))+ .(za(j2,j4)*za(j4,j5)*zb(j1,j2)* .(-(za(j1,j3)*zb(j1,j4)*zb(j3,j6))- .za(j2,j3)*zb(j2,j4)*zb(j3,j6))*(t(j1,j2,j3)-t(j1,j2,j4)))/ .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3)))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)+ .(za(j3,j5)*zb(j1,j4)*(-(za(j1,j4)*zb(j1,j6))- .za(j2,j4)*zb(j2,j6)))/t(j1,j2,j4)))/ .(za(j1,j4)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))) Fsc6=Fsc6+ .(L1(-s(j5,j6),-t(j1,j2,j4))*za(j2,j3)*za(j3,j5)*zb(j3,j6))/ .(za(j1,j4)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))*t(j1,j2,j4 .))+ .(L0(-t(j1,j2,j4),-s(j5,j6))*za(j2,j3)* .(-(za(j1,j4)*zb(j1,j6))-za(j2,j4)*zb(j2,j6))*zb(j3,j6)*t(j1,j2,j4) .)/ .(s(j5,j6)*za(j1,j4)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))** .2* .zb(j5,j6))-(Lnrat(-s(j3,j4),-s(j5,j6))* .((6d0*za(j1,j2)*za(j3,j4)* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))* .((-s(j1,j2)-s(j3,j4)+s(j5,j6))*zb(j1,j3)+ .2d0*za(j2,j4)*zb(j1,j2)*zb(j3,j4)))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)+ .(za(j3,j5)*(-(za(j1,j2)*za(j3,j4)*zb(j1,j3)* .(-(za(j1,j5)*zb(j1,j4))-za(j2,j5)*zb(j2,j4)))- .3d0*(-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j2,j4)*za(j3,j5)* .zb(j3,j4)-za(j2,j4)*za(j3,j4)* .(-(za(j1,j5)*zb(j1,j4))-za(j2,j5)*zb(j2,j4))*zb(j3,j4)- .2d0*za(j1,j2)*zb(j1,j4)* .((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j4,j5)- .2d0*za(j3,j4)*za(j5,j6)*zb(j3,j6))))/za(j5,j6)- .(za(j3,j4)*zb(j4,j6)*((-(za(j1,j3)*zb(j1,j6))- .za(j2,j3)*zb(j2,j6))* .(-(za(j1,j2)*zb(j1,j3))-za(j2,j4)*zb(j3,j4))+ .2d0*((-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j2,j3)*zb(j3,j6)- .(-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j2,j5)*zb(j5,j6))))/ .zb(j5,j6)+(za(j2,j4)* .(-((za(j3,j5)*za(j4,j5)*zb(j3,j4))/za(j5,j6))- .(za(j3,j4)*zb(j3,j6)*zb(j4,j6))/zb(j5,j6))* .(-s(j1,j2)**2+2d0*s(j1,j2)*s(j3,j4)-s(j3,j4)**2+ .2d0*s(j1,j2)*s(j5,j6)+2d0*s(j3,j4)*s(j5,j6)-s(j5,j6)**2+ .(-s(j1,j2)+s(j3,j4)-s(j5,j6))*(-t(j1,j2,j3)+t(j1,j2,j4))))/ .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))))/ .(2d0*(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s .(j5,j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*za(j1,j4)* .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))) Fsc6=Fsc6+ .(I3m(s(j1,j2),s(j3,j4),s(j5,j6))*zb(j1,j2)* .((3d0*za(j1,j2)*((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j2,j4)*zb(j1,j2)- .(-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j3,j4)*zb(j1,j3))* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))*zb(j3,j4))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)**2+ .(za(j2,j3)*(-(za(j1,j5)*zb(j1,j3))-za(j2,j5)*zb(j2,j3))* .zb(j4,j6))/t(j1,j2,j3)- .(za(j2,j3)*((-s(j1,j2)+s(j3,j4)-s(j5,j6))*za(j4,j5)*zb(j3,j4)+ .(-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j5,j6)*zb(j3,j6))* .zb(j4,j6)+(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j2,j4)*(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))* .zb(j3,j4))+3d0*za(j1,j2)*za(j4,j5)*zb(j1,j3)*zb(j4,j6))+ .za(j1,j2)*za(j3,j5)*zb(j1,j4)* .((s(j1,j2)-s(j3,j4)-s(j5,j6))*zb(j3,j6)- .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3))*zb(j4,j6)- .2d0*za(j4,j5)*zb(j3,j4)*zb(j5,j6))- .(za(j1,j2)*zb(j1,j3)* .(-(za(j1,j3)*za(j4,j5)*zb(j1,j4))- .za(j2,j3)*za(j4,j5)*zb(j2,j4))*zb(j3,j6)* .(-t(j1,j2,j3)+t(j1,j2,j4)))/ .(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3)))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)))/ .(zb(j2,j3)*(-(za(j1,j4)*zb(j1,j3))-za(j2,j4)*zb(j2,j3)))+ .(za(j1,j3)*zb(j1,j4)*(za(j1,j3)*zb(j3,j6)+za(j1,j4)*zb(j4,j6))* .((L0(-t(j1,j3,j4),-s(j3,j4))* .((-2d0*za(j1,j3)*zb(j2,j6))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))+ .(-(za(j1,j3)*zb(j1,j6))+za(j3,j4)*zb(j4,j6))/t(j1,j3,j4)))/ .s(j3,j4)-(L1(-s(j3,j4),-t(j1,j3,j4))*za(j1,j3)*zb(j1,j6))/t(j1,j3, .j4)**2 .))/(2d0*za(j1,j4)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*zb( .j5,j6)) Fsc6=Fsc6+ .(za(j1,j3)*za(j2,j3)*zb(j2,j6)**2* .((-2d0*L0(-s(j5,j6),-t(j1,j3,j4))*za(j1,j3))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))+ .(L1(-s(j5,j6),-t(j1,j3,j4))*za(j2,j3))/t(j1,j3,j4)))/ .(2d0*za(j1,j4)*za(j3,j4)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j .4))* .zb(j5,j6))-(za(j1,j3)*zb(j4,j6)* .(-(zb(j4,j6)/zb(j3,j4))+ .(-(za(j1,j3)*zb(j1,j6))+za(j3,j4)*zb(j4,j6))/t(j1,j3,j4)))/ .(2d0*za(j1,j4)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*zb(j5, .j6))- .(Lsm1_2mh(s(j1,j2),t(j1,j3,j4),s(j3,j4),s(j5,j6))*za(j1,j3)**3*zb( .j2,j6)**2* .t(j1,j3,j4))/ .(za(j1,j4)*za(j3,j4)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))* .*3* .zb(j5,j6))+(I3m(s(j1,j2),s(j3,j4),s(j5,j6))*zb(j2,j4)* .((s(j1,j2)*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6)))/2- .za(j1,j2)*zb(j1,j4)*((-(za(j1,j3)*zb(j1,j2))- .za(j3,j4)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))- .za(j1,j5)*za(j3,j4)*zb(j1,j2)*zb(j4,j6))- .(za(j1,j3)*za(j1,j5)*(-(za(j2,j3)*zb(j1,j3))- .za(j2,j4)*zb(j1,j4))*zb(j2,j4)* .(-(zb(j2,j3)*(-(za(j1,j3)*zb(j1,j6))-za(j2,j3)*zb(j2,j6)))- .zb(j2,j4)*(-(za(j1,j4)*zb(j1,j6))-za(j2,j4)*zb(j2,j6))- .zb(j1,j2)*(za(j1,j3)*zb(j3,j6)+za(j1,j4)*zb(j4,j6))))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))+ .za(j1,j5)*za(j2,j3)*(zb(j2,j4)* .(-(zb(j1,j3)*(-(za(j1,j3)*zb(j1,j6))- .za(j2,j3)*zb(j2,j6)))- .zb(j1,j4)*(-(za(j1,j4)*zb(j1,j6))-za(j2,j4)*zb(j2,j6))+ .zb(j1,j2)*(za(j2,j3)*zb(j3,j6)+za(j2,j4)*zb(j4,j6)))- .zb(j1,j2)*(-(zb(j1,j4)* .(za(j1,j3)*zb(j3,j6)+za(j1,j4)*zb(j4,j6)))- .zb(j2,j4)*(za(j2,j3)*zb(j3,j6)+za(j2,j4)*zb(j4,j6))))+ .(3d0*s(j1,j2)*(s(j1,j2)-s(j3,j4)-s(j5,j6))* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))* .(t(j1,j3,j4)-t(j2,j3,j4)))/ .(2d0*(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2- .2d0*s(j1,j2)*s(j5,j6)-2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2))))/ .((s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5, .j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*zb(j2,j3)* .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))) Fsc6=Fsc6+ .(Lnrat(-s(j3,j4),-s(j5,j6))*zb(j2,j4)* .(-2d0*(-(za(j1,j3)*zb(j1,j4))+za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))+ .4d0*(-(za(j1,j2)*za(j3,j5)*zb(j1,j4)*zb(j2,j6))+ .za(j1,j3)*za(j2,j5)*zb(j1,j2)*zb(j4,j6))+ .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4)+ .(4d0*s(j5,j6)*za(j1,j3)*zb(j2,j4))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4)))* .(-((za(j1,j5)*za(j2,j5)*zb(j1,j2))/za(j5,j6))+ .(za(j1,j2)*zb(j1,j6)*zb(j2,j6))/zb(j5,j6))- .(-s(j1,j2)+s(j3,j4)-s(j5,j6))* .((za(j1,j5)*za(j3,j5)*zb(j1,j4))/za(j5,j6)+ .(2d0*za(j1,j5)*zb(j2,j4)* .(za(j1,j5)*za(j2,j3)*zb(j1,j2)- .za(j1,j5)*za(j3,j4)*zb(j1,j4)- .za(j2,j3)*za(j5,j6)*zb(j2,j6)))/ .(za(j5,j6)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4)))- .(za(j2,j3)*zb(j2,j6)*zb(j4,j6))/zb(j5,j6))+ .(8d0*s(j1,j2)*za(j1,j5)*za(j3,j5)*zb(j2,j4)*zb(j5,j6))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))+ .(3d0*(-s(j1,j2)-s(j3,j4)+s(j5,j6))* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))* .(t(j1,j3,j4)-t(j2,j3,j4)))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)))/ .(2d0*(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s .(j5,j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*zb(j2,j3)* .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4)))- .(za(j2,j3)*zb(j2,j4)*(za(j3,j5)*zb(j2,j3)+za(j4,j5)*zb(j2,j4))* .((L0(-t(j2,j3,j4),-s(j3,j4))* .((-2d0*za(j1,j5)*zb(j2,j4))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))+ .(-(za(j2,j5)*zb(j2,j4))-za(j3,j5)*zb(j3,j4))/t(j2,j3,j4)))/ .s(j3,j4)-(L1(-s(j3,j4),-t(j2,j3,j4))*za(j2,j5)*zb(j2,j4))/t(j2,j3, .j4)**2 .))/(2d0*za(j5,j6)*zb(j2,j3)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j .2,j4))) Fsc6=Fsc6+ .(za(j1,j5)**2*zb(j1,j4)*zb(j2,j4)* .((-2d0*L0(-s(j5,j6),-t(j2,j3,j4))*zb(j2,j4))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))+ .(L1(-s(j5,j6),-t(j2,j3,j4))*zb(j1,j4))/t(j2,j3,j4)))/ .(2d0*za(j5,j6)*zb(j2,j3)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j .4))* .zb(j3,j4))+(za(j3,j5)*zb(j2,j4)* .(za(j3,j5)/za(j3,j4)+(-(za(j2,j5)*zb(j2,j4))-za(j3,j5)*zb(j3,j4))/ .t(j2,j3,j4)))/ .(2d0*za(j5,j6)*zb(j2,j3)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j .4)))- .(Lsm1_2mh(s(j1,j2),t(j2,j3,j4),s(j3,j4),s(j5,j6))*za(j1,j5)**2*zb( .j2,j4)**3* .t(j2,j3,j4))/ .(za(j5,j6)*zb(j2,j3)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))* .*3* .zb(j3,j4))+(Lnrat(-s(j1,j2),-s(j5,j6))*zb(j2,j4)* .(((s(j1,j2)-s(j3,j4)-s(j5,j6))*zb(j1,j2)* .(-(za(j2,j5)*(-(za(j1,j5)*za(j2,j3))+za(j1,j2)*za(j3,j5))* .zb(j2,j4))+(za(j1,j3)*za(j1,j5)**2* .(-(za(j2,j3)*zb(j1,j3))-za(j2,j4)*zb(j1,j4))*zb(j2,j4) .)/(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))- .za(j1,j5)*za(j2,j5)* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))))/ .((s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2- .2d0*s(j1,j2)*s(j5,j6)-2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*za(j5,j6) .) .-(2d0*(-s(j1,j2)-s(j3,j4)+s(j5,j6))*za(j1,j2)*za(j3,j5)* .(zb(j1,j6)*zb(j2,j4)+zb(j1,j4)*zb(j2,j6)))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)- .(za(j1,j3)*zb(j2,j4)*(-((za(j1,j5)*za(j2,j5)*zb(j1,j2))/ .za(j5,j6))-(za(j1,j2)*zb(j1,j6)*zb(j2,j6))/zb(j5,j6)))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))- .(za(j2,j3)*zb(j2,j6)*zb(j4,j6))/zb(j5,j6)+ .((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j1,j2)* .(zb(j1,j4)*zb(j2,j6)* .(-(za(j1,j3)*zb(j1,j6))-za(j2,j3)*zb(j2,j6))- .za(j1,j3)*zb(j2,j4)* .(zb(j1,j6)**2+((-(za(j2,j3)*zb(j1,j3))- .za(j2,j4)*zb(j1,j4))*zb(j2,j6)**2)/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4)))+ .zb(j1,j2)*zb(j4,j6)* .(-(za(j2,j3)*zb(j2,j6))+za(j3,j4)*zb(j4,j6))))/ .((s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2- .2d0*s(j1,j2)*s(j5,j6)-2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*zb(j5,j6) .) .+(6d0*s(j1,j2)*s(j3,j4)*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4 .))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))* .(t(j1,j3,j4)-t(j2,j3,j4)))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)**2- .(s(j1,j2)*((za(j3,j5)**2*zb(j3,j4))/za(j5,j6)+ .(za(j3,j4)*zb(j4,j6)**2)/zb(j5,j6))*t(j2,j3,j4))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)))/ .(2d0*s(j3,j4)*zb(j2,j3)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4 .))) Fsc6=Fsc6+ .(Lnrat(-s(j1,j2),-s(j5,j6))*za(j1,j3)* .((za(j1,j5)*za(j3,j5)*zb(j1,j4))/za(j5,j6)+ .((s(j1,j2)-s(j3,j4)-s(j5,j6))*zb(j1,j2)* .(za(j1,j5)*za(j2,j3)* .(-(za(j1,j5)*zb(j1,j4))-za(j2,j5)*zb(j2,j4))- .za(j1,j3)*zb(j2,j4)* .(za(j2,j5)**2+(za(j1,j5)**2* .(-(za(j2,j3)*zb(j1,j3))-za(j2,j4)*zb(j1,j4)))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4)))- .za(j1,j2)*za(j3,j5)* .(-(za(j1,j5)*zb(j1,j4))-za(j3,j5)*zb(j3,j4))))/ .((s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2- .2d0*s(j1,j2)*s(j5,j6)-2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*za(j5,j6) .) .+(2d0*(-s(j1,j2)-s(j3,j4)+s(j5,j6))* .(za(j1,j5)*za(j2,j3)+za(j1,j3)*za(j2,j5))*zb(j1,j2)*zb(j4,j6) .)/(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5 .,j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)- .(za(j1,j3)*zb(j2,j4)*(-((za(j1,j5)*za(j2,j5)*zb(j1,j2))/ .za(j5,j6))-(za(j1,j2)*zb(j1,j6)*zb(j2,j6))/zb(j5,j6)))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))+ .((s(j1,j2)-s(j3,j4)-s(j5,j6))*za(j1,j2)* .(-(zb(j1,j6)*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .zb(j2,j6))+(za(j1,j3)* .(-(za(j2,j3)*zb(j1,j3))-za(j2,j4)*zb(j1,j4))* .zb(j2,j4)*zb(j2,j6)**2)/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))- .za(j1,j3)*zb(j1,j6)* .(-(zb(j1,j4)*zb(j2,j6))-zb(j1,j2)*zb(j4,j6))))/ .((s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2- .2d0*s(j1,j2)*s(j5,j6)-2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*zb(j5,j6) .) .-(s(j1,j2)*((za(j3,j5)**2*zb(j3,j4))/za(j5,j6)+ .(za(j3,j4)*zb(j4,j6)**2)/zb(j5,j6))*t(j1,j3,j4))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)+ .(6d0*s(j1,j2)*s(j3,j4)*(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4) .)* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))* .(-t(j1,j3,j4)+t(j2,j3,j4)))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)**2))/ .(2d0*s(j3,j4)*za(j1,j4)*(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4 .))) Fsc6=Fsc6+ .(I3m(s(j1,j2),s(j3,j4),s(j5,j6))*za(j1,j3)* .(-((za(j1,j3)*(-(za(j2,j3)*zb(j1,j3))-za(j2,j4)*zb(j1,j4))* .zb(j2,j4)*(-(za(j1,j3)* .(-(za(j1,j5)*zb(j1,j3))-za(j2,j5)*zb(j2,j3)))- .za(j1,j4)*(-(za(j1,j5)*zb(j1,j4))- .za(j2,j5)*zb(j2,j4))+ .za(j1,j2)*(za(j3,j5)*zb(j2,j3)+za(j4,j5)*zb(j2,j4)))* .zb(j2,j6))/(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4)))+ .zb(j1,j4)*(za(j1,j3)*(-(za(j1,j2)* .(za(j3,j5)*zb(j1,j3)+za(j4,j5)*zb(j1,j4)))- .za(j2,j3)*(-(za(j1,j5)*zb(j1,j3))-za(j2,j5)*zb(j2,j3))- .za(j2,j4)*(-(za(j1,j5)*zb(j1,j4))-za(j2,j5)*zb(j2,j4)))+ .za(j1,j2)*(-(za(j1,j3)* .(za(j3,j5)*zb(j1,j3)+za(j4,j5)*zb(j1,j4)))- .za(j2,j3)*(za(j3,j5)*zb(j2,j3)+za(j4,j5)*zb(j2,j4))))* .zb(j2,j6)+(s(j1,j2)* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6)))/2+ .za(j2,j3)*zb(j1,j2)*(-(za(j1,j2)*za(j3,j5)*zb(j2,j6)*zb(j3,j4))+ .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))* .(za(j1,j2)*zb(j2,j4)+za(j1,j3)*zb(j3,j4)))+ .(3d0*s(j1,j2)*(s(j1,j2)-s(j3,j4)-s(j5,j6))* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))* .(-t(j1,j3,j4)+t(j2,j3,j4)))/ .(2d0*(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2- .2d0*s(j1,j2)*s(j5,j6)-2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2))))/ .((s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5, .j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*za(j1,j4)* .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))) Fsc6=Fsc6+ .(Lnrat(-s(j3,j4),-s(j5,j6))*za(j1,j3)* .(-2d0*(za(j1,j3)*zb(j1,j4)-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))- .(8d0*s(j1,j2)*za(j1,j3)*za(j5,j6)*zb(j2,j6)*zb(j4,j6))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))+ .4d0*(-(za(j1,j2)*za(j3,j5)*zb(j1,j6)*zb(j2,j4))+ .za(j1,j5)*za(j2,j3)*zb(j1,j2)*zb(j4,j6))+ .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4)+ .(4d0*s(j5,j6)*za(j1,j3)*zb(j2,j4))/ .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4)))* .((za(j1,j5)*za(j2,j5)*zb(j1,j2))/za(j5,j6)- .(za(j1,j2)*zb(j1,j6)*zb(j2,j6))/zb(j5,j6))- .(-s(j1,j2)+s(j3,j4)-s(j5,j6))* .((za(j1,j5)*za(j3,j5)*zb(j1,j4))/za(j5,j6)- .(za(j2,j3)*zb(j2,j6)*zb(j4,j6))/zb(j5,j6)- .(2d0*za(j1,j3)*zb(j2,j6)* .(-(za(j1,j2)*zb(j1,j4)*zb(j2,j6))+ .za(j2,j3)*zb(j2,j6)*zb(j3,j4)+ .za(j1,j5)*zb(j1,j4)*zb(j5,j6)))/ .((-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))*zb(j5,j6)))+ .(3d0*(-s(j1,j2)-s(j3,j4)+s(j5,j6))* .(-(za(j1,j3)*zb(j1,j4))-za(j2,j3)*zb(j2,j4))* .(-(za(j1,j5)*zb(j1,j6))-za(j2,j5)*zb(j2,j6))* .(-t(j1,j3,j4)+t(j2,j3,j4)))/ .(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s(j5,j .6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)))/ .(2d0*(s(j1,j2)**2-2d0*s(j1,j2)*s(j3,j4)+s(j3,j4)**2-2d0*s(j1,j2)*s .(j5,j6)- .2d0*s(j3,j4)*s(j5,j6)+s(j5,j6)**2)*za(j1,j4)* .(-(za(j1,j3)*zb(j2,j3))-za(j1,j4)*zb(j2,j4))) return end