This is the mail archive of the
gsl-discuss@sourceware.cygnus.com
mailing list for the GSL project.
Re: bug in bessel_j.c
- To: Thomas Walter <walter at pctc dot chemie dot uni-erlangen dot de>
- Subject: Re: bug in bessel_j.c
- From: Gerard Jungman <jungman at lanl dot gov>
- Date: Fri, 28 Apr 2000 16:31:59 -0600
- CC: gsl-discuss at sourceware dot cygnus dot com
- Organization: LANL T-8
- References: <200004280946.LAA18051@server2.chemie.uni-erlangen.de>
Thomas Walter wrote:
>
> I think there is a bug in 'specfunc/bessel_j.c' and friends.
> Calling 'bessel_jl_impl (0,0)' should return '1' like 'bessel_j0_impl
> (0)' does.
Oops. I just fixed it in CVS.
> PS:
> What about adding incomplete beta function?
Not likely to happen soon. In the meantime, you
could try using the relation to the confluent
hypergeometric function. That might work for you
if the arguments are not in one of the very
difficult limits.
--
G. Jungman