I have just wasted about three hours of sunday afternoon trying to compute a modified bessel function of second kind of fractional order (also known as hyperbolic Bessel function). I hope that this tiny C++ snippet will be googlable, and it will save a few minutes of someone’s live:

#include <boost/math/special_functions/bessel.hpp> #include

(…)

double v = 5.0/3.0; /* an order of Bessel function */

double x = 2.5; /* an argument of bessel function */

double r; /* will contain a value */

r = boost::math::cyl_bessel_k(v, x);

std::cout<<r<<std::endl;

Just in case someone is interested what is a modified Bessel function, this is an integral representation of Bessel funtion of order *v*: