function relbw3(ecm) c common /pawpar/ g1,em1,w1,p2,g2,em2,w2,p3,g3,em3,w3 c complex fpi,drho c parameter mpi=139.57 c s=ecm**2 c gamfact=(em1/ecm)*sqrt((s-4*mpi**2)/(em1**2-4*mpi**2))**3 gamma=w1*gamfact c drho=cmplx(s-em1**2,em1*gamma) fpi=cmplx(g1,0.0)/drho c gamfact=(em2/ecm)*sqrt((s-4*mpi**2)/(em2**2-4*mpi**2))**3 gamma=w2*gamfact c drho=cmplx(s-em2**2,em2*gamma) fpi=fpi+g2*cmplx(cos(p2),sin(p2))/drho c gamfact=(em3/ecm)*sqrt((s-4*mpi**2)/(em3**2-4*mpi**2))**3 gamma=w3*gamfact c drho=cmplx(s-em3**2,em3*gamma) fpi=fpi+g3*cmplx(cos(p3),sin(p3))/drho c relbw3=abs(fpi)**2 return end