#! /usr/bin/awk -f ### binomial coeffecient # https://rosettacode.org/wiki/Evaluate_binomial_coefficients function binom(n, k) { b = 1 for (i=1; i<(k+1); i++) { b *= (n - i + 1) / i } return b } BEGIN { ARGV[1] ? N = ARGV[1] : N = 1 ARGV[2] ? K = ARGV[2] : K = 1 ARGV[3] ? OFMT = "%." ARGV[3] "g" : OFMT = "%g" printf(OFMT ORS, binom(N, K)) }