gsl_sf_psi_1piy

Hiroshi Imamura hima@cmt.is.tohoku.ac.jp
Tue May 28 05:53:00 GMT 2002


Hi,

gsl_sf_psi_1piy() is an even function.
(see 6.3.10 of Abramowitz and Stegun)

The following changes should be made in psi.c:

line 491
    const double lny = log(y); /* FIXME: y < 0 ?? */
--->
    const double lny = log(ay);

line 500
    const double lny = log(y);
--->
    const double lny = log(ay);

=====================================
Dr. Hiroshi Imamura (hima@cmt.ims.tohoku.ac.jp)

Graduate School of Information Sciences,
Tohoku University, Aobaku-aramaki
Sendai 980-8579 Japan



More information about the Gsl-discuss mailing list