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, 1 insertions, 1 deletions

diff --git a/gaussian.awk b/gaussian.awk
index dc62d37..402ea10 100644
--- a/gaussian.awk
+++ b/gaussian.awk
@@ -27,7 +27,7 @@ BEGIN {
         if (N > 1 && M > 0) {

             # sigma <= 0.5

             sigma = 0.4

-            en[n] = (n - M)/(sigma*M)

+            e[n] = (n - M)/(sigma*M)

             w[n] = exp(-0.5*e[n]*e[n])

             print n, w[n]/M

         }