comb.m (view raw)
1function r = comb(n, k) 2 r = gamma(n + 1) / (gamma(k + 1) * gamma(n - k + 1)); 3end