function [hyge] = numHygeCdf(x,t,m,n)
  x = 0:x;
  hyge = sum(exp(gammaln(m+1)-gammaln(m-x+1)-gammaln(x+1)+gammaln(t-m+1)-gammaln(t-m-n+x+1)-gammaln(n-x+1)-gammaln(t+1)+gammaln(t-n+1)+gammaln(n+1)));
end