PS Looks like Boost has Bessel functions: http://www.boost.org/doc/libs/1_60_0/libs/math/doc/html/math_toolkit/bessel/bessel_over.html So maybe try that. We should have access to that in the BH package. Dirk -- http://dirk.eddelbuettel.com | @eddelbuettel | edd at debian.org