# This procedure takes in n and d and computes the Ehrhart polynomial of the uniform matroid on n elements with rand d. Uses the equation in Corollary 2.2 in Katzman

uniehr := proc(n,d)
local totalSum;
with (combinat);

totalSum := 0;
for s from 0 to d-1 do
    totalSum := totalSum + expand( (-1)^s*binomial(n,s)*binomial(t*(d-s)-s+n-1,n-1  )):
od;
return totalSum;
end proc;

# This procedure tests if the uniform matroids have positive coefficients up to l

testuniehr := proc(l)
local i,j;

for i from 2 to l do
    for j from 1 to floor(l/2) do
        postest(uniehr(i,j));
    od;
    if (i mod 10 = 0) then
        printf("i=%d\n",i);
    fi;
od;

end proc;