samer@4: function S=hsphere_surfln(N), S=log(2)+(N/2)*log(pi)-gammaln(N/2); end
samer@4: % hsphere_surfln - log of surface area of unit hypersphere
samer@4: %
samer@4: % hsphere_surfln :: natural -> real.
samer@4: %
samer@4: % hsphere_surfln(N)=log(hsphere_surf(N))
samer@4: