| 1 | greg | 3.1 | #ifndef lint | 
| 2 | schorsch | 3.2 | static const char       RCSid[] = "$Id: biggerlib.c,v 3.1 2003/02/22 02:07:21 greg Exp $"; | 
| 3 | greg | 3.1 | #endif | 
| 4 |  |  | /* | 
| 5 |  |  | *  biggerlib.c - functions for an even bigger library. | 
| 6 |  |  | * | 
| 7 |  |  | *   10/2/86 | 
| 8 |  |  | */ | 
| 9 |  |  |  | 
| 10 |  |  | #include <math.h> | 
| 11 | schorsch | 3.2 |  | 
| 12 |  |  | #include "calcomp.h" | 
| 13 | greg | 3.1 |  | 
| 14 |  |  | double  argument(); | 
| 15 |  |  | static double  l_j0(), l_j1(), l_jn(), l_y0(), l_y1(), l_yn(); | 
| 16 |  |  | static double  l_erf(), l_erfc(); | 
| 17 |  |  |  | 
| 18 |  |  |  | 
| 19 |  |  | biggerlib()                     /* expand the library */ | 
| 20 |  |  | { | 
| 21 |  |  | /* the Bessel functions */ | 
| 22 |  |  | funset("j0", 1, ':', l_j0); | 
| 23 |  |  | funset("j1", 1, ':', l_j1); | 
| 24 |  |  | funset("jn", 2, ':', l_jn); | 
| 25 |  |  | funset("y0", 1, ':', l_y0); | 
| 26 |  |  | funset("y1", 1, ':', l_y1); | 
| 27 |  |  | funset("yn", 2, ':', l_yn); | 
| 28 |  |  | funset("erf", 1, ':', l_erf); | 
| 29 |  |  | funset("erfc", 1, ':', l_erfc); | 
| 30 |  |  | } | 
| 31 |  |  |  | 
| 32 |  |  |  | 
| 33 |  |  | static double | 
| 34 |  |  | l_j0() | 
| 35 |  |  | { | 
| 36 |  |  | return(j0(argument(1))); | 
| 37 |  |  | } | 
| 38 |  |  |  | 
| 39 |  |  |  | 
| 40 |  |  | static double | 
| 41 |  |  | l_j1() | 
| 42 |  |  | { | 
| 43 |  |  | return(j1(argument(1))); | 
| 44 |  |  | } | 
| 45 |  |  |  | 
| 46 |  |  |  | 
| 47 |  |  | static double | 
| 48 |  |  | l_jn() | 
| 49 |  |  | { | 
| 50 |  |  | return(jn((int)(argument(1)+.5), argument(2))); | 
| 51 |  |  | } | 
| 52 |  |  |  | 
| 53 |  |  |  | 
| 54 |  |  | static double | 
| 55 |  |  | l_y0() | 
| 56 |  |  | { | 
| 57 |  |  | return(y0(argument(1))); | 
| 58 |  |  | } | 
| 59 |  |  |  | 
| 60 |  |  |  | 
| 61 |  |  | static double | 
| 62 |  |  | l_y1() | 
| 63 |  |  | { | 
| 64 |  |  | return(y1(argument(1))); | 
| 65 |  |  | } | 
| 66 |  |  |  | 
| 67 |  |  |  | 
| 68 |  |  | static double | 
| 69 |  |  | l_yn() | 
| 70 |  |  | { | 
| 71 |  |  | return(yn((int)(argument(1)+.5), argument(2))); | 
| 72 |  |  | } | 
| 73 |  |  |  | 
| 74 |  |  |  | 
| 75 |  |  | static double | 
| 76 |  |  | l_erf() | 
| 77 |  |  | { | 
| 78 |  |  | extern double  erf(); | 
| 79 |  |  |  | 
| 80 |  |  | return(erf(argument(1))); | 
| 81 |  |  | } | 
| 82 |  |  |  | 
| 83 |  |  |  | 
| 84 |  |  | static double | 
| 85 |  |  | l_erfc() | 
| 86 |  |  | { | 
| 87 |  |  | extern double  erfc(); | 
| 88 |  |  |  | 
| 89 |  |  | return(erfc(argument(1))); | 
| 90 |  |  | } |