A numerical scheme for computing approximations to the inf-sup constant of the divergence operator in bounded Lipschitz polytopes in ℝ𝑛 is proposed. The method is based on a conforming approximation ...