1 files changed, 14 insertions, 5 deletions
diff --git a/gaussian.awk b/gaussian.awk
index 48c3f9a..dfb5fb4 100644
--- a/gaussian.awk
+++ b/gaussian.awk
@@ -4,19 +4,28 @@
 # generate a Gaussian window
 # https://en.wikipedia.org/wiki/Window_function
 
+
+### Gaussian
+function gaussian(n) {
+    # sigma <= 0.5
+    sigma = 0.4
+    en[n] = (n - M)/(sigma*M)
+    return exp(-0.5*e[n]*e[n])
+}
+
+
 BEGIN {
 
     ARGV[1] ? N = ARGV[1] : N = 0
     ARGV[2] ? OFMT = "%." ARGV[2] "g" : OFMT = "%g"
-    M = 0.5*(N - 1)
 
-    # sigma <= 0.5
-    sigma = 0.4
+    # window interval goes from -M to M
+    M = 0.5*(N - 1)
+    pi = 4*atan2(1,1)
 
     for (n=0; n<N; n++) {
         if (N > 1 && M > 0) {
-            e[n] = (n - M)/(sigma*M)
-            w[n] = exp(-0.5*e[n]*e[n])
+            w[n] = gaussian(n)
             print n, w[n]/M
         }
         else {