Displaying 1 result from an estimated 1 matches for "bessel_y_old".
Did you mean:
bessel_i_old
2007 Jun 18
1
two bessel function bugs for nu<0
....c bessel_i_new.c
#57c57
#< bessel_k(x, -alpha, expo) * ((ize == 1)? 2. : 2.*exp(-x))/M_PI
#---
#> bessel_k(x, -alpha, expo) * ((ize == 1)? 2. : 2.*exp(-2.0*x))/M_PI
#bug 2: besselY() for nu<0
#don't know how to check in R; a few random checks against mathematica 5.2
#fix:
#$ diff bessel_y_old.c bessel_y_new.c
#55c55
#< return(bessel_y(x, -alpha) + bessel_j(x, -alpha) * sin(-M_PI * alpha));
#---
#> return(bessel_y(x, -alpha) * cos(M_PI * alpha) - bessel_j(x,
-alpha) * sin(M_PI * alpha));
h.
--
----------------------------------
Hiroyuki Kawakatsu
Business School
Dublin City Uni...