Hahn-Exton Bessel function