Re: [R] Bessel function with large index value

2009-11-20 Thread Martin Maechler
> "DS" == David Scott > on Sat, 21 Nov 2009 02:29:38 +1300 writes: DS> This is a reply to my own question. I thought I had found an answer but DS> it seems not so (some analysis follows below). Maybe Martin Maechler or DS> Robin Hankin or Duncan Murdoch may have some ide

Re: [R] Bessel function with large index value

2009-11-20 Thread David Scott
This is a reply to my own question. I thought I had found an answer but it seems not so (some analysis follows below). Maybe Martin Maechler or Robin Hankin or Duncan Murdoch may have some ideas---I know the question is a bit specialized. David Scott wrote: I am looking for a method of dealing

[R] Bessel function with large index value

2009-11-20 Thread David Scott
I am looking for a method of dealing with the modified Bessel function K_\nu(x) for large \nu. The besselK function implementation of this allows for dealing with large values of x by allowing for exponential scaling, but there is no facility for dealing with large \nu. What would work for