N = 1:100; idx = 0; for nn = N for kk = 1:nn; idx = idx + 1; gold = nchoosek(nn, kk); a1 = exp(nchoosekln(nn, kk)); a2 = exp(nchoosekln2(nn, kk)); err1(idx) = a1/gold - 1; err2(idx) = a2/gold - 1; end end [max(err1), max(err2)] [median(err1), median(err2)] [mean(err1), mean(err2)] test_exit