The Ehrhart series of magic squares of orders seven and eight
Let $\mathrm{IMS}_n(m)$ count the $n\times n$ nonnegative integer matrices whose row sums, column sums, main-diagonal sum, and antidiagonal sum are all $m$. We determine the Ehrhart series $F_n(q)=\sum_{m\geq0}\mathrm{IMS}_n(m)q^m$ as reduced rational functions for $n=7$ and $n=8$. Their numerator--denominator degrees...