subroutine qqb_dm_monophot_v_Vamps(p,i1,i2,i3,i4,i5,qgqb) implicit none !----- combined colour and Born intefered amplitudes as a function of quark line helicity !------ default is q(i1)+g(i2)+qb(i3)+x(i4)+x(i5) include 'constants.f' include 'qcdcouple.f' ! include 'zprods_decl.f' double precision qgqb(2) integer i1,i2,i3,i4,i5 double precision p(mxpart,4) integer h1,h2,h3,h4 double complex amp_tree(2,2,2,2),amp_slc(2,2,2,2) !------ tree call qqb_dm_monojet_Vamps(p,i1,i2,i3,i4,i5,amp_tree) !------ subleading colour (swap i2 and i3) call qqb_dm_monophot_Vamps_fill(p,i1,i3,i2,i4,i5,amp_slc) qgqb(:)=0d0 do h1=1,2 do h2=1,2 do h3=1,2 do h4=1,2 qgqb(h1)=qgqb(h1) & -ason2pi*Dble(Dconjg(amp_tree(h1,h2,h3,h4))* & amp_slc(h1,h2,h3,h4)) enddo enddo enddo enddo return end subroutine qqb_dm_monophot_Vamps_fill(p,i1,i2,i3,i4,i5,amp) implicit none include 'constants.f' include 'dm_params.f' include 'zprods_decl.f' include 'sprods_com.f' include 'scale.f' include 'epinv.f' !----- SUB LEADING COLOR !----- fills amplitude for q qb g chi,chib double complex amp(2,2,2,2),Vamp(2,2,2,2),Famp(2,2,2,2) double complex amp_tree(2,2,2,2) double precision p(mxpart,4),q(mxpart,4) integer i1,i2,i3,i4,i5 integer h1,h2,h3,h4 double complex l12,l45,L0,Lsm1,L1 double complex lnrat double complex vfac amp(:,:,:,:)=czip if(xmass.gt.1d-8) then !--------- generate massless phase space call gen_masslessvecs(p,q,i4,i5) !--------- generate spinors call spinoru(5,q,za,zb) else !-------- massless dm can use usual spinoru call spinoru(5,p,za,zb) endif !------ basis integrals l12=lnrat(musq,-s(i1,i2)) l45=lnrat(musq,-s(i4,i5)) !----- use same notation as Born !----- split into Vamp and Famp with amp=Vamp+Famp !------ Vamp piece universial vfac= &-(epinv**2+epinv*l12+0.5d0*l12**2) & -2d0*(epinv+l45)-4d0 & +0.5d0*(epinv+l45)+0.5d0 !------ note have swapped order of i3 and i2 to calculate proper tree !------ need minus sign too call qqb_dm_monojet_Vamps(p,i1,i3,i2,i4,i5,amp_tree) do h1=1,2 do h2=1,2 do h3=1,2 do h4=1,2 Vamp(h1,h2,h3,h4)=-amp_tree(h1,h2,h3,h4)*vfac enddo enddo enddo enddo ! write(6,*) 'amp_tree ',amp_tree(2,2,1,2)/s(i4,i5) !------ Fpieces !====== Helicity conserving Famp(1,1,1,2)= -(za(i1,i3)*za(i1,i4)*zb(i5,i1))/ - (2.*za(i1,i2)*zb(i3,i1)) - - (za(i1,i3)**2*za(i2,i4)*zb(i5,i1))/ - (2.*za(i1,i2)*za(i2,i3)*zb(i3,i1)) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i1,i4)*zb(i2,i1)*zb(i5,i1))/ - (s(i1,i2)*zb(i3,i1)) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i3,i4)*zb(i3,i2)*zb(i5,i1))/ - (2.*s(i2,i3)**2) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i2,i4)*zb(i2,i1)*zb(i3,i2)*zb(i5,i1) - )/(s(i1,i2)*zb(i3,i1)**2) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i2,i4)*zb(i2,i1)*zb(i3,i2)*zb(i5,i1) - )/(s(i2,i3)*zb(i3,i1)**2) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i3,i4)*zb(i3,i2)*zb(i5,i1))/ - (s(i2,i3)*zb(i3,i1)) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i2,i4)*zb(i2,i1)*zb(i3,i2)* - zb(i5,i1))/(2.*s(i2,i3)**2*zb(i3,i1)) - - (2d0*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i3,i4)*zb(i5,i2))/ - (s(i1,i2) + s(i1,i3) + s(i2,i3)) - - (za(i1,i4)*za(i2,i3)*zb(i5,i2))/ - (2.*za(i1,i2)*zb(i3,i1)) - - (za(i1,i3)*za(i2,i4)*zb(i5,i2))/ - (2.*za(i1,i2)*zb(i3,i1)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i4)*zb(i5,i2))/zb(i3,i1) - - (2d0*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i2,i4)*zb(i2,i1)*zb(i5,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i3,i1)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i2,i1)*zb(i5,i2))/ - (zb(i3,i1)*zb(i3,i2)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i4)*zb(i2,i1)**2*zb(i5,i3))/zb(i3,i1)**3 + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)* - za(i3,i4)*zb(i2,i1)*zb(i5,i3))/ - (s(i1,i2)**2*zb(i3,i1)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i2,i1)**2*zb(i5,i3))/ - (zb(i3,i1)**2*zb(i3,i2)) Famp(2,1,1,2)=(-2*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i3)*za(i3,i4)*zb(i5,i1))/ - (s(i1,i2) + s(i1,i3) + s(i2,i3)) + - (za(i1,i4)*za(i2,i3)*zb(i5,i1))/ - (2.*za(i1,i2)*zb(i3,i2)) + - (za(i1,i3)*za(i2,i4)*zb(i5,i1))/ - (2.*za(i1,i2)*zb(i3,i2)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i4)*zb(i5,i1))/zb(i3,i2) + - (2d0*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*za(i2,i3)*zb(i2,i1)*zb(i5,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i3,i2)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i4)*zb(i2,i1)*zb(i5,i1))/ - (zb(i3,i1)*zb(i3,i2)) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i2,i3)**2*za(i3,i4)*zb(i3,i1)*zb(i5,i2))/ - (2.*s(i1,i3)**2) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i4)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i5,i2) - )/(s(i1,i2)*zb(i3,i2)**2) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i4)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i5,i2) - )/(s(i1,i3)*zb(i3,i2)**2) + - (za(i1,i4)*za(i2,i3)**2*zb(i5,i2))/ - (2.*za(i1,i2)*za(i1,i3)*zb(i3,i2)) + - (za(i2,i3)*za(i2,i4)*zb(i5,i2))/ - (2.*za(i1,i2)*zb(i3,i2)) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i2,i3)*za(i2,i4)*zb(i2,i1)*zb(i5,i2))/ - (s(i1,i2)*zb(i3,i2)) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i2,i3)*za(i3,i4)*zb(i3,i1)*zb(i5,i2))/ - (s(i1,i3)*zb(i3,i2)) + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i4)*za(i2,i3)**2*zb(i2,i1)*zb(i3,i1)* - zb(i5,i2))/(2.*s(i1,i3)**2*zb(i3,i2)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i2,i1)**2*zb(i5,i3))/zb(i3,i2)**3 - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i4)*zb(i2,i1)**2*zb(i5,i3))/ - (zb(i3,i1)*zb(i3,i2)**2) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)* - za(i3,i4)*zb(i2,i1)*zb(i5,i3))/ - (s(i1,i2)**2*zb(i3,i2)) Famp(1,2,1,2)= (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i4)*zb(i5,i1))/za(i2,i3)**3 - - (2d0*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i4)*zb(i3,i2)*zb(i5,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i3)*za(i2,i4)*zb(i3,i2)*zb(i5,i1) - )/(s(i1,i2)*za(i2,i3)**2) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i4)*zb(i3,i2)*zb(i5,i1) - )/(s(i1,i3)*za(i2,i3)**2) - - (za(i1,i4)*zb(i3,i2)*zb(i5,i1))/ - (2.*za(i2,i3)*zb(i2,i1)) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i4)*zb(i3,i2)**2* - zb(i5,i1))/(2.*s(i1,i3)**2*za(i2,i3)) - - (za(i2,i4)*zb(i3,i2)**2*zb(i5,i1))/ - (2.*za(i2,i3)*zb(i2,i1)*zb(i3,i1)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i4)*zb(i5,i2))/ - (za(i1,i3)*za(i2,i3)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i4)*zb(i5,i2))/ - (za(i1,i3)*za(i2,i3)**2) - - (za(i1,i4)*zb(i3,i1)*zb(i5,i2))/ - (2.*za(i2,i3)*zb(i2,i1)) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i2,i4)*zb(i3,i2)*zb(i5,i2))/ - (s(i1,i2)*za(i2,i3)) - - (za(i2,i4)*zb(i3,i2)*zb(i5,i2))/ - (2.*za(i2,i3)*zb(i2,i1)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i5,i3))/za(i2,i3) + - (2d0*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i3,i2)*zb(i5,i3))/ - (s(i1,i2) + s(i1,i3) + s(i2,i3)) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i4)*zb(i3,i2)*zb(i5,i3))/ - (s(i1,i3)*za(i2,i3)) + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i4)*zb(i3,i2)*zb(i5,i3))/ - (s(i1,i2)**2*za(i2,i3)) + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i4)*zb(i3,i2)**2*zb(i5,i3))/ - (2.*s(i1,i3)**2) Famp(2,2,1,2)=(Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i2,i4)*zb(i5,i1))/ - (za(i1,i3)*za(i2,i3)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i4)*zb(i5,i1))/ - (za(i1,i3)**2*za(i2,i3)) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i4)*zb(i3,i1)*zb(i5,i1))/ - (s(i1,i2)*za(i1,i3)) + - (za(i1,i4)*zb(i3,i1)*zb(i5,i1))/ - (2.*za(i1,i3)*zb(i2,i1)) + - (za(i2,i4)*zb(i3,i2)*zb(i5,i1))/ - (2.*za(i1,i3)*zb(i2,i1)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i4)*zb(i5,i2))/za(i1,i3)**3 - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i4)*za(i2,i3)*zb(i3,i1)*zb(i5,i2) - )/(s(i1,i2)*za(i1,i3)**2) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i4)*za(i2,i3)*zb(i3,i1)*zb(i5,i2) - )/(s(i2,i3)*za(i1,i3)**2) + - (2d0*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i2,i4)*zb(i3,i1)*zb(i5,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)) + - (za(i2,i4)*zb(i3,i1)*zb(i5,i2))/ - (2.*za(i1,i3)*zb(i2,i1)) + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i4)*za(i2,i3)*zb(i3,i1)**2* - zb(i5,i2))/(2.*s(i2,i3)**2*za(i1,i3)) + - (za(i1,i4)*zb(i3,i1)**2*zb(i5,i2))/ - (2.*za(i1,i3)*zb(i2,i1)*zb(i3,i2)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i4)*zb(i5,i3))/za(i1,i3) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i4)*za(i2,i3)*zb(i3,i1)*zb(i5,i3))/ - (s(i2,i3)*za(i1,i3)) + - (2d0*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i4)*zb(i3,i1)*zb(i5,i3))/ - (s(i1,i2) + s(i1,i3) + s(i2,i3)) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i4)*zb(i3,i1)*zb(i5,i3))/ - (s(i1,i2)**2*za(i1,i3)) + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i4)*za(i2,i3)*zb(i3,i1)**2*zb(i5,i3))/ - (2.*s(i2,i3)**2) Famp(1,1,2,1)=-(za(i1,i3)*za(i1,i5)*zb(i4,i1))/ - (2.*za(i1,i2)*zb(i3,i1)) - - (za(i1,i3)**2*za(i2,i5)*zb(i4,i1))/ - (2.*za(i1,i2)*za(i2,i3)*zb(i3,i1)) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i1,i5)*zb(i2,i1)*zb(i4,i1))/ - (s(i1,i2)*zb(i3,i1)) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i3,i5)*zb(i3,i2)*zb(i4,i1))/ - (2.*s(i2,i3)**2) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i2,i5)*zb(i2,i1)*zb(i3,i2)*zb(i4,i1) - )/(s(i1,i2)*zb(i3,i1)**2) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i2,i5)*zb(i2,i1)*zb(i3,i2)*zb(i4,i1) - )/(s(i2,i3)*zb(i3,i1)**2) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i3,i5)*zb(i3,i2)*zb(i4,i1))/ - (s(i2,i3)*zb(i3,i1)) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i2,i5)*zb(i2,i1)*zb(i3,i2)* - zb(i4,i1))/(2.*s(i2,i3)**2*zb(i3,i1)) - - (2d0*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i3,i5)*zb(i4,i2))/ - (s(i1,i2) + s(i1,i3) + s(i2,i3)) - - (za(i1,i5)*za(i2,i3)*zb(i4,i2))/ - (2.*za(i1,i2)*zb(i3,i1)) - - (za(i1,i3)*za(i2,i5)*zb(i4,i2))/ - (2.*za(i1,i2)*zb(i3,i1)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i5)*zb(i4,i2))/zb(i3,i1) - - (2d0*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i2,i5)*zb(i2,i1)*zb(i4,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i3,i1)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i2,i1)*zb(i4,i2))/ - (zb(i3,i1)*zb(i3,i2)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i5)*zb(i2,i1)**2*zb(i4,i3))/zb(i3,i1)**3 + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)* - za(i3,i5)*zb(i2,i1)*zb(i4,i3))/ - (s(i1,i2)**2*zb(i3,i1)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i2,i1)**2*zb(i4,i3))/ - (zb(i3,i1)**2*zb(i3,i2)) Famp(2,1,2,1)=(-2*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i3)*za(i3,i5)*zb(i4,i1))/ - (s(i1,i2) + s(i1,i3) + s(i2,i3)) + - (za(i1,i5)*za(i2,i3)*zb(i4,i1))/ - (2.*za(i1,i2)*zb(i3,i2)) + - (za(i1,i3)*za(i2,i5)*zb(i4,i1))/ - (2.*za(i1,i2)*zb(i3,i2)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i5)*zb(i4,i1))/zb(i3,i2) + - (2d0*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*za(i2,i3)*zb(i2,i1)*zb(i4,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i3,i2)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i5)*zb(i2,i1)*zb(i4,i1))/ - (zb(i3,i1)*zb(i3,i2)) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i2,i3)**2*za(i3,i5)*zb(i3,i1)*zb(i4,i2))/ - (2.*s(i1,i3)**2) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i5)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i4,i2) - )/(s(i1,i2)*zb(i3,i2)**2) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i5)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i4,i2) - )/(s(i1,i3)*zb(i3,i2)**2) + - (za(i1,i5)*za(i2,i3)**2*zb(i4,i2))/ - (2.*za(i1,i2)*za(i1,i3)*zb(i3,i2)) + - (za(i2,i3)*za(i2,i5)*zb(i4,i2))/ - (2.*za(i1,i2)*zb(i3,i2)) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i2,i3)*za(i2,i5)*zb(i2,i1)*zb(i4,i2))/ - (s(i1,i2)*zb(i3,i2)) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i2,i3)*za(i3,i5)*zb(i3,i1)*zb(i4,i2))/ - (s(i1,i3)*zb(i3,i2)) + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i5)*za(i2,i3)**2*zb(i2,i1)*zb(i3,i1)* - zb(i4,i2))/(2.*s(i1,i3)**2*zb(i3,i2)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i2,i1)**2*zb(i4,i3))/zb(i3,i2)**3 - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i5)*zb(i2,i1)**2*zb(i4,i3))/ - (zb(i3,i1)*zb(i3,i2)**2) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)* - za(i3,i5)*zb(i2,i1)*zb(i4,i3))/ - (s(i1,i2)**2*zb(i3,i2)) Famp(1,2,2,1)=(Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i4,i1))/za(i2,i3)**3 - - (2d0*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i5)*zb(i3,i2)*zb(i4,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i3)*za(i2,i5)*zb(i3,i2)*zb(i4,i1) - )/(s(i1,i2)*za(i2,i3)**2) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i5)*zb(i3,i2)*zb(i4,i1) - )/(s(i1,i3)*za(i2,i3)**2) - - (za(i1,i5)*zb(i3,i2)*zb(i4,i1))/ - (2.*za(i2,i3)*zb(i2,i1)) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i5)*zb(i3,i2)**2* - zb(i4,i1))/(2.*s(i1,i3)**2*za(i2,i3)) - - (za(i2,i5)*zb(i3,i2)**2*zb(i4,i1))/ - (2.*za(i2,i3)*zb(i2,i1)*zb(i3,i1)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i5)*zb(i4,i2))/ - (za(i1,i3)*za(i2,i3)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i4,i2))/ - (za(i1,i3)*za(i2,i3)**2) - - (za(i1,i5)*zb(i3,i1)*zb(i4,i2))/ - (2.*za(i2,i3)*zb(i2,i1)) + - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i2,i5)*zb(i3,i2)*zb(i4,i2))/ - (s(i1,i2)*za(i2,i3)) - - (za(i2,i5)*zb(i3,i2)*zb(i4,i2))/ - (2.*za(i2,i3)*zb(i2,i1)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i4,i3))/za(i2,i3) + - (2d0*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i3,i2)*zb(i4,i3))/ - (s(i1,i2) + s(i1,i3) + s(i2,i3)) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i5)*zb(i3,i2)*zb(i4,i3))/ - (s(i1,i3)*za(i2,i3)) + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i5)*zb(i3,i2)*zb(i4,i3))/ - (s(i1,i2)**2*za(i2,i3)) + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i5)*zb(i3,i2)**2*zb(i4,i3))/ - (2.*s(i1,i3)**2) Famp(2,2,2,1)=(Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i2,i5)*zb(i4,i1))/ - (za(i1,i3)*za(i2,i3)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i4,i1))/ - (za(i1,i3)**2*za(i2,i3)) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i5)*zb(i3,i1)*zb(i4,i1))/ - (s(i1,i2)*za(i1,i3)) + - (za(i1,i5)*zb(i3,i1)*zb(i4,i1))/ - (2.*za(i1,i3)*zb(i2,i1)) + - (za(i2,i5)*zb(i3,i2)*zb(i4,i1))/ - (2.*za(i1,i3)*zb(i2,i1)) + - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i4,i2))/za(i1,i3)**3 - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i5)*za(i2,i3)*zb(i3,i1)*zb(i4,i2) - )/(s(i1,i2)*za(i1,i3)**2) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i5)*za(i2,i3)*zb(i3,i1)*zb(i4,i2) - )/(s(i2,i3)*za(i1,i3)**2) + - (2d0*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i2,i5)*zb(i3,i1)*zb(i4,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)) + - (za(i2,i5)*zb(i3,i1)*zb(i4,i2))/ - (2.*za(i1,i3)*zb(i2,i1)) + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i5)*za(i2,i3)*zb(i3,i1)**2* - zb(i4,i2))/(2.*s(i2,i3)**2*za(i1,i3)) + - (za(i1,i5)*zb(i3,i1)**2*zb(i4,i2))/ - (2.*za(i1,i3)*zb(i2,i1)*zb(i3,i2)) - - (Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i5)*zb(i4,i3))/za(i1,i3) - - (L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i5)*za(i2,i3)*zb(i3,i1)*zb(i4,i3))/ - (s(i2,i3)*za(i1,i3)) + - (2d0*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i5)*zb(i3,i1)*zb(i4,i3))/ - (s(i1,i2) + s(i1,i3) + s(i2,i3)) - - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i5)*zb(i3,i1)*zb(i4,i3))/ - (s(i1,i2)**2*za(i1,i3)) + - (L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i5)*za(i2,i3)*zb(i3,i1)**2*zb(i4,i3))/ - (2.*s(i2,i3)**2) !======= Helicity violating Famp(1,1,1,1)=(xmass*za(i1,i3)*za(i1,i4)*zb(i4,i1))/ - (2.*za(i1,i2)*zb(i3,i1)*zb(i5,i4)) + - (xmass*za(i1,i3)**2*za(i2,i4)*zb(i4,i1))/ - (2.*za(i1,i2)*za(i2,i3)*zb(i3,i1)*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i1,i4)*zb(i2,i1)*zb(i4,i1))/ - (s(i1,i2)*zb(i3,i1)*zb(i5,i4)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i3,i4)*zb(i3,i2)*zb(i4,i1))/ - (2.*s(i2,i3)**2*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i2,i4)*zb(i2,i1)*zb(i3,i2)*zb(i4,i1))/ - (s(i1,i2)*zb(i3,i1)**2*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i2,i4)*zb(i2,i1)*zb(i3,i2)*zb(i4,i1))/ - (s(i2,i3)*zb(i3,i1)**2*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i3,i4)*zb(i3,i2)*zb(i4,i1))/ - (s(i2,i3)*zb(i3,i1)*zb(i5,i4)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i2,i4)*zb(i2,i1)*zb(i3,i2)*zb(i4,i1) - )/(2.*s(i2,i3)**2*zb(i3,i1)*zb(i5,i4)) + - (2d0*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i3,i4)*zb(i4,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i5,i4)) + - (xmass*za(i1,i4)*za(i2,i3)*zb(i4,i2))/ - (2.*za(i1,i2)*zb(i3,i1)*zb(i5,i4)) + - (xmass*za(i1,i3)*za(i2,i4)*zb(i4,i2))/ - (2.*za(i1,i2)*zb(i3,i1)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i4)*zb(i4,i2))/(zb(i3,i1)*zb(i5,i4)) + - (2d0*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i2,i4)*zb(i2,i1)*zb(i4,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i3,i1)*zb(i5,i4)) - + (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i1,i4)* - zb(i2,i1)*zb(i4,i2))/(zb(i3,i1)*zb(i3,i2)*zb(i5,i4)) - + (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i2,i4)* - zb(i2,i1)**2*zb(i4,i3))/(zb(i3,i1)**3*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)* - za(i3,i4)*zb(i2,i1)*zb(i4,i3))/ - (s(i1,i2)**2*zb(i3,i1)*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i2,i1)**2*zb(i4,i3))/ - (zb(i3,i1)**2*zb(i3,i2)*zb(i5,i4)) - - (xmass*za(i1,i3)*za(i1,i5)*zb(i5,i1))/ - (2.*za(i1,i2)*zb(i3,i1)*zb(i5,i4)) - - (xmass*za(i1,i3)**2*za(i2,i5)*zb(i5,i1))/ - (2.*za(i1,i2)*za(i2,i3)*zb(i3,i1)*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i1,i5)*zb(i2,i1)*zb(i5,i1))/ - (s(i1,i2)*zb(i3,i1)*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i3,i5)*zb(i3,i2)*zb(i5,i1))/ - (2.*s(i2,i3)**2*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i2,i5)*zb(i2,i1)*zb(i3,i2)*zb(i5,i1))/ - (s(i1,i2)*zb(i3,i1)**2*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i2,i5)*zb(i2,i1)*zb(i3,i2)*zb(i5,i1))/ - (s(i2,i3)*zb(i3,i1)**2*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i3,i5)*zb(i3,i2)*zb(i5,i1))/ - (s(i2,i3)*zb(i3,i1)*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i2,i5)*zb(i2,i1)*zb(i3,i2)*zb(i5,i1) - )/(2.*s(i2,i3)**2*zb(i3,i1)*zb(i5,i4)) - - (2d0*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i3,i5)*zb(i5,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i5,i4)) - - (xmass*za(i1,i5)*za(i2,i3)*zb(i5,i2))/ - (2.*za(i1,i2)*zb(i3,i1)*zb(i5,i4)) - - (xmass*za(i1,i3)*za(i2,i5)*zb(i5,i2))/ - (2.*za(i1,i2)*zb(i3,i1)*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i5)*zb(i5,i2))/(zb(i3,i1)*zb(i5,i4)) - - (2d0*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i2,i5)*zb(i2,i1)*zb(i5,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i3,i1)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i1,i5)* - zb(i2,i1)*zb(i5,i2))/(zb(i3,i1)*zb(i3,i2)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i2,i5)* - zb(i2,i1)**2*zb(i5,i3))/(zb(i3,i1)**3*zb(i5,i4)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)* - za(i3,i5)*zb(i2,i1)*zb(i5,i3))/ - (s(i1,i2)**2*zb(i3,i1)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i2,i1)**2*zb(i5,i3))/ - (zb(i3,i1)**2*zb(i3,i2)*zb(i5,i4)) Famp(2,1,1,1)=(2d0*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) & - s(i2,i3))* - za(i2,i3)*za(i3,i4)*zb(i4,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i5,i4)) - - (xmass*za(i1,i4)*za(i2,i3)*zb(i4,i1))/ - (2.*za(i1,i2)*zb(i3,i2)*zb(i5,i4)) - - (xmass*za(i1,i3)*za(i2,i4)*zb(i4,i1))/ - (2.*za(i1,i2)*zb(i3,i2)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i4)*zb(i4,i1))/(zb(i3,i2)*zb(i5,i4)) - - (2d0*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*za(i2,i3)*zb(i2,i1)*zb(i4,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i3,i2)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i2,i4)* - zb(i2,i1)*zb(i4,i1))/(zb(i3,i1)*zb(i3,i2)*zb(i5,i4)) - + (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3))*za(i2,i3)**2*za(i3,i4)*zb(i3,i1)* - zb(i4,i2))/(2.*s(i1,i3)**2*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i4)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i4,i2))/ - (s(i1,i2)*zb(i3,i2)**2*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i4)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i4,i2))/ - (s(i1,i3)*zb(i3,i2)**2*zb(i5,i4)) - - (xmass*za(i1,i4)*za(i2,i3)**2*zb(i4,i2))/ - (2.*za(i1,i2)*za(i1,i3)*zb(i3,i2)*zb(i5,i4)) - - (xmass*za(i2,i3)*za(i2,i4)*zb(i4,i2))/ - (2.*za(i1,i2)*zb(i3,i2)*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i2,i3)*za(i2,i4)*zb(i2,i1)*zb(i4,i2))/ - (s(i1,i2)*zb(i3,i2)*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i2,i3)*za(i3,i4)*zb(i3,i1)*zb(i4,i2))/ - (s(i1,i3)*zb(i3,i2)*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i4)*za(i2,i3)**2*zb(i2,i1)*zb(i3,i1)*zb(i4,i2) - )/(2.*s(i1,i3)**2*zb(i3,i2)*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i2,i1)**2*zb(i4,i3))/ - (zb(i3,i2)**3*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i4)*zb(i2,i1)**2*zb(i4,i3))/ - (zb(i3,i1)*zb(i3,i2)**2*zb(i5,i4)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)* - za(i3,i4)*zb(i2,i1)*zb(i4,i3))/ - (s(i1,i2)**2*zb(i3,i2)*zb(i5,i4)) - - (2d0*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i3)*za(i3,i5)*zb(i5,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i5,i4)) + - (xmass*za(i1,i5)*za(i2,i3)*zb(i5,i1))/ - (2.*za(i1,i2)*zb(i3,i2)*zb(i5,i4)) + - (xmass*za(i1,i3)*za(i2,i5)*zb(i5,i1))/ - (2.*za(i1,i2)*zb(i3,i2)*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i5)*zb(i5,i1))/(zb(i3,i2)*zb(i5,i4)) + - (2d0*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*za(i2,i3)*zb(i2,i1)*zb(i5,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i3,i2)*zb(i5,i4)) - + (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i2,i5)* - zb(i2,i1)*zb(i5,i1))/(zb(i3,i1)*zb(i3,i2)*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3))*za(i2,i3)**2*za(i3,i5)*zb(i3,i1)* - zb(i5,i2))/(2.*s(i1,i3)**2*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i5)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i5,i2))/ - (s(i1,i2)*zb(i3,i2)**2*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i5)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i5,i2))/ - (s(i1,i3)*zb(i3,i2)**2*zb(i5,i4)) + - (xmass*za(i1,i5)*za(i2,i3)**2*zb(i5,i2))/ - (2.*za(i1,i2)*za(i1,i3)*zb(i3,i2)*zb(i5,i4)) + - (xmass*za(i2,i3)*za(i2,i5)*zb(i5,i2))/ - (2.*za(i1,i2)*zb(i3,i2)*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i2,i3)*za(i2,i5)*zb(i2,i1)*zb(i5,i2))/ - (s(i1,i2)*zb(i3,i2)*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i2,i3)*za(i3,i5)*zb(i3,i1)*zb(i5,i2))/ - (s(i1,i3)*zb(i3,i2)*zb(i5,i4)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i5)*za(i2,i3)**2*zb(i2,i1)*zb(i3,i1)*zb(i5,i2) - )/(2.*s(i1,i3)**2*zb(i3,i2)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i2,i1)**2*zb(i5,i3))/ - (zb(i3,i2)**3*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i5)*zb(i2,i1)**2*zb(i5,i3))/ - (zb(i3,i1)*zb(i3,i2)**2*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)* - za(i3,i5)*zb(i2,i1)*zb(i5,i3))/ - (s(i1,i2)**2*zb(i3,i2)*zb(i5,i4)) Famp(1,2,1,1)=-((xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i1,i2)**2* - za(i3,i4)*zb(i4,i1))/(za(i2,i3)**3*zb(i5,i4))) + - (2d0*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i4)*zb(i3,i2)*zb(i4,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i2))*za(i1,i2)*za(i1,i3)*za(i2,i4)*zb(i3,i2)* - zb(i4,i1))/(s(i1,i2)*za(i2,i3)**2*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i4)*zb(i3,i2)*zb(i4,i1))/ - (s(i1,i3)*za(i2,i3)**2*zb(i5,i4)) + - (xmass*za(i1,i4)*zb(i3,i2)*zb(i4,i1))/ - (2.*za(i2,i3)*zb(i2,i1)*zb(i5,i4)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i4)*zb(i3,i2)**2*zb(i4,i1) - )/(2.*s(i1,i3)**2*za(i2,i3)*zb(i5,i4)) + - (xmass*za(i2,i4)*zb(i3,i2)**2*zb(i4,i1))/ - (2.*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i4)*zb(i4,i2))/ - (za(i1,i3)*za(i2,i3)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i4)*zb(i4,i2))/ - (za(i1,i3)*za(i2,i3)**2*zb(i5,i4)) + - (xmass*za(i1,i4)*zb(i3,i1)*zb(i4,i2))/ - (2.*za(i2,i3)*zb(i2,i1)*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i2,i4)*zb(i3,i2)*zb(i4,i2))/ - (s(i1,i2)*za(i2,i3)*zb(i5,i4)) + - (xmass*za(i2,i4)*zb(i3,i2)*zb(i4,i2))/ - (2.*za(i2,i3)*zb(i2,i1)*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i4,i3))/(za(i2,i3)*zb(i5,i4)) - - (2d0*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i3,i2)*zb(i4,i3))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i4)*zb(i3,i2)*zb(i4,i3))/ - (s(i1,i3)*za(i2,i3)*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i4)*zb(i3,i2)*zb(i4,i3))/ - (s(i1,i2)**2*za(i2,i3)*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i4)*zb(i3,i2)**2*zb(i4,i3))/ - (2.*s(i1,i3)**2*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i5,i1))/ - (za(i2,i3)**3*zb(i5,i4)) - - (2d0*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i5)*zb(i3,i2)*zb(i5,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)*zb(i5,i4)) - + (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i2))*za(i1,i2)*za(i1,i3)*za(i2,i5)*zb(i3,i2)* - zb(i5,i1))/(s(i1,i2)*za(i2,i3)**2*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i5)*zb(i3,i2)*zb(i5,i1))/ - (s(i1,i3)*za(i2,i3)**2*zb(i5,i4)) - - (xmass*za(i1,i5)*zb(i3,i2)*zb(i5,i1))/ - (2.*za(i2,i3)*zb(i2,i1)*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i5)*zb(i3,i2)**2*zb(i5,i1) - )/(2.*s(i1,i3)**2*za(i2,i3)*zb(i5,i4)) - - (xmass*za(i2,i5)*zb(i3,i2)**2*zb(i5,i1))/ - (2.*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i5)*zb(i5,i2))/ - (za(i1,i3)*za(i2,i3)*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i5,i2))/ - (za(i1,i3)*za(i2,i3)**2*zb(i5,i4)) - - (xmass*za(i1,i5)*zb(i3,i1)*zb(i5,i2))/ - (2.*za(i2,i3)*zb(i2,i1)*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i2,i5)*zb(i3,i2)*zb(i5,i2))/ - (s(i1,i2)*za(i2,i3)*zb(i5,i4)) - - (xmass*za(i2,i5)*zb(i3,i2)*zb(i5,i2))/ - (2.*za(i2,i3)*zb(i2,i1)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i5,i3))/(za(i2,i3)*zb(i5,i4)) + - (2d0*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i3,i2)*zb(i5,i3))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i5)*zb(i3,i2)*zb(i5,i3))/ - (s(i1,i3)*za(i2,i3)*zb(i5,i4)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i5)*zb(i3,i2)*zb(i5,i3))/ - (s(i1,i2)**2*za(i2,i3)*zb(i5,i4)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i5)*zb(i3,i2)**2*zb(i5,i3))/ - (2.*s(i1,i3)**2*zb(i5,i4)) Famp(2,2,1,1)=-((xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i1,i2)* - za(i2,i4)*zb(i4,i1))/ - (za(i1,i3)*za(i2,i3)*zb(i5,i4))) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i4)*zb(i4,i1))/ - (za(i1,i3)**2*za(i2,i3)*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i4)*zb(i3,i1)*zb(i4,i1))/ - (s(i1,i2)*za(i1,i3)*zb(i5,i4)) - - (xmass*za(i1,i4)*zb(i3,i1)*zb(i4,i1))/ - (2.*za(i1,i3)*zb(i2,i1)*zb(i5,i4)) - - (xmass*za(i2,i4)*zb(i3,i2)*zb(i4,i1))/ - (2.*za(i1,i3)*zb(i2,i1)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i4)*zb(i4,i2))/ - (za(i1,i3)**3*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i4)*za(i2,i3)*zb(i3,i1)*zb(i4,i2))/ - (s(i1,i2)*za(i1,i3)**2*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i4)*za(i2,i3)*zb(i3,i1)*zb(i4,i2))/ - (s(i2,i3)*za(i1,i3)**2*zb(i5,i4)) - - (2d0*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i2,i4)*zb(i3,i1)*zb(i4,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)*zb(i5,i4)) - - (xmass*za(i2,i4)*zb(i3,i1)*zb(i4,i2))/ - (2.*za(i1,i3)*zb(i2,i1)*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i4)*za(i2,i3)*zb(i3,i1)**2*zb(i4,i2) - )/(2.*s(i2,i3)**2*za(i1,i3)*zb(i5,i4)) - - (xmass*za(i1,i4)*zb(i3,i1)**2*zb(i4,i2))/ - (2.*za(i1,i3)*zb(i2,i1)*zb(i3,i2)*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i4)*zb(i4,i3))/(za(i1,i3)*zb(i5,i4)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i4)*za(i2,i3)*zb(i3,i1)*zb(i4,i3))/ - (s(i2,i3)*za(i1,i3)*zb(i5,i4)) - - (2d0*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i4)*zb(i3,i1)*zb(i4,i3))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i5,i4)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i4)*zb(i3,i1)*zb(i4,i3))/ - (s(i1,i2)**2*za(i1,i3)*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i4)*za(i2,i3)*zb(i3,i1)**2*zb(i4,i3))/ - (2.*s(i2,i3)**2*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i2,i5)*zb(i5,i1))/ - (za(i1,i3)*za(i2,i3)*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i5,i1))/ - (za(i1,i3)**2*za(i2,i3)*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i5)*zb(i3,i1)*zb(i5,i1))/ - (s(i1,i2)*za(i1,i3)*zb(i5,i4)) + - (xmass*za(i1,i5)*zb(i3,i1)*zb(i5,i1))/ - (2.*za(i1,i3)*zb(i2,i1)*zb(i5,i4)) + - (xmass*za(i2,i5)*zb(i3,i2)*zb(i5,i1))/ - (2.*za(i1,i3)*zb(i2,i1)*zb(i5,i4)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i5,i2))/ - (za(i1,i3)**3*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i5)*za(i2,i3)*zb(i3,i1)*zb(i5,i2))/ - (s(i1,i2)*za(i1,i3)**2*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i5)*za(i2,i3)*zb(i3,i1)*zb(i5,i2))/ - (s(i2,i3)*za(i1,i3)**2*zb(i5,i4)) + - (2*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i2,i5)*zb(i3,i1)*zb(i5,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)*zb(i5,i4)) - + (xmass*za(i2,i5)*zb(i3,i1)*zb(i5,i2))/ - (2.*za(i1,i3)*zb(i2,i1)*zb(i5,i4)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i5)*za(i2,i3)*zb(i3,i1)**2*zb(i5,i2) - )/(2.*s(i2,i3)**2*za(i1,i3)*zb(i5,i4)) + - (xmass*za(i1,i5)*zb(i3,i1)**2*zb(i5,i2))/ - (2.*za(i1,i3)*zb(i2,i1)*zb(i3,i2)*zb(i5,i4)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i5)*zb(i5,i3))/(za(i1,i3)*zb(i5,i4)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i5)*za(i2,i3)*zb(i3,i1)*zb(i5,i3))/ - (s(i2,i3)*za(i1,i3)*zb(i5,i4)) + - (2*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i5)*zb(i3,i1)*zb(i5,i3))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*zb(i5,i4)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i5)*zb(i3,i1)*zb(i5,i3))/ - (s(i1,i2)**2*za(i1,i3)*zb(i5,i4)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i5)*za(i2,i3)*zb(i3,i1)**2*zb(i5,i3))/ - (2.*s(i2,i3)**2*zb(i5,i4)) Famp(1,1,2,2)=-(xmass*za(i1,i3)*za(i1,i4)*zb(i4,i1))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i1)) - - (xmass*za(i1,i3)**2*za(i2,i4)*zb(i4,i1))/ - (2.*za(i1,i2)*za(i2,i3)*za(i4,i5)*zb(i3,i1)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i1,i4)*zb(i2,i1)*zb(i4,i1))/ - (s(i1,i2)*za(i4,i5)*zb(i3,i1)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i3,i4)*zb(i3,i2)*zb(i4,i1))/ - (2.*s(i2,i3)**2*za(i4,i5)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i2,i4)*zb(i2,i1)*zb(i3,i2)*zb(i4,i1))/ - (s(i1,i2)*za(i4,i5)*zb(i3,i1)**2) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i2,i4)*zb(i2,i1)*zb(i3,i2)*zb(i4,i1))/ - (s(i2,i3)*za(i4,i5)*zb(i3,i1)**2) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i3,i4)*zb(i3,i2)*zb(i4,i1))/ - (s(i2,i3)*za(i4,i5)*zb(i3,i1)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i2,i4)*zb(i2,i1)*zb(i3,i2)*zb(i4,i1) - )/(2.*s(i2,i3)**2*za(i4,i5)*zb(i3,i1)) - - (2*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i3,i4)*zb(i4,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)) - - (xmass*za(i1,i4)*za(i2,i3)*zb(i4,i2))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i1)) - - (xmass*za(i1,i3)*za(i2,i4)*zb(i4,i2))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i1)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i4)*zb(i4,i2))/(za(i4,i5)*zb(i3,i1)) - - (2*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i2,i4)*zb(i2,i1)*zb(i4,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)*zb(i3,i1)) - - (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i1,i4)* - zb(i2,i1)*zb(i4,i2))/(za(i4,i5)*zb(i3,i1)*zb(i3,i2)) - - (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i2,i4)* - zb(i2,i1)**2*zb(i4,i3))/(za(i4,i5)*zb(i3,i1)**3) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)* - za(i3,i4)*zb(i2,i1)*zb(i4,i3))/ - (s(i1,i2)**2*za(i4,i5)*zb(i3,i1)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i2,i1)**2*zb(i4,i3))/ - (za(i4,i5)*zb(i3,i1)**2*zb(i3,i2)) + - (xmass*za(i1,i3)*za(i1,i5)*zb(i5,i1))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i1)) + - (xmass*za(i1,i3)**2*za(i2,i5)*zb(i5,i1))/ - (2.*za(i1,i2)*za(i2,i3)*za(i4,i5)*zb(i3,i1)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i1,i5)*zb(i2,i1)*zb(i5,i1))/ - (s(i1,i2)*za(i4,i5)*zb(i3,i1)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i3,i5)*zb(i3,i2)*zb(i5,i1))/ - (2.*s(i2,i3)**2*za(i4,i5)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i3)*za(i2,i5)*zb(i2,i1)*zb(i3,i2)*zb(i5,i1))/ - (s(i1,i2)*za(i4,i5)*zb(i3,i1)**2) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i2,i5)*zb(i2,i1)*zb(i3,i2)*zb(i5,i1))/ - (s(i2,i3)*za(i4,i5)*zb(i3,i1)**2) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)*za(i3,i5)*zb(i3,i2)*zb(i5,i1))/ - (s(i2,i3)*za(i4,i5)*zb(i3,i1)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i3)**2*za(i2,i5)*zb(i2,i1)*zb(i3,i2)*zb(i5,i1) - )/(2.*s(i2,i3)**2*za(i4,i5)*zb(i3,i1)) + - (2*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i3,i5)*zb(i5,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)) + - (xmass*za(i1,i5)*za(i2,i3)*zb(i5,i2))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i1)) + - (xmass*za(i1,i3)*za(i2,i5)*zb(i5,i2))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i1)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i5)*zb(i5,i2))/(za(i4,i5)*zb(i3,i1)) + - (2*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i3)*za(i2,i5)*zb(i2,i1)*zb(i5,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)*zb(i3,i1)) - + (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i1,i5)* - zb(i2,i1)*zb(i5,i2))/(za(i4,i5)*zb(i3,i1)*zb(i3,i2)) - + (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i2,i5)* - zb(i2,i1)**2*zb(i5,i3))/(za(i4,i5)*zb(i3,i1)**3) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)* - za(i3,i5)*zb(i2,i1)*zb(i5,i3))/ - (s(i1,i2)**2*za(i4,i5)*zb(i3,i1)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i2,i1)**2*zb(i5,i3))/ - (za(i4,i5)*zb(i3,i1)**2*zb(i3,i2)) Famp(2,1,2,2)=(-2*xmass*L0(-s(i1,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i2,i3)* - za(i3,i4)*zb(i4,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)) + - (xmass*za(i1,i4)*za(i2,i3)*zb(i4,i1))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i2)) + - (xmass*za(i1,i3)*za(i2,i4)*zb(i4,i1))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i2)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i4)*zb(i4,i1))/(za(i4,i5)*zb(i3,i2)) + - (2*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*za(i2,i3)*zb(i2,i1)*zb(i4,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)*zb(i3,i2)) - + (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i2,i4)* - zb(i2,i1)*zb(i4,i1))/(za(i4,i5)*zb(i3,i1)*zb(i3,i2)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3))*za(i2,i3)**2*za(i3,i4)*zb(i3,i1)* - zb(i4,i2))/(2.*s(i1,i3)**2*za(i4,i5)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i4)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i4,i2))/ - (s(i1,i2)*za(i4,i5)*zb(i3,i2)**2) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i4)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i4,i2))/ - (s(i1,i3)*za(i4,i5)*zb(i3,i2)**2) + - (xmass*za(i1,i4)*za(i2,i3)**2*zb(i4,i2))/ - (2.*za(i1,i2)*za(i1,i3)*za(i4,i5)*zb(i3,i2)) + - (xmass*za(i2,i3)*za(i2,i4)*zb(i4,i2))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i2)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i2,i3)*za(i2,i4)*zb(i2,i1)*zb(i4,i2))/ - (s(i1,i2)*za(i4,i5)*zb(i3,i2)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i2,i3)*za(i3,i4)*zb(i3,i1)*zb(i4,i2))/ - (s(i1,i3)*za(i4,i5)*zb(i3,i2)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i4)*za(i2,i3)**2*zb(i2,i1)*zb(i3,i1)*zb(i4,i2) - )/(2.*s(i1,i3)**2*za(i4,i5)*zb(i3,i2)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i2,i1)**2*zb(i4,i3))/ - (za(i4,i5)*zb(i3,i2)**3) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i4)*zb(i2,i1)**2*zb(i4,i3))/ - (za(i4,i5)*zb(i3,i1)*zb(i3,i2)**2) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)* - za(i3,i4)*zb(i2,i1)*zb(i4,i3))/ - (s(i1,i2)**2*za(i4,i5)*zb(i3,i2)) + - (2*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i3)*za(i3,i5)*zb(i5,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)) - - (xmass*za(i1,i5)*za(i2,i3)*zb(i5,i1))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i2)) - - (xmass*za(i1,i3)*za(i2,i5)*zb(i5,i1))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i2)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i3,i5)*zb(i5,i1))/(za(i4,i5)*zb(i3,i2)) - - (2*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*za(i2,i3)*zb(i2,i1)*zb(i5,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)*zb(i3,i2)) - - (xmass*Lsm1(-s(i1,i2), - -s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3), - -s(i1,i2) - s(i1,i3) - s(i2,i3))*za(i2,i5)* - zb(i2,i1)*zb(i5,i1))/(za(i4,i5)*zb(i3,i1)*zb(i3,i2)) - + (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3))*za(i2,i3)**2*za(i3,i5)*zb(i3,i1)* - zb(i5,i2))/(2.*s(i1,i3)**2*za(i4,i5)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i5)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i5,i2))/ - (s(i1,i2)*za(i4,i5)*zb(i3,i2)**2) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i5)*za(i2,i3)*zb(i2,i1)*zb(i3,i1)*zb(i5,i2))/ - (s(i1,i3)*za(i4,i5)*zb(i3,i2)**2) - - (xmass*za(i1,i5)*za(i2,i3)**2*zb(i5,i2))/ - (2.*za(i1,i2)*za(i1,i3)*za(i4,i5)*zb(i3,i2)) - - (xmass*za(i2,i3)*za(i2,i5)*zb(i5,i2))/ - (2.*za(i1,i2)*za(i4,i5)*zb(i3,i2)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i2,i3)*za(i2,i5)*zb(i2,i1)*zb(i5,i2))/ - (s(i1,i2)*za(i4,i5)*zb(i3,i2)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i2,i3)*za(i3,i5)*zb(i3,i1)*zb(i5,i2))/ - (s(i1,i3)*za(i4,i5)*zb(i3,i2)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i5)*za(i2,i3)**2*zb(i2,i1)*zb(i3,i1)*zb(i5,i2) - )/(2.*s(i1,i3)**2*za(i4,i5)*zb(i3,i2)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i2,i1)**2*zb(i5,i3))/ - (za(i4,i5)*zb(i3,i2)**3) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i5)*zb(i2,i1)**2*zb(i5,i3))/ - (za(i4,i5)*zb(i3,i1)*zb(i3,i2)**2) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)* - za(i3,i5)*zb(i2,i1)*zb(i5,i3))/ - (s(i1,i2)**2*za(i4,i5)*zb(i3,i2)) Famp(1,2,2,2)=(xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) & - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i4)*zb(i4,i1))/ - (za(i2,i3)**3*za(i4,i5)) - - (2*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i4)*zb(i3,i2)*zb(i4,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)*za(i4,i5)) - + (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i2))*za(i1,i2)*za(i1,i3)*za(i2,i4)*zb(i3,i2)* - zb(i4,i1))/(s(i1,i2)*za(i2,i3)**2*za(i4,i5)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i4)*zb(i3,i2)*zb(i4,i1))/ - (s(i1,i3)*za(i2,i3)**2*za(i4,i5)) - - (xmass*za(i1,i4)*zb(i3,i2)*zb(i4,i1))/ - (2.*za(i2,i3)*za(i4,i5)*zb(i2,i1)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i4)*zb(i3,i2)**2*zb(i4,i1) - )/(2.*s(i1,i3)**2*za(i2,i3)*za(i4,i5)) - - (xmass*za(i2,i4)*zb(i3,i2)**2*zb(i4,i1))/ - (2.*za(i2,i3)*za(i4,i5)*zb(i2,i1)*zb(i3,i1)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i4)*zb(i4,i2))/ - (za(i1,i3)*za(i2,i3)*za(i4,i5)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i4)*zb(i4,i2))/ - (za(i1,i3)*za(i2,i3)**2*za(i4,i5)) - - (xmass*za(i1,i4)*zb(i3,i1)*zb(i4,i2))/ - (2.*za(i2,i3)*za(i4,i5)*zb(i2,i1)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i2,i4)*zb(i3,i2)*zb(i4,i2))/ - (s(i1,i2)*za(i2,i3)*za(i4,i5)) - - (xmass*za(i2,i4)*zb(i3,i2)*zb(i4,i2))/ - (2.*za(i2,i3)*za(i4,i5)*zb(i2,i1)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i4,i3))/(za(i2,i3)*za(i4,i5)) + - (2*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i4)*zb(i3,i2)*zb(i4,i3))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i4)*zb(i3,i2)*zb(i4,i3))/ - (s(i1,i3)*za(i2,i3)*za(i4,i5)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i4)*zb(i3,i2)*zb(i4,i3))/ - (s(i1,i2)**2*za(i2,i3)*za(i4,i5)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i4)*zb(i3,i2)**2*zb(i4,i3))/ - (2.*s(i1,i3)**2*za(i4,i5)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i5,i1))/ - (za(i2,i3)**3*za(i4,i5)) + - (2*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i5)*zb(i3,i2)*zb(i5,i1))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i2,i3)*za(i4,i5)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i2))*za(i1,i2)*za(i1,i3)*za(i2,i5)*zb(i3,i2)* - zb(i5,i1))/(s(i1,i2)*za(i2,i3)**2*za(i4,i5)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i5)*zb(i3,i2)*zb(i5,i1))/ - (s(i1,i3)*za(i2,i3)**2*za(i4,i5)) + - (xmass*za(i1,i5)*zb(i3,i2)*zb(i5,i1))/ - (2.*za(i2,i3)*za(i4,i5)*zb(i2,i1)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i2)*za(i1,i3)*za(i2,i5)*zb(i3,i2)**2*zb(i5,i1) - )/(2.*s(i1,i3)**2*za(i2,i3)*za(i4,i5)) + - (xmass*za(i2,i5)*zb(i3,i2)**2*zb(i5,i1))/ - (2.*za(i2,i3)*za(i4,i5)*zb(i2,i1)*zb(i3,i1)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i1,i5)*zb(i5,i2))/ - (za(i1,i3)*za(i2,i3)*za(i4,i5)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i5,i2))/ - (za(i1,i3)*za(i2,i3)**2*za(i4,i5)) + - (xmass*za(i1,i5)*zb(i3,i1)*zb(i5,i2))/ - (2.*za(i2,i3)*za(i4,i5)*zb(i2,i1)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i2,i5)*zb(i3,i2)*zb(i5,i2))/ - (s(i1,i2)*za(i2,i3)*za(i4,i5)) + - (xmass*za(i2,i5)*zb(i3,i2)*zb(i5,i2))/ - (2.*za(i2,i3)*za(i4,i5)*zb(i2,i1)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i5,i3))/(za(i2,i3)*za(i4,i5)) - - (2*xmass*L0(-s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i5)*zb(i3,i2)*zb(i5,i3))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i5)*zb(i3,i2)*zb(i5,i3))/ - (s(i1,i3)*za(i2,i3)*za(i4,i5)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i5)*zb(i3,i2)*zb(i5,i3))/ - (s(i1,i2)**2*za(i2,i3)*za(i4,i5)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i3))* - za(i1,i3)*za(i2,i5)*zb(i3,i2)**2*zb(i5,i3))/ - (2.*s(i1,i3)**2*za(i4,i5)) Famp(2,2,2,2)=(xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) & - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i2,i4)*zb(i4,i1))/ - (za(i1,i3)*za(i2,i3)*za(i4,i5)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i4)*zb(i4,i1))/ - (za(i1,i3)**2*za(i2,i3)*za(i4,i5)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i4)*zb(i3,i1)*zb(i4,i1))/ - (s(i1,i2)*za(i1,i3)*za(i4,i5)) + - (xmass*za(i1,i4)*zb(i3,i1)*zb(i4,i1))/ - (2.*za(i1,i3)*za(i4,i5)*zb(i2,i1)) + - (xmass*za(i2,i4)*zb(i3,i2)*zb(i4,i1))/ - (2.*za(i1,i3)*za(i4,i5)*zb(i2,i1)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i4)*zb(i4,i2))/ - (za(i1,i3)**3*za(i4,i5)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i4)*za(i2,i3)*zb(i3,i1)*zb(i4,i2))/ - (s(i1,i2)*za(i1,i3)**2*za(i4,i5)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i4)*za(i2,i3)*zb(i3,i1)*zb(i4,i2))/ - (s(i2,i3)*za(i1,i3)**2*za(i4,i5)) + - (2*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i2,i4)*zb(i3,i1)*zb(i4,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)*za(i4,i5)) - + (xmass*za(i2,i4)*zb(i3,i1)*zb(i4,i2))/ - (2.*za(i1,i3)*za(i4,i5)*zb(i2,i1)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i4)*za(i2,i3)*zb(i3,i1)**2*zb(i4,i2) - )/(2.*s(i2,i3)**2*za(i1,i3)*za(i4,i5)) + - (xmass*za(i1,i4)*zb(i3,i1)**2*zb(i4,i2))/ - (2.*za(i1,i3)*za(i4,i5)*zb(i2,i1)*zb(i3,i2)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i4)*zb(i4,i3))/(za(i1,i3)*za(i4,i5)) - - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i4)*za(i2,i3)*zb(i3,i1)*zb(i4,i3))/ - (s(i2,i3)*za(i1,i3)*za(i4,i5)) + - (2*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i4)*zb(i3,i1)*zb(i4,i3))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i4)*zb(i3,i1)*zb(i4,i3))/ - (s(i1,i2)**2*za(i1,i3)*za(i4,i5)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i4)*za(i2,i3)*zb(i3,i1)**2*zb(i4,i3))/ - (2.*s(i2,i3)**2*za(i4,i5)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i2,i5)*zb(i5,i1))/ - (za(i1,i3)*za(i2,i3)*za(i4,i5)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i5,i1))/ - (za(i1,i3)**2*za(i2,i3)*za(i4,i5)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i5)*zb(i3,i1)*zb(i5,i1))/ - (s(i1,i2)*za(i1,i3)*za(i4,i5)) - - (xmass*za(i1,i5)*zb(i3,i1)*zb(i5,i1))/ - (2.*za(i1,i3)*za(i4,i5)*zb(i2,i1)) - - (xmass*za(i2,i5)*zb(i3,i2)*zb(i5,i1))/ - (2.*za(i1,i3)*za(i4,i5)*zb(i2,i1)) - - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)**2*za(i3,i5)*zb(i5,i2))/ - (za(i1,i3)**3*za(i4,i5)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - za(i1,i2)*za(i1,i5)*za(i2,i3)*zb(i3,i1)*zb(i5,i2))/ - (s(i1,i2)*za(i1,i3)**2*za(i4,i5)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i5)*za(i2,i3)*zb(i3,i1)*zb(i5,i2))/ - (s(i2,i3)*za(i1,i3)**2*za(i4,i5)) - - (2*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i1,i2)*za(i2,i5)*zb(i3,i1)*zb(i5,i2))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i3)*za(i4,i5)) - - (xmass*za(i2,i5)*zb(i3,i1)*zb(i5,i2))/ - (2.*za(i1,i3)*za(i4,i5)*zb(i2,i1)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i2)*za(i1,i5)*za(i2,i3)*zb(i3,i1)**2*zb(i5,i2) - )/(2.*s(i2,i3)**2*za(i1,i3)*za(i4,i5)) - - (xmass*za(i1,i5)*zb(i3,i1)**2*zb(i5,i2))/ - (2.*za(i1,i3)*za(i4,i5)*zb(i2,i1)*zb(i3,i2)) + - (xmass*Lsm1(-s(i1,i2),-s(i1,i2) - s(i1,i3) - s(i2,i3), - -s(i1,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i5)*zb(i5,i3))/(za(i1,i3)*za(i4,i5)) + - (xmass*L0(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i5)*za(i2,i3)*zb(i3,i1)*zb(i5,i3))/ - (s(i2,i3)*za(i1,i3)*za(i4,i5)) - - (2*xmass*L0(-s(i2,i3),-s(i1,i2) - s(i1,i3) - s(i2,i3))* - za(i2,i5)*zb(i3,i1)*zb(i5,i3))/ - ((s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i4,i5)) + - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i1,i2))* - (s(i1,i2) + s(i1,i3) + s(i2,i3))*za(i1,i2)* - za(i3,i5)*zb(i3,i1)*zb(i5,i3))/ - (s(i1,i2)**2*za(i1,i3)*za(i4,i5)) - - (xmass*L1(-s(i1,i2) - s(i1,i3) - s(i2,i3),-s(i2,i3))* - za(i1,i5)*za(i2,i3)*zb(i3,i1)**2*zb(i5,i3))/ - (2.*s(i2,i3)**2*za(i4,i5)) do h1=1,2 do h2=1,2 do h3=1,2 do h4=1,2 amp(h1,h2,h3,h4)=Vamp(h1,h2,h3,h4)+Famp(h1,h2,h3,h4) enddo enddo enddo enddo ! write(6,*) 'vec' ! write(6,*) ! do h1=1,2 ! do h2=1,2 ! do h3=1,2 ! do h4=1,2 ! write(6,*) h1,h2,h3,h4,Vamp(h1,h2,h3,h4),Famp(h1,h2,h3,h4) ! enddo ! enddo ! enddo ! enddo ! call writeout(p) ! write(6,*) 'my tree ',amp_tree(2,2,1,2)/s(i4,i5) ! call A52(i1,i2,i3,i4,i5,za,zb) ! write(6,*) ' my imp ', amp(2,2,1,2)/s(i4,i5) ! pause return end