diff options
| author | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
|---|---|---|
| committer | Joseph Hunkeler <jhunkeler@gmail.com> | 2015-07-08 20:46:52 -0400 |
| commit | fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4 (patch) | |
| tree | bdda434976bc09c864f2e4fa6f16ba1952b1e555 /noao/imred/vtel/gauss.x | |
| download | iraf-linux-fa080de7afc95aa1c19a6e6fc0e0708ced2eadc4.tar.gz | |
Initial commit
Diffstat (limited to 'noao/imred/vtel/gauss.x')
| -rw-r--r-- | noao/imred/vtel/gauss.x | 16 |
1 files changed, 16 insertions, 0 deletions
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 |