function rhomint(ecm)
c
      common /pawpar/ grho,gomega,phase
c
      complex fpi,drho
c
      parameter (mrho=775.9)
      parameter (momega=781.9)
      parameter (wrho=150.5)
      parameter (womega=8.43)
c
      parameter (mpi=139.57)
c
      s=ecm**2
c
      gamfact=(mrho/ecm)*sqrt((s-4*mpi**2)/(mrho**2-4*mpi**2))**3
      gamma=wrho*gamfact
c
      drho=cmplx(s-mrho**2,mrho*gamma)
      fpi=cmplx(grho,0.0)/drho
c
      gamfact=(momega/ecm)*sqrt((s-4*mpi**2)/(momega**2-4*mpi**2))**3
      gamma=womega*gamfact
c
      drho=cmplx(s-momega**2,momega*gamma)
      fpi=fpi+gomega*cmplx(cos(phase),sin(phase))/drho
c
      rhomint=abs(fpi)**2
      return
      end