Exact integration for height probabilities in the abelian sandpile model

The height probabilities for the recurrent configurations in the Abelian Sandpile model on the square lattice have analytic expressions, in terms of multidimensional quadratures. At first, these quantities were evaluated numerically with high accuracy and conjectured to be certain cubic rational- coefficient polynomials in 1/π. Later their values were determined by different methods. We revert to the direct derivation of these probabilities, by computing analytically the corresponding integrals. Once again, we confirm the predictions on the probabilities, and thus, as a corollary, the conjecture on the average height, ⟨ρ⟩ = 17/8.