I would like to compute integral of polynomial functions on polytopes. I'm using Sagemath in a Jupyter notebook and I've also installed the latte-integrale package with conda. In the notebook, I manage to import "latte", define the polytope and the function, but the result of the integral is always to 0. My script is :
from sage.features.latte import Latte
Latte().is_present()
FeatureTestResult('latte_int', True)
P=polytopes.cube()
x,y,z = polygens(QQ,names='x,y,z')
f=x^2*y^2*z^2
P.integrate(f)
0
Thanks in advance!