      FUNCTION BETA(Z,W)
      BETA=EXP(GAMMLN(Z)+GAMMLN(W)-GAMMLN(Z+W))
      RETURN
      END
