This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
gsl_sf_psi_1piy
- From: Hiroshi Imamura <hima at cmt dot is dot tohoku dot ac dot jp>
- To: gsl-discuss at sources dot redhat dot com
- Date: Tue, 28 May 2002 18:47:35 +0900
- Subject: gsl_sf_psi_1piy
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