gsl_sf_psi_1piy

[Brian Gough]

Hiroshi Imamura writes:
 > 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);

Thanks.  I'll add those changes for the next release.
regards
Brian Gough