Print this page

        

*** 70,80 **** 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){
*** 210,227 **** 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; }