Hi, thanks for the report, when I tried it locally I got the following output: polytope > $m = new Matrix([1,0,0],[1,1/2,0],[1,0,1/2]); polytope > $p = new Polytope<Rational>(VERTICES=>$m); polytope > print join("\n",@{$p->EHRHART_QUASI_POLYNOMIAL}); 1/8*x^2 + 3/4*x + 1 1/8*x^2 + 1/2*x + 3...