Bessel_I problem



On Sun, 4 May 2008, Raymond E. Rogers wrote:

> Richard Hennessy wrote:
>> Maxima cannot evaluate bessel_i(-1/6,1).  It seems like it does not evaluate for any negative first argument.  Is there a work around for this?
>>
>> Rich
>>
>>
>>
>>
> Are you willing to use the recurrence from A&S?
> I have the book open or I can send a copy of the page.
>
> RayR
> _______________________________________________


For what it is worth, the R implementation of Bessel functions (which 
comes from netlib), claims to use 9.1.2 and 9.6.2 from A&S and gives:

> besselI(1,-1/6)
[1] 1.313664

(Index/order and argument are in the reverse order in the R 
implementation)


David Scott

_________________________________________________________________
David Scott	Department of Statistics, Tamaki Campus
 		The University of Auckland, PB 92019
 		Auckland 1142,    NEW ZEALAND
Phone: +64 9 373 7599 ext 86830		Fax: +64 9 373 7000
Email:	d.scott at auckland.ac.nz

Graduate Officer, Department of Statistics
Director of Consulting, Department of Statistics