From 309c5d8f7ff8c36235222d079955cd3783bb7ad0 Mon Sep 17 00:00:00 2001 From: wukong Date: Mon, 18 Dec 2023 12:37:31 -0800 Subject: 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. --- fib.awk | 13 +++++++++---- 1 file changed, 9 insertions(+), 4 deletions(-) (limited to 'fib.awk') 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 } -- cgit v1.2.3