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, 9 insertions, 4 deletions

diff --git a/fib.awk b/fib.awk
index 0d46479..7baef76 100644
--- a/fib.awk
+++ b/fib.awk
@@ -1,13 +1,18 @@
 #!/usr/bin/awk -f
 
 ### fib.awk
-# fib sol'n from Hamming
 
 BEGIN {
     ARGV[1] ? n = ARGV[1] : n = 0
     ARGV[2] ? OFMT = "%." ARGV[2] "g" : OFMT = "%g"
-    C1 = 1.0/sqrt(5)
-    C2 = -1.0/sqrt(5)
-    print C1*(0.5*(1 + sqrt(5)))^n + C2*(0.5*(1 - sqrt(5)))^n
+
+    # fib sol'n from Hamming
+    C1 = 1.0/sqrt(5.0)
+    C2 = -1.0/sqrt(5.0)
+    print C1*(0.5*(1 + sqrt(5.0)))^n + C2*(0.5*(1 - sqrt(5.0)))^n
+
+    phi = 0.5*(1 + sqrt(5.0))
+    binet = (phi^n - (-1.0/phi)^n)/sqrt(5.0)
+    print binet
 }