Hi,
the attached patch makes the bessel functions return NA if one of the arguments of the
bessel function is NA. this is needed because not all cephes functions handle NAs
properly. Examples are:
scalar d=bessel(K,0,g)
and
scalar d=bessel(K,g,0)
Christoph