This is the mail archive of the
gsl-discuss@sources.redhat.com
mailing list for the GSL project.
Re: [suggestion] for gsl_sf_Ci_e (in specfunc/sinint.c)
- To: Arvid Requate <arvid at Physik dot Uni-Bielefeld dot DE>
- Subject: Re: [suggestion] for gsl_sf_Ci_e (in specfunc/sinint.c)
- From: Gerard Jungman <jungman at lanl dot gov>
- Date: Tue, 29 May 2001 15:29:57 -0600
- CC: gsl-discuss at sources dot redhat dot com
- Organization: LANL T-3
- References: <20010529131222.A6062@Physik.Uni-Bielefeld.DE>
Arvid Requate wrote:
>
> while using the special function routine gsl_sf_Ci_e, I noticed
> that given an argument of zero it aborts with an error, while
> mathematically it should simply return zero.
Ci(x) has a logarithmic singularity at x=0. There may
be varying definitions, but the standard one has this
behaviour, as does the implementation in GSL.
Thanks.
--
G. Jungman