function gv=GammaConst(n)
  gv=sum(1./[1:n])-log(n);
endfunction