  We present a simple and convenient analytical formula for efficient exact
computation of the hafnian of Toeplitz matrices of a special type. An
interpretation of the obtained results is given in the language of perfect
matchings and Bessel polynomials.
