Intrinsic volumes of the quantum state space and mutually unbiased bases

This program estimates, via Monte Carlo (MC) simulation, the coefficients of the polynomial describing the volume of the ε-neighborhood of the complementarity polytope constructed from $(N+1)$ mutually orthogonal N-vertex simplices in RN²-1. Note that, in our paper, we use d in place of N, since it denotes the dimension of the original complex vector space Cᵈ. The main purpose of this program is to confirm our derived analytic formulas (Theorem 2 in the paper https://arxiv.org/abs/2511.11523) to several decimal places, and the code may be useful to other researchers for estimating further intrinsic volumes.