*** 56,65 ****
--- 56,66 ----
   *      values of n>1.
   *      
   */
  
  #include "libm.h"
+ #include "longdouble.h"
  #include <float.h>      /* LDBL_MAX */
  
  #define GENERIC long double
  
  static const GENERIC
*** 69,79 ****
  one  = 1.0L;
  
  GENERIC
  jnl(n,x) int n; GENERIC x;{
          int i, sgn;
!         GENERIC a, b, temp, 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){                
--- 70,80 ----
  one  = 1.0L;
  
  GENERIC
  jnl(n,x) int n; GENERIC x;{
          int i, sgn;
!         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,226 ****
  
  GENERIC ynl(n,x) 
  int n; GENERIC x;{
          int i;
          int sign;
!         GENERIC a, b, temp;
  
!         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;
          }
--- 210,229 ----
  
  GENERIC ynl(n,x) 
  int n; GENERIC x;{
          int i;
          int sign;
!         GENERIC a, b, temp = 0;
  
!         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;
          }