function emrat(x)
c
      common /pawpar/ xnorm,ratio,bckg
c
      xe1p=(1+x**2)/2.0
      xm1p=(5-3*x**2)/2.0
c
      emrat=xnorm*xm1p*(1.0+ratio*xe1p/xm1p)+bckg
      return
      end