AW: [Gretl-devel] bessel K function

thank you *so* much!!! this function is much more accurate than my version. Additionaly, I propose to dump my bessel_I function, too. With the attached patch bessel_I for negative orders is calculated using the new besselK function. bessel_I for negative orders can't be calculated with the cephes function, because cephes assumes symmetry around 0 for the order (see code), which is false. Christoph ________________________________________ Von: gretl-devel-bounces(a)lists.wfu.edu [gretl-devel-bounces(a)lists.wfu.edu] im Auftrag von Allin Cottrell [cottrell(a)wfu.edu] Gesendet: Donnerstag, 18. Juni 2009 21:49 An: Gretl development Betreff: [Gretl-devel] bessel K function On Wed, 17 Jun 2009, Allin Cottrell wrote: Now fixed, courtesy of the netlib/R code. Allin. _______________________________________________ Gretl-devel mailing list Gretl-devel(a)lists.wfu.edu http://lists.wfu.edu/mailman/listinfo/gretl-devel