From fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 Mon Sep 17 00:00:00 2001
From: Joseph Hunkeler <jhunkeler@gmail.com>
Date: Wed, 8 Jul 2015 20:46:52 -0400
Subject: Initial commit

---
 noao/imred/vtel/gauss.x | 16 ++++++++++++++++
 1 file changed, 16 insertions(+)
 create mode 100644 noao/imred/vtel/gauss.x

(limited to 'noao/imred/vtel/gauss.x')

diff --git a/noao/imred/vtel/gauss.x b/noao/imred/vtel/gauss.x
new file mode 100644
index 00000000..fc5f9211
--- /dev/null
+++ b/noao/imred/vtel/gauss.x
@@ -0,0 +1,16 @@
+procedure gauss (x, a, ymod, dyda, ma)
+
+real	x, a[ma], ymod, dyda[ma]
+int	ma
+
+real	arg, ex, fac
+
+begin
+	arg = (x - a(2))/a(3)
+	ex = exp(-arg**2)
+	fac = a(1)*ex*2.0*arg
+	ymod = a(1)*ex
+	dyda(1) = ex
+	dyda(2) = fac/a(3)
+	dyda(3) = fac*arg/a(3)
+end
-- 
cgit