1 #====================================================================== |
1 #====================================================================== |
2 # L I B F U N S . P L |
2 # L I B F U N S . P L |
3 # doc: Wed Mar 24 11:49:13 1999 |
3 # doc: Wed Mar 24 11:49:13 1999 |
4 # dlm: Fri Apr 16 15:58:47 2010 |
4 # dlm: Fri Sep 7 11:11:09 2012 |
5 # (c) 1999 A.M. Thurnherr |
5 # (c) 1999 A.M. Thurnherr |
6 # uE-Info: 269 45 NIL 0 0 72 2 2 4 NIL ofnI |
6 # uE-Info: 286 38 NIL 0 0 72 2 2 4 NIL ofnI |
7 #====================================================================== |
7 #====================================================================== |
8 |
8 |
9 # HISTORY: |
9 # HISTORY: |
10 # Mar 24, 1999: - copied from the c-version of NR |
10 # Mar 24, 1999: - copied from the c-version of NR |
11 # Mar 26, 1999: - added stuff for better [./fit] |
11 # Mar 26, 1999: - added stuff for better [./fit] |
12 # Sep 18, 1999: - argument typechecking |
12 # Sep 18, 1999: - argument typechecking |
13 # Oct 04, 1999: - added gauss(), normal() |
13 # Oct 04, 1999: - added gauss(), normal() |
14 # Jan 25, 2001: - added f(), sgn() |
14 # Jan 25, 2001: - added f(), sgn() |
15 # Apr 16, 2010: - added sinc() |
15 # Apr 16, 2010: - added sinc() |
|
16 # Sep 7, 2012: - added acosh() |
16 |
17 |
17 require "$ANTS/libvec.pl"; # rad() |
18 require "$ANTS/libvec.pl"; # rad() |
18 |
19 |
19 #---------------------------------------------------------------------- |
20 #---------------------------------------------------------------------- |
20 |
21 |
274 my($piX) = 3.14159265358979 * $_[0]; |
275 my($piX) = 3.14159265358979 * $_[0]; |
275 return $piX==0 ? 1 : sin($piX)/$piX; |
276 return $piX==0 ? 1 : sin($piX)/$piX; |
276 } |
277 } |
277 |
278 |
278 #---------------------------------------------------------------------- |
279 #---------------------------------------------------------------------- |
|
280 # inverse hyperbolic cosine; mathworld |
|
281 # - requires argument >= 1 |
|
282 #---------------------------------------------------------------------- |
|
283 |
|
284 sub acosh($) |
|
285 { |
|
286 return log($_[0] + sqrt($_[0]**2-1)); |
|
287 } |
|
288 |
|
289 #---------------------------------------------------------------------- |
279 |
290 |
280 1; |
291 1; |