[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

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

Re: [R] Bessel function with large index value

2009-11-20 Thread Martin Maechler
DS == David Scott d.sc...@auckland.ac.nz 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