BESSELINT
Introduction
BESSELINT is a Matlab program to integrate arbitrary products of Bessel functions over the real halfline. Only Bessel functions of the first kind are currently supported.
Implementation
Our program was written to work in Matlab 7.0 and Octave 2.1.69. However, we specifically implemented it to be backward compatible with earlier versions of Matlab. We extensively tested our code under
 Matlab 7.2.0.294 (R2006a) for Linux
 Matlab 7.0.1.24074 (R14) Service Pack 1 for Linux
 Matlab 6.5.0.180913a (R13) for Linux and Windows
 Matlab 6.1.0.450 (R12.1) for Windows
 Octave 2.1.69
Availability
Once the main reference given below is published, it will be available at the ACM Digital Library and Netlib. In the mean time, you can obtain it by asking one of the authors.
Publications
The main reference for this software project.

J. Van Deun and R. Cools.
Algorithm 858: Computing infinite range integrals of an arbitrary product
of Bessel functions.
ACM Trans. Math. Software, Vol. 32(4) pp. 580596.
2006.
http://doi.acm.org/10.1145/1186785.1186790
Some implementation issues and related work.
 J. Van Deun and R. Cools. Integrating products of Bessel functions using the incomplete gamma function. In T.E. Simos, G. Psihoyios, and Ch. Tsitouras, editors, International Conference on Numerical Analysis and Applied Mathematics 2005, pp. 668671. WileyVCH. 2005.

J. Van Deun and R. Cools.
A Matlab implementation of an algorithm for computing integrals of products
of Bessel functions.
In Second International Congress on Mathematical Software, Lecture Notes in Computer Science, pp. 284295.
Springer, ISBN: 9783540380849.
2006.
http://dx.doi.org/10.1007/11832225_29 
J. Van Deun and R. Cools.
Note on "Electromagnetic Response of a Large Circular Loop Source on a Layered Earth:
A New Computation Method".
Pure and Applied Geophysics, Vol. 164(5) pp. 11071111.
2007.
http://dx.doi.org/10.1007/s000240070202y 
J. Van Deun and R. Cools.
A stable recurrence for the incomplete gamma function with imaginary second argument.
Numerische Mathematik, Vol. 104(4) pp. 445456.
2006.
http://dx.doi.org/10.1007/s0021100600261 
J. Van Deun and R. Cools.
Integrating products of Bessel functions with an additional exponential or rational factor.
Computer Physics Communications, Vol. 178(8) pp. 578590.
2008.
http://dx.doi.org/10.1016/j.cpc.2007.11.010