“…A fast and reliable numerical evaluation of J n (x, y) would also speed up calculation of wave packet evolution in laser fields [17,18]. In nonrelativistic calculations, the generalized Bessel function has been employed mainly for strong-field ionization [19,20,21,22], but also for high-harmonic generation [23,24]. On the mathematical side, a thorough study of J n (x, y) has been initiated in a series of papers [25,26,27], and even further generalizations of the Bessel function to multiple arguments and indices have been considered [28,29,30,31].…”