On Wed, 17 Jun 2009, Allin Cottrell wrote: > On an initial scan it seems that most things are fine but the > Bessel 'K' function goes a bit strange for large values of X (as > you say, values < 20 are OK). Now fixed, courtesy of the netlib/R code. Allin.