Added a partially working function performing LU decomposition on a square

matrix, ludcmp.awk. This is providing correct answers on _some_ test data,
further testing debugging required here.

1 files changed, 18 insertions, 0 deletions

diff --git a/binom_coeff.awk b/binom_coeff.awk
new file mode 100644
index 0000000..a19f464
--- /dev/null
+++ b/binom_coeff.awk
@@ -0,0 +1,18 @@
+#! /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))
+}