7 form1 = fpikm(sqrt(s1),ampi,ampi)
8 form1 = form1*wigfor(qq,xm2,gam2)
22 form1 = bwigk1(qq,xm2,gam2,amres,gamres,am3,am2,am1,l1,l2)
23 form1 = form1*bwigkst(s1,amres,gamres,am2,am3,l2)
49 COMPLEX FUNCTION bwigk1(S,M,G,XMRES0,XGAM0,XMRES1,XMRES2,XMLEFT
54 * ,xmres1,xmres2,xmleft
69 COMMON /k1mass/ w,mres0,gam0,mres1,mres2,mleft,n1,n2
72 wigner(a,b,c)= cmplx(1.0,0.0)/cmplx(a-b**2,b*c)
84 rhok1 = gauss(func,xmres1+xmres2,m-xmleft,eps)
86 rhos = gauss(func,xmres1+xmres2,sqrt(s)-xmleft,eps)
88 bwigk1 = (m**2)*wigner(s,m,gam)
103 * ,mres0,mres1,mres2,mleft
109 COMMON /k1mass/ w,mres0,gam0,mres1,mres2,mleft,l1,l2
114 IF (x.GT.(mres1+mres2))
THEN
115 q0 = (mres0**2-(mres1+mres2)**2)*(mres0**2-(mres1-mres2)**2)
116 q0 = sqrt(q0)/(2.0*mres0)
118 q = (x**2-(mres1+mres2)**2)*(x**2-(mres1-mres2)**2)
121 gam = gam0*(q/q0)**(2*l2+1)
122 bwig = (x-mres0)**2 + (gam/2.0)**2
123 bwig = (gam/(2.0*pival))/bwig
129 q = (w**2-(mleft+x)**2)*(w**2-(mleft-x)**2)
132 func = bwig*(q)**(2*l1+1)/w
143 COMPLEX FUNCTION bwigkst(S,M,G,XM1,XM2,L)
152 wigner(a,b,c)= cmplx(1.0,0.0)/cmplx(a-b**2,b*c)
154 IF (s.GT.(xm1+xm2)**2.AND.m.GT.(xm1+xm2))
THEN
155 q0 = (m**2-(xm1+xm2)**2)*(m**2-(xm1-xm2)**2)
156 q0 = sqrt(q0)/(2.0*m)
158 q = (s-(xm1+xm2)**2)*(s-(xm1-xm2)**2)
159 q = sqrt(q)/(2.0*sqrt(s))
161 gam = g*(q/q0)**(2*l+1)
162 bwigkst = (m**2)*wigner(s,m,gam)
164 bwigkst = cmplx(0.0,0.0)