Print this page

        

@@ -56,10 +56,11 @@
  *      values of n>1.
  *      
  */
 
 #include "libm.h"
+#include "longdouble.h"
 #include <float.h>      /* LDBL_MAX */
 
 #define GENERIC long double
 
 static const GENERIC

@@ -69,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){                

@@ -209,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;
         }