This blog is updated daily.
A general description is here.
besselI(x,.., expon.scaled=TRUE) now gives accurate results for x < 1e5 whereas it previous underflowed for x >~= 1500. Further, besselJ(x,..) also computes up to 1e5 (instead of 1e4).
permanent link
