We study the joint probability generating function for k occupancy numbers on disjoint intervals in the Bessel point process. This generating function can be expressed as a Fredholm determinant. We obtain an expression for it in terms of a system of coupled Painlevé V equations, which are derived from a Lax pair of a Riemann-Hilbert problem. This generalizes a result of Tracy and Widom [24], which corresponds to the case k = 1. We also provide some examples and applications. In particular, several relevant quantities can be expressed in terms of the generating function, like the gap probability on a union of disjoint bounded intervals, the gap between the two smallest particles, and large n asymptotics for n × n Hankel determinants with a Laguerre weight possessing several jumps discontinuities near the hard edge.
2(x − y), α > −1, (1.1)where J α stands for the Bessel function of the first kind of order α (see [22, formula 10.2.2] for a definition of J α ).