[Bug-gsl] bug report

[Yingxiang XU]

hi,
when I try to compute the modified bessel function of the second type, by using 
of the function gsl_sf_bessel_Knu(), I got some problems: for example
gsl_sf_bessel_Knu(pi,0.5) gives "domain error"
gsl_sf_bessel_Kn_scaled(3.5, 1) and gsl_sf_bessel_Kn_scaled(3.1, 1) give the 
same results (of gsl_sf_bessel_Kn_scaled(3, 1)), it seems it does not work for 
fractional order
 
How can we solve these problems?
 
Looking forward to hearing from you.
 
Yingxiang XU