Print this page
@@ -70,11 +70,11 @@
one = 1.0L;
GENERIC
jnl(n,x) int n; GENERIC x;{
int i, sgn;
- GENERIC a, b, temp, z, w;
+ GENERIC a, b, temp = 0, z, w;
/* J(-n,x) = (-1)^n * J(n, x), J(n, -x) = (-1)^n * J(n, x)
* Thus, J(-n,x) = J(n,-x)
*/
if(n<0){
@@ -210,18 +210,20 @@
GENERIC ynl(n,x)
int n; GENERIC x;{
int i;
int sign;
- GENERIC a, b, temp;
+ GENERIC a, b, temp = 0;
- if(x!=x) return x+x;
- if (x <= zero)
+ if(x!=x)
+ return x+x;
+ if (x <= zero) {
if(x==zero)
return -one/zero;
else
return zero/zero;
+ }
sign = 1;
if(n<0){
n = -n;
if((n&1) == 1) sign = -1;
}