Automatic computation of Bessel function integrals