12 |
|
|
13 |
|
#include "calcomp.h" |
14 |
|
|
15 |
< |
double argument(); |
16 |
< |
static double l_j0(), l_j1(), l_jn(), l_y0(), l_y1(), l_yn(); |
17 |
< |
static double l_erf(), l_erfc(); |
15 |
> |
double argument(int); |
16 |
> |
static double l_j0(char *), l_j1(char *), l_jn(char *); |
17 |
> |
static double l_y0(char *), l_y1(char *), l_yn(char *); |
18 |
> |
static double l_erf(char *), l_erfc(char *); |
19 |
|
|
20 |
|
|
21 |
+ |
void |
22 |
|
biggerlib() /* expand the library */ |
23 |
|
{ |
24 |
|
/* the Bessel functions */ |
34 |
|
|
35 |
|
|
36 |
|
static double |
37 |
< |
l_j0() |
37 |
> |
l_j0(char *nm) |
38 |
|
{ |
39 |
|
return(j0(argument(1))); |
40 |
|
} |
41 |
|
|
42 |
|
|
43 |
|
static double |
44 |
< |
l_j1() |
44 |
> |
l_j1(char *nm) |
45 |
|
{ |
46 |
|
return(j1(argument(1))); |
47 |
|
} |
48 |
|
|
49 |
|
|
50 |
|
static double |
51 |
< |
l_jn() |
51 |
> |
l_jn(char *nm) |
52 |
|
{ |
53 |
|
return(jn((int)(argument(1)+.5), argument(2))); |
54 |
|
} |
55 |
|
|
56 |
|
|
57 |
|
static double |
58 |
< |
l_y0() |
58 |
> |
l_y0(char *nm) |
59 |
|
{ |
60 |
|
return(y0(argument(1))); |
61 |
|
} |
62 |
|
|
63 |
|
|
64 |
|
static double |
65 |
< |
l_y1() |
65 |
> |
l_y1(char *nm) |
66 |
|
{ |
67 |
|
return(y1(argument(1))); |
68 |
|
} |
69 |
|
|
70 |
|
|
71 |
|
static double |
72 |
< |
l_yn() |
72 |
> |
l_yn(char *nm) |
73 |
|
{ |
74 |
|
return(yn((int)(argument(1)+.5), argument(2))); |
75 |
|
} |
76 |
|
|
77 |
|
|
78 |
|
static double |
79 |
< |
l_erf() |
79 |
> |
l_erf(char *nm) |
80 |
|
{ |
81 |
< |
extern double erf(); |
81 |
> |
extern double erf(double x); |
82 |
|
|
83 |
|
return(erf(argument(1))); |
84 |
|
} |
85 |
|
|
86 |
|
|
87 |
|
static double |
88 |
< |
l_erfc() |
88 |
> |
l_erfc(char *nm) |
89 |
|
{ |
90 |
< |
extern double erfc(); |
90 |
> |
extern double erfc(double x); |
91 |
|
|
92 |
|
return(erfc(argument(1))); |
93 |
|
} |