The R Project SVN R

Rev

Rev 59202 | Rev 64376 | Go to most recent revision | Show entire file | Ignore whitespace | Details | Blame | Last modification | View Log | RSS feed

Rev 59202 Rev 59206
Line 58... Line 58...
58 
	return(bessel_y(x, -alpha) * cos(M_PI * alpha) -
 58 
	return(bessel_y(x, -alpha) * cos(M_PI * alpha) -
59 
	       ((alpha == na) ? 0 :
 59 
	       ((alpha == na) ? 0 :
60 
		bessel_j(x, -alpha) * sin(M_PI * alpha)));
 60 
		bessel_j(x, -alpha) * sin(M_PI * alpha)));
61 
    }
 61 
    }
62 
    nb = 1+ (long)na;/* nb-1 <= alpha < nb */
 62 
    nb = 1+ (long)na;/* nb-1 <= alpha < nb */
63 
    alpha -= (double)nb-1);
 63 
    alpha -= (double)(nb-1);
64 
#ifdef MATHLIB_STANDALONE
 64 
#ifdef MATHLIB_STANDALONE
65 
    by = (double *) calloc(nb, sizeof(double));
 65 
    by = (double *) calloc(nb, sizeof(double));
66 
    if (!by) MATHLIB_ERROR("%s", _("bessel_y allocation error"));
 66 
    if (!by) MATHLIB_ERROR("%s", _("bessel_y allocation error"));
67 
#else
 67 
#else
68 
    vmax = vmaxget();
 68 
    vmax = vmaxget();